The Chinese white pear genome database was obtained from the Pear Genome Project.¹ To identify putative BGLU genes from Chinese white pear, several approaches were employed. We downloaded the HMM profile (PF00232) from the Pfam website² (Finn et al., 2010) and searched for candidate genes using the HMM with E-values < 1e-10 from the Chinese white pear genome (Zhou et al., 2016). After performing multiple sequence alignment, the retrieved domain sequences of the pear are aligned, and the HMM of the BGLU genes family of the pear is constructed again. Using the constructed new hidden Markov model with E-values < 1e-10 to blasted against the Chinese pear genome for carryout PbBGLU genes. The online site SMART³ (Letunic et al., 2012) was used to determine the domain PF00232 of BGLU. The online analysis tool EXPASY⁴ was used to predict and analyze the physical and chemical properties of the obtained BGLU amino acids sequences.

Sequence alignment of PbBGLUs proteins was done by using the MUSCLE method in MEGA 7.0 software. The phylogenetic tree was constructed with MEGA 7.0 software using the neighbor-joining (NJ) (bootstrap = 2,500) (Kumar et al., 2016). The conservative motifs were analyzed by MEME⁵ software the maximum value of the motif was set to 20, and the motif length was set between 6 and 50 (Bailey et al., 2015). The Gene Structure Server⁶ was used to generate exon-intron maps (Hu et al., 2015).

Phylogenetic analysis of BGLU gene family in Arabidopsis, P. bretschneideri, Oryza sative, and Populus. The software MEGA 7.0 was used to construct the phylogenetic tree. The amino-acid sequences of Arabidopsis, O. sative, and Populus were obtained from phytozome database.⁷

The materials were from “Dangshan Su” pear (P. bretschneideri cv. Dangshan Su) with the same growing environment for 40 years in Dangshan County, Anhui Province, China. We collected the mature leaves, buds, stem segments, and flowers on March 30, 2021; and collected the fruit on April 15 (15 days after flowering), 2021; May 14 (39 DAF), 2021; May 22 (47 DAF), 2021; June 4 (55 DAF), 2021; June 12 (63 DAF), 2021; June 28 (79 DAF), 2021; and August 30 (145 DAF), 2021. Each sample was frozen in liquid nitrogen and brought back to the laboratory for use. All samples were stored at–80°C for subsequent RNA extraction.

Arabidopsis thaliana L. wild type was purchased from the American Arabidopsis Biological Resources Center. Arabidopsis At1g61810 bglu45-2 mutant seeds corresponding to the Salk_117269 flanking sequence tag were ordered from the Salk Institute Genomic Analysis Laboratory. Homozygous BGLU45-2 mutant plants were isolated by genotype analysis with Arabidopsis BGLU45 gene and T-DNA specific primers, as described by Chapelle et al. (2012). The Arabidopsis plants were cultivated in a growth chamber at 23°C under 16/8 h light/dark cycles (constant illumination 100 μEm-2s-1). Ten replicates of each line were planted, and three similarly growing plants were collected for further analysis.

The total RNA of all the used “Dangshan Su” pear fruits in this study was extracted using a Plant RNA extraction kit from Tiangen (Beijing, China). Total RNA (1 μg) from each sample was used in reverse transcription. First-strand cDNA was synthesized with a PrimeScript™ RT reagent kit with a gDNA Eraser kit (TaKaRa, Kyoto, Japan). All RT-qPCR primers were designed using Primer Premier 5.0 software (Supplementary Table 3). We chose the tubulin (No. AB239680.1) served as the reference gene for RT-qPCR analysis (Wu T. et al., 2012). The RT-qPCR was performed with a CFX96 Touch™ Real-Time PCR Detection System (BIO-RAD), each sample was subjected to three biological replicates. The relative changes in gene expression levels were calculated using the 2^{–ΔΔ CT} method (Livak and Schmittgen, 2001).

