Genome-Wide Analysis of C-Repeat Binding Factor Gene Family in Capsicum baccatum and Functional Exploration in Low-Temperature Response

Capsicum baccatum is a close relative of edible chili peppers (Capsicum annuum) with high economic value. The CBF gene family plays an important role in plant stress resistance physiology. We detected a total of five CBF genes in the C. baccatum genome-wide sequencing data. These genes were scattered irregularly across four chromosomes. The genes were categorized into three groupings according to their evolutionary relationships, with genes in the same category showing comparable principles for motif composition. The 2000 bp upstream of CbCBF contains many resistance-responsive elements, hormone-responsive elements, and transcription factor binding sites. These findings emphasize the crucial functions of these genes in responding to challenging conditions and physiological regulation. Analysis of tissue-specific expression revealed that CbCBF3 exhibited the greatest level of expression among all tissues. Under conditions of low-temperature stress, all CbCBF genes exhibited different levels of responsiveness, with CbCBF3 showing a considerable up-regulation after 0.25 h of cold stress, indicating a high sensitivity to low-temperature response. The importance of the CbCBF3 gene in the cold response of C. baccatum was confirmed by the use of virus-induced gene silencing (VIGS) technology, as well as the prediction of its protein interaction network. To summarize, this study conducts a thorough bioinformatics investigation of the CbCBF gene family, showcases the practicality of employing VIGS technology in C. baccatum, and confirms the significance of the CbCBF3 gene in response to low temperatures. These findings provide significant references for future research on the adaptation of C. baccatum to low temperatures.


Introduction
Plants are often subjected to a variety of abiotic stresses during growth, such as cold, heat, drought, heavy metals, and salt stress, which can limit the geographic range of plants and reduce their productivity [1].Low-temperature stress has a huge impact on plants, destroying their plasma membrane structure, reducing their photosynthesis capacity, and over-accumulating reactive oxygen species (ROS), which can retard or even stagnate their growth and development [2][3][4].During the long evolutionary process, plants have developed a series of complex and efficient regulatory mechanisms to cope with the threat of low-temperature.With the continuous development of molecular biology technology, the mechanisms by which plants cope with cold stress have been gradually clarified [5].

Identification and Characterization Analysis of the CBF Genes in C. baccatum
Based on the HMM file of CBF, C. baccatum was searched genome-wide, and the protein sequences of Arabidopsis and C. annuum CBF were used for alignment in C. baccatum genome-wide.After manual screening and removal of redundant genes, a total of five CbCBFs were identified.Physicochemical analysis revealed diverse amino acid structure for CbCBF, the length ranging from 162 amino acids (CbCBF1) to 332 amino acids (CbCBF2).The predicted molecular weights spanned from 18.4 kDa (CbCBF1) to 36.3 kDa (CbCBF2).The theoretical pI of all CbCBF proteins was less than 7, indicating their potential negative charge in an alkaline pH environment.The predicted instability indexes of all CbCBF proteins except CbCBF4 were greater than 40, signifying their predominantly unstable nature.Moreover, the grand average of hydropathicity (GRAVY) values for all CbCBF proteins were less than 0, indicating their hydrophilic nature, with CbCBF1 being the most hydrophilic.Subcellular predictions suggest that localization of CbCBF proteins in the nucleus and cytoplasm (Table 1).

Evolutionary Analysis of CbCBF Proteins
To investigate the evolutionary relationship of the CBF, an evolutionary tree of CBF proteins in C. baccatum, Arabidopsis, C. annuum, eggplant, tomato, and C. pubescens was constructed (Figure 1) and the evolutionary distances between these proteins were calculated (Table S1).The results showed that the evolutionary tree could be divided into four subgroups.CbCBF genes were distributed in three of these subgroups, with one CbCBF gene in the S1 subgroup, one CbCBF gene in the S2 subgroup, and three CbCBF genes in the S4 subgroup.All of the S3 subgroups were Arabidopsis CBF genes, and there were no Arabidopsis members distributed in the remaining subgroups.In addition, no eggplant CBF members were distributed in subgroup S4.CbCBF genes are evolutionary distant from Arabidopsis CBF genes.Among Solanaceae, CbCBFs are even more distant from eggplant CBFs.It can be concluded that the high stability of CBF genes during the evolution of Solanaceae played an indispensable and important role in life activities.
gene in the S1 subgroup, one CbCBF gene in the S2 subgroup, and three CbCBF genes in the S4 subgroup.All of the S3 subgroups were Arabidopsis CBF genes, and there were no Arabidopsis members distributed in the remaining subgroups.In addition, no eggplant CBF members were distributed in subgroup S4.CbCBF genes are evolutionary distant from Arabidopsis CBF genes.Among Solanaceae, CbCBFs are even more distant from eggplant CBFs.It can be concluded that the high stability of CBF genes during the evolution of Solanaceae played an indispensable and important role in life activities.

