Family-1 UDP glycosyltransferases in pear (Pyrus bretschneideri): Molecular identification, phylogenomic characterization and expression profiling during stone cell formation

Identification and Functional Verification of the Glycosyltransferase Gene Family Involved in Flavonoid Synthesis in Rubus chingii Hu

Glycosylation is catalyzed by UDP-glycosyltransferase (UGT) and plays an important role in enriching the diversity of flavonoids. Rubus plants contain a lot of natural flavonoid glycosides, which are important plants with a homology of medicine and food. However, information about the Rubus UGT gene family is very limited. In this study, we carried out genome-wide analysis and identified the 172, 121, 130, 121 UGT genes in R. chingii, R. corchorifolius, R. idaeus, and R. occidentalis, respectively, and divided them into 18 groups. The analysis of the protein motif and gene structure showed that there were structural and functional conservations in the same group, but there were differences among different groups. Gene replication analysis showed that raspberry and dicotyledons had a higher homology. The expansion of the UGTs gene family was mainly driven by tandem replication events, and experienced purified selection during the long evolution of the raspberry. Cis-acting element analysis showed that they were related to plant growth and development, hormone regulation, and stress response. In addition, according to a comprehensive analysis of the co-expression network constructed by transcriptome data and phylogenetic homology, RchUGT169 was identified as a flavonoid glucosyltransferase. Through the transient expression in tobacco, it was verified that RchUGT169 could catalyze the conversion of kaempferol and quercetin to the corresponding flavonoid glycosides. In conclusion, this research enriched the understanding of the diversity of UGTs in Rubus and determined that RcUGT169 can catalyze flavonoids.

Genome-Wide Analysis of UGT Genes in Petunia and Identification of PhUGT51 Involved in the Regulation of Salt Resistance

UDP-glycosyltransferase (UGT) plays an essential role in regulating the synthesis of hormones and secondary metabolites in plants. In this study, 129 members of the Petunia UGT family were identified and classified into 16 groups (A-P) based on phylogenetic analysis. The same subgroups have conserved motif compositions and intron/exon arrangement. In the promoters of the Petunia UGT genes, several cis-elements associated with plant hormones, growth and development, and abiotic stress have been discovered. Their expression profiles in five tissues were revealed by tissue expression based on RNA-seq data. Subcellular localization analysis showed that PhUGT51 was located in the nucleus and cell membrane. Salt stress caused an increase in the expression level of PhUGT51, but the expression level remained stable with the growth over time. In addition, the overexpression of PhUGT51 caused a significant increase in salt resistance. Our study systematically analyses the UGT gene family in Petunia for the first time and provides some valuable clues for the further functional studies of UGT genes.

PbUGT72AJ2-Mediated Glycosylation Plays an Important Role in Lignin Formation and Stone Cell Development in Pears (Pyrus bretschneideri)

Glycosylation is necessary for many processes of plant secondary metabolism. It can maintain plant homeostasis and is of great significance to normal plant growth and development. At present, the significance of glycosylation for lignin biosynthesis has been proven in some plants, but it has not yet been reported in pears. We used in situ hybridization, in vitro expression, substrate catalysis, transgenic Arabidopsisthaliana, and transient transformation of pear fruit in our investigation, which was predicated on the identification of a gene PbUGT72AJ2 that may be involved in lignin monolignol glycosylation according to our previous work. These results revealed that PbUGT72AJ2 transcripts were localized to some pulp cell walls, lignin deposition, and stone cell areas of pear fruit. The recombinant PbUGT72AJ2-pGEX4T-1 protein had activity against coniferyl alcohol and sinapyl alcohol, and its catalytic efficiency against coniferyl alcohol was higher than that against sinapyl alcohol. When PbUGT72AJ2 was transferred into Arabidopsisthaliana mutants, it was found that some characteristics of Arabidopsisthalianaugt72e3 mutants were restored. In Arabidopsisthaliana, overexpression of PbUGT72AJ2 enhanced the contents of coniferin and syringin, whereas lignification did not change significantly. Transient transformation of pear fruit showed that when PbUGT72AJ2 in pear fruit was silenced by RNA interference, the content of lignin and stone cells in pear fruit increased, whereas the gene PbUGT72AJ2 was overexpressed in pear fruit, and there was almost no change in the pear fruit compared with the control. Lignin deposition in pear fruit was closely related to stone cell development. In this study, we proved that PbUGT72AJ2 plays an important role in lignin deposition and stone cell development in pear fruit, which provides a molecular biological basis for improving pear fruit quality at the molecular level.