The three genes RT-PCR primers were designed using Primer Premier 5.0 software (Supplementary Table 4). The full-length CDS of PbBGLU1, PbBGLU15, and PbBGLU16 without stop codon were cloned based on genomic information and constructed into pCAMBIA1305 vectors (Clontech, Beijing, country-region China) between the XbaI and BamHI (NEB) sites which have both CaMV35S promoter and GFP gene. After electroporation of these constructions into Agrobacterium tumefaciens EHA105, using pCAMBIA1305 vector as a negative control. The infection solution was injected into the epidermis of Nicotiana tabacum leaves, after culturing in the dark for 3 days, the glass slide was made and placed under the FV1200 laser confocal microscope (Olympus Corporation, Tokyo, Japan) to observe the distribution of fusion protein.

Segments of the 39 DAF fresh pear fruit were fixed in in situ hybridization fixed solution overnight. Pear pulp segments fixed: Take out the segment, wash, clean, put in the fixed fluid (DEPC) immediately to fix above 12 h. Dehydration: The segment is dehydrated by gradient alcohol, paraffin, embedding, and vacuum pumping in the dehydration process. Section: The paraffin is sliced through the slicer, the piece of the slicing machine, and the 62–degree oven roast for 2 h. Dewaxing and dehydration: Soak sections in 2 changes of xylene, 15 min each. Dehydrate in 2 changes of pure ethanol for 5 min each. Then, followed respectively by dehydrating in gradient ethanol of 85 and 75% ethanol 5 min each. Wash in DEPC dilution. Digestion: Mark the objective tissue with a liquid blocker pen, according to the characteristics of tissues, add proteinase K (20 ug/ml) working solution to cover objectives and incubate at 37°C for 22 min. Washing with pure water, then wash three times with PBS (pH 7.4) in a rocker device, 5 min each. Hybridization: Discard the pre-hybridization solution, add the probe hybridization solution, concentration, and incubate the section in a humidity chamber and hybridize overnight. Add the Alkaline Phosphatase-conjugated IgG Fraction Monoclonal Mouse Anti-Digoxin Antibody: (anti-DIG-AP): Remove the blocking solution and add anti-DIG-HRP. Incubate at 37°C for 40 min, then wash sections in TBS four times for 5 min each. NBT/BCIP developing: dry sections slightly, add fresh prepared NBT/BCIP (Thermo Scientific 1-Step NBT/BCIP) chromogenic reagent to marked tissue. Manage reaction time by observing in microscopy until positive expression appears shows blue-purple. Then stop developing reaction by a wash in running tap water. The probe sequence was unique to the PbBGLU1 and PbBGLU16 locus (Supplementary Table 4) and resulted in a single hit when used as a quarry to BLAST the Chinese white pear genome. In situ hybridization probes were synthesized by Sangon Bio Shang hai and labeled with Digoxigenin.

The RT-PCR primers were designed using Primer Premier 5.0 software (Supplementary Table 4). The full-length coding sequence of PbBGLU1 and PbBGLU16 were cloned and then expressed in pGEX4T-1 vector (GE Healthcare, Chicago, IL, United States) between the SamI and NotI (NEB) sites. Two recombinant pGEX4T-1-PbBGLUs were transformed into Escherichia coli BL21 (DE3) Competent Cells. Add Escherichia coli BL21 (DE3) cells carrying pGEX4T-1-PbBGLU1 and pGEX4T-1-PbBGLU16 plasmid to 100 ml of LB liquid containing ampicillin for shaking. Expand the culture to an OD value of about 0.6, add isopropyl-β-D-thioacetamide IPTG to a final concentration of 1 mM/L to induce the fusion expression of the recombinant protein, culture with shaking at 16°C for 24 h, and collect the cells by centrifugation at 4°C. The centrifuged cells were suspended in PBS (pH 7.2–7.4) buffer, using an ultrasonic cell disruptor for disruption, then centrifuge to take the supernatant and discard the pellet. Using affinity chromatography resin with GST tag GST⋅ Bind™ resin (Novagen) to Purify protein. The Elution Buffer (pH 8.0) for GST-Sefinose (TM) Resin (containing reducing glutathione reagent at a concentration of 5 mM) Eluent was used for elution.

