Xin Jia1†
Xing Han1†Yuan Cheng1Xiaoli Ren1Guiyan Fan1Xiaocong Jiao1Yueyue Cai1
Lu Li1
Chenwei Zhang1*Hongguang Pang1,2*- 1Modern Agricultural Science and Technology Laboratory, Department of Agriculture and Food Science, Shijiazhuang University, Shijiazhuang, China
- 2Hebei International Joint Research Center for Green Agricultural Biological Agents, Department of Agriculture and Food Science, Shijiazhuang University, Shijiazhuang, China
Introduction: The BASIC PENTACYSTEINE (BPC) family comprises plant-specific transcription factors that regulate diverse developmental programs and stress responses. Pear (Pyrus bretschneideri), an economically significant fruit crop, often experiences marked declines in fruit yield and quality under drought stress. Although BPC genes have been identified in several plant species, a comprehensive characterization of this family in pear is lacking.
Methods: In this study, we systematically characterized PbBPC genes in the pear genome using various bioinformatic approaches. We examined their expression profiles across diverse tissues and under dehydration conditions and further validated the role of PbBPC5 in drought tolerance using virus-induced gene silencing (VIGS).
Results: This study identified seven PbBPC genes in the pear genome, which were subsequently classified into three distinct groups through phylogenetic analysis. Comprehensive bioinformatics analyses were performed, examining their phylogenetic relationships, gene structures, conserved motifs, protein domains, chromosomal locations, and gene duplication events. Promoter analyses showed that all PbBPC genes contained various cis-acting elements associated with growth and development, stress response, and phytohormone signaling. Quantitative real-time PCR (qRT-PCR) showed that most PbBPC transcripts were upregulated by dehydration, with PbBPC5 exhibiting the strongest upregulation. Furthermore, subcellular localization experiments indicated that PbBPC5 was localized to the nucleus. Silencing PbBPC5 reduced drought tolerance, as indicated by more severe wilting under water deficit, lower relative water content, higher electrolyte leakage, and elevated malondialdehyde levels. PbBPC5 silencing also weakened antioxidant defenses during drought by reducing antioxidant enzyme activities. These results suggest that PbBPC5 functions on drought tolerance regulation in pear mainly by influencing reactive oxygen species scavenging.
Discussion: This study provides a genome-wide characterization of the PbBPC family and reveals PbBPC5 as a key regulator of the drought response, offering a foundation for improving pear drought tolerance through genetic approaches.
1 Introduction
Pear (Pyrus bretschneideri) is an economically, nutritionally, and medicinally important fruit crop. However, climate change and groundwater scarcity increasingly expose pear trees to drought stress, which negatively affects fruit yield and quality (Lin et al., 2025). Identifying drought-responsive genes is essential for breeding drought-resistant cultivars through genetic improvement and for supporting the sustainable development of the pear industry.
The BASIC PENTACYSTEINE (BPC)/BARLEY B RECOMBINANT (BBR) gene family comprises plant-specific transcription factors (TFs) that regulate diverse developmental processes (Ma et al., 2022; Sahu et al., 2023). BPC TFs, also known as GAGA-binding proteins, recognize and bind to GAGA repeat sequences or C-box elements (RGARAGRRAA) in gene promoter regions (Meister et al., 2004; Kooiker et al., 2005; Hecker et al., 2015). In Arabidopsis, approximately 7% of the 3-kb upstream regions of annotated genes contain at least one (GA)6 motif, and nearly 80% contain a C-box (Hecker et al., 2015). The widespread distribution of these motifs suggests that BPC TFs regulate the transcription of large numbers of genes. To date, seven BPC members have been identified in Arabidopsis. All of these proteins share a conserved C-terminal region containing five cysteine residues that recognize GA-rich boxes or C-boxes in target gene promoters (Monfared et al., 2011; Theune et al., 2019). The Arabidopsis BPC proteins are classified into three groups based on sequence similarity and N-terminal structural features: class I (AtBPC1–3), class II (AtBPC4–6), and class III (AtBPC7) (Monfared et al., 2011; Petrella et al., 2020).
BPC TFs regulate a broad range of developmental processes, including seed and ovule development, bud dormancy, endosperm growth, vegetative-to-reproductive transitions, sex determination, stomatal and leaf morphogenesis, root development, inflorescence formation, fertility, and fruit development (Hecker et al., 2015; Ezquer et al., 2016; Charlesworth, 2021; Iwasaki et al., 2021; Lloret et al., 2021; Sahu et al., 2023). The quadruple atbpc1 atbpc2 atbpc4 atbpc6 mutant exhibits extensive defects in both vegetative and reproductive growth (Monfared et al., 2011). In rice, knockout of OsBPC1 promotes seedling growth and increases grain length, highlighting its negative role in regulating yield (Gong et al., 2018). In Arabidopsis, BPCs modulate lateral root development by repressing ABSCISIC ACID INSENSITIVE4 (ABI4) expression, linking their function to the auxin pathway (Mu et al., 2017a). In cucumber, CsBPC proteins bind to the ABI3 promoter to inhibit its expression, thereby influencing seed germination (Mu et al., 2017b). Knockout of CsBPC2 further impairs root growth by disrupting gibberellin biosynthesis (Feng et al., 2023). In addition, AtBPC3 contributes to circadian rhythm regulation and leaf edge formation in Arabidopsis thaliana (Lee et al., 2022). In apple, MdBPC2 interacts with LIKE HETEROCHROMATIN PROTEIN 1 (LHP1) to repress the expression of two auxin synthesis genes (MdYUC2a and MdYUC6b), thereby modulating plant height, leaf morphology, and root growth (Zhao et al., 2024). In Marchantia, BPC TFs are also critical regulators of sex determination (Iwasaki et al., 2021). Many of these developmental roles are mediated through phytohormone signaling, including ethylene, abscisic acid (ABA), and cytokinin (Monfared et al., 2011; Simonini and Kater, 2014; Mu et al., 2017a; Shanks et al., 2018).
The functions of BPC TFs in stress responses are less understood compared with their developmental roles. Previous studies have reported that BPC gene expression can be induced by diverse stresses in flowering Chinese cabbage, cucumber, and oilseed rape (Li et al., 2019; Zhang et al., 2023; Wang et al., 2024). In Arabidopsis, BPC1 and BPC2 enhance salt tolerance by reducing β-1,4-galactan accumulation in the cell wall (Yan et al., 2021). Conversely, AtBPC2 negatively regulates tolerance to mannitol during the seedling stage by repressing a late embryogenesis abundant gene (LEA4-5) (Li et al., 2021). In tobacco, overexpression of CsBPC2 inhibits seed germination under NaCl and PEG treatment (Li et al., 2019). More recently, CsBPC2 was shown to promote salt tolerance by enhancing osmoprotectant biosynthesis, reactive oxygen species (ROS) scavenging, ion homeostasis, and ABA signaling (Li et al., 2023). Similarly, BcBPC9 has been identified as a positive regulator of cadmium stress resistance in flowering Chinese cabbage (Zhang et al., 2023). Collectively, these studies suggest that BPC TFs are involved in the responses of plants to multiple types of stress, though additional research is needed to clarify their mechanisms.
Members of the BPC gene family have been identified in several species, such as coconut palm, grape, and oilseed rape, broadening our understanding of this TF family (Lao et al., 2024; Hu et al., 2025; Zhang et al., 2025). However, little is known about BPCs in pear. In this study, we identified seven PbBPC genes from the Chinese white pear genome and analyzed their sequence features, phylogenetic relationships, protein motifs, gene structures, promoter cis-elements, chromosomal distributions, and collinearity. We further examined their expression patterns in different tissues and under dehydration treatment. PbBPC5 silencing via virus-induced gene silencing (VIGS) reduced drought tolerance in pear. Overall, we provide the first comprehensive characterization of the BPC family in pear, uncover a role for PbBPC5 in the drought response, and offer candidate genetic resources for breeding drought-resistant cultivars through genetic engineering.
2 Materials and methods
2.1 Identification of BPC gene family members in pear
To identify BPC gene family members in pear (Pyrus bretschneideri Rehd.), we first retrieved genome and proteome sequences of pear, apple (Malus domestica), and peach (Prunus persica) from the Rosaceae Genome Database (https://www.rosaceae.org/). BPC protein sequences from Arabidopsis thaliana were downloaded from the TAIR database (https://www.arabidopsis.org/) and used as queries to identify all putative BPC homologs in the genomes of pear, apple, and peach using BLASTp with an e-value of ≤ 1e−5. Subsequently, the hidden Markov model (HMM) of the GAGA_bind domain (PF06217) was downloaded from the Pfam database (http://pfam.xfam.org/) and employed to identify BPC genes using HMMER 3.0 software with default parameters (http://hmmer.org/). Finally, all putative BPC genes were verified using SMART databases. Protein lengths and genomic information for PbBPC genes were obtained from the NCBI database (https://www.ncbi.nlm.nih.gov/). Fundamental physicochemical properties, including isoelectric point (pI), molecular weight, and instability index, were calculated using the ExPASy online tool (https://www.expasy.org/). Subcellular localization of the predicted proteins was assessed using WoLF PSORT (https://wolfpsort.hgc.jp/).
2.2 Analysis of phylogenetic relationships, gene structures, and conserved motifs
Multiple sequence alignments of BPC proteins were carried out using MAFFT v7.526 with default settings (Katoh and Standley, 2013). A maximum-likelihood (ML) phylogenetic tree was then constructed using the resulting alignment in IQ-TREE2 v2.3.4 with the parameters “-m MFP -B 1000 –bnni”. The final phylogenetic tree was visualized using the Interactive Tree of Life (iTOL) web server (https://itol.embl.de). The exon–intron structures of PbBPC genes were illustrated using the Gene Structure Display Server (https://gsds.gao-lab.org/) based on an annotation file of all PbBPC genes in pear. Conserved motifs within PbBPC proteins were identified using MEME v5.0.3 with default settings (Bailey et al., 2015).
2.3 Chromosomal locations and collinearity analysis
To measure the chromosomal distribution of PbBPC members, genome annotation files were obtained from the NCBI database, and their positions were mapped using TBtools (Chen et al., 2020). Collinearity analysis was carried out using the MCScanX toolkit with default parameters to analyze gene duplication events (Wang et al., 2012). The rates of synonymous (Ks) and nonsynonymous (Ka) substitutions were estimated using DnaSP v5.10.01 (Librado and Rozas, 2009).
2.4 Cis-acting elements in the promoters of PbBPC genes
Promoter regions of PbBPC genes were analyzed for cis-acting regulatory elements utilizing the PlantCARE database (https://bioinformatics.psb.ugent.be/webtools/plantcare/html/). For each gene, a 2,000 base pair (bp) sequence upstream of the initiation codon was extracted and submitted to the PlantCARE database for cis-acting element prediction. The distribution of these cis-acting elements in the promoter regions was subsequently visualized using TBtools.
2.5 Expression patterns of PbBPC genes in various tissues of pear
The expression profiles of PbBPC genes across various pear tissues were analyzed based on the transcriptome data retrieved from the NCBI database (accession number PRJNA498777). The tissues included leaf, stem, bud, petal, sepal, ovary, and fruit. Bud samples were collected during the flower bud differentiation stage; leaf, stem, petal, sepal, and ovary samples were collected at full flowering; and fruit samples were obtained at commercial maturity (Li et al., 2019). Gene expression levels were quantified using featureCounts v2.0.6 and normalized as fragments per kilobase per million (FPKM) values (Liao et al., 2014).
2.6 Expression patterns of PbBPC genes under dehydration treatment
P. betulaefolia seedlings were cultivated in a growth chamber maintained at 23°C, with a 14 h light/10 h dark photoperiod and a light intensity of 100 μmol m−2s−1. Following a cultivation period of three months, seedlings of P. betulaefolia with uniform size were subjected to dehydration by placing them on dry filter paper at a temperature of 26°C. Leaves were collected at intervals of 0, 1, 3, 6, and 12 h, with three independent biological replicates per time point. Total RNA was isolated from the leaves using the Wolact Plant RNA Isolation Kit (Wolact, Hong Kong, China), and first-strand cDNA synthesis was performed using the Thermo RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, USA). qRT-PCR was performed using a Roche LightCycler 96 System (Roche Diagnostics, Basel, Switzerland) via the methodology described by Huo et al. (2020). Each experiment was conducted with three biological replicates, and each biological replicate was analyzed in triplicate. Tubulin was used as the internal reference gene. The 2-ΔΔCt method was used to calculate the relative expression levels of PbBPC genes. Primers for the target genes and the internal reference gene are shown in Supplementary Table S1.
2.7 Subcellular localization of PbBPC5
Mature leaves of P. bretschneideri were collected for cloning PbBPC5. The open reading frame (ORF) of PbBPC5 was amplified and inserted into the pRI101-GFP vector to construct a subcellular localization vector. The resulting recombinant plasmid (pRI101-PbBPC5) was inserted into Agrobacterium GV3101 using electroporation. Healthy tobacco leaves were co-infiltrated with GV3101 carrying the recombinant plasmid and a membrane marker (CBL1n::OFP). Nuclear DNA was stained with 4′,6-diamidino-2-phenylindole (DAPI). Subcellular localization of the PbBPC5 protein was examined using a laser scanning confocal microscope (FV1200, Olympus, Japan).
2.8 Plant transformation and drought treatment
To construct a virus-induced gene silencing (VIGS) vector, a 300-bp fragment of the PbBPC5 ORF was cloned into the tobacco rattle virus (TRV)-based vector pTRV2. The recombinant pTRV2-PbBPC5 vector and the empty pTRV1 vector were introduced separately into GV3101 by electroporation. Co-infiltration of GV3101 carrying pTRV2-PbBPC5 and pTRV1 was performed on 45-day-old pear seedlings following a previously described protocol (Han et al., 2023). Control plants were generated by co-infiltrating GV3101 containing empty pTRV2 and pTRV1 vectors. Silencing efficiency was confirmed by quantifying PbBPC5 transcript levels using qRT-PCR. Seedlings of similar size from both the silenced group (pTRV2-PbBPC5 plants) and the control group were subjected to 20 days of water deprivation following adequate irrigation. Leaves were sampled and stored at −80 °C for physiological analyses.
2.9 Determination of stress-related physiological parameters
Relative electrolyte leakage was assessed in accordance with the protocol outlined by Lutts et al. (1996), while relative water content (RWC) was determined using the methodology described by Gaxiola et al. (2001). Quantification of malondialdehyde (MDA), hydrogen peroxide (H2O2), and superoxide radical (O2–), along with the enzymatic activities of superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT), was conducted utilizing commercial assay kits provided by Suzhou Comin Biotechnology Co., Ltd., Suzhou, China. For the histochemical detection of ROS, leaves from both silenced and control plants were harvested at the end of the drought treatment and subjected to staining with nitro blue tetrazolium (NBT) and 3,3-diaminobenzidine (DAB), following established protocols (Fryer et al., 2002).
2.10 Statistical analysis
Statistical analyses were performed utilizing SPSS v27.0. Variations among means were assessed through a one-way analysis of variance (ANOVA), and significant differences were identified using Tukey’s multiple range test, with a significance threshold set at P < 0.05.
3 Results
3.1 Genome-wide identification of BPC genes in pear
Using the seven Arabidopsis BPC proteins as queries in BLASTp searches against the pear genome, we identified seven PbBPC family members. These were verified based on the presence of the conserved BPC domain and designated as PbBPC1, PbBPC2, PbBPC3, PbBPC4, PbBPC5, PbBPC6.1, and PbBPC6.2, according to their evolutionary relationships with homologs from other species. As summarized in Table 1, the PbBPC proteins ranged in size from 279 amino acids (PbBPC1 and PbBPC2) to 340 amino acids (PbBPC6.1 and PbBPC6.2). Their predicted molecular weights spanned from 30,802.48 Da (PbBPC2) to 38,137.06 Da (PbBPC5). All members exhibited isoelectric points (pI) greater than 9, indicating enrichment in alkaline amino acids. Instability indices ranged from 46.22 (PbBPC1) to 65.70 (PbBPC5), suggesting that all PbBPC proteins are relatively unstable (values exceeding 40 were considered unstable) (Guruprasad et al., 1990). Additionally, all proteins had negative GRAVY values, reflecting their hydrophilic nature. Subcellular localization predictions placed all PbBPC proteins in the nucleus, consistent with their roles as transcription factors (Table 1).
Table 1. The basic physicochemical properties of BPC gene family members in Pyrus bretschneideri.
3.2 Phylogenetic analysis of BPC family members
To investigate the evolutionary relationships of BPC proteins, a phylogenetic tree was constructed using BPC sequences from pear, apple, peach, and Arabidopsis. As shown in Figure 1, 23 BPC proteins from the four species were categorized into three distinct classes. Class I comprised eight members, including two pear genes (PbBPC1 and PbBPC2), two apple genes (MdBPC1 and MdBPC2), one peach gene (PrBPC1), and three Arabidopsis genes (AtBPC1, AtBPC2, and AtBPC3). Class II contained 14 members, including five pear genes (PbBPC3, PbBPC4, PbBPC5, PbBPC6.1, and PbBPC6.2), four apple genes (MdBPC3–6), two peach genes (PrBPC2 and PrBPC3), and three Arabidopsis genes (AtBPC4, AtBPC5, and AtBPC6). Class III consisted of a single gene, AtBPC7, from Arabidopsis. This pattern suggests that Class III BPC genes are more prone to evolutionary loss compared with those in Classes I and II. Within the phylogenetic tree, each pear BPC gene clustered most closely with its apple homolog, reflecting their shared ancestry within the Rosaceae family.
Figure 1. Phylogenetic tree of BPCs in Pyrus bretschneideri (Pb), Arabidopsis thaliana (At), Malus domestica (Md), and Prunus persica (Pp).
3.3 Conserved motif and gene structure analysis of PbBPC members in pear
To gain deeper insight into PbBPC genes, we analyzed their conserved motifs, domains, and gene structures. Motif analysis identified 10 motifs among the seven PbBPC proteins, and their compositions were more similar within each phylogenetic class than between classes (Figure 2A). All PbBPC members possessed motifs 1, 2, 3, and 9. Two Class I members (PbBPC1 and PbBPC2) shared an identical composition and order (motifs 6, 9, 3, 7, 10, 1, and 2). In Class II, PbBPC3 and PbBPC4 also had a similar motif order (motifs 3, 9, 4, 7, 5, 1, and 2), whereas PbBPC5, PbBPC6.1, and PbBPC6.2 exhibited highly similar motif compositions and arrangements (Figure 2A). All seven proteins contained the GAGA-binding superfamily domain (Figure 2B). Gene structure analysis further showed strong conservation within classes. With the exception of PbBPC4, all genes possessed untranslated regions (UTRs). PbBPC1 and PbBPC2 contained a single exon, whereas most Class II genes had three exons, except PbBPC3, which had two (Figure 2C). Collectively, these results demonstrate that the phylogenetic classification of PbBPC genes is well supported by their motif compositions and gene structures.
Figure 2. Schematics diagrams of conserved motifs, gene structure, and conserved protein domains of BPC TFs in pear. (A) Conserved motifs in PbBPC proteins. (B) Conserved domains in PbBPC proteins. (C) Gene structures of PbBPC members.
3.4 Chromosomal locations and synteny analysis of PbBPC members in pear
The chromosomal distribution of PbBPC members was examined to clarify their genomic organization. As shown in Figure 3A, five PbBPC genes were mapped to distinct chromosomes: PbBPC3 on Chr5, PbBPC5 on Chr8, PbBPC4 on Chr10, PbBPC6.2 on Chr15, and PbBPC2 on Chr16; PbBPC1 and PbBPC6.1 were located on unassembled scaffolds. Collinearity analysis within the pear genome revealed two pairs of segmentally duplicated genes (PbBPC5/PbBPC6.2 and PbBPC3/PbBPC4) (Figure 3B). Both duplicated pairs exhibited Ka/Ks ratios less than 1, indicating that they mainly underwent purifying selection (Supplementary Table S2). To further explore evolutionary relationships across species, collinearity analysis was conducted among P. bretschneideri, A. thaliana, and M. domestica. PbBPC2, PbBPC3, PbBPC4, and PbBPC6.2 showed strong collinearity with BPC genes in Arabidopsis, while 13 collinear gene pairs were identified between pear and apple, reflecting high homology and strong conservation within the Rosaceae family (Figure 3C).
Figure 3. Chromosomal distribution and collinearity analysis of PbBPC genes. (A) Locations of PbBPC genes on five pear chromosomes. (B) Intragenomic collinearity of PbBPC genes in pear, with duplicated gene pairs connected by red lines. (C) Intergenomic collinearity of BPC transcription factors among Pyrus bretschneideri, Arabidopsis thaliana, and Malus domestica, with collinear gene pairs connected by blue lines.
3.5 Analysis of cis-elements in the upstream sequences of BPC genes in pear
To better understand the potential functions and regulatory mechanisms of the PbBPC gene family, cis-acting elements were analyzed in the 2,000-bp promoter regions upstream of each PbBPC gene. This analysis identified the presence of development-related cis-elements, including the CAT-box, AT-rich element, and O2-site, underscoring the potential involvement of PbBPC genes in the regulation of plant development (Figure 4). Furthermore, the promoters were found to contain abundant stress-responsive elements, such as ARE (anaerobic induction), MBS (MYB-binding site involved in drought inducibility), and TC-rich elements (defense and stress response), among others. Additionally, several hormone-responsive elements with known links to stress regulation were identified, including ABRE (abscisic acid responsiveness), TCA-element (salicylic acid responsiveness), and TGACG-motif (MeJA responsiveness) (Figure 4). The presence of drought-related cis-elements indicates that specific PbBPC genes may play important roles in modulating drought tolerance in pear.
Figure 4. Analysis of cis-acting elements in the promoter regions of PbBPC genes.
3.6 Expression of BPC genes in different tissues of pear and under dehydration treatment
To investigate the expression patterns of PbBPC genes across pear tissues, transcriptome data from the bud, fruit, leaf, ovary, petal, sepal, and stem of Chinese white pear were analyzed. Most PbBPC genes showed high expression in bud and ovary, particularly PbBPC1, PbBPC2, PbBPC3, and PbBPC4, suggesting that they play a role in early growth and reproductive development (Figure 5A). PbBPC5 exhibited relatively high expression in leaf and fruit, whereas other members displayed low expression in these tissues. Notably, PbBPC6.1/6.2 were more highly expressed in the stem than in other tissues (Figure 5A). Because promoter analysis revealed that most PbBPC genes contained cis-elements associated with the drought response, their expression under dehydration treatment was further examined. Under dehydration treatment, the transcript levels of PbBPC genes were generally increased. PbBPC1/2, PbBPC3, PbBPC4, and PbBPC5 showed an initial rise followed by a decline, while PbBPC6.1/6.2 steadily increased, reaching 2.6-fold higher expression at 12 h compared with 0 h. Among them, PbBPC5 responded most rapidly and strongly, peaking at 1 h after dehydration with a 7.8-fold increase relative to the expression level at 0 h (Figures 5B-F). PbBPC5 was thus selected for functional studies of the role of BPC transcription factors in drought tolerance.
Figure 5. Expression patterns of PbBPC genes in pear. (A) Heatmap of RNA-seq expression of PbBPC genes across different pear tissues. (B-F) Expression patterns of PbBPC genes in leaves under dehydration treatment. Data represent mean ± SD (n=3). Different letters denote statistically significant differences (P < 0.05, Tukey’s multiple range test).
3.7 Subcellular localization of PbBPC5
As a transcription factor, PbBPC5 was predicted to function primarily within the nucleus. To verify its subcellular localization, the PbBPC5::GFP fusion construct was co-expressed with the membrane marker CBL1n::OFP in Nicotiana benthamiana leaf cells. Previous studies have shown that CBL1n::OFP is localized to the plasma membrane and was therefore used as a membrane marker in this experiment (Batistic et al., 2010). As shown in Figure 6, fluorescence from PbBPC5::GFP was detected in the nucleus, and DAPI staining further confirmed the nuclear localization. These results demonstrate that PbBPC5 is indeed a nuclear protein, consistent with its predicted role as a transcription factor.
Figure 6. Subcellular localization of PbBPC5.
3.8 Silencing of PbBPC5 reduced the drought tolerance of pear
To elucidate the function of PbBPC5 in regulating drought tolerance in pear, a VIGS assay was performed to transiently silence PbBPC5 in P. betulaefolia. PbBPC5 transcript levels were approximately 70% lower in pTRV-PbBPC5 plants than in control plants (Figure 7B). Under normal conditions, no phenotypic differences in growth were observed between the control and pTRV-PbBPC5 plants. However, after 20 days of drought stress, pTRV-PbBPC5 plants displayed more pronounced wilting and leaf curling than controls (Figure 7A). Physiological indices commonly used to assess drought damage, including RWC, electrolyte leakage, and malondialdehyde (MDA) content, were measured. Prior to drought treatment, no significant differences were detected between groups. Following drought, RWC was significantly reduced in pTRV-PbBPC5 plants relative to control plants, while both electrolyte leakage and the MDA content increased markedly in all plants, with greater increases in silenced lines compared with controls (Figures 7C-E). These results suggest that silencing PbBPC5 reduced drought tolerance in pear.
Figure 7. Effect of PbBPC5 silencing on the drought tolerance of pear. (A) Comparative phenotypic analysis of control and PbBPC5-silenced pear plants prior to and following a 20-day drought treatment. (B) The expression level of PbBPC5 was measured in both control and PbBPC5-silenced pear plants after 7 days of drought exposure. (C) RWC, (D) relative electrolyte leakage, and (E) MDA content in leaves of control and PbBPC5-silenced pear plants before and after drought treatment. Scale bar: 3cm. Data represent mean ± SD (n = 3). Different letters indicate statistically significant differences (P < 0.05, Tukey’s multiple range test).
3.9 PbBPC5 affects ROS metabolism in pear under drought stress
Drought is known to cause excessive ROS accumulation, resulting in oxidative damage to plant cells. After 20 days of drought treatment, leaves from both control and pTRV-PbBPC5 plants were subjected to NBT and DAB staining to visualize O2– and H2O2 accumulation, respectively. As shown in Figure 8A, pTRV-PbBPC5 plants exhibited more intense blue staining with NBT and stronger brown coloration with DAB compared with control plants under drought stress, indicating higher ROS levels (Figure 8A). Consistent with these results, quantitative assays revealed that O2− and H2O2 concentrations were significantly elevated in silenced plants compared with control plants under drought stress, whereas no differences were observed before treatment (Figures 8B, C). To assess antioxidant capacity, the activities of SOD, POD, and CAT were measured. Prior to drought, enzyme activities were comparable between groups; however, under stress, all three enzyme activities were markedly lower in pTRV-PbBPC5 plants (Figures 8D–F). In addition, the expression levels of genes related to antioxidant enzymes, including PbFe-SOD, PbMn-SOD, PbPOD3, and PbCAT3, were significantly reduced in silenced lines relative to controls under drought stress (Figures 8G–J). Together, these results suggest that silencing PbBPC5 decreases antioxidant capacity in pear, thereby reducing drought tolerance.
Figure 8. Effect of PbBPC5 silencing on the antioxidant defense of pear under drought stress. (A) Visualization of O2− and H2O2 accumulation in the leaves of control and PbBPC5-silenced pear plants before and after drought treatment by NBT and DAB staining, respectively. (B) O2– content, (C) H2O2 content, (D) SOD activity, (E) POD activity, and (F) CAT activity in leaves of control and PbBPC5-silenced plants before and after drought treatment. Expression levels of (G) PbFe-SOD, (H) PbMn-SOD, (I) PbPOD3, and (J) PbCAT3 in leaves of control and PbBPC5-silenced plants after drought treatment. Data represent mean ± SD (n = 3). Different letters indicate statistically significant differences (P < 0.05, Tukey’s multiple range test).
4 Discussion
In recent years, the combined effects of global warming and groundwater scarcity have intensified the adverse effects of drought on agricultural productivity (Jury and Vaux, 2005). Transcription factors are widely recognized as key regulators of plant responses to multiple environmental stresses (Xiao et al., 2017). BPC proteins comprise a plant-specific transcription factor family implicated in the regulation of plant growth and development, as well as in responses to unfavorable environments (Mu et al., 2017a, Mu et al., 2017b; Theune et al., 2017; Sahu et al., 2023). BPC family members have been identified in multiple plant species, and the number of BPC genes varies considerably among species. For example, 4, 7, 12, 25, and 8 BPC members have been reported in Arabidopsis (Meister et al., 2004), Cucumis sativus (Li et al., 2019), Brassica campestris (Zhang et al., 2023), Brassica napus (Wang et al., 2024), and Cocos nucifera (Lao et al., 2024), respectively. However, studies examining BPC family members in pear and their roles in stress responses remain limited. We systematically characterized PbBPC genes in the pear genome, examined their expression in various tissues and under dehydration, and further validated the role of PbBPC5 in drought tolerance using VIGS technology. Our results expand our knowledge of the BPC family in pear and provide valuable genetic resources for drought-resistance breeding.
Seven PbBPC members were identified in P. bretschneideri and compared with homologs from apple, peach, and Arabidopsis. Phylogenetic analysis grouped these 23 proteins into three classes, consistent with the established classification in Arabidopsis (Meister et al., 2004). All PbBPC proteins were assigned to Class I and Class II, consistent with their classification in apple and peach. Similar patterns have been reported in cucumber, where four BPC proteins were divided into two groups (Li et al., 2019). This suggests that Class III members may have been lost during evolution. Supporting this classification, analyses of gene structures and conserved motifs revealed strong similarities in gene structure and motif organization among PbBPC members within the same subgroup, indicating close evolutionary relationships among members within the same subgroup. Comparable findings in other species indicate that BPC genes are highly conserved within subgroups but divergent across groups (Zhu et al., 2025; Pang et al., 2025). Previous studies suggest that the variable N-terminal domains of BPC proteins form coiled-coil structures involved in dimerization, protein–protein interactions, and nucleolar localization (Kooiker et al., 2005; Gong et al., 2018). In our analysis, N-terminal motifs varied between subgroups, while motifs 1 and 2 at the C-terminus were highly conserved across all PbBPC proteins. These motifs contained five conserved cysteine residues predicted to form a zinc finger-like structure for recognizing GAGA motifs, although alternative mechanisms have also been proposed (Monfared et al., 2011; Theune et al., 2017; Ma et al., 2022).
Gene family expansion in plants is primarily driven by tandem and segmental duplications (Cannon et al., 2004; Song et al., 2016). To evaluate the contribution of these duplication events to PbBPC gene expansion, we conducted a collinearity analysis of the PbBPC gene family. Two segmental duplications involving PbBPC genes were identified, suggesting that the family did not expand extensively during pear evolution. Similarly, limited expansions were observed in grape, coconut, cucumber, and camellia (Li et al., 2019; Zhang et al., 2025; Lao et al., 2024; Yu et al., 2025). Comparative collinearity analysis showed stronger homology of PbBPC genes with apple than with Arabidopsis, reflecting their closer evolutionary relationship within the Rosaceae family and suggesting potentially conserved functions among related crops.
Gene expression patterns are typically correlated with the types and numbers of cis-regulatory elements present in promoter regions (Eerapagula et al., 2025). Previous studies have shown that BPC TFs regulate diverse plant growth and development processes (Sahu et al., 2023). Mutations in BPC genes can cause diverse vegetative and reproductive abnormalities, such as dwarfism, small leaves, and abnormal ovules (Monfared et al., 2011; Simonini and Kater, 2014). BPC1 and BPC2 affect ovule development and the growth of endosperm and embryos in Arabidopsis by suppressing FUSCA3 expression (Wu et al., 2020). Consistent with these findings, the promoter regions of PbBPC members comprised several cis-regulatory elements linked to plant growth and development, including CAT-box, AT-rich element, and O2-site (Wang et al., 2024). Studies of tissue-specific expression provide valuable insights into gene functions during plant growth and development (Li et al., 2019). We found that most PbBPC members were highly expressed in buds and ovaries, indicating that they might play a role in the early growth and reproductive stages of pear.
In addition, the promoter regions of PbBPC contained numerous stress-related cis-acting elements, especially in PbBPC5, which harbored multiple ABRE and MBS motifs. Previous studies showed contrasting roles of BPC proteins in stress responses; for example, CsBPC2 enhanced salt tolerance in cucumber by regulating ROS metabolism, osmotic adjustment, and ion balance (Li et al., 2023), whereas AtBPC2 suppressed osmotic stress resistance in Arabidopsis by downregulating LEA4-5 (Li et al., 2021). Although numerous studies have revealed that BPC TFs respond to diverse environmental stresses across multiple plant species, their roles in regulating stress tolerance and the underlying mechanisms remain poorly understood (Hu et al., 2020; Xian et al., 2020; Zhang et al., 2023; Sahu et al., 2023; Wang et al., 2024). We analyzed the transcriptional responses of PbBPC genes during dehydration. All PbBPC genes were induced to some degree, with PbBPC5 displaying the most pronounced response. Given the enrichment of drought-responsive elements in its promoter and its strong induction during dehydration, PbBPC5 was selected for functional analysis to further clarify its role in drought resistance in pear.
Drought stress disrupts multiple physiological and biochemical processes and is a major cause of yield reduction in crops (Zhu, 2016). To explore the role of PbBPC5 in regulating the drought response, both control plants and PbBPC5-silenced plants were subjected to drought stress. Silencing PbBPC5 significantly weakened pear drought tolerance, as evidenced by reduced RWC, increased electrolyte leakage, and elevated MDA levels compared with control plants. These indicators reflect compromised cell membrane stability, a hallmark of stress-induced cellular injury (Bajji et al., 2004; Jia et al., 2021a). These results demonstrated that PbBPC5 silencing in pear could exacerbate the cellular damage inflicted by drought stress. ROS overaccumulation is another key factor contributing to drought-induced damage (Wang et al., 2016), and enhanced ROS metabolism has been strongly associated with improved drought tolerance in many plants (Sun et al., 2018; Wang et al., 2019; Zhao et al., 2020; Huang et al., 2021; Jia et al., 2021a, Jia et al., 2021b). For example, MdATG8i enhanced apple drought tolerance by promoting ROS scavenging (Jia et al., 2021a). In our study, PbBPC5 silencing increased ROS accumulation, reduced antioxidant enzyme activities (SOD, POD, CAT), and downregulated antioxidant enzyme-related enzymes (PbFe-SOD, PbMn-SOD, PbPOD3, PbCAT3), indicating that compromised antioxidant defenses in pTRV-PbBPC5 plants led to greater cellular damage, increasing drought susceptibility. This is consistent with recent research suggesting that BPC TFs regulate ROS homeostasis under stress, such as CsBPC2, which positively regulates salt stress resistance in cucumber by modulating ROS scavenging (Li et al., 2023), and BcBPC9, which promotes the expression of antioxidant enzyme genes under Cd stress (Zhang et al., 2023). Together, these findings demonstrate that PbBPC5 silencing reduces drought resistance in pear by weakening ROS scavenging and disrupting cellular homeostasis.
5 Conclusions
In this study, we conducted a systematic analysis of the BPC gene family in pear, identifying seven PbBPC members. Using comprehensive bioinformatics approaches, we examined their phylogenetic relationships, conserved motifs, gene structures, protein domains, chromosomal distributions, collinearity, and promoter cis-elements. Given the enrichment of drought-responsive cis-elements in their promoters, we further investigated their expression under dehydration treatment. Most PbBPC genes were upregulated, with PbBPC5 showing the strongest induction. Functional validation using VIGS revealed that silencing PbBPC5 reduced drought tolerance, as pTRV-PbBPC5 plants accumulated more ROS and experienced greater cellular damage than controls. Together, these findings offer novel insights into the evolutionary and functional characteristics of the PbBPC gene family and identify PbBPC5 as a potential positive regulator of drought tolerance in pear. This work will aid future research on the roles of BPC transcription factors in stress adaptation and provide valuable genetic resources for enhancing drought resistance in pear and related species.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
XJ: Methodology, Writing – original draft, Investigation, Funding acquisition, Data curation, Conceptualization, Validation. XH: Methodology, Conceptualization, Writing – review & editing, Investigation, Data curation, Software. YC: Software, Investigation, Funding acquisition, Writing – review & editing, Data curation. XR: Funding acquisition, Investigation, Data curation, Writing – review & editing. GF: Writing – review & editing, Software, Methodology, Data curation. XCJ: Investigation, Methodology, Data curation, Writing – review & editing. YYC: Writing – review & editing, Visualization, Data curation, Validation. LL: Validation, Data curation, Visualization, Writing – review & editing. CZ: Investigation, Resources, Writing – review & editing, Supervision, Conceptualization, Project administration. HP: Writing – review & editing, Supervision, Resources, Project administration, Investigation, Conceptualization.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This work was funded by the National Natural Science Foundation of China (32402535), the Natural Science Foundation of Hebei Province (C2024106013), the Natural Science Foundation of Hebei Province (C2024106010), the Shijiazhuang University Doctoral Scientific Research Starting Foundation (23BS029), and the Shijiazhuang University Doctoral Scientific Research Starting Foundation (23BS030).
Conflict of interest
The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1752990/full#supplementary-material
References
Bailey, T. L., Johnson, J., Grant, C. E., and Noble, W. S. (2015). The MEME suite. Nucleic Acids Res. 43, 39–49. doi: 10.1093/nar/gkv416
Bajji, M., Kinet, J., and Lutts, S. (2004). The use of the electrolyte leakage method for assessing cell membrane stability as a water stress tolerance test in durum wheat. Plant Growth Regul. 36, 61–70. doi: 10.1023/A:1014732714549
Batistic, O., Waadt, R., Steinhorst, L., Held, K., and Kudla, J. (2010). CBL-mediated targeting of CIPKs facilitates the decoding of calcium signals emanating from distinct cellular stores. Plant J. 61, 211–222. doi: 10.1111/j.1365-313X.2009.04045.x
Cannon, S. B., Mitra, A., Baumgarten., A., Young., N. D., and May, G. (2004). The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 4, 10. doi: 10.1186/1471-2229-4-10
Charlesworth, D. (2021). Evolution: The oldest sex chromosomes. Curr. Biol. 31, R1585–R1588. doi: 10.1016/j.cub.2021.10.062
Chen, C., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y., et al. (2020). TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 13, 1194e1202. doi: 10.1016/j.molp.2020.06.009
Eerapagula, R., Singh, A., Maurya, A., Singh, R., Singh, N. K., and Mahato, A. K. (2025). Genome-wide analysis of NAC transcription factors in grain amaranth reveals structural diversity and regulatory features. Sci. Rep. 15, 39968. doi: 10.1038/s41598-025-23630-7
Ezquer, I., Mizzotti, C., Nguema-Ona, E., Gotté, M., Beauzamy, L., Viana, V. E., et al. (2016). The developmental regulator SEEDSTICK controls structural and mechanical properties of the Arabidopsis seed coat. Plant Cell 28, 2478–2492. doi: 10.1105/tpc.16.00454
Feng, X., Li, S., Meng, D., Di, Q., Zhou, M., Yu, X., et al. (2023). CsBPC2 is a key regulator of root growth and development. Physiol. Plant 175, e13977. doi: 10.1111/ppl.13977
Fryer, M. J., Oxborough, K., Mullineaux, P. M., and Baker, N. R. (2002). Imaging of photo-oxidative stress responses in leaves. J. Exp. Bot. 53, 1249–1254. doi: 10.1093/jexbot/53.372.1249
Gaxiola, R. A., Li, J., Undurraga, S., Dang, L. M., Allen, G. J., Alper, S. L., et al. (2001). Drought- and salt-tolerant plants result from overexpression of the AVP1 H+-pump. Proc. Natl. Acad. Sci. 98, 11444–11449. doi: 10.1073/pnas.191389398
Gong, R., Cao, H., Zhang, J., Xie, K., Wang, D., and Yu, S. (2018). Divergent functions of the GAGA-binding transcription factor family in rice. Plant J. 94, 32–47. doi: 10.1111/tpj.13837
Guruprasad, K., Reddy, B. V., and Pandit, M. W. (1990). Correlation between stability of a protein and its dipeptide composition: a novel approach for predicting in vivo stability of a protein from its primary sequence. Protein Eng. 4, 155–161. doi: 10.1093/protein/4.2.155
Han, C., Dong, H., Qiao, Q., Dai, Y., Huang, X., and Zhang, S. (2023). Comparative genomic analysis of N6-methyladenosine regulators in nine rosaceae species and functional characterization in response to drought stress in pear. Hortic. Plant J. 9, 693–704. doi: 10.1016/j.hpj.2022.09.008
Hecker, A., Brand, L. H., Peter, S., Simoncello, N., Kilian, J., Harter, K., et al. (2015). The Arabidopsis GAGA-binding factor basic pentacysteine6 recruits the polycomb-repressive complex1 component like heterochromatin protein1 to GAGA DNA motifs. Plant Physiol. 168, 1013–1024. doi: 10.1104/pp.15.00409
Hu, H., Jiang, Y., Liu, C., Zhang, Y., Chen, M., and Liu, Z. (2025). Genome-wide identification and characterization of basic pentacysteine transcription factors in. Brassica. napus. Plants 14, 1136. doi: 10.3390/plants14071136
Hu, Z., Fu, Q., Zheng, J., Zhang, A., and Wang, H. (2020). Transcriptomic and metabolomic analyses reveal that melatonin promotes melon root development under copper stress by inhibiting jasmonic acid biosynthesis. Hortic. Res. 7, 79. doi: 10.1038/s41438-020-0293-5
Huang, D., Wang, Q., Jing, G. Q., Ma, M. N., Li, C., and Ma, F. W. (2021). Overexpression of MdIAA24 improves apple drought resistance by positively regulating strigolactone biosynthesis and mycorrhization. Tree Physiol. 41, 134–146. doi: 10.1093/treephys/tpaa109
Huo, L. Q., Guo, Z. J., Zhang, Z. J., Jia, X., Sun, Y. M., Sun, X., et al. (2020). The apple autophagy-related gene MdATG9 confers tolerance to low nitrogen in transgenic apple callus. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.00423
Iwasaki, M., Kajiwara, T., Yasui, Y., Yoshitake, Y., Miyazaki, M., Kawamura, S., et al. (2021). Identification of the sex-determining factor in the liverwort Marchantia polymorpha reveals unique evolution of sex chromosomes in a haploid system. Curr. Biol. 31, 5522–5532. doi: 10.1016/j.cub.2021.10.023
Jia, X., Gong, X. Q., Jia, X. M., Li, X. P., Wang, Y., Wang, P., et al. (2021a). Overexpression of MdATG8i enhances drought tolerance by alleviating oxidative damage and promoting water uptake in transgenic apple. Int. J. Mol. Sci. 22, 5517. doi: 10.3390/ijms22115517
Jia, X., Mao, K., Wang, P., Wang, Y., Jia, X. M., Huo, L. Q., et al. (2021b). Overexpression of MdATG8i improves water use efficiency in transgenic apple by modulating photosynthesis, osmotic balance, and autophagic activity under moderate water deficit. Horticult Res. 8, 81. doi: 10.1038/s41438-021-00521-2
Jury, W. and Vaux, H. (2005). The role of science in solving the world’s emerging water problems. Proc. Natl. Acad. Sci. 102, 15715–15720. doi: 10.1073/pnas.0506467102
Katoh, K. and Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010
Kooiker, M., Airoldi, C. A., Losa, A., Manzotti, P. S., Finzi, L., Kater, M. M., et al. (2005). BASIC PENTACYSTEINE1, a GA binding protein that induces confor-mational changes in the regulatory region of the homeotic Arabidopsis gene SEEDSTICK. Plant Cell 17, 722–729. doi: 10.1105/tpc.104.030130
Lao, Z., Mao, J., Chen, R., Xu, R., Yang, Z., Wang, Y., et al. (2024). Genome-wide identification and characterization of BASIC PENTACYSTEINE transcription factors and their binding motifs in coconut palm. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1491139
Lee, Y. C., Tsai, P. T., Huang, X. X., and Tsai, H. L. (2022). Family members additively repress the ectopic expression of BASIC PENTA CYSTEINE3 to prevent disorders in Arabidopsis circadian vegetative development. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.919946
Li, S., Miao, L., Huang, B., Gao, L., He, C., Yan, Y., et al. (2019). Genome-wide identification and characterization of cucumber BPC transcription factors and their responses to abiotic stresses and exogenous phytohormones. Int. J. Mol. Sci. 2, 5048. doi: 10.3390/ijms20205048
Li, Q., Qiao, X., Yin, H., Zhou, Y., Dong, H., Qi, K., et al. (2019). Unbiased subgenome evolution following a recent whole-genome duplication in pear (Pyrus bretschneideri Rehd.). Horticult Res. 6, 34. doi: 10.1038/s41438-018-0110-6
Li, S., Sun, M., Miao, L., Di, Q., Lv, L., Yu, X., et al. (2023). Multifaceted regulatory functions of CsBPC2 in cucumber under salt stress conditions. Hortic. Res. 10, uhad051. doi: 10.1093/hr/uhad051
Li, Q., Wang, M., and Fang, L. (2021). BASIC PENTACYSTEINE2 negatively regulates osmotic stress tolerance by modulating LEA4–5 expression in Arabidopsis thaliana. Plant Physiol. Biochem. 168, 373–380. doi: 10.1016/j.plaphy.2021.10.030
Liao, Y., Smyth, G. K., and Shi, W. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930. doi: 10.1093/bioinformatics/btt656
Librado, P. and Rozas, J. (2009). DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25, 1451–1452. doi: 10.1093/bioinformatics/btp187
Lin, L., Yuan., K., Huang., X., and Zhang., S. (2025). Genome-wide identification of the Phospholipase D (PLD) gene family in Chinese white pear (Pyrus bretschneideri) and the role of PbrPLD2 in drought resistance. Plant Sci. 350, 112286. doi: 10.1016/j.plantsci.2024.112286
Lloret, A., Quesada-Traver, C., Conejero, A., Arbona, V., Gómez-Mena, C., Petri, C., et al. (2021). Regulatory circuits involving bud dormancy factor PpeDAM6. Horticult Res. 8, 261. doi: 10.1038/s41438-021-00706-9
Lutts, S., Kinet, J., and Bouharmont, J. (1996). NaCl-induced senescence in leaves of rice (Oryza sativa L.) cultivars differing in salinity resistance. Ann. Bot. 78, 389–398. doi: 10.1006/anbo.1996.0134
Ma, X. J., Yu, Y. F., Hu, Z. K., Huang, H., Li, S. J., and Yin, H. F. (2022). Characterizations of a Class-I BASIC PENTACYSTEINE gene reveal conserved roles in the transcriptional repression of genes involved in seed development. Curr. Issues Mol. Biol. 44, 4059–4069. doi: 10.3390/cimb44090278
Meister, R. J., Williams, L. A., Monfared, M. M., Gallagher, T. L., Kraft, E. A., Nelson, C. G., et al. (2004). Definition and interactions of a positive regulatory element of the Arabidopsis INNER NO OUTER promoter. Plant J. 37, 426–438. doi: 10.1046/j.1365-313x.2003.01971.x
Monfared, M. M., Simon, M. K., Meister, R. J., Roig-Villanova, I., Kooiker, M., Colombo, L., et al. (2011). Overlapping and antagonistic activities of BASIC PENTACYSTEINE genes affect a range of developmental processes in Arabidopsis. Plant J. 66, 1020–1031. doi: 10.1111/j.1365-313X.2011.04562.x
Mu, Y., Liu, Y., Bai, L., Li, S., He, C., Yan, Y., et al. (2017b). Cucumber CsBPCs regulate the expression of CsABI3 during seed germination. Front. Plant Sci. 8. doi: 10.3389/fpls.2017.00459
Mu, Y., Zou, M., Sun, X., He, B., Xu, X., Liu, Y., et al. (2017a). BASIC PENTACYSTEINE proteins repress ABSCISIC ACID INSENSITIVE4 expression via direct recruitment of the polycomb-repressive complex 2 in Arabidopsis root development. Plant Cell Physiol. 58, 607–621. doi: 10.1093/pcp/pcx006
Pang, D., Chen, C., Tian, Y., Jiang, H., Deng, S., Liu, Y., et al. (2025). Functional characterization of CsSTK in floral organ development and its positive regulation by CsBPC5 in tea plants. BMC Plant Biol. 25, 1153. doi: 10.1186/s12870-025-07266-3
Petrella, R., Caselli, F., Roig-Villanova, I., Vignati, V., Chiara, M., Ezquer, I., et al. (2020). BPC transcription factors and a Polycomb Group protein confine the expression of the ovule identity gene SEEDSTICK in Arabidopsis. Plant J. 102, 582–599. doi: 10.1111/tpj.14673
Sahu, A., Singh, R., and Verma, P. K. (2023). Plant BBR/BPC transcription factors: unlocking multilayered regulation in development, stress and immunity. Planta 258, 31. doi: 10.1007/s00425-023-04188-y
Shanks, C. M., Heckerm, A., Cheng, C. Y., Brand, L., Collani, S., Schmid, M., et al. (2018). Role of BASIC PENTACYSTEINE transcription factors in a subset of cytokinin signaling responses. Plant J. 95, 458–473. doi: 10.1111/tpj.13962
Simonini, S. and Kater, M. M. (2014). Class I BASIC PENTACYSTEINE factors regulate HOMEOBOX genes involved in meristem size maintenance. J. Exp. Bot. 65, 1455–1465. doi: 10.1093/jxb/eru003
Song, L., Nguyen, N., Deshmukh, R. K., Patil, G. B., Prince, S. J., Valliyodan, B., et al. (2016). Soybean TIP gene family analysis and characterization of GmTIP1;5 and GmTIP2;5 water transport activity. Front. Plant Sci. 7. doi: 10.3389/fpls.2016.01564
Sun, X., Wang, P., Jia, X., Huo, L. Q., Che, R. M., and Ma, F. W. (2018). Improvement of drought tolerance by overexpressing MdATG18a is mediated by modified antioxidant system and activated autophagy in transgenic apple. Plant Biotechnol. J. 16, 545–557. doi: 10.1111/pbi.12794
Theune, M. L., Bloss, U., Brand, L. H., Ladwig, F., and Wanke, D. (2019). Phylogenetic analyses and GAGA-motif binding studies of BBR/BPC proteins lend to clues in GAGA-motif recognition and a regulatory role in brassinosteroid signaling. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.00466
Theune, M. L., Hummel, S., Jaspert, N., Lafos, M., and Wanke, D. (2017). Dimerization of the BASIC PENTACYSTEINE domain in plant GAGA-factors is mediated by disulfide bonds and required for DNA-binding. J. Adv. Plant Biol. 1, 26–39. doi: 10.14302/issn.2638-4469.japb-17-1563
Wang, L., Chen, W., Zhao, Z., Li, H., Pei, D., Huang, Z., et al. (2024). Genome-wide identification, conservation, and expression pattern analyses of the BBR-BPC gene family under abiotic stress in Brassica napus L. Genes 16, 36. doi: 10.3390/genes16010036
Wang, Q., Liu, C. H., Huang, D., Dong, Q. L., Li, P. M., and Ma, F. W. (2019). High-efficient utilization and uptake of N contribute to higher NUE of ‘Qinguan’ apple under drought and N-deficient conditions compared with ‘Honeycrisp’. Tree Physiol. 39, 1880–1895. doi: 10.1093/treephys/tpz093
Wang, Q. J., Sun, H., Dong, Q. L., Sun, T. Y., Jin, Z. X., Hao, Y. J., et al. (2016). The enhancement of tolerance to salt and cold stresses by modifying the redox state and salicylic acid content via the cytosolic malate dehydrogenase gene in transgenic apple plants. Plant Biotechnol. J. 14, 1986–1997. doi: 10.1111/pbi.12556
Wu, J., Mohamed, D., Dowhanik, S., Petrella, R., Gregis, V., Li, J., et al. (2020). Spatiotemporal restriction of FUSCA3 expression by class I BPCs promotes ovule development and coordinates embryo and endosperm growth. Plant Cell 32, 1886–1904. doi: 10.1105/tpc.19.00764
Xian, J., Wang, Y., Niu, K., Ma, H., and Ma, X. (2020). Transcriptional regulation and expression network responding to cadmium stress in a Cd-tolerant perennial grass Poa Pratensis. Chemosphere 250, 126158. doi: 10.1016/j.chemosphere.2020.126158
Xiao, J., Jin, R., Yu, X., Shen, M., Wagner, J. D., Pai, A., et al. (2017). Cis and trans determinants of epigenetic silencing by polycomb repressive complex 2 in Arabidopsis. Nat. Genet. 49, 1546–1552. doi: 10.1038/ng.3937
Yan, J., Liu, Y., Yang, L., He, H., Huang, Y., Fang, L., et al. (2021). Cell wall β-1,4-galactan regulated by the BPC1/BPC2-GALS1 module aggravates salt sensitivity in Arabidopsis thaliana. Mol. Plant 14, 411–425. doi: 10.1016/j.molp.2020.11.023
Yu, Y., Chu, X., Ma, X., Huang, M., Hu, Z., Li, S., et al. (2025). Diverse roles for a class II BPC gene in Camellia japonica through tissue-specific regulation of gene expression. Int. J. Biol. Macromol. 311, 144035. doi: 10.1016/j.ijbiomac.2025.144035
Zhang, S., Wang, J., Feng, Y., Xue, Y., Wang, Y., Zhao, M., et al. (2023). Genome-wide identification of flowering Chinese cabbage BPC family genes and BcBPC9 functional analysis in Cd stress tolerance. Plant Stress 10, 100220. doi: 10.1016/j.stress.2023.100220
Zhang, S., Zhong, H., Zhang, F., Zheng, J., Zhang, C., Yadav, V., et al. (2025). Identification of grapevine BASIC PENTACYSTEINE transcription factors and functional characterization of VvBPC1 in ovule d-evelopment. Plant Sci. 356, 112491. doi: 10.1016/j.plantsci.2025.112491
Zhao, S., Gao, H. B., Luo, J. W., Wang, H. B., Dong, Q. L., Wang, Y., et al. (2020). Genome-wide analysis of the light-harvesting chlorophyll a/b-binding gene family in apple (Malus domestica) and functional characterization of MdLhcb4.3, which confers tolerance to drought and osmotic stress. Plant Physiol. Bioch. 154, 517–529. doi: 10.1016/j.plaphy.2020.06.022
Zhao, H. Y., Wan, S., Huang, Y., Li, X., Jiao, T., Zhang, Z., et al. (2024). The transcription factor MdBPC2 alters apple growth and promotes dwarfing by regulating auxin biosynthesis. Plant Cell 36, 585–604. doi: 10.1093/plcell/koad297
Zhu, J. K. (2016). Abiotic stress signaling and responses in plants. Cell 167, 313–324. doi: 10.1016/j.cell.2016.08.029
Keywords: BPC transcription factor, drought stress, genome-wide analysis, pear, ROS
Citation: Jia X, Han X, Cheng Y, Ren X, Fan G, Jiao X, Cai Y, Li L, Zhang C and Pang H (2026) Genome-wide characterization of BPC transcription factors in pear and functional validation of PbBPC5 in drought tolerance regulation. Front. Plant Sci. 16:1752990. doi: 10.3389/fpls.2025.1752990
Received: 24 November 2025; Accepted: 29 December 2025; Revised: 19 December 2025;
Published: 21 January 2026.
Edited by:
Wei Qian, Chinese Academy of Agricultural Sciences, ChinaReviewed by:
Qian Wang, Guizhou University, ChinaQinglong Dong, Hebei Agricultural University, China
Copyright © 2026 Jia, Han, Cheng, Ren, Fan, Jiao, Cai, Li, Zhang and Pang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chenwei Zhang, Y3d6aGFuZzQxOTJAZm94bWFpbC5jb20=; Hongguang Pang, aGdwYW5nQHNqemMuZWR1LmNu
†These authors have contributed equally to this work