Evolutionary and Integrative Analysis of Gibberellin-Dioxygenase Gene Family and Their Expression Profile in Three Rosaceae Genomes (F. vesca, P. mume, and P. avium) Under Phytohormone Stress

The gibberellin-dioxygenase (GAox) gene family plays a crucial role in regulating plant growth and development. GAoxs, which are encoded by many gene subfamilies, are extremely critical in regulating bioactive GA levels by catalyzing the subsequent stages in the biosynthesis process. Moreover, GAoxs are important enzymes in the GA synthesis pathway, and the GAox gene family has not yet been identified in Rosaceae species (Prunus avium L., F. vesca, and P. mume), especially in response to gibberellin and PCa (prohexadione calcium; reduce biologically active GAs). In the current investigation, 399 GAox members were identified in sweet cherry, Japanese apricot, and strawberry. Moreover, they were further classified into six (A-F) subgroups based on phylogeny. According to motif analysis and gene structure, the majority of the PavGAox genes have a remarkably well-maintained exon-intron and motif arrangement within the same subgroup, which may lead to functional divergence. In the systematic investigation, PavGAox genes have several duplication events, but segmental duplication occurs frequently. A calculative analysis of orthologous gene pairs in Prunus avium L., F. vesca, and P. mume revealed that GAox genes are subjected to purifying selection during the evolutionary process, resulting in functional divergence. The analysis of cis-regulatory elements in the upstream region of the 140 PavGAox members suggests a possible relationship between genes and specific functions of hormone response-related elements. Moreover, the PavGAox genes display a variety of tissue expression patterns in diverse tissues, with most of the PavGAox genes displaying tissue-specific expression patterns. Furthermore, most of the PavGAox genes express significant expression in buds under phytohormonal stresses. Phytohormones stress analysis demonstrated that some of PavGAox genes are responsible for maintaining the GA level in plant-like Pav co4017001.1 g010.1.br, Pav sc0000024.1 g340.1.br, and Pav sc0000024.1 g270.1.mk. The subcellular localization of PavGAox protein utilizing a tobacco transient transformation system into the tobacco epidermal cells predicted that GFP signals were mostly found in the cytoplasm. These findings will contribute to a better understanding of the GAox gene family's interaction with prohexadione calcium and GA, as well as provide a strong framework for future functional characterization of GAox genes in sweet cherry.

Analysis of the β-Glucosidase Family Reveals Genes Involved in the Lignification of Stone Cells in Chinese White Pear (Pyrus bretschneideri Rehd.)

BGLU β-glucosidases in glycoside hydrolase family 1 (GH1) are involved in many processes of plant secondary metabolism. In particular, its de-glycosylation function plays an important role in the transport of lignin monolignols. No comprehensive study of the BGLU family in Chinese pear (Pyrus bretschneideri Rehd.) has been reported yet. In this study, the 50 BGLU family members from Chinese white pear were identified. Three candidate genes, PbBGLU1, PbBGLU15, and PbBGLU16, that may be involved in lignin synthesis were screened by bioinformatics analysis and qRT-PCR. Subcellular localization showed that all three of these candidate genes were expressed in the extracellular region. Then, we analyzed the functions of PbBGLU1 and PbBGLU16. In situ hybridization analysis showed that PbBGLU1 transcripts were not only localized to some pulp cell walls, lignin deposition, and stone cell areas of a pear fruit, but that was also a small amount of enrichment in normal pear flesh cells. PbBGLU16 transcripts were only enriched in lignin deposition and stone cell areas of pear fruit. Enzyme activity analysis revealed that GST-PbBGLU1 and GST-PbBGLU16 had a stronger activity and higher catalytic efficiency for coniferin than syringin. In addition, GST-PbBGLU16 exhibited the higher activity and catalytic efficiency for the two substrates compared with GST-PbBGLU1. The transformation of PbBGLU1 and PbBGLU16 into Arabidopsis identified that the lignin contents of Arabidopsis BGLU-45 mutant, PbBGLU1-RE, and PbBGLU16-RE were not changed than that of wild-type. However, compared with wild-type Arabidopsis, the overexpression of the plant's lignin increased in varying degrees. The effect of PbBGLU16 on the lignin increment was greater than that of PbBGLU1 in Arabidopsis. In pear fruits, with transient overexpression of PbBGLU1, the contents of lignin and stone cells were significantly higher (0.01 < P < 0.05) than those with empty vector injection pear fruits. After transient expression of PbBGLU16, lignin in pear fruit increased significantly (0.01 < P < 0.05) and stone cells showed a very significant difference (P < 0.01) compared with the control group. However, RNA interference silenced these two genes in pear fruit, which seemed to have no impression on lignin and stone cells. This study provides a molecular biological basis for improving pear fruit quality at the molecular level.