Conserved Motif Analysis of CbCBF Genes
Motif analysis, facilitated by the MEME tool, was conducted, and the resulting motifs were visualized in conjunction with evolutionary relationships, as explicated in Figure 2A.Five conserved motifs were identified in total, with their sequences displayed in Fig- ure 2B.Motifs 3 and 2 were found in all CbCBFs, while members of subgroup S1 and CbCBFs in subgroup S4 exhibited similar motif patterns, including motifs 5 exclusive to these two subgroups.This suggested a potential similarity in underlying functions for these two subgroups.

Conserved Motif Analysis of CbCBF Genes
Motif analysis, facilitated by the MEME tool, was conducted, and the resulting motifs were visualized in conjunction with evolutionary relationships, as explicated in Figure 2A.Five conserved motifs were identified in total, with their sequences displayed in Figure 2B.Motifs 3 and 2 were found in all CbCBFs, while members of subgroup S1 and CbCBFs in subgroup S4 exhibited similar motif patterns, including motifs 5 exclusive to these two subgroups.This suggested a potential similarity in underlying functions for these two subgroups.

Cis-Acting Element Analysis of CbCBF Genes
To unravel the regulatory mechanisms of CBF genes in C. baccatum, the cis-acting elements within the 2000 bp promoter region upstream of CbCBF genes were analyzed.The findings revealed an abundance of stress-responsive, light-responsive, hormone-responsive, and transcription factor binding sites in the promoter region of CbCBF genes.The presence of light-responsive elements upstream of all CbCBF genes suggested potential regulation by light signals.All members except CbCBF1 possessed low-temperature response elements, underscoring the critical role of the CbCBF gene family in response of C. baccatum to low-temperature.The hormone-response elements encompassed abscisic acid, auxin, gibberellin, MeJA, and salicylic acid response elements, indicating that CbCBF genes may

Cis-Acting Element Analysis of CbCBF Genes
To unravel the regulatory mechanisms of CBF genes in C. baccatum, the cis-acting elements within the 2000 bp promoter region upstream of CbCBF genes were analyzed.The findings revealed an abundance of stress-responsive, light-responsive, hormone-responsive, and transcription factor binding sites in the promoter region of CbCBF genes.The presence of light-responsive elements upstream of all CbCBF genes suggested potential regulation by light signals.All members except CbCBF1 possessed low-temperature response elements, underscoring the critical role of the CbCBF gene family in response of C. baccatum to low-temperature.The hormone-response elements encompassed abscisic acid, auxin, gibberellin, MeJA, and salicylic acid response elements, indicating that CbCBF genes may be regulated by a complex network of hormones with implications for plant growth, development, and metabolism.Additionally, the drought-associated MYB transcription factor binding sites were identified in the promoter regions of CbCBF2, CbCBF3, and CbCBF5.It meant that these CbCBF genes may be regulated by MYB transcription factors in modulating the drought physiology of C. baccatum, as demonstrated in Figure 3.

Collinearity Analysis and Chromosomal Localization of CbCBF Genes
The chromosomal localization analysis revealed that CbCBFs are situated on four chromosomes: Chr01, Chr03, Chr09, and Chr10.The nomenclature of these members, namely CbCBF1-5, was based on their positional arrangement on the chromosomes.Notably, CbCBF2 and CbCBF3 were found to be tandemly arranged gene clusters.To delve into the origin and evolutionary relationships of CbCBF genes, the collinearity of CbCBF genes was analyzed, as illustrated in Figure 4.As it unveiled, in the C. baccatum genome, there were two pairs of CbCBF collinear genes: CbCBF4/CbCBF5 and CbCBF1/CbCBF2.As a result, there were multiple duplication events within the CbCBF gene family, which contributed to the expansion of the CbCBF family.In addition, we calculated the nonsynony-