To analyze the activity of recombinant protein GST-PbBGLU1 and GST-PbBGLU16, 10 μg of purified protein was incubated at 35°C with 20 μL buffer (50 mM MgSO4, 200 mM KCL, 100 mM PBS pH 7.2–7.4) and 1 mM substrates (coniferin and syringin). The water was added to a final volume of 50 μL. After 1 h of reaction, add 50 μL methanol termination reaction. All reactions were supplemented with 0.1% (V/V) β- Mercaptoethanol. Then the reaction results were analyzed by HPLC.

To calculate the K_m, V_max, and k_cat value of the recombinant protein, under the same original reaction conditions, the substrates (coniferin and syringin) were set to 10 concentrations (0.025, 0.05, 0.075, 0.1, 0.15, 0.2, 0.3, 0.4, 0.5, and 0.6 mM). Then the product concentration C is calculated by an external standard method, and then the reaction rate V is further calculated. Finally, use the double reciprocal method to plot to calculate K_m and V_max.

The RT-PCR primers were designed using Primer Premier 5.0 software (Supplementary Table 4). The full-length coding sequence of PbBGLU1 and PbBGLU16 were cloned from “Dangshan Su” pear first-strand complementary DNA (cDNA) with PrimeSTAR^® Max DNA Polymerase (Takara, Japan). The fragment was cloned into pCAMBIA1301-35S binary vectors between the BamHI and XbaI (NEB) sites and used in plant transformation. Wild type and Arabidopsis At1g61810 BGLU45-2 mutant were used for transformation with A. tumefaciens GV3101 carrying the above binary plasmid using the floral dip method.

The full-length coding fragment was cloned into pCAMBIA1301-35S binary vectors and used in the transient injection of pear fruit overexpression. Hairpin constructs were based on the pRNAiDE001 vector (Heneghan et al., 2007; Eastwood et al., 2008) where self-complementary sense (PbBGLU1: 905-1220bP; PbBGLU16: 832-1183bP) and antisense sequences are separated by a non-functional sequence (loop), then construct this fragment into the pCAMBIA1301 vector with the 35S promoter (Supplementary Figures 1, 2). Pear trees (40-year-old) and 39 days after flowering young fruits were selected for transient injection. The constructed pCAMBIA1301-35s binary vectors of PbBGLUs overexpression and PbBGLUs RNA interference gene were injected on the left of the fruit, and the empty pCAMBIA1301-35s binary vectors were injected on the right for control. The method followed a previously reported protocol (Zhang et al., 2021) nine fruits were injected with each construct in an experiment that was repeated three times independently. After 10 days, the transient injected fruits were picked up.

We collected the Arabidopsis plants for approximately 6 weeks, removed the leaves, and dried them in an oven at 65°C for 48 h. The lignin content of the Arabidopsis plants was estimated following the method of Anderson et al. (Anderson et al., 2015).

We obtained the inflorescence stem segments (young stem regions) of the transgenic and wild-type Arabidopsis plants for approximately 6 weeks. The samples were placed in a solution containing 95% methanol, 70% (v/v) ethanol and glacial acetic acid for 12 h and embedded in paraffin for sectioning with the pathology slicer (RM 2018). Plant tissue sections were stained following standard phloroglucinol staining protocols (Pradhan and Loqué, 2014).

After the transiently injected pear fruits was brought back to the laboratory, samples 5 g pulp from 2.0 mm under the peel to 0.5 mm outside the core, was collected and frozen at –80°C for 24 h, homogenized with a high-speed homogenizer for 3 min, rotating speed was 20,000 RPM/min, added water and stood for a while, then poured out the upper suspension, repeated several times until the upper layer was clear, then filtered and dried stone cells were weighed, repeated for 3 times, stone cell content = measured stone cell dry weight/5 × 100%. The lignin content was measured using the Klason method (Raiskila et al., 2007). A small amount (0.2 g) of stone cells was extracted with 15 ml of 72% H₂SO4 at 30°C for 1 h, combined with 115 ml of distilled water, and boiled for 1 h. The volume was kept constant during boiling. The liquid mixture was filtered and the residue was rinsed with 500 ml of hot water, air-dried and weighed. The lignin content was shown as a percentage (calculated lignin content/dry weight of stone cells × 100%) (Zhang et al., 2021).