UGT72, a Major Glycosyltransferase Family for Flavonoid and Monolignol Homeostasis in Plants

Simple Summary

Phenylpropanoids are specialized metabolites playing crucial roles in plant developmental processes and in plant defense towards pathogens. The attachment of sugar moieties to these small hydrophobic molecules renders them more hydrophilic and increases their solubility. The UDP-glycosyltransferase 72 family (UGT72) of plants has been shown to glycosylate mainly two classes of phenylpropanoids, (i) the monolignols that are the building blocks of lignin, the second most abundant polymer after cellulose, and (ii) the flavonoids, which play determinant roles in plant interactions with other organisms and in response to stress. The purpose of this review is to bring an overview of the current knowledge of the UGT72 family and to highlight its role in the homeostasis of these molecules. Potential applications in pharmacology and in wood, paper pulp, and bioethanol production are given within the perspectives.

Abstract

Plants have developed the capacity to produce a diversified range of specialized metabolites. The glycosylation of those metabolites potentially decreases their toxicity while increasing their stability and their solubility, modifying their transport and their storage. The UGT, forming the largest glycosyltransferase superfamily in plants, combine enzymes that glycosylate mainly hormones and phenylpropanoids by using UDP-sugar as a sugar donor. Particularly, members of the UGT72 family have been shown to glycosylate the monolignols and the flavonoids, thereby being involved in their homeostasis. First, we explore primitive UGTs in algae and liverworts that are related to the angiosperm UGT72 family and their role in flavonoid homeostasis. Second, we describe the role of several UGT72s glycosylating monolignols, some of which have been associated with lignification. In addition, the role of other UGT72 members that glycosylate flavonoids and are involved in the development and/or stress response is depicted. Finally, the importance to explore the subcellular localization of UGTs to study their roles in planta is discussed.

Genomic-Wide Identification and Characterization of the Uridine Diphosphate Glycosyltransferase Family in Eucommia ulmoides Oliver

Eucommia ulmoides Oliver is a woody plant with great economic and medicinal value. Its dried bark has a long history of use as a traditional medicinal material in East Asia, which led to many glycosides, such as aucubin, geniposide, hyperoside, astragalin, and pinoresinol diglucoside, being recognized as pharmacologically active ingredients. Uridine diphosphate glycosyltransferases (UGTs) catalyze a glycosyl-transferring reaction from the donor molecule uridine-5'-diphosphate-glucose (UDPG) to the substrate, which plays an important role in many biological processes, such as plant growth and development, secondary metabolism, and environmental adaptation. In order to explore the biosynthetic pathways of glycosides in E. ulmoides, 91 putative EuUGT genes were identified throughout the complete genome of E. ulmoides through function annotation and an UDPGT domain search. Phylogenetic analysis categorized them into 14 groups. We also performed GO annotations on all the EuUGTs to gain insights into their functions in E. ulmoides. In addition, transcriptomic analysis indicated that most EuUGTs showed different expression patterns across diverse organs and various growing seasons. By protein-protein interaction predication, a biosynthetic routine of flavonoids and their glycosides was also proposed. Undoubtedly, these results will help in future research into the biosynthetic pathways of glycoside compounds in E. ulmoides.

Genome-wide identification and characterization of bZIP transcription factors and their expression profile under abiotic stresses in Chinese pear (Pyrus bretschneideri)

BACKGROUND

In plants, basic leucine zipper transcription factors (TFs) play important roles in multiple biological processes such as anthesis, fruit growth & development and stress responses. However, systematic investigation and characterization of bZIP-TFs remain unclear in Chinese white pear. Chinese white pear is a fruit crop that has important nutritional and medicinal values.

RESULTS