Collinearity Analysis and Chromosomal Localization of CbCBF Genes
The chromosomal localization analysis revealed that CbCBFs are situated on four chromosomes: Chr01, Chr03, Chr09, and Chr10.The nomenclature of these members, namely CbCBF1-5, was based on their positional arrangement on the chromosomes.Notably, CbCBF2 and CbCBF3 were found to be tandemly arranged gene clusters.To delve Plants 2024, 13, 549 6 of 19 into the origin and evolutionary relationships of CbCBF genes, the collinearity of CbCBF genes was analyzed, as illustrated in Figure 4.As it unveiled, in the C. baccatum genome, there were two pairs of CbCBF collinear genes: CbCBF4/CbCBF5 and CbCBF1/CbCBF2.As a result, there were multiple duplication events within the CbCBF gene family, which contributed to the expansion of the CbCBF family.In addition, we calculated the nonsynonymous/synonymous substitution rates (Ka/Ks) in the covariate gene pairs, and the results showed that the Ka/Ks values were 0.127 in CbCBF4/CbCBF5 and 0.281 in CbCBF1/CbCBF2, which is much smaller than 1.The Ka/Ks values were all much smaller than 1, so purifying selection played a major role in the evolutionary process.

Tissue-Specific Expression Patterns of CbCBF Genes
To understand the regulatory ability of CbCBF in distinct plant parts, qRT-PCR analysis was employed to examine the expression levels in the in roots, stems, and leaves of C. baccatum.The findings, illustrated in Figure 5, revealed intricate expression patterns within the CbCBF family across diverse tissues.Specifically, CbCBF2 and CbCBF3 exhibited remarkably higher expressions in contrast to other CbCBF genes in leaves and roots, the expression level of CbCBF3 gene in stems was significantly higher than that of other CbCBF genes.The CbCBF1 and CbCBF5 consistently demonstrated lower expressions compared to the other family members in various tissues.Notably, the relative expression of CbCBF3 was the highest in all the tissues, indicating its pivotal role in the environmental adaptation of C. baccatum.

Tissue-Specific Expression Patterns of CbCBF Genes
To understand the regulatory ability of CbCBF in distinct plant parts, qRT-PCR analysis was employed to examine the expression levels in the in roots, stems, and leaves of C. baccatum.The findings, illustrated in Figure 5, revealed intricate expression patterns within the CbCBF family across diverse tissues.Specifically, CbCBF2 and CbCBF3 exhibited remarkably higher expressions in contrast to other CbCBF genes in leaves and roots, the expression level of CbCBF3 gene in stems was significantly higher than that of other CbCBF genes.The CbCBF1 and CbCBF5 consistently demonstrated lower expressions compared to the other family members in various tissues.Notably, the relative expression of CbCBF3 was the highest in all the tissues, indicating its pivotal role in the environmental adaptation of C. baccatum.
ited remarkably higher expressions in contrast to other CbCBF genes in leaves and roots, the expression level of CbCBF3 gene in stems was significantly higher than that of other CbCBF genes.The CbCBF1 and CbCBF5 consistently demonstrated lower expressions compared to the other family members in various tissues.Notably, the relative expression of CbCBF3 was the highest in all the tissues, indicating its pivotal role in the environmental adaptation of C. baccatum.

Expression Patterns of CbCBF Genes in Cold Stress
To investigate the response of CbCBF genes to low temperature, we examined the expression levels of C. baccatum under low-temperature stress (4 • C) 0 h, 0.25 h, 0.5 h, 1 h, 2 h, and 4 h by qRT-PCR to understand the expression pattern of CbCBF genes under low temperature.The results showed that all CbCBF genes responded to low temperature to different degrees.Under low-temperature stress, CbCBF1 and CbCBF5 showed a significant response at 0.25 h, and the relative expression reached the highest at 0.5 h, after which the relative expression decreased significantly.CbCBF2 increased significantly at 0.5 h of the low-temperature treatment.CbCBF3 demonstrated a pronounced elevation in expression at 0.25 h, followed by a gradual decrease, ultimately returning to pre-low-temperature treatment levels at 2-4 h.In contrast, CbCBF4 experienced a noticeable decrease in its expression at 0-1 h post-low-temperature treatment and elevated greatly at 2 h.It can be seen that different CbCBF genes respond to low temperature with different mechanisms.It is noteworthy that the expression level of CbCBF3 gene increased significantly and reached the highest level at 0.25 h of low-temperature treatment, and then declined gradually after sustained cold stress, and its sensitivity to response to low-temperature was the highest within the CbCBF family, implying that it has an important role in cold stress and the early sensing of low temperature (Figure 6).