HPLC (Thermo Scientific UltiMate 3000 HPLC) analysis was carried out using a Columbus (Thermo Scientific 5-μm C18 column; 250×4.60 mm). Acetonitrile (solvent A) and H₂O (solvent B) with gradient of 10–30% acetonitrile in water (all solutions contained 0.1% trifluoroacetic acid) with a flow rate at 1 ml min^–1 over 35 min was used. Each peak on the chromatogram was scanned between 200 and 400 nm (photodiode array profile) and was integrated at 264, 280, and 324 nm. The data were acquired and analyzed using the software ChromQuest version. Extraction of tissue was carried out as described in Lanot et al. (2006).

The BGLU family conserved domain (Pfam: PF00232) was used as a target query sequence, and the corresponding hidden Markov model was obtained from the Pfam website (see text footnote 2) which was used to identify the BGLU members. SMART (see text footnote 3) was used to check whether there are characteristic domains, with the redundant and repeated sequences ultimately removed. As a result, we identified 50 BGLU genes from Chinese white pear (P. bretschneideri). The detailed information on the gene ID, gene name, chromosomal location, protein structure, and characteristics of the corresponding BGLU proteins are listed in Supplementary Table 1.

Inclusive of all PbBGLU genes, a phylogenetic tree was constructed with MEGA 7.0 software using the neighbor-joining (NJ) method (bootstraps = 1,000), which divided them into 8 groups (I, II, III, IV, V-A, V-B, V-C, V-D) with supported bootstrap values (Figure 1). We used MEME online software to analyze the amino acid sequence of PbBGLU genes and identified 20 motifs among all the sequences (Supplementary Table 2 and Figure 1). All the PbBGLUs contain motif1 (CTGGBSATEPYJVAHHQLLAHAAAVKLYREKYQA), and only two genes do not contain motif9 (WFEPASES KEDKAAALRALDF), suggesting that motif1 and motif9 are conserved gene motifs of PbBGLUs family. Only group I, group II, and group III genes contain motif19, whereas the other groups do not, and no members of group V-D contain motif2, it is speculated that these motifs may have special significance. Previous studies implied that gene structural diversity can lead to the evolution of multi-gene families (Dong et al., 2019). To better characterize and understand the structural diversity of the PbBGLUs, gene exon-intron analysis was carried out (Figure 1). Notably, PbBGLU23 and PbBGLU24 have only one exon, whereas the others contain 7 or more exons; PbBGLU34 contains the most exons (21), and most genes had between 10 and 15 exons. These results indicate the possible occurrence of alternative splicing during gene evolution, which in turn leads to functional differences between these closely related PbBGLUs (Cao et al., 2017; Dong et al., 2019).