In this study, 62 bZIP genes were comprehensively identified from Chinese Pear, and 54 genes were distributed among 17 chromosomes. Frequent whole-genome duplication (WGD) and dispersed duplication (DSD) were the major driving forces underlying the bZIP gene family in Chinese white pear. bZIP-TFs are classified into 13 subfamilies according to the phylogenetic tree. Subsequently, purifying selection plays an important role in the evolution process of PbbZIPs. Synteny analysis of bZIP genes revealed that 196 orthologous gene pairs were identified between Pyrus bretschneideri, Fragaria vesca, Prunus mume, and Prunus persica. Moreover, cis-elements that respond to various stresses and hormones were found on the promoter regions of PbbZIP, which were induced by stimuli. Gene structure (intron/exon) and different compositions of motifs revealed that functional divergence among subfamilies. Expression pattern of PbbZIP genes differential expressed under hormonal treatment abscisic acid, salicylic acid, and methyl jasmonate  in pear fruits by real-time qRT-PCR.

CONCLUSIONS

Collectively, a systematic analysis of gene structure, motif composition, subcellular localization, synteny analysis, and calculation of synonymous (Ks) and non-synonymous (Ka) was performed in Chinese white pear. Sixty-two bZIP-TFs in Chinese pear were identified, and their expression profiles were comprehensively analyzed under ABA, SA, and MeJa hormones, which respond to multiple abiotic stresses and fruit growth and development. PbbZIP gene occurred through Whole-genome duplication and dispersed duplication events. These results provide a basic framework for further elucidating the biological function characterizations under multiple developmental stages and abiotic stress responses.

Identification and Characterization of Glucosyltransferase That Forms 1-Galloyl-β-d-Glucogallin in Canarium album L., a Functional Fruit Rich in Hydrolysable Tannins

Hydrolysable tannins (HTs) are useful secondary metabolites that are responsible for pharmacological activities and astringent taste, flavor, and quality in fruits. They are also the main polyphenols in Canarium album L. (Chinese olive) fruit, an interesting and functional fruit that has been cultivated for over 2000 years. The HT content of C. album fruit was 2.3-13 times higher than that of berries with a higher content of HT. 1-galloyl-β-d-glucose (βG) is the first intermediate and the key metabolite in the HT biosynthesis pathway. It is catalyzed by UDP-glucosyltransferases (UGTs), which are responsible for the glycosylation of gallic acid (GA) to form βG. Here, we first reported 140 UGTs in C. album. Phylogenetic analysis clustered them into 14 phylogenetic groups (A, B, D-M, P, and Q), which are different from the 14 typical major groups (A~N) of Arabidopsis thaliana. Expression pattern and correlation analysis showed that UGT84A77 (Isoform0117852) was highly expressed and had a positive correlation with GA and βG content. Prokaryotic expression showed that UGT84A77 could catalyze GA to form βG. These results provide a theoretical basis on UGTs in C. album, which will be helpful for further functional research and availability on HTs and polyphenols.

Tandem UGT71B5s Catalyze Lignan Glycosylation in Isatis indigotica With Substrates Promiscuity

Lignans are a class of chemicals formed by the combination of two molecules of phenylpropanoids with promising nutritional and pharmacological activities. Lignans glucosides, which are converted from aglycones catalyzed by uridine diphosphate (UDP) glycosyltransferases (UGTs), have abundant bioactivities. In the present study, two UGTs from Isatis indigotica Fort., namely IiUGT71B5a and IiUGT71B5b, were characterized to catalyze the glycosylation of lignans with promiscuities toward various sugar acceptors and sugar donors, and pinoresinol was the preferred substrate. IiUGT71B5a was capable of efficiently producing both pinoresinol monoglycoside and diglycoside. However, IiUGT71B5b only produced monoglycoside, and exhibited considerably lower activity than IiUGT71B5a. Substrate screening indicated that ditetrahydrofuran is the essential structural characteristic for sugar acceptors. The transcription of IiUGT71B5s was highly consistent with the spatial distribution of pinoresinol glucosides, suggesting that IiUGT71B5s may play biological roles in the modification of pinoresinol in I. indigotica roots. This study not only provides insights into lignan biosynthesis, but also elucidates the functional diversity of the UGT family.

Genome-wide characterization of the cellulose synthase gene superfamily in Pyrus bretschneideri and reveal its potential role in stone cell formation