Expression Patterns of CbCBF Genes in Cold Stress
To investigate the response of CbCBF genes to low temperature, we examined the expression levels of C. baccatum under low-temperature stress (4 °C) 0 h, 0.25 h, 0.5 h, 1 h, 2 h, and 4 h by qRT-PCR to understand the expression pattern of CbCBF genes under low temperature.The results showed that all CbCBF genes responded to low temperature to different degrees.Under low-temperature stress, CbCBF1 and CbCBF5 showed a significant response at 0.25 h, and the relative expression reached the highest at 0.5 h, after which the relative expression decreased significantly.CbCBF2 increased significantly at 0.5 h of the low-temperature treatment.CbCBF3 demonstrated a pronounced elevation in expression at 0.25 h, followed by a gradual decrease, ultimately returning to pre-low-temperature treatment levels at 2-4 h.In contrast, CbCBF4 experienced a noticeable decrease in its expression at 0-1 h post-low-temperature treatment and elevated greatly at 2 h.It can be seen that different CbCBF genes respond to low temperature with different mechanisms.It is noteworthy that the expression level of CbCBF3 gene increased significantly and reached the highest level at 0.25 h of low-temperature treatment, and then declined gradually after sustained cold stress, and its sensitivity to response to low-temperature was the highest within the CbCBF family, implying that it has an important role in cold stress and the early sensing of low temperature (Figure 6).

Virus-Induced Gene Silencing in C. baccatum
Combining the outcomes from the tissue-specific expression analysis and low-temperature stress expression modeling, the potential involvement of the CbCBF3 gene was identified in the response of C. baccatum to low temperatures.To delve deeper into its role in the cold stress response of C. baccatum, the VIGS experiment was employed to silence the CbCBF3 gene.To ensure experimental rigor, four distinct treatment groups were set up, namely H2O, Empty-TRV2, TRV2-PDS, and TRV2-CbCBF3.The results on the 18th day of silencing revealed that the leaves in the TRV2-PDS treatment group showed distinct albino traits (Figure 7B).PCR assay using TRV2 primers confirmed successful virus trans- cold stress response of C. baccatum, the VIGS experiment was employed to silence the CbCBF3 gene.To ensure experimental rigor, four distinct treatment groups were set up, namely H 2 O, Empty-TRV2, TRV2-PDS, and TRV2-CbCBF3.The results on the 18th day of silencing revealed that the leaves in the TRV2-PDS treatment group showed distinct albino traits (Figure 7B).PCR assay using TRV2 primers confirmed successful virus transmission in the Empty-TRV2, TRV2-PDS, and TRV2-CbCBF3 treatment groups, with amplified bands indicating the presence of the virus.Notably, the amplified fragments of TRV2-PDS and TRV2-CbCBF3 were approximately 300 bp longer than the amplified fragment of Empty TRV2 (Figure 7A), signifying the successful transmission of silenced fragments with the virus.Gene expressions in different groups were detected by qRT-PCR, with results displayed in Figure 7C, the VIGS significantly reduced the expression levels of PDS and CbCBF3.

Silencing the CbCBF3 Gene Reduces Cold Resistance in C. baccatum
To examine the impacts of CbCBF3 silencing on the low-temperature tolerance of C. baccatum, silenced plants were subjected to low-temperature treatment at 4 • C, with the plant survival rate monitored at 0.5 h intervals (Figure 8A).It was evident that survival rate in each treatment group experienced a decreasing trend with the increase of treatment duration under low-temperature stress.Noticeably, the survival rates of the H 2 O and the Empty TRV2 treatment groups did not exhibit sharp differences at each time point, suggesting that the transmission of tobacco rattle virus did not markedly influence the ability of C. baccatum to tolerate low temperature.In addition, plants silencing the CbCBF3 gene began to die after 1.0 h of low-temperature treatment.The survival rate of CbCBF3 silenced plants at 1.5 h was notably lower in comparison to that of the control group, and this trend persisted at all subsequent time points.All CbCBF3 silenced plants died at 5.0 h of low-temperature treatment, while the survival rate of control plants remained above 60%.Additionally, CbCBF3-silenced plants exhibited more pronounced wilting than control plants starting from 0.5 h of low-temperature treatment (Figure 8B).Physiological indices were measured for control plants and CbCBF3-silenced plants under low-temperature treatment.The results showed that at 3 h of treatment at 4 • C, the leaves of CbCBF3-silenced plants accumulated more proline (PRO) and malondialdehyde (MDA), while the content of photosynthetic pigments was significantly reduced (Figure 8C).

Protein Interaction Network of CbCBF3 in C. baccatum
To further understand the function of CbCBF3 protein, its protein interaction network was predicted and analyzed, as demonstrated in Figure 9. CbCBF3 protein possessed potential interactions with abscisic acid receptor (Cbp03g18280.1),which can interact with several genes like Cbp12g06220.1,Cbp07g06420.1,Cbp05g01240.1, and Cbp08g10830.Further exploration through KEGG pathway prediction highlighted the involvement of these interacting proteins in the resistance-associated MAPK signaling pathway and the plant hormone signaling pathway.This observation suggests that CbCBF3 protein may play a role in regulating plant resistance physiology by interacting with proteins associated with abscisic acid regulation.act with several genes like Cbp12g06220.1,Cbp07g06420.1,Cbp05g01240.1, and Cbp08g10830.Further exploration through KEGG pathway prediction highlighted the involvement of these interacting proteins in the resistance-associated MAPK signaling pathway and the plant hormone signaling pathway.This observation suggests that CbCBF3 protein may play a role in regulating plant resistance physiology by interacting with proteins associated with abscisic acid regulation.