To further understand the evolutionary relationship and related functions of PbBGLUs, a phylogenetic tree comprising BGLUs from A. thaliana (Xu et al., 2004), O. sative (Opassiri et al., 2006), poplar (Tsuyama and Takabe, 2015), and P. bretschneideri was constructed (Figure 2). We marked the position of the eight groups in pear on this phylogenetic tree and found that group I, group II, and group IV contained only PbBGLUs but not AtBGLUs, Osbglus, and PtrBGLUs, indicating that some changes occurred among the BGLUs of different species during the evolutionary process (Cao et al., 2016a,b). Arabidopsis AtBGLU17 (Roepke et al., 2017) was related to the use of Arabidopsis flavonoids which was clustered together with PbBGLU19 (group III), and it is speculated that PbBGLU19 might be involved in the hydrolysis of flavonoid glycosides in pear. PbBGLU group V-D clustered together with AtBGLU42, and which encodes an MYB72-dependent key regulator of rhizobacterium-induced systemic resistance and modulates iron deficiency responses in Arabidopsis roots (Zamioudis et al., 2014). AtBGLU1-6, which plays a role in the accumulation of flavonoids in Arabidopsis (Ishihara et al., 2016), clustered together with PbBGLU group V-A genes, speculating that PbBGLU29 and PbBGLU30 may have similar functions in pear. The functions of AtBGLU45 and AtBGLU46 in Arabidopsis (Escamilla-Treviño et al., 2006; Chapelle et al., 2012); Os4BGlu14, Os4BGlu16, and Os4BGlu18 in O. sative (Baiya et al., 2014; Baiya et al., 2018); PtrBGLU6 in poplar (Tsuyama and Takabe, 2015) have been clarified, and their products catalyze the hydrolysis of monolignol glucosides during lignification. The three genes (PbBGLU1, PbBGLU15, and PbBGLU16) in PbBGLU group V-C clustered together with these functional genes, speculating that PbBGLU1, PbBGLU15, and PbBGLU16 may be able to affect lignification by catalyzing the hydrolysis of monolignol glucosides in pear fruit, which in turn affects the development of pear stone cells.

We constructed a phylogenetic tree comprising BGLUs genes from A. thaliana (Xu et al., 2004), O. sative (Opassiri et al., 2006), poplar (Tsuyama and Takabe, 2015), and P. bretschneideri, and found that PbBGLU1, PbBGLU15, and PbBGLU16 in the group V-C PbBGLUs of pear clustered together with AtBGLU45, AtBGLU46, Os4BGlu14, Os4BGlu16, Os4BGlu18, and PtrBGLU6 functional genes, which can hydrolyze monolignol glucosides during lignification (Escamilla-Treviño et al., 2006; Chapelle et al., 2012; Baiya et al., 2014; Tsuyama and Takabe, 2015; Baiya et al., 2018). To verify whether PbBGLU1, PbBGLU15, and PbBGLU16 have similar functions, we selected all PbBGLU genes in group V (a total of 11 PbBGLUs) for qRT-PCR analysis. The expression profiles of these genes in different tissues of pear (leaves, buds, stems, flowers) and pear fruit at different developmental stages (15, 39, 47, 55, 63, 79, and 145 days after flowering) were also investigated (Figure 3). It found that PbBGLU43 was not expressed in different tissue parts or at different developmental stages of the fruit. The expression levels of the other ten PbBGLUs were generally high in leaves, buds, flowers, and young fruits (15, 39, 47, and 55 days after flowering). This result was consistent with those of previous research showing that β-glucosidase can be found mainly in young vegetative parts (Kristoffersen et al., 2000; Biely et al., 2003). The expression levels of the PbBGLU1, PbBGLU15, and PbBGLU16 genes were relatively high in the leaves and flowers as compared to buds and stems. All three genes were expressed at different stages of fruit development, and their expression showed a trend of first increasing then decreasing, with a similar trend of the content for lignin and stone cells in “Dangshan Su” pear fruit (Cai et al., 2010; Jin et al., 2013). The expression patterns of these three genes were also similar to those of the key genes involved in the regulation of the lignin synthesis pathway (Xie et al., 2013). These results suggested that the three genes (PbBGLU1, PbBGLU15, and PbBGLU16) in the group V-C PbBGLUs might be functional genes involved in the hydrolysis of monolignol glucosides in the pear lignin synthesis pathway.

As early as 1995, it found that Pinus contorta β-glucosidase activity occurred in the cell wall of secondary xylem tissue (Dharmawardhana et al., 1995). AtBGLU45 and AtBGLU46 of A. thaliana were also expressed in the cell wall (Chapelle et al., 2012). The rice Os4BGlu16 tagged with a C-terminal GFP was transiently expressed in tobacco leaf epithelial cells, which revealed that the GFP signal was exclusively localized to the plasma membrane and extracellular space in both intact and plasmolyzed cells (Baiya et al., 2018). It is consistent with apoplastic or cell wall localization. In this study, to explore the subcellular localization of PbBGLU1, PbBGLU15, and PbBGLU16, PbBGLU-GFP expression vectors were constructed and transformed into N. tabacum. As shown in Figure 4, green fluorescence signals from the expressed PbBGLU1-GFP, PbBGLU15-GFP, and PbBGLU16-GFP fusion constructs were specifically distributed in the extracellular region.