Members of the cellulose synthase (CesA) and cellulose synthase-like (Csl) families from the cellulose synthase gene superfamily participate in cellulose and hemicellulose synthesis in the plasma membrane. The members of this superfamily are vital for cell wall construction during plant growth and development. However, little is known about their function in pear fruit, a model for Rosaceae species and for fleshy fruit development. In our research, a total of 36 CesA/Csl family members were identified from the pear and were grouped into six subfamilies (CesA, CslB, CslC, CslD, CslE, and CslG) according to phylogenetic relationships. We performed a protein sequence physicochemical analysis, phylogenetic tree construction, a gene structure, a conserved domain, and chromosomal localization analysis. The results indicated that most of the CesA/Csl genes from pear are closely related to genes in Arabidopsis, but these families have unique characteristics in terms of their gene structure, chromosomal localization, phylogeny, and deduced protein sequences, suggesting that they have evolved through different processes. Tissue expression analysis results showed that most of the CesA/Csl genes were constitutively expressed at different levels in different organs. Furthermore, the expression levels of four genes (Pbr032894.2, Pbr016107.1, Pbr00518.1, and Pbr034218.1) tended to first increase and then decrease during fruit development, implying that these four genes may be involved in the development of stone cells of pear fruit. Our results may help elucidate the evolutionary history and functional differences of the CesA/Csl genes in pear and lay a foundation for further investigation of the CesA/Csl genes in pear and other Rosaceae species.

Functional analysis of a novel C-glycosyltransferase in the orchid Dendrobium catenatum

Flavonoids, which are a diverse class of phytonutrients, are used by organisms to respond to nearly all abiotic stresses and are beneficial for human health. Glycosyltransferase, used during the last step of flavonoid biosynthesis, is important in flavonoid enrichment. However, little is known about glycosyltransferase in the orchid Dendrobium catenatum (D. officinale). In this study, we isolated a novel C-glycosyltransferase (designated DcaCGT) from the orchid D. catenatum by identifying and analyzing 82 putative genes in the GT1 family. DcaCGT could specifically catalyze not only di-C-glycosylation but also O-glycosylation. Apart from the normal function of catalyzing 2-hydroxynaringenin and phloretin to the respective di-C-glycosides, DcaCGT also catalyzes apigenin to cosmosiin. Targeted metabolic profiling of the substrates (2-hydroxynaringenin, phloretin, and apigenin) and products (vitexin, isovitexin, vicenin-2, nothofagin, 3',5'-di-C-glucosylphloretin, and cosmosiin) in different tissues showed that vicenin-2 was the most abundant product of this novel enzyme. Cosmosiin was detected in flowers and flower buds. We also established that DcaCGT functions expanded throughout the evolution of D. catenatum. Residual OGT activity may help D. catenatum resist drought stress. Our study illustrates the function, origin, and differentiation of DcaCGT and provides insights into glycosylation and molecular propagation processes, which can be used to improve the production of flavonoids by the cultivated medicinal plant D. catenatum.

Effects of Metaxenia on Stone Cell Formation in Pear (Pyrus bretschneideri) Based on Transcriptomic Analysis and Functional Characterization of the Lignin-Related Gene PbC4H2

The deposition of lignin in flesh parenchyma cells for pear stone cells, and excessive stone cells reduce the taste and quality of the fruit. The effect of metaxenia on the quality of fruit has been heavily studied, but the effect of metaxenia on stone cell formation has not been fully elucidated to date. This study used P. bretschneideri (Chinese white pear) cv. ‘Yali’ (high-stone cell content) and P. pyrifolia (Sand pear) cv. ‘Cuiguan’ (low-stone cell content) as pollination trees to pollinate P. bretschneideri cv. ‘Lianglizaosu’ separately to fill this gap in the literature. The results of quantitative determination, histochemical staining and electron microscopy indicated that the content of stone cells and lignin in YL fruit (‘Yali’ (pollen parent) × ‘Lianglizaosu’ (seed parent)) was significantly higher than that in CL fruit (‘Cuiguan’ (pollen parent) × ‘Lianglizaosu’ (seed parent)). The transcriptome sequencing results that were obtained from the three developmental stages of the two types of hybrid fruits indicated that a large number of differentially expressed genes (DEGs) related to auxin signal transduction (AUX/IAAs and ARFs), lignin biosynthesis, and lignin metabolism regulation (MYBs, LIMs, and KNOXs) between the CL and YL fruits at the early stage of fruit development. Therefore, metaxenia might change the signal transduction process of auxin in pear fruit, thereby regulating the expression of transcription factors (TFs) related to lignin metabolism, and ultimately affecting lignin deposition and stone cell development. In addition, we performed functional verification of a differentially expressed gene, PbC4H2 (cinnamate 4-hydroxylase). Heterologous expression of PbC4H2 in the c4h mutant not only restored its collapsed cell wall, but also significantly increased the lignin content in the inflorescence stem. The results of our research help to elucidate the metaxenia-mediated regulation of pear stone cell development and clarify the function of PbC4H2 in cell wall development and lignin synthesis, which establishes a foundation for subsequent molecular breeding.