Discussion
Drought, salinity, low temperatures, and nutrient deficits are the primary environmental factors that hinder plant growth, development, and crop yields [53][54][55].To cope with abiotic stresses, plants have developed stress response mechanisms at various levels-including morphological, physiological, biochemical, cellular, and molecular-over a long period of evolution.Among them, transcription factors play important roles in plant growth and development as well as in the response process to abiotic adversity stress, where a number of transcription factors are activated and subsequently a large number of defense-related genes are transcriptionally regulated, and these changes are

Discussion
Drought, salinity, low temperatures, and nutrient deficits are the primary environmental factors that hinder plant growth, development, and crop yields [53][54][55].To cope with abiotic stresses, plants have developed stress response mechanisms at various levelsincluding morphological, physiological, biochemical, cellular, and molecular-over a long period of evolution.Among them, transcription factors play important roles in plant growth and development as well as in the response process to abiotic adversity stress, where a number of transcription factors are activated and subsequently a large number of defense-related genes are transcriptionally regulated, and these changes are essential for the establishment of defense mechanisms [56].CBF transcription factors are important members of the AP2/ERF transcription factor family and are widely involved in a variety of plant responses to adversity [57].In recent years, with the continuous progress of sequencing technology, the reference genomes of various plants have been increasing, and a large number of transcription factor gene families have been identified in different plants.C. baccatum is a plant with high potential value and low tolerance to low temperature [58].However, the CBF transcription factors in C. baccatum is not well known, and the mechanism of their regulation in response to low temperature is poorly understood.Therefore, it is necessary to identify and systematically analyze the CBF gene family in C. baccatum.
In this study, bioinformatics analysis revealed the identification of five CBF genes in the genome of C. baccatum, exhibiting irregular distribution across four chromosomes.Previous reports have shown that the number of CBF genes varies among plants.For example, there are 8 CBF genes in potato [59], 3 CBF genes in cucumber [60], 16 CBF genes in walnut [61], and at least 15 CBF genes in cotton [62].C. baccatum has fewer CBF genes compared to other plants, which may explain its weaker cold tolerance.The CbCBF2 and CbCBF3 genes formed tandemly arranged gene clusters, a phenomenon observed in the CBF family of various plants, such as Arabidopsis [63], cotton [62], and wheat [64].Analysis of predicted physicochemical properties indicated that all CbCBF proteins are acidic, hydrophilic, and unstable proteins localized in the nucleus and cytoplasm.Evolutionary analyses showed that the Arabidopsis CBF family is evolutionarily more distantly related to CbCBF, and in Solanaceae, the CbCBF family is evolutionarily more distantly related to eggplant CBF.Previous studies have shown that tandem arrangement of CBF genes play important roles in family expansion and evolution [56], and their roles in the evolution of C. baccatum need to be further investigated.Collinearity analysis unveiled two pairs of segmental duplication genes within CbCBF members, with purifying selection being important in their evolution.In previous studies, it was noted that purifying selection eliminates deleterious mutations during the evolution of a species, thereby preserving ancestral functions in duplicated genes and therefore increasing the amount of ancestral gene products through dosage effects [57,58].The CbCBF genes were all responsive to low temperature under low-temperature treatment, which may be influenced by a dose effect.This inference is consistent with previous studies in Arabidopsis [65].Analysis of promoter cis-acting elements identified numerous light-responsive, hormone-responsive, and adversity-responsive elements upstream of the CbCBF genes.Previous studies have indicated the regulatory effect of light signaling on the CBF genes and the regular changes in expression of the CBF genes with light intensity throughout the day [66].In addition, variations in light-quality affected CBF genes expression in plants, with the lower-red/farred light ratio promoting its expression, enhancing the cold hardiness of plants at low temperatures.Therefore, it is inferred that the CbCBF genes may be regulated by light, thus influencing the cold hardiness of C. baccatum, and that its specific mechanisms need to be verified through further experiments.
The CBF gene functions in the early stages of cold stress in plants and enables plants to respond rapidly to cold stress [67].Although the CBF1/2/3 genes are functionally similar, the expression patterns are quite different [68].In this study, we comprehensively analyzed the tissue-specific expression of CbCBF genes and the expression patterns under low-temperature stress, and the results showed that CbCBF genes had similar expression patterns in roots, stems, and leaves.The relative expression levels of CbCBF1 and CbCBF5 were low in all organs, and we believe that their expression must be induced by external factors, while CbCBF3 had the highest expression level in all organs, and therefore may have an important regulatory function.CbCBF genes demonstrated complex expression patterns under low-temperature stress.Specifically, CbCBF1 and CbCBF5 were similar in this aspect, with a gradual increase under low-temperature treatment, peaking at 0.5 h, followed by a slow decline.In contrast, CbCBF2 and CbCBF3 rapidly increased and reached the highest level under low-temperature conditions.CbCBF4 gene displayed a continuous decrease within 1 h upon this treatment but transiently increased at 2 h.As a result, CbCBF genes exhibited varying expression times under low-temperature stress, suggesting a feedback regulatory mechanism among CbCBF genes.Such a result is consistent with observations in Arabidopsis [63].Notably, CbCBF3 demonstrated the most significant increase in its expression at 0.25 h of low-temperature stress, aligning with the results of the tissue-specific expression analysis.It indicates that the CbCBF3 gene greatly involves cold tolerance in C. baccatum.
VIGS serves as an approach to rapidly validate the gene function [69] and has been extensively applied in diverse plants, like tomato, tobacco, Arabidopsis, maize, loofah, cotton, and more [70][71][72][73].The effectiveness of the VIGS system varies among plants due to different factors, such as virus-host interaction, plant growth status, and virus invasion [74].Despite this, no previous studies have reported the application of VIGS technology in C. baccatum.In this research, VIGS analysis revealed a significant expression reduction in PDS gene in C. baccatum and obvious albino trait on the 18th day of silencing.These results demonstrated the feasibility of VIGS technology in C. baccatum for the first time.Subsequently, CbCBF3 was silenced by VIGS technology to compute the survival rate of silenced plants at low-temperature conditions, verifying its important role in cold tolerance.
In addition, a predictive analysis was conducted on the CbCBF proteins interaction network, revealing interactions between the CbCBF protein and several abscisic-acid-Plants 2024, 13, 549 13 of 19 regulated genes.Abscisic acid exerts a pivotal role in plant responses to environmental stress [75].Previous studies have shown that abscisic acid and CBF are closely linked in signaling pathways and that abscisic acid can affect the CRT cis-acting element on the promoter of the CBF genes, thereby activating the expression of the CBF genes [38,76].Under stress conditions, abscisic acid accumulation triggers the activation of the open stomata1 (OST1) gene, enhancing CBF genes expression [77].This evidence suggests that the CbCBF3 gene may engage in a complex and tight interaction with the abscisic acid pathway, regulating the stress tolerance physiology of C. baccatum.