To confirm the PbBGLU1 and PbBGLU16 gene expression associated with those stone cells undergoing lignin deposition and secondary wall thickening, we carried out RNA in situ hybridization (Wu and Wagner, 2012) on 39 DAF (days after flowering) pear fruit. The results of PbBGLU1 transcripts in situ hybridization with PbBGLU1 anti-sense probe showed that PbBGLU1 transcripts were not only localized to some pulp cell walls, lignin deposition and stone cell areas of a pear fruit, there was also a small amount of enrichment in normal pear flesh cells (Figure 5A). This result shows that PbBGLU1 is not only involved in lignin biosynthesis but also may be involved in other metabolic activities during pear fruit development. However, PbBGLU16 transcripts were only enriched in lignin deposition and stone cell areas of pear fruit (Figure 5C), indicating that it may be a specific gene for lignin synthesis in pear fruit. Hybridization with control PbBGLU1 and PbBGLU16 sense probe showed no labeling as expected (Figures 5B,D).

PbBGLU1 and PbBGLU16 were cloned into pGEX4T-1 vector, and then the constructed vector was induced and expressed in BL21 (ED3) E. coli. A target protein GST-PbBGLU1 with a weight of 78.07 kDa and another target protein GST-PbBGLU16 with a weight of 87.23 kDa obtained (GST: 26kDa; PbBGLU1: 52.07 kDa; PbBGLU16: 61.23 kDa). The sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) was used to analysis and verification (Supplementary Figure 3). To verify whether the recombinant protein has the ability to glycosylate the monolignol glucosides coniferin, syringin. We added recombinant protein, coniferin, and syringin and buffer in the certain reaction system. After the reaction for a while, methanol was added to terminate the reaction and then put into the HPLC for detection. The standard of coniferyl alcohol, sinapyl alcohol, coniferin and syringin were put into the HPLC for reference (Figures 6A,B). The results showed that after the recombinant protein was added, two peaks of the substrate (coniferin and syringin) and product (coniferyl alcohol and sinapyl alcohol) were detected by HPLC. This phenomenon indicates that the recombinant protein of GST-PbBGLU1 and GST-PbBGLU16 were active against both substrates (coniferin and syringin). To calculate the kinetic parameters of the recombinant protein for coniferin and syringin enzyme, 10 substrates with different concentrations were set in this study. The K_m values, V_max and k_cat of GST-PbBGLU1 and GST-PbBGLU16 were calculated through the double reciprocal method. The results are shown in Supplementary Figure 4 and Table 1. Compared with syringin, GST-PbBGLU1 and GST-PbBGLU16 had a stronger activity and faster catalytic efficiency for coniferin. In addition, compared with GST-PbBGLU1, GST-PbBGLU16 shows higher catalytic activity and catalytic efficiency for the two substrates.