Phylogenomic analysis of UDP-dependent glycosyltransferases provides insights into the evolutionary landscape of glycosylation in plant metabolism

Glycosylated metabolites generated by UDP-dependent glycosyltransferases (UGTs) play critical roles in plant interactions with the environment as well as human and animal nutrition. The evolution of plant UGTs has previously been explored, but with a limited taxon sampling. In this study, 65 fully sequenced plant genomes were analyzed, and stringent criteria for selection of candidate UGTs were applied to ensure a more comprehensive taxon sampling and reliable sequence inclusion. In addition to revealing the overall evolutionary landscape of plant UGTs, the phylogenomic analysis also resolved the phylogenetic association of UGTs from free-sporing plants and gymnosperms, and identified an additional UGT group (group R) in seed plants. Furthermore, lineage-specific expansions and contractions of UGT groups were detected in angiosperms, with the total number of UGTs per genome remaining constant generally. The loss of group Q UGTs in Poales and Brassicales, rather than functional convergence in the group Q containing species, was supported by a gene tree of group Q UGTs sampled from many species, and further corroborated by the absence of group Q homologs on the syntenic chromosomal regions in Arabidopsis thaliana (Brassicales). Branch-site analyses of the group Q UGT gene tree allowed for identification of branches and amino acid sites that experienced episodic positive selection. The positively selected sites are located on the surface of a representative group Q UGT (PgUGT95B2), away from the active site, suggesting their role in protein folding/stability or protein-protein interactions.

In Silico Genome-Wide Analysis of Respiratory Burst Oxidase Homolog (RBOH) Family Genes in Five Fruit-Producing Trees, and Potential Functional Analysis on Lignification of Stone Cells in Chinese White Pear

: The accumulation of lignin in fruit has a significant negative impact on the quality of fruit-producing trees, and in particular the lignin formation stimulates the development of stone cells in pear fruit. Reactive oxygen species (ROS) are essential for lignin polymerization. However, knowledge of the RBOH family, a key enzyme in ROS metabolism, remains unknown in most fruit trees. In this study, a total of 40 RBOHs were identified from five fruit-producing trees (Pyrusbretschneideri, Prunuspersica, Citrussinensis, Vitisvinifera, and Prunusmume), and 10 of these sequences came from Pyrusbretschneideri. Multiple sequence alignments revealed that all 10 PbRBOHs contained the NADPH_Ox domain and the six alpha-helical transmembrane domains (TM-I to TM-VI). Chromosome localization and interspecies phylogenetic tree analysis showed that 10 PbRBOHs irregularly distributed on 8 chromosomes and 3 PbRBOHs (PbRBOHA, PbRBOHB, and PbRBOHD) are closely related to known lignification-related RBOHs. Furthermore, hormone response pattern analysis showed that the transcription of PbRBOHs is regulated by SA, ABA and MeJA. Reverse transcription-quantitative real-time polymerase chain reaction (qRT-PCR) and transcriptome sequencing analysis showed that PbRBOHA, PbRBOHB, and PbRBOHD accumulated high transcript abundance in pear fruit, and the transcriptional trends of PbRBOHA and PbRBOHD was consistent with the change of stone cell content during fruit development. In addition, subcellular localization revealed that PbRBOHA and PbRBOHD are distributed on the plasma membrane. Combining the changes of apoplastic superoxide (O2.-) content and spatio-temporal expression analysis, these results indicate that PbRBOHA and PbRBOHD, which are candidate genes, may play an important role in ROS metabolism during the lignification of pear stone cells. This study not only provided insight into the molecular characteristics of the RBOH family in fruit-producing trees, but also lays the foundation for studying the role of ROS in plant lignification.