Genes Conserved Motifs Analysis
The amino acid sequences of CbCBF were analyzed utilizing the procedures mentioned in the MEME website (https://meme-suite.org/) (accessed on 5 July 2023) to obtain the motif annotation file.The genomic annotation file (.gff) was obtained from the C. baccatum genome database.Information on gene length and position was extracted from the genomic annotation file based on TBtools v2.052, with the analysis results visualized by the motif annotation file [78].

Cis-Acting Element Prediction
The genomic annotation file was analyzed using TBtools v2.052 for an extraction of 2000 bp sequence of the upstream promoter region of the CbCBF family.The promoter sequences were analyzed by referring to the procedures given by PlantCARE (https: //bioinformatics.psb.ugent.be/webtools/plantcare/html/)(accessed on 9 July 2023), yielding a promoter annotation file.Hormone-responsive, adversity-regulated, environment-responsive, and transcription-factor-binding-site-related regulatory elements were retained.Finally, cis-acting elements were visualized using TBtools v2.052 software [78].

Collinearity Analysis
C. baccatum genome-wide collinearity relationships were analyzed using the MCScanX integration plug-in for TBtools v2.052 software [78].The collinearity relationships were visualized using the Advanced Circos plug-in in TBtools v2.052 software to highlight the CbCBF gene collinear relationships.Then the CbCBF collinear gene pair nonsynonymous to synonymous substitution ratio was calculated using TBtools v2.052 software Simple Ka/Ks Calculator software [78].

Plant Materials and Treatments
C. baccatum (AjCB0071: A strain from C. baccatum cv."Bishop Crown") variety used in this experiment was planted in the experimental greenhouse of Southwest Forestry University (102.454308• E, 25.035229 • N).The seeds were subjected to a 24 h soaking period in pure water and subsequently germinated on moist filter paper.After germination, the seeds were transplanted into seedling trays and then transferred to 10 cm diameter planting pots once they had grown mature leaves.The culture substrate used was a mixture of peat and vermiculite in a 3:1 ratio.Planting conditions were maintained at a temperature of 25 • C and a light cycle of 16 h per day.The plants were used for subsequent experiments upon reaching the stage of 4-6 leaves.The low-temperature stress treatment program was determined by referring to the experimental method of low-temperature treatment of peppers by previous researchers [79][80][81] and the expression pattern of CBF genes in other plants under a low-temperature environment [18,58].Using an artificial climate chamber, the whole plant was subjected to low-temperature treatment.The temperature was maintained at 4 • C for varying lengths of time (0, 0.25, 0.5, 1, 2 and 4 h).Following the treatment, samples were carefully rinsed with ddH 2 O, soaked with liquid nitrogen, and kept at −80 • C for prolonged preservation.

Functional Analysis of CbCBF3 Based on VIGS
The VIGS experiment was modified based on the C. annuum VIGS system [42].Four treatment groups were established-involving the injection of H 2 O, Empty-TRV2, TRV2-PDS, and TRV2-CbCBF3-to reduce the error of the VIGS experiment on the results.Agrobacterium tumefaciens (GV3101) carrying Empty-TRV2, TRV2-PDS, TRV2-CbCBF3, and TRV1 were cultured and resuspended to OD600 = 0.8.Equal volumes of TRV2 and TRV1 were mixed to create an infiltration solution.After being protected from light for 1 h, the seedling cotyledons were injected using a needleless syringe.The injected seedlings were then incubated in the dark at 24 • C for 24 h before being returned to the normal culture environment.Upon observing the albino trait in the TRV2-PDS treatment group, virus replication was assessed using TRV2 primers, and expression levels of target genes were determined by qRT-PCR.Subsequently, the silenced plants were subjected to cold stress treatment at 4 • C, and three parallel experiments were conducted for each treatment, consisting of 10 plants in each parallel group, and the survival rate was counted and recorded every 0.5 h.

Physiological Parameter Determination
The photosynthetic pigment content was determined using the spectrophotometric method [83,84].To do this, 0.2 g of fresh leaf tissue was weighed and ground with 12 mL of 95% ethanol and a small amount of quartz sand until the tissue turned white.The mixture was then left for 5 min.The ground solution was filtered into a 25 mL brown volumetric flask using filter paper.The filter paper was rinsed with 95% ethanol, and the volume was adjusted to 25 mL with 95% ethanol.The solution was then shaken well.The absorbance of each sample was measured at 470 nm, 649 nm, and 665 nm, and the concentration of each photosynthetic pigment was calculated according to the following formula: Ca = 13.95A665− 6.99A649, Cb = 24.96A649− 7.32A665, Cx = (1000A470 − 2.05Ca − 114.8Cb)/245.Proline and malondialdehyde content were determined using kits provided by Comin Biotechnology Co. (Comin Biotechnology Co., Ltd., Suzhou, China).

Interaction Network Analysis of CbCBF3
Predictive analysis of the CbCBF3 protein interacting network was performed using the STRING website (https://cn.string-db.org)(accessed on 9 September 2023), using CbCBF3 as a lookup sequence for interacting protein prediction in the whole genome of C. baccatum.The potential regulatory mechanisms were analyzed based on the analysis results provided by the STRING website.

Conclusions
In this study, we identified five CBF gene family members in C. baccatum.The CbCBFs were subjected to evolutionary analysis, physical and chemical property prediction analysis, conserved motifs, gene structure analysis, and collinearity analysis.A large number of resistance-associated cis-acting elements were found in the promoter region of CbCBF genes.Expression pattern analysis showed that all CbCBF genes responded to low temperature.The cold-responsive function of CbCBF3 was analyzed by VIGS technology, and its proteininteraction network was analyzed.The results indicate that CbCBFs play a significant role in C. baccatum's response to adversity.This study provides a theoretical basis for breeding for improved cold tolerance in peppers.

Figure 1 .
Figure 1.Evolutionary tree of the CBF proteins from six species.Different colors on the outside of the circle represent different subgroups (S1-S4).Different color and shape icons represent different species.The numbers in the figure represent bootstrap.

Figure 1 .
Figure 1.Evolutionary tree of the CBF proteins from six species.Different colors on the outside of the circle represent different subgroups (S1-S4).Different color and shape icons represent different species.The numbers in the figure represent bootstrap.

Figure 3 .
Figure 3. Analysis of cis-acting elements of CbCBF promoters.(A) The number of cis-acting elements of CbCBF promoters.(B) The position of cis-acting elements of CbCBF promoters.

Figure 3 .
Figure 3. Analysis of cis-acting elements of CbCBF promoters.(A) The number of cis-acting elements of CbCBF promoters.(B) The position of cis-acting elements of CbCBF promoters.

Figure 4 .
Figure 4. Collinearity analysis and chromosomal localization of CbCBF genes.The red line represents collinear genes.The color of the bars represents the gene density on the chromosome, with red indicating high density and blue indicating low density.

Figure 4 .
Figure 4. Collinearity analysis and chromosomal localization of CbCBF genes.The red line represents collinear genes.The color of the bars represents the gene density on the chromosome, with red indicating high density and blue indicating low density.

Figure 5 .
Figure 5.The expression patterns of CbCBF genes in different tissues.The y-axis represents the relative expression levels, the x-axis represents the different CbCBF genes.According to one-way ANOVA and Tukey's tests (p < 0.01), different capital letters between samples indicated significant differences.The error bars represent the standard deviation.

Figure 5 .
Figure 5.The expression patterns of CbCBF genes in different tissues.The y-axis represents the relative expression levels, the x-axis represents the different CbCBF genes.According to one-way ANOVA and Tukey's tests (p < 0.01), different capital letters between samples indicated significant differences.The error bars represent the standard deviation.

Figure 6 .
Figure 6.Expression patterns of CbCBF genes under different time points of 4 °C cold stress.The yaxis represents the relative expression levels, and the x-axis represents the time points of cold stress.According to one-way ANOVA and Tukey's test (p < 0.01), different capital letters between samples indicated significant differences.The error bars represent the standard deviation.

Figure 6 .
Figure 6.Expression patterns of CbCBF genes under different time points of 4 • C cold stress.The y-axis represents the relative expression levels, and the x-axis represents the time points of cold stress.According to one-way ANOVA and Tukey's test (p < 0.01), different capital letters between samples indicated significant differences.The error bars represent the standard deviation.2.8.Virus-Induced Gene Silencing in C. baccatumCombining the outcomes from the tissue-specific expression analysis and low-temperature stress expression modeling, the potential involvement of the CbCBF3 gene was identified in the response of C. baccatum to low temperatures.To delve deeper into its role in the

Figure 7 .
Figure 7. Virus-induced gene silencing in C. baccatum.(A) Gel electrophoresis of different silencing treatment groups.M: 5000 bp marker, 1-3: Empty TRV2, 5-7: H 2 O, 9-11: TRV2-CbCBF3, 13-15: TRV2-PDS.(B) Silencing of the PDS gene resulted in leaves with the distinct albino trait, indicating the success of the silencing system.The planting pots used in the picture are directly 10 cm.(C) qRT-PCR demonstrated that silencing significantly decreased relative gene expression levels.According to one-way ANOVA and Tukey's test (p < 0.01), different capital letters between samples indicated significant differences.The error bars represent the standard deviation.

Figure 8 .Figure 8 .
Figure 8.Effect of silencing of CbCBF3 on C. baccatum under cold stress.(A) Plant survival was measured at different times of treatment at 4 °C.Each treatment group consisted of three replicate groups, Figure 8.Effect of silencing of CbCBF3 on C. baccatum under cold stress.(A) Plant survival was measured at different times of treatment at 4 • C. Each treatment group consisted of three replicate groups, with 10 plants per replicate group.(B) Phenotype of CbCBF3-silenced and control plants challenged with cold stress at 0.5 h post-treatment.The planting pots used in the picture are directly 10 cm.(C) Physiological indicators of CbCBF3-silenced and control plants challenged with cold stress at 3 h post-treatment.Chl-a: chlorophyll a; Chl-b: chlorophyll b; PRO: proline; MDA: malondialdehyde; FW: fresh weight.According to one-way ANOVA and Tukey's test (p < 0.01), different capital letters between samples indicated significant differences.The error bars represent the standard deviation.

Figure 9 .
Figure 9. Interaction network analysis of CbCBF3 protein.Nodes represent different proteins, with colored nodes indicating proteins that may have a first shell of interactions with the CbCBF3 protein and uncolored nodes indicating proteins that may have a second shell of interactions with the CbCBF3 protein.Lines represent known or predicted protein-protein interactions.The red arrow is used to identify the CbCBF3 protein.

Figure 9 .
Figure 9. Interaction network analysis of CbCBF3 protein.Nodes represent different proteins, with colored nodes indicating proteins that may have a first shell of interactions with the CbCBF3 protein and uncolored nodes indicating proteins that may have a second shell of interactions with the CbCBF3 protein.Lines represent known or predicted protein-protein interactions.The red arrow is used to identify the CbCBF3 protein.

Table 2 .
The primer sequences used in this research.