The PbBGLU1 and PbBGLU16 were transferred into Arabidopsis BGLU-45 mutants (PbBGLU1-RE and PbBGLU16-RE) and wild-type Arabidopsis (PbBGLU1-OE and PbBGLU16-OE), reverse transcription PCR validation of Arabidopsis is shown in Figures 7A,B. To directly observe the in situ distribution of lignin in inflorescence stems of these Arabidopsis plants, we performed the Wiesner (phloroglucinol-HCl) histochemical staining (Pradhan and Loqué, 2014) of inflorescence stems of Arabidopsis plants. The Wiesner staining results indicated that the xylem and intravascular fibers in the stems of Arabidopsis BGLU-45 mutants showed no increase to any extent compared with wild-type Arabidopsis (Figures 7C-2,D-2), this result is consistent with the results of previous studies (Chapelle et al., 2012). Then, we observed the stained sections of PbBGLU1-OE and PbBGLU16-OE Arabidopsis and the Wiesner staining results indicated that the xylem and intravascular fibers in the stems of Arabidopsis PbBGLU1-OE and PbBGLU16-OE showed stronger phloroglucinol staining compared to the wild-type plants (Figures 7C-4,D-4). According to the staining results, we can also find that PbBGLU16-OE xylem and intravascular fibers Arabidopsis increase more strongly than PbBGLU1-OE Arabidopsis. Subsequently, we measured the lignin content of wild-type, Arabidopsis BGLU-45 mutant, PbBGLU1-RE, and PbBGLU16-RE, PbBGLU1-OE, and PbBGLU16-OE Arabidopsis by acetyl bromide method (Anderson et al., 2015; Figures 7E,F). It found that the lignin content of Arabidopsis BGLU-45 mutant, PbBGLU1-RE and PbBGLU16-RE was no change than that of wild-type. However, compared with wild-type Arabidopsis the overexpression plants lignin increased in varying degrees, the effect of PbBGLU16 on the increase of lignin in Arabidopsis is greater than that of PbBGLU1.

In order to further explore how BGLU-45 mutants and overexpression of PbBGLU1 and PbBGLU16 affect the lignin of A. thaliana, we measured the A. thaliana contents of coniferin and syringin by HPLC. The result showed that coniferin and syringin increased significantly in BGLU-45 mutant compared with the wild type (Table 2), but the lignin did not increase. This phenomenon occurred probably because coniferin and syringin are not precursors of lignin synthesis and usually exist in vacuoles as transport or storage forms. This is also probably the reason why there is no significant change in Arabidopsis lignin after knocking out BGLU45 in previous research (Chapelle et al., 2012; Baiya et al., 2018). Then, the contents of coniferin and syringin rising phenomenon have been restored through transferred PbBGLU1 and PbBGLU16 into Arabidopsis BGLU-45 mutant (Table 2). However, the lignin monolignol glycosides of transgenic Arabidopsis after overexpression of PbBGLU1 and PbBGLU16 decreased to some extent, but the decrease was not obvious, it may be because the decrease of lignin monolignol glycosides will increase the expression of UGT.

Since the stable genetic transformation system of pears is in its infancy, it is difficult to perform a functional analysis of the PbBGLU1 and PbBGLU16 in pears. Therefore, we chose the method of transient expression by injection for analysis. We injected transiently the PbBGLU1, PbBGLU16, and PbBGLU1 and PbBGLU16 RNA interference in the 39 DAF (days after flowering) fruit of the pear tree. By Wiesner (phloroglucinol-HCl) histochemical staining (Pradhan and Loqué, 2014), the results show (Figures 8B,D) that when PbBGLU1 and PbBGLU16 were silenced by RNA interference, the change of pear lignin and stone cells was not obvious compared with the empty vector injection. But after PbBGLU1 and PbBGLU16 gene were overexpressed in the 39 DAF pear fruit, it found that the lignin and stone cells increased compared with the empty vector injection control in the pear fruit (Figures 8A,C). Moreover, compared with the transient expression of PbBGLU1, PbBGLU16 had a stronger effect on the lignin and expression cells of pear fruit. To explore the specific effects of overexpression and RNAi of PbBGLU1 and PbBGLU16 on stone cells and lignin of pear fruit, we measured the contents of stone cells and lignin in pear fruit after empty vector, overexpression, and RNAi injection. The results are shown in Figures 8B,D which indicate that when PbBGLU1 and PbBGLU16 were silenced, the contents of stone cells and lignin in pear fruit did not change compared with empty vector injection. However, in pear fruits with transient overexpression of PbBGLU1, the contents of lignin and stone cells in pear fruits were significantly higher (0.01 < P < 0.05) than those in pear fruits with empty vector injection. After transient expression of PbBGLU16, lignin in pear fruit increased significantly compared with the control group (0.01 < P < 0.05), while the content of stone cells in pear fruit injected with PbBGLU16 showed a very significant difference (P < 0.01) compared with the control group.