Next Article in Journal
A Hybrid Prediction Method for Plant lncRNA-Protein Interaction
Next Article in Special Issue
In Silico Analysis of Bioactive Peptides in Invasive Sea Grass Halophila stipulacea
Previous Article in Journal
Obesity as an Inflammatory Agent Can Cause Cellular Changes in Human Milk due to the Actions of the Adipokines Leptin and Adiponectin
Previous Article in Special Issue
Multi-Path Dilated Residual Network for Nuclei Segmentation and Detection
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 8
Issue 6
10.3390/cells8060520
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Xi Cheng
1,†,
Guohui Li
1,†,
Muhammad Aamir Manzoor
1,
Han Wang
1,
Muhammad Abdullah
1,
Xueqiang Su
1,
Jingyun Zhang
2,
Taoshan Jiang
1,
Qing Jin
1,
Yongping Cai
1,* and
Yi Lin
1,*
1
School of Life Science, Anhui Agricultural University, No. 130, Changjiang West Road, Hefei 230036, China
2
Horticultural Institute, Anhui Academy of Agricultural Sciences, Hefei, Anhui 230031, China
*
Authors to whom correspondence should be addressed.
†
These authors have contributed equally.
Cells 2019, 8(6), 520; https://doi.org/10.3390/cells8060520
Submission received: 8 April 2019 / Revised: 22 May 2019 / Accepted: 28 May 2019 / Published: 29 May 2019
(This article belongs to the Special Issue Bioinformatics and Computational Biology 2019)
Browse Figures
Versions Notes

Abstract

:
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 (Pyrus bretschneideri, Prunus persica, Citrus sinensis, Vitis vinifera, and Prunus mume), and 10 of these sequences came from Pyrus bretschneideri. 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.

1. Introduction

Pyrus spp. is one of the world’s most important deciduous fruit trees. Asiatic pears or Oriental pears are mainly distributed in East Asia, and the main cultivars include P. bretschneideri, P. sinkiangensis, P. sinkiangensis, P. ussuriensis and P. pyrifolia [1,2,3]. Among them, Pyrus bretschneideri cv. Dangshan Su is the largest diploid (2n = 34) pear species in China. However, ‘Dangshan Su’ has the defect of a high content of large-diameter stone cell cluster in the fruit, which seriously affects its quality and economic value [4,5,6]. In addition, the large diameter and high content of stone cell clusters are also the “common problems” of some P. ussuriensis and P. pyrifolia varieties [1]. Therefore, reducing the size and content of the stone cell cluster of pear fruit and improving the fruit quality is an urgent problem to be solved in the pear industry.
Stone cells (also known as sclereids) exist in the xylem, phloem, leaves, fruits and seeds of trees, with support and protection functions [7,8]. The stone cells of pears accumulate specifically in the fruit, which is formed by the thickening and lignification of the secondary cell wall of the parenchyma cells [9,10,11]. The multiple aggregated stone cells are composed of stone cell clusters. At present, a large number of studies have shown that the size, content and density of stone cell clusters not only have a significant negative impact on the texture of pear fruit, but also have a close relationship with its sucrose content [6,12,13,14].
The average content of lignin in pear stone cells is approximately 37%, and the highest is approximately 44% [15]. Pear stone cells can not only undergo lignin-specific color reaction with phloroglucinol-HCl (Wiesner staining), but also have autofluorescence [16,17]. Transmission electron microscopy (TEM) observations revealed a large amount of lignin deposition in the cell corners, compound middle lamellas and secondary cell walls of pear stone cells [18,19]. It can be seen that lignin is one of the key components of pear stone cells, and the biosynthesis of lignin is closely related to stone cell development [20,21]. It is worth noting that lignin synthesis not only has a negative impact on the texture of pear fruit, but also has a significant effect on the intrinsic quality of peach (Prunus persica), loquat (Eriobotrya japonica), pomegranate (Punica granatum), sweet orange (Citrus sinensis), and other fruit trees [22,23,24,25,26,27]. Therefore, studying lignin metabolism is important for fruit-producing trees.
The polymerization mechanism of lignin monomers is still unclear. Previous studies have speculated that lignin monomer-coupled polymeric proteins mainly include dirigent (DIR), peroxidase (POD), laccase (LAC), superoxide dismutase (SOD) and nicotinamide adenine dinucleotide phosphate (NADPH) oxidases [8]. Notably, reactive oxygen species (ROS) play an important role in the lignification of plant cells [8,28,29,30,31]. Plant NADPH oxidase, also known as respiratory burst oxidase homolog (RBOH) [28], is a homolog of the mammalian macrophage NADPH oxidase catalytic subunit gp91phox [8,32]. It catalyzes the production of O2.−, NADP+ and H+ by O2 and NADPH. Subsequent O2.− can be catalyzed as a substrate for SOD (EC 1.15.1.1) to form hydrogen peroxide (H2O2) and molecular oxygen (O2). H2O2 and O2 participate in the polymerization of lignin monomers as a second substrate for POD and LAC, respectively [8,28,29]. Finally, POD and LAC activate the lignin monomer radicals, resulting in a terminal polymerization of the free radical oligo/monomers with the extension units, and finally the polymerization is coupled to form a long chain of lignin polymer [8,28,33]. Currently, although the functional study of RBOH in biotic and abiotic stress is not uncommon in the field of stress biology, the role of RBOH in plant lignin synthesis and cell wall development has not been well-investigated until now. In addition, the DIR, POD, LAC, and SOD families related to lignin polymerization in pears have been identified [4,34,35,36], but the study of the RBOH family of pears has not been reported.
To this end, we identified the RBOH family at the genome-wide level of five fruit-producing trees, Pyrus bretschneideri, Prunus persica, Citrus sinensis, Vitis vinifera, and Prunus mume, and systematically analyzed molecular evolution, chromosomal localization, microsynteny, conserved domains, gene structures and other aspects. Secondly, through phylogenetic tree clustering, gene expression profiling and determination of reactive oxygen species, the lignification-related RBOHs in pears were screened and subcellular localization analysis was performed. This study not only clarified the molecular characteristics of the five fruit-producing tree RBOH families, but also laid the foundation for the later study of lignification of pear stone cells and improvement of fruit texture.

2. Materials and Methods

2.1. Collection and Identification of Respiratory Burst Oxidase Homologs (RBOHs) in Five Fruit-Producing Trees

Pyrus bretschneideri (Chinese white pear) genome files were downloaded from the GigaDB data (http://gigadb.org/dataset/100083). The genome databases of Prunus persica (peach) and Vitis vinifera (common grape vine) were obtained from the Phytozome database (https://phytozome.jgi.doe.gov/pz/portal.html). Prunus mume (mei) genome files were downloaded from the Genome Database for Rosaceae (GDR) (http://www.rosaceae.org/). Citrus sinensis (sweet orange) genome files were downloaded from the Citrus sinensis Annotation Project (CAP) (http://citrus.hzau.edu.cn/orange/download/index.php).
Using the amino acid sequences of ten Arabidopsis RBOH family members [37] as the query sequences, candidate RBOHs were screened in the genome database using the BLASTp program (p-value = 10−10) in BioEdit v7.0.9.0 [38]. After removing repeating and redundant sequences, Pfam protein families database 32.0 (http://pfam.sanger.ac.uk/search) [39] and the Simple Modular Architecture Research Tool (http://smart.embl-heidelberg.de/) [40] were used to identify whether the candidate sequence contained the NADPH_Ox (PF08414) domain. The properties of the sequence of each RBOH [including amino acid sequence length (aa), molecular weight (MW), isoelectric point (pI) and subcellular localization] were completed by the ExPASy-Compute pI/Mw tool (http://web.expasy.org/compute_pi/) and WoLF PSORT (http://www.genscript.com/wolf-psort.html). Amino acid multiple sequence alignments were performed using BioEdit v7.0.9.0. The three-dimensional structures of RBOHs were confirmed using the PHYRE2 Protein Fold Recognition Server [41].

2.2. Conserved Motif, Gene Structure and Evolution Analysis

The conserved motifs and gene structures of the RBOH families were performed using Gene Structure Display Server (http://gsds.cbi.pku.edu.cn/) and Multiple Em for Motif Elicitation (http://meme-suite.org/), respectively [42,43,44,45].
The species phylogenetic tree was obtained from Common Taxonomy Tree (https://www.ncbi.nlm.nih.gov/Taxonomy/CommonTree/wwwcmt.cgi). All phylogenetic trees in this study were completed using MEGA v. 5.1 software. The construction parameters are as follows: Statistical Method: Neighbor-joining; Test of phylogeny: Bootstrap method; No. of Bootstrap Replications: 1000; Model: p-distance; Gaps: Pairwise deletion.

2.3. RBOH Family Genes’ Physical Localization and Gene Duplications

Visualization of the chromosomal location of the RBOHs was done using MapInspect v.1 software [46,47]. The identification of tandem duplication and segmental duplication events is based on the method used as previously described [16].

2.4. Microsynteny and Cis-Acting Elements Analysis of RBOH Family Genes

Microsynteny analysis was completed by using the Plant Genome Duplication Database and Circos v. 0.69 software (http://circos.ca/software/download/circos/) [4,48].
A sequence of 1500 bp upstream of the start codon of each PbRBOH coding region (CDS) was extracted from the pear genome as a putative promoter. The identification of the type and distribution of cis-acting elements was carried out using the software PLANT CARE v. 1 (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/).

2.5. Plant Materials and Treatments

The annual stem segments, buds, flowers, mature leaves and fruits were sampled from 40-year-old Pyrus bretschneideri cv. Dangshan Su trees were grown in the orchard of Dangshan County, Anhui Province, China. We collected fruits at eight developmental stages: 15 days after flowering (DAF), 39 DAF, 47 DAF, 63 DAF, 79 DAF, 102 DAF and 145 DAF (mature period). For exogenous hormone treatment of pear fruit, pear trees with no pests and diseases, with the same age and plant height, were selected. Hormones were sprayed on the entire surface of the pear fruit using a sprayer at 39 DAF. The concentration of the hormone treatment [0.5 mmol/L abscisic acid (ABA), 0.5 mmol/L methyl jasmonate (MeJA), or 0.2 mmol/L salicylic acid (SA)] was determined according to Cheng et al. (2018 and 2019) [4,5]. Samples were collected at 0, 1, 2, and 3 h post-treatment (HPT), and pear fruits of uniform size were harvested at each time point and immediately frozen in liquid nitrogen.

2.6. Determination of Apoplastic Superoxide (O2.−) Content in Pear Fruit

Histochemical staining with nitroblue tetrazolium (NBT) was used to analyse the in situ accumulation of O2.− according to Huang et al. (2011) [49]. Apoplastic superoxide was measured using superoxide anion free radical detection kit (sulfonamide colorimetry) following the manufacturer’s instructions (Beijing Leagene Biotech Co. Ltd., Beijing, China).

2.7. RNA Extraction and Reverse Transcription-Quantitative Real-Time Polymerase Chain Reaction (RT-qPCR) Analysis

RNA extraction and reverse transcription were performed according to the instructions in the Plant RNA maxi kit (Biomiga, San Diego, CA, USA) and the EasyScript One-Step gDNA Removal and cDNA Synthesis SuperMix (TransGen Biotech Co., Beijing, China), respectively. Reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR) analysis was performed using TransStart Probe qPCR SuperMix (TransGen Biotech Co.) and detected using a CFX96 Touch™ Real-Time PCR Detection System (Bio-rad Laboratories, Inc., Singapore, Singapore). Relative expression levels were calculated using the 2−△△Ct method [50]. All RT-qPCR primers are listed in Table S1. In this study, Tubulin (accession No. AB239680.1) was used as an internal reference [4]. The transcriptome sequencing data have been deposited in the Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra) under accession numbers SUB2967341 [6].

2.8. Subcellular Localization of PbRBOHs

The specific primers (Table S1) were designed to amplify the PbRBOH (PbRBOHA and PbRBOHD) CDS, and two restriction sites (Xba I, Sma I) were introduced at both ends according to the information provided for the eukaryotic expression vector pCAMBIA1305-GFP. By enzyme digestion and ligation, the plant expression vectors pCAMBIA1304-PbRBOHA/D-GFP were obtained. The vector was transformed into Agrobacterium tumefaciens EHA105 by electroporation, and 4-6-week-old Nicotiana benthamiana were used for the subcellular localization analysis. The method of subcellular localization was based on Cheng et al. (2017) [16].

3. Results

3.1. Identification, Characterization and Genomic Distribution of RBOH Family Genes in Five Fruit-Producing Trees

BLASTp was performed in the genome database of five fruit-producing trees using the amino acid sequence of AtRBOHA~J as the query sequence [51]. After deleting the repeat sequence and the redundant sequence, it was finally determined that there were 10 RBOHs in Chinese white pear (PbRBOH), and 7 RBOHs in peach (PpRBOH), 8 RBOHs in mei (PmRBOH), 7 RBOHs in common grape vine (VvRBOH), and 8 RBOHs in sweet orange (CsRBOH). We named the 40 RBOHs separately and analyzed their amino acid properties (Table 1). The 40 predicted full-length RBOH proteins varied from 116 (PmRBOHA) to 1067 (PmRBOHE) amino acid residues and the relative molecular mass ranged from 13.74 (PmRBOHA) to 121.09 (PmRBOHE) kDa, with isoelectric points in the range of 4.58–9.37.
Across PbRBOH family members, the MW varied from 79.09 to 111.34 kDa (average MW is 100.41 kDa), and the theoretical pI of these members ranged from 8.07 to 9.37 (average theoretical pI is 8.86). For PpRBOHs, VvRBOHs, PmRBOHs, and CsRBOHs, the average MW are 102.76 kDa, 99.64 kDa, 93.39 kDa, and 93.20 kDa, and the average theoretical pI are 9.00, 8.82, 8.38, and 8.75 respectively. Subcellular localization predictions indicated that 40 RBOHs derived from pear, peach, mei, grape and sweet orange were mainly located in chloroplast (22.5%), cytoplasmic (30%), nucleus (15%) and plasmodesmata (32.5%) (Table 1).
To expand knowledge of the RBOH family of the five fruit-producing trees, we compared the genome size, chromosome number and RBOH family distribution of the five fruit-producing trees from the common hexaploid ancestor [26]. We constructed a species phylogenetic tree and performed statistics on the number of RBOHs (Figure 1). The results showed that the range of the number (7–10) of RBOH family members among the five species was small. The percentage of putative RBOHs in the total number of predicted genes in five fruit-producing trees was also similar. The number of genes per megabase (Mb) ranged from 0.015 (common grape vine) to 0.034 (mei). This indicates that the number of RBOH family members has not increased due to the size of the species genome. Therefore, the expansion of these five fruit-producing tree RBOH families was not associated with an increase in the genome size and in the total number of predicted genes.

3.2. Chromosome Distribution and Duplication Events of RBOHs in Five Fruit-Producing Trees

In order to clarify the distribution of RBOH on the chromosomes of five fruit trees, we used MapInspect v.1 software to map the location of RBOH family members (Figure S1). The RBOH family members of the five fruit trees showed an irregular distribution on the chromosome and did not form a large number of gene clusters. According to the phylogenetic tree clustering of each RBOH family (Figure S2) and the rule of determining gene duplication events [16], only two segmental duplication events were identified in the pear RBOH family. The segmental duplicated gene pairs are PbRBOHA/PbRBOHB and PbRBOHH/PbRBOHG, respectively (Figure S1).

3.3. Phylogenetic Relationship, Conserved Motifs and Intron/Exon Structures of RBOHs in Five Fruit-Producing Trees

Through systematic phylogenetic tree analysis, 40 non-repetitive RBOHs derived from five fruit-producing trees can be divided into four subfamilies (I, II, III and IV). In addition, CsRBOF forms an independent clade. Except for subfamily IV, members of subfamily I–III are composed of RBOHs of five fruit-producing trees. Eleven gene pairs were formed among RBOHs with higher bootstrap values (≥99), except for the lower bootstrap values of two pairs: CsRBOHB/VvRBOHD and CsRBOHH/VvRBOHF (Figure 2A).
To investigate the structural diversity of RBOHs, a total of 20 conserved motifs in the RBOHs were captured by MEME v. 5.05 software (Figure 2B and Table S2). It is worth noting that the type and distribution of the C-terminal domains (motifs 3, 15, 2, 6, 10, 1, and 12) of most RBOHs are similar. However, the type and distribution of the N-terminal domains exhibit a certain subfamily specificity. For example, motif 16 is only present in members of subfamilies I and II. Members of subfamily IV do not contain motif 13.
In addition, the N-terminal domain of the members of subfamily I is generally composed of motifs 16, 19, and 13. Members of subfamily II can have two major clades, one of which has the same N-terminal domain composition as the members of subfamily I, and the N-terminal domain of another clade member consists only of motifs 16 and 13. The N-terminal domain of the subfamily IV members contains motifs 17 and 14, which lack a motif 13 compared to members of subfamily III (Figure 2B). These results suggest that differences in the composition of the N-terminal domain may confer different biological functions to members of each subfamily.
Gene structure analysis showed that the number of introns of RBOHs varied from 3 to 16 (Figure 2C). The exons of the RBOHs of the five fruit trees are all shorter, which may be characteristic of the RBOH family. Generally, most of RBOH family members in the same subfamilies showed similar motif characteristic and exon-intron structure, which supports their close evolutionary relationship and the classification of subfamilies.

3.4. Amino Acid Sequence and Characteristic Domain Analysis of PbRBOHs

We further analyzed the analysis of the characteristic domains of the RBOH family of pear (Figure 3). The results showed that the characteristic domains of the PbRBOH family members mainly include: NADPH_Ox (PF08414), EFh (IPR002048), Ferric_reduct (PF01794), FAD_binding_8 (PF08022), and NAD_binding_6 (PF08030). As can be seen from Figure 3, NADPH_Ox and EFh are mainly distributed at the N-terminus, while Ferric_reduct, FAD_binding_8, and NAD_binding_6 are distributed at the C-terminus. According to the annotated results of SMART (http://smart.embl-heidelberg.de/smart/show_motifs.pl), Ferric_reduct domain is a common region in the transmembrane proteins mammalian cytochrome B-245 heavy chain (gp91-phox). Plant RBOH encode homologs of the mammalian gp91phox [52], therefore, each PbRBOH contains this domain. NADPH_Ox is a characteristic domain of plant respiratory burst NADPH oxidase proteins, and PbRBOHA-J all have this domain, further demonstrating that the PbRBOH family members we identified are reliable. Interestingly, only PbRBOHJ lacks the NAD_binding_6 domain, which may cause changes in its function.
In addition, multiple sequence alignments of 10 PbRBOHs and 10 AtRBOHs were performed according to previously reported methods [32,53] (Figure 4). It can be found that the amino acid sequences of AtRBOH family members and PbRBOH family members are highly conserved, especially the typical domain of the RBOH family, such as two calcium-binding motifs (EF-hands) and six α-helical transmembrane domains (TM-I to TM-VI), and flavin adenine dinucleotide (FAD), nicotinamide adenine dinucleotide phosphate [NAD(P)H]-ribose, and NAD(P)H-adenine conserved binding sites in the C-terminal region. Notably, two amino acid residues (Pro-415 and Asp-500) in gp91phox are critical for catalytic activity [53]. However, PbRBOHJ lacks these two sites and the tertiary structure changes significantly, suggesting that the catalytic activity of this member will be affected (Figure 4 and Figure S3).

3.5. Microsynteny Analysis of PbRBOHs

To explore the evolutionary relationship of RBOHs, we performed microsynteny analysis of PbRBOHs and the coding genes in their surrounding chromosomal regions (Figure 5). It was found that there are 6 pairs of collinear gene pairs in the pear genome, including Pbr003403.1/Pbr036006.1, Pbr003403.1/Pbr037399.1, Pbr018609.1/Pbr037815.1, Pbr023445.1/Pbr033955.1, Pbr036006.1/ Pbr037399.1, and Pbr038641.1/Pbr038667.1. These 6 pairs of collinear gene pairs did not undergo chromosomal inversion during duplication (Figure 5).
Pear, sweet orange, mei, peach and grape share a common ancestor [4,26]. To further clarify the evolutionary origin and orthologous relationship of the RBOH family among the five fruit-producing trees, we identified the interspecific collinear relationship of these five species (Spreadsheets S1). Our results suggest that 18 collinear gene pairs (7 collinear gene pairs between pear and grape, 1 pair between pear and sweet orange, 3 pairs between sweet orange and grape, 3 pairs between mei and grape, 3 pairs between peach and grape, 1 gene pair is within the peach) were found among pear, sweet orange, mei, peach and grape (Spreadsheets S1). Overall, the degree of interspecies collinearity of CsRBOHs, PpRBOHs, VvRBOHs, and PmRBOHs is lower than intraspecies collinearity. These results indicate that the RBOH family of four species is relatively independent in evolution, and that chromosome duplication, deletion and rearrangement are later than species differentiation.

3.6. Comparative Phylogenetic Analysis and Functional Prediction of PbRBOHs

To characterize the evolutionary relationship and possible functions between RBOHs from five fruit trees and other known RBOHs from Arabidopsis, rice and poplar, a neighbor-joining tree was created (Figure 6).
The results showed that the 70 RBOHs from eight species (pear, peach, grape, mei, sweet orange, Arabidopsis, rice and poplar) could be assigned to 5 phylogenetic groups (I-V). With the exception of Group V, the other four phylogenetic groups are composed of RBOHs from eight species. Ten PbRBOHs were distributed in five phylogenetic groups, suggesting the existence of a diversified RBOH family in pear with diverse functions. Generally, RBOHs from pear have closer relationships with the PpRBOHs or PmRBOHs than that from VvRBOHs and CsRBOHs, which is in accordance with the current understanding of plant evolutionary history. In addition, evolutionary analysis also identified some closely related orthologous RBOHs between peach and mei, indicating that an ancestral set of RBOHs existed prior to the divergence of peach and mei.
For members of Group I and III, AtRBOHD and AtRBOHF are considered to be closely related to lignification, and they are also involved in the plant’s stress response process, which has a significant impact on the production of ROS in plants [8,28,31,54,55]. Phylogenetic tree clustering results indicated that RBOHA, RBOHB, and RBOHD in pears belonged to the same clades as AtRBOHD and AtRBOHF, suggesting that they may have similar biological functions (Figure 6).
It has now been demonstrated that AtRBOHC mutations result in inhibition of Arabidopsis root and root hair growth, and the level of ROS in plants is significantly reduced [54]. In addition, AtRBOHA, -B, and -G are all specifically expressed in the root [37]. It is speculated that RBOHs that are closely related to them may play a role mainly in the roots.
The functions of members in Group II and IV are currently poorly understood. AtRBOHH and -J are specifically expressed in pollen, AtRBOHE and OsRbohD are expressed in root and seeds [37]. Therefore, whether PbRBOHs of the same group has similar expression patterns and functions remains to be further verified. Group V contains only 5 RBOHs, while other groups contain at least 11 RBOHs. The five RBOHs in this group are derived from pear, sweet orange, mei and peach, respectively. It is speculated that these five RBOHs (CsRBOHD, CsRBOHF, PbRBOHI, PmRBOHG, and PpRBOHD) are a unique class of these four species, which may have new unknown functions.

3.7. Analysis of Cis-Acting Elements in Putative PbRBOH Promoters

Previous studies have shown that RBOH plays an important role in plant stress response, and their expression is regulated by a variety of plant hormones [37,56,57]. To this end, cis-acting elements of 10 PbRBOH promoter regions were identified in this study to explore their possible expression regulation mechanisms (Spreadsheets S2).
Analysis of the results showed that except for PbRBOHG, the other nine members contained G-box elements, which could respond to a variety of hormones, including ABA and MeJA [58]. In addition, W-box elements that can respond to SA, MeJA and ethylene were identified in the promoter regions of PbRBOHG, PbRBOHA, PbRBOHH, and PbRBOHC [5]. We found that 70% of the PbRBOH family members have ABRE (ABA-responsive) in their promoters and 50% of members contain TGACG-motif (MeJA-responsive). In addition to the hormone response elements, we also identified some abiotic stress response elements in the PbRBOH promoter region, such as LTR (low-temperature responsiveness), MBS (drought-inducibility), and ARE (anaerobic induction). Therefore, the transcription of PbRBOH is likely to be multi-regulated by hormones and the environment.

3.8. Differentially Expressed PbRBOHs under Hormonal Treatments

Based on the analysis of the cis-acting elements in the promoters of the PbRBOH family members, we found that most of the promoters of PbRBOHs contain a variety of plant hormone response-related elements (Spreadsheets S2). Therefore, we used three exogenous hormone (SA, ABA, and MeJA)-treated pear fruits as materials to analyze the effect of hormone treatment on PbRBOH expression patterns.
For SA-treated pear fruits, it was found that the transcription levels of PbRBOHC and PbRBOHJ were inhibited, and the remaining members peaked at 1 to 3 HPT. Among them, PbRBOHA, PbRBOHE, PbRBOHH and PbRBOHI showed a rise–fall tendency. PbRBOHB, PbRBOHF, PbRBOHG and PbRBOHD peaked at 3 HPT (Figure 7).
In the MeJA-treated pear fruits, the transcription levels of PbRBOHE and PbRBOHJ were inhibited, and the expression of other members was induced to varying degrees. Moreover, the transcriptional abundance of most members peaked at 1 HPT and 2 HPT, and decreased at 3 HPT (Figure 7).
ABA plays an inducing role in the transcription of 10 PbRBOHs. Among them, the expression levels of PbRBOHA and PbRBOHB reached a peak at 1 HPT, and then basically returned to the pre-treatment level. PbRBOHC, PbRBOHI and PbRBOJ have the highest expression at 2 HPT. The inducing effect of ABA on PbRBOHD lasted for a long time, and the transcript abundance of PbRBOHD from 1 to 3HPT was maintained at a high level. In addition, the transcript levels of PbRBOHG, PbRBOHF, PbRBOHE, and PbRBOHH were gradually increased after ABA treatment, reaching a peak at 3 HPT (Figure 7).

3.9. Histochemical Staining and Determination of Apoplastic Superoxide (O2.−) in Pear Fruit

To obtain some clues on the roles of PbRBOHs in the development of pear stone cells, this study determined the content of RBOH product (O2.−) in different developmental stages of pear fruit. The accumulation of O2.− in pear fruit was also analyzed by NBT staining. It can be seen from Figure 8 that the content of O2.− firstly increased and then decreased during different developmental stages of pear fruit, and reached a peak at 39 DAF. NBT-specific staining results also showed that O2.− accumulated at high levels in 39 DAF and 55 DAF pear fruits, and lower in 79 DAF and mature fruits. Previous studies have shown that high-density stone cells are distributed in pear flesh at 39 DAF and 55 DAF, and stone cell development and lignification have basically stopped at 79 DAF [10,16]. This indicates that the O2.− content change pattern is consistent with the development of stone cells and lignification, further confirming that the secondary cell wall lignification of stone cells requires the participation of O2.−.

3.10. Tissue Specificity and Temporal Expression Patterns in Pear Fruit during Development

Our initial interest in PbRBOHs came from the proposed role of PbRBOHs in cell wall lignification. In order to screen out RBOHs involved in the production of reactive oxygen species during the lignification of stone cells, we analyzed the expression patterns of 10 PbRBOHs in different developmental stages of pear fruit. As can be seen from Figure 9, the transcriptional levels of PbRBOHA, PbRBOHD, and PbRBOHF showed a rise-fall tendency, consistent with changes in the content of stone cells and O2.−. This suggests that they may be involved in the production of ROS during the lignification of stone cells.
Subsequently, we used the fruit transcriptome database of two Chinese white pear (Pyrus bretschneideri) cultivars that were previously completed [35] to analyze the transcript abundance of 10 PbRBOHs in different developmental stages of pear fruit (Figure 10). The fragments per kilobase of exon per million fragments mapped (FPKM) values for each PbRBOH are listed in Table S3. The results showed that although the expression changes of some members in pear fruit were consistent with the content of stone cells and O2.−, only PbRBOHA, PbRBOHB, and PbRBOHD had higher transcript abundance in fruits. The remaining members are weakly expressed in the fruit and may not play a major role in the fruit, such as PbRBOHF.
Therefore, combined with the results of phylogenetic tree clustering, trend of expression changes and transcript abundance, this study speculated that PbRBOHA and PbRBOHD are mainly responsible for the supply of ROS during the lignification of pear stone cells.
For buds, stems, leaves and flowers, except for PbRBOHC, the other 9 PbRBOHs expressed the highest in leaves. PbRBOHB and PbRBOHD also have higher levels of transcription in stems (Figure 11).

3.11. Subcellular Localization Analysis of Candidate PbRBOHs for Pear Stone Cell Lignification

Previous studies have shown that the RBOH is mainly localized on the plasma membrane and is responsible for catalyzing the formation of O2.− and thus participating in the polymerization of lignin monomers [8,37]. In order to clarify the subcellular distribution of the candidate genes, the eukaryotic expression vectors of PbRBOHA and PbRBOHD were constructed and transiently transformed into tobacco leaves. Green fluorescent protein (GFP) fluorescence was observed by laser scanning confocal microscopy (LSCM) (Olympus, Tokyo, Japan). As can be seen from Figure 12, both PbRBOHA-GFP and PbRBOHD-GFP are localized on the plasma membrane.

4. Discussion

There are five main types of lignified cells currently known in plants, including tracheary elements (TEs), sclerenchyma cells, endodermal cells, and seed coat cells [8]. Lignified cells are critical to plants, including resistance to external environmental conditions, seed protection, mechanical support, water and nutrient transport [7,8]. However, excessive lignification of the cells can have many negative effects on the fruit quality of trees [18,23,25,59]. Therefore, studying the molecular mechanism of lignification is of great significance for regulating fruit quality.
Numerous studies have shown that ROS play an important role in the lignification process [28,31]. However, the research on the RBOH family of catalytically lignification-related ROS production in fruit trees has rarely been reported. For this reason, 40 RBOHs were screened from the genome database of five fruit trees. Among them, the number of PbRBOH family members was the highest, and the number of PpRBOH family members was the least (Table 1). Phylogenetic tree analysis showed that the RBOH family members of pear, mei and peach were closely related, but they were far from sweet orange and grape, probably because the first three species belonged to Rosaceae (Figure 2). In particular, pear, mei, peach and sweet orange all differentiated into a single group (Group V) RBOH family members. The members of this subfamily are very distantly related to the RBOH members of poplar, rice, grape and Arabidopsis (Figure 6). It may already have completely different biological functions.
The formation of pear stone cells is a typical lignification process [19]. It provides a good model for the formation of lignification-related ROS. Therefore, this study focused on the content of O2.− in pear fruit and the expression pattern of PbRBOHs. The distribution of stone cell clusters in different developmental stages of pear fruit has been revealed in our previous studies [16]. The results show that in the fruit of 39 DAF to 55 DAF, the stone cell cluster accumulates in the pulp between the fruit core and the peels. As expected, NBT staining is the deepest in the pear fruit at this stage of development, suggesting that a large amount of O2.− is produced at this time to promote stone cell development and lignification (Figure 8). Peroxidase and laccase catalyzed the formation of monolignol radicals must have H2O2 and O2 as additional substrates, while O2.− is H2O2 and O2 direct precursor [8,29], so the content of O2.− directly affects the lignification polymerization process.
Analysis of cis-acting elements revealed that most of the PbRBOH promoters contained a large number of elements related to hormones and stress response (Spreadsheets S2). This is consistent with previous studies [57]. We selected three hormones (SA, ABA, and MeJA) to treat pear fruit and found that it has a significant effect on the expression of PbRBOHs, and the transcription of most PbRBOHs is induced in a short time (Figure 7). Notably, not all of the 10 PbRBOH promoters have elements that respond to the three hormones. For example, the SA response element has not been found in the PbRBOHI promoter. However, the expression level of PbRBOHI still changed in the fruits treated with SA (Spreadsheets S2 and Figure 7). This may be due to interactions between various plant hormones in plants, which can promote synergy and induce one another [60]. Therefore, spraying a hormone is likely to induce an increase in the content of another hormone, thereby affecting the transcription of the gene. In addition, photoresponsive elements such as Box4, GATA-motif, GT1-motif are also present in the RBOH promoters. This suggests that bagging treatment changes the lignin and stone cell content of pear fruit may be related to the transcription of RBOH [61,62].
Importantly, in order to identify the RBOHs responsible for lignification-related ROS synthesis in pear fruit, we compared the expression patterns of 10 PbRBOHs in pear fruit at 7 developmental stages (Figure 9). Phylogenetic tree clustering indicated that PbRBOHA, PbRBOHB, and PbRBOHD were most likely related to lignification, and transcriptome sequencing also showed that the transcript abundance of PbRBOHA, PbRBOHB and PbRBOHD was the highest among the 10 members in different pear varieties (Figure 10). It is indicated that these three RBOHs play a role in pear fruit. However, the expression level of PbRBOHB showed a downward trend during the large accumulation period of stone cell (47 DAF and 55 DAF) (Figure 9) [10,16], and the expression level was up-regulated after the stone cell development was basically stopped (63 DAF) [10], which was inconsistent with the process of stone cell development. It is speculated that although this gene has a high expression level in fruits, it is not the major gene involved in the synthesis of ROS during the development of stone cells. In addition, the expression levels of PbRBOHA, PbRBOHD and PbRBOHF in pear fruit showed a rise-fall tendency, which was consistent with the trend of stone cell content (Figure 9) [10,61,63]. However, transcriptome sequencing showed that PbRBOHF has a weak transcription level in fruits (Figure 10), so it is speculated that this gene has little correlation with stone cell development. Therefore, PbRBOHA and PbRBOHD are considered to be lignification-related RBOH in pears. Notably, the up-regulation of PbRBOHA is most significant at 55 DAF (Figure 10 and Table S3), suggesting that the gene may be the major gene responsible for lignification-related ROS production.
More strikingly, we summarized the results of this study with the reported pear uridine diphosphate-glycosyltransferases (UGTs), PODs, LACs and DIRs in an ideograph of the polymerization of lignin monomer in pear fruit (Figure 13) [4,19,35,36]. The monolignols in pear fruit may be catalyzed by PbUGT72AJ1, PbUGT72AJ2 and PbUGT72AK1 to form lignin glycoside, which is then transported across the membrane to cell wall deposition. Subsequently, the lignin oligomers are formed by the catalysis of PbDIR4. Finally, POD and LAC use ROS and O2 produced by RBOH (PbRBOHA and PbRBOHD) and SOD to catalyze the elongation of lignin polymer chains.

5. Conclusions

In this study, we identified 10, 7, 8, 8 and 7 members of the RBOH family from the pear, peach, mei, sweet orange and grape genomes, respectively. RBOHs derived from fruit-producing trees can be divided into five subfamilies, and the distribution of conserved motifs and gene structures have a certain subfamily specificity. The 10 RBOH family members of the pear all have six α-helical transmembrane domains, and there are a large number of cis-acting elements related to biotic and abiotic stress in the promoter. Hormone treatment revealed that the transcription of most PbRBOHs was induced by SA, ABA and MeJA. Analysis of spatio-temporal expression patterns showed that PbRBOHA, PbRBOHB, and PbRBOHD had higher transcript abundance in pear fruit, and the expression trends of PbRBOHA and PbRBOHD were consistent with the content of stone cells and O2.−. Combined with the results of phylogenetic tree clustering and expression pattern analysis, it is speculated that PbRBOHA and PbRBOHD participate in the production of reactive oxygen species during the formation of stone cells. Subcellular localization analysis indicated that both RBOHs were localized to the plasma membrane.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4409/8/6/520/s1.

Author Contributions

Conceptualization, X.C, G.L., Y.C. and Y.L.; Data curation, X.C., G.L., M.A.M., H.W., M.A. and X.S.; Formal analysis, X.C., G.L., M.A.M., H.W., M.A., X.S. and J.Z.; Funding acquisition, T.J., Q.J., Y.C. and Y.L.; Investigation, X.C., G.L., M.A.M., H.W., M.A. and X.S.; Methodology, X.C., G.L., M.A.M., H.W., M.A., X.S., J.Z. and T.J.; Project administration, T.J., Q.J., Y.C. and Y.L.; Resources, X.C., G.L., M.A.M., J.Z. and T.J.; Software, X.C., H.W., M.A. and J.Z.; Validation, G.L., M.A.M., H.W., M.A., X.S. and Q.J.; Visualization, X.C. and H.W.; Writing—original draft, X.C., G.L., Y.C. and Y.L.; Writing—review and editing, X.C., M.A.M., Y.C. and Y.L.

Funding

This study was supported by the National Natural Science Foundation of China (grant #31640068) and Anhui Provincial Natural Science Foundation (grant #1808085QC79).

Acknowledgments

We thank the graduate student Yanming Sun for her assistance in the experimental process. We also extend our thanks to the reviewers and editors for their careful reading and helpful comments on this manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Wu, J.; Wang, Y.; Xu, J.; Korban, S.S.; Fei, Z.; Tao, S.; Ming, R.; Tai, S.; Khan, A.M.; Postman, J.D.; et al. Diversification and independent domestication of Asian and European pears. Genome Biol. 2018, 19, 1–16. [Google Scholar] [CrossRef] [PubMed]
  2. Tian, L.; Cao, Y.; Gao, Y.; Dong, X. Effect of stone cells size and flesh texture in pear cultivars. Acta Hortic. Sin. 2011, 38, 1225–1234. [Google Scholar]
  3. Bao, L.; Chen, K.; Zhang, D.; Cao, Y.; Yamamoto, T.; Teng, Y. Genetic diversity and similarity of pear (Pyrus L.) cultivars native to East Asia revealed by SSR (simple sequence repeat) markers. Genet. Resour. Crop Evol. 2007, 54, 959–971. [Google Scholar] [CrossRef]
  4. Cheng, X.; Su, X.; Muhammad, A.; Li, M.; Zhang, J.; Sun, Y.; Li, G.; Jin, Q.; Cai, Y.; Lin, Y. Molecular characterization, evolution, and expression profiling of the dirigent (DIR) family genes in Chinese white pear (Pyrus bretschneideri). Front. Genet. 2018, 9, 1–15. [Google Scholar] [CrossRef]
  5. Cheng, X.; Li, G.; Muhammad, A.; Zhang, J.; Jiang, T.; Jin, Q.; Zhao, H.; Cai, Y.; Lin, Y. Molecular identification, phylogenomic characterization and expression patterns analysis of the LIM (LIN-11, Isl1 and MEC-3 domains) gene family in pear (Pyrus bretschneideri) reveal its potential role in lignin metabolism. Gene 2019, 686, 237–249. [Google Scholar] [CrossRef]
  6. Zhang, J.; Cheng, X.; Jin, Q.; Su, X.; Li, M.; Yan, C.; Jiao, X.; Li, D.; Lin, Y.; Cai, Y. Comparison of the transcriptomic analysis between two Chinese white pear (Pyrus bretschneideri Rehd.) genotypes of different stone cells contents. PLoS ONE 2017, 12, 1–22. [Google Scholar] [CrossRef]
  7. Zhao, M.; Zhu, H. Development and morphology of stone cells in phloem of Toxicodendron vernicifluum. Trees Struct. Funct. 2014, 28, 1553–1558. [Google Scholar] [CrossRef]
  8. Barros, J.; Serk, H.; Granlund, I.; Pesquet, E. The cell biology of lignification in higher plants. Ann. Bot. 2015, 115, 1053–1074. [Google Scholar] [CrossRef] [Green Version]
  9. Xue, C.; Yao, J.L.; Qin, M.F.; Zhang, M.Y.; Allan, A.C.; Wang, D.F.; Wu, J. PbrmiR397a regulates lignification during stone cell development in pear fruit. Plant Biotechnol. J. 2018, 1–15. [Google Scholar] [CrossRef]
  10. Cai, Y.; Li, G.; Nie, J.; Lin, Y.; Nie, F.; Zhang, J.; Xu, Y. Study of the structure and biosynthetic pathway of lignin in stone cells of pear. Sci. Hortic. 2010, 125, 374–379. [Google Scholar] [CrossRef]
  11. Cheng, X.; Yan, C.; Zhang, J.; Ma, C.; Li, S.; Jin, Q.; Zhang, N.; Cao, Y.; Lin, Y.; Cai, Y. The effect of different pollination on the expression of Dangshan Su pear microRNA. Biomed. Res. Int. 2017, 2017, 2794040. [Google Scholar] [CrossRef]
  12. Yan, C.; Yin, M.; Zhang, N.; Jin, Q.; Fang, Z.; Lin, Y.; Cai, Y. Stone cell distribution and lignin structure in various pear varieties. Sci. Hortic. 2014, 174, 142–150. [Google Scholar] [CrossRef]
  13. Choi, J.; Choi, J.; Hong, K.; Kim, W. Cultivar differences of stone cells in pear flesh and their effects on fruit quality. Hortic. Environ. Biotechnol. 2007, 48, 17–31. [Google Scholar]
  14. Choi, J.; Lee, S. Distribution of stone cell in Asian, Chinese, and European pear fruit and its morphological changes. J. Appl. Bot. Food Qual. 2013, 189, 185–189. [Google Scholar] [CrossRef]
  15. Tian, L.; Dong, X.; Cao, Y.; Zhang, Y.; Qi, D. Correlation of flesh in Pyrus fruit with its stone cells lignin. Southwest China J. Agric. Sci. 2017, 30, 2091–2096. [Google Scholar] [CrossRef]
  16. Cheng, X.; Li, M.; Li, D.; Zhang, J.; Jin, Q.; Sheng, L.; Cai, Y.; Lin, Y. Characterization and analysis of CCR and CAD gene families at the whole-genome level for lignin synthesis of stone cells in pear (Pyrus bretschneideri) fruit. Biol. Open 2017, 6, 1602–1613. [Google Scholar] [CrossRef]
  17. Tao, S.; Khanizadeh, S.; Zhang, H.; Zhang, S. Anatomy, ultrastructure and lignin distribution of stone cells in two Pyrus species. Plant Sci. 2009, 176, 413–419. [Google Scholar] [CrossRef]
  18. Jin, Q.; Yan, C.; Qiu, J.; Zhang, N.; Lin, Y.; Cai, Y. Structural characterization and deposition of stone cell lignin in Dangshan Su pear. Sci. Hortic. 2013, 155, 123–130. [Google Scholar] [CrossRef]
  19. Cheng, X.; Muhammad, A.; Li, G.; Zhang, J.; Cheng, J.; Qiu, J.; Jiang, T. Family-1 UDP glycosyltransferases in pear (Pyrus bretschneideri): Molecular identification, phylogenomic characterization and expression profiling during stone cell formation. Mol. Biol. Rep. 2019, 46, 2153–2175. [Google Scholar] [CrossRef]
  20. Brahem, M.; Renard, C.M.G.C.; Gouble, B.; Bureau, S.; Le Bourvellec, C. Characterization of tissue specific differences in cell wall polysaccharides of ripe and overripe pear fruit. Carbohydr. Polym. 2017, 156, 152–164. [Google Scholar] [CrossRef]
  21. Li, S.; Su, X.; Abdullah, M.; Sun, Y.; Li, G.; Cheng, X.; Lin, Y.; Cai, Y.; Jin, Q. Effects of different pollens on primary metabolism and lignin biosynthesis in pear. Int. J. Mol. Sci. 2018, 19, 2273. [Google Scholar] [CrossRef]
  22. Zarei, A.; Zamani, Z.; Fatahi, R.; Mousavi, A.; Salami, S.A.; Avila, C.; Cánovas, F.M. Differential expression of cell wall related genes in the seeds of soft- and hard-seeded pomegranate genotypes. Sci. Hortic. 2016, 205, 7–16. [Google Scholar] [CrossRef]
  23. Wang, Y.; Zhang, X.; Yang, S.; Yuan, Y. Metabolite and Transcriptome analyses indicate the involvement of lignin in programmed changes in peach fruit texture. J. Agric. Food Chem. 2018, 66, acs.jafc.8b04284. [Google Scholar] [CrossRef]
  24. Jia, N.; Liu, J.; Sun, Y.; Tan, P.; Cao, H.; Xie, Y.; Wen, B.; Gu, T.; Liu, J.; Li, M.; et al. Citrus sinensis MYB transcription factors CsMYB330 and CsMYB308 regulate fruit juice sac lignification through fine-tuning expression of the Cs4CL1 gene. Plant Sci. 2018, 277, 334–343. [Google Scholar] [CrossRef] [PubMed]
  25. Xu, Q.; Yin, X.R.; Zeng, J.K.; Ge, H.; Song, M.; Xu, C.J.; Li, X.; Ferguson, I.B.; Chen, K.S. Activator-and repressor-type MYB transcription factors are involved in chilling injury induced flesh lignification in loquat via their interactions with the phenylpropanoid pathway. J. Exp. Bot. 2014, 65, 4349–4359. [Google Scholar] [CrossRef]
  26. Wu, J.; Wang, Z.; Shi, Z.; Zhang, S.; Ming, R.; Zhu, S.; Khan, M.A.; Tao, S.; Korban, S.S.; Wang, H.; et al. The genome of the pear (Pyrus bretschneideri Rehd.). Genome Res. 2013, 23, 396–408. [Google Scholar] [CrossRef] [PubMed]
  27. Wang, W.; Wang, C.; Xie, Z. Function analysis of VvmiR397a and its target genes VvLACs in grape berry development. Acta Hortic. Sin. 2018, 45, 1441–1455. [Google Scholar] [CrossRef]
  28. Lee, Y.; Rubio, M.C.; Alassimone, J.; Geldner, N. A mechanism for localized lignin deposition in the endodermis. Cell 2013, 153, 402–412. [Google Scholar] [CrossRef] [PubMed]
  29. Koutaniemi, S. Lignin biosynthesis studies in plant tissue cultures. J. Integr. Plant Biol. 2010, 52, 176–185. [Google Scholar] [CrossRef]
  30. Shafi, A.; Gill, T.; Zahoor, I.; Singh, P.; Yelam, A.; Sanjay, S. Ectopic expression of SOD and APX genes in Arabidopsis alters metabolic pools and genes related to secondary cell wall cellulose biosynthesis and improve salt tolerance. Mol. Biol. Rep. 2019, 46, 1985–2002. [Google Scholar] [CrossRef] [PubMed]
  31. Lee, Y.; Yoon, T.H.; Lee, J.; Cho, H.K.; Mang, H.; Kwak, J.M.; Lee, Y.; Yoon, T.H.; Lee, J.; Jeon, S.Y.; et al. A lignin molecular brace controls precision processing of cell walls critical for surface integrity in Arabidopsis. Cell 2018, 173, 1468–1480. [Google Scholar] [CrossRef]
  32. Zhang, Y.; Li, Y.; He, Y.; Hu, W.; Zhang, Y.; Wang, X.; Tang, H. Identification of NADPH oxidase family members associated with cold stress in strawberry. FEBS Open Bio 2018, 8, 593–605. [Google Scholar] [CrossRef]
  33. Meents, M.J.; Watanabe, Y.; Samuels, A.L. The cell biology of secondary cell wall biosynthesis. Ann. Bot. 2018, 1107–1125. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, L.; Wang, L.; Zhang, Z.; Ma, M.; Wang, R.; Qian, M.; Zhang, S. Genome-wide identification and comparative analysis of the superoxide dismutase gene family in pear and their functions during fruit ripening. Postharvest Biol. Technol. 2018, 143, 68–77. [Google Scholar] [CrossRef]
  35. Xi, C.; Li, G.; Ma, C.; Abdullah, M.; Zhang, J.; Zhao, H.; Jin, Q.; Cai, Y.; Lin, Y. Comprehensive genome-wide analysis of the pear (Pyrus bretschneideri) laccase gene (PbLAC) family and functional identification of PbLAC1 involved in lignin biosynthesis. PLoS ONE 2019, 14, e0210892. [Google Scholar]
  36. Cao, Y.; Han, Y.; Meng, D.; Li, D.; Jin, Q.; Lin, Y. Structural, evolutionary, and functional analysis of the Class III peroxidase gene family in Chinese pear (Pyrus bretschneideri). Front. Plant Sci. 2016, 7, 1874. [Google Scholar] [CrossRef] [PubMed]
  37. Sagi, M. Production of reactive oxygen species by plant NADPH oxidases. Plant Physiol. 2006, 141, 336–340. [Google Scholar] [CrossRef]
  38. Hall, T.A. BioEdit: A user-friendly biological sequence alignment program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  39. El-Gebali, S.; Mistry, J.; Bateman, A.; Eddy, S.R.; Luciani, A.; Potter, S.C.; Qureshi, M.; Richardson, L.J.; Salazar, G.A.; Smart, A.; et al. The Pfam protein families database in 2019. Nucleic Acids Res. 2018, 47, 427–432. [Google Scholar] [CrossRef]
  40. Letunic, I.; Bork, P. 20 years of the SMART protein domain annotation resource. Nucleic Acids Res. 2018, 46, D493–D496. [Google Scholar] [CrossRef]
  41. Kelley, L.A.; Mezulis, S.; Yates, C.M.; Wass, M.N.; Sternberg, M.J.E. The Phyre2 web portal for protein modeling, prediction and analysis. Nat. Protoc. 2015, 10, 845–858. [Google Scholar] [CrossRef] [Green Version]
  42. Liu, F.; Xu, Y.; Jiang, H.; Jiang, C.; Du, Y.; Gong, C.; Wang, W.; Zhu, S.; Han, G.; Cheng, B. Systematic identification, evolution and expression analysis of the Zea mays PHT1 gene family reveals several new members involved in root colonization by arbuscular mycorrhizal fungi. Int. J. Mol. Sci. 2016, 17, 1. [Google Scholar] [CrossRef] [PubMed]
  43. Lin, Y.X.; Jiang, H.Y.; Chu, Z.X.; Tang, X.L.; Zhu, S.W.; Cheng, B.J. Genome-wide identification, classification and analysis of heat shock transcription factor family in maize. BMC Genom. 2011, 12, 76. [Google Scholar] [CrossRef] [PubMed]
  44. Hu, B.; Jin, J.; Guo, A.; Zhang, H.; Luo, J. GSDS 2.0: An upgraded gene feature visualization server. Bioinformatics 2015, 31, 1296–1297. [Google Scholar] [CrossRef] [PubMed]
  45. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME SUITE: Tools for motif discovery and searching. Nucleic Acids Res. 2009, 37, 202–208. [Google Scholar] [CrossRef] [PubMed]
  46. Cheng, X.; Xiong, Y.; Li, D.H.; Cheng, J.; Cao, Y.P.; Yan, C.C.; Jin, Q.; Sun, N.; Cai, Y.P.; Lin, Y. Bioinformatic and expression analysis of the OMT gene family in Pyrus bretschneideri cv. Dangshan Su. Genet. Mol. Res. 2016, 15, 1–17. [Google Scholar] [CrossRef]
  47. Abdullah, M.; Cao, Y.; Cheng, X.; Shakoor, A.; Su, X.; Gao, J.; Cai, Y. Genome-wide analysis characterization and evolution of SBP Genes in Fragaria vesca, Pyrus bretschneideri, Prunus persica and Prunus mume. Front. Genet. 2018, 9, 1–12. [Google Scholar] [CrossRef]
  48. Krzywinski, M.; Schein, J.; Birol, I.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Huang, X.; Luo, T.; Fu, X.; Fan, Q.; Liu, J. Cloning and molecular characterization of a mitogen-activated protein kinase gene from Poncirus trifoliata whose ectopic expression confers dehydration/drought tolerance in transgenic tobacco. J. Exp. Bot. 2011, 62, 5191–5206. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  51. Kaur, G.; Pati, P.K. In silico insights on diverse interacting partners and phosphorylation sites of respiratory burst oxidase homolog (Rbohs) gene families from Arabidopsis and rice. BMC Plant Biol. 2018, 18, 161. [Google Scholar] [CrossRef]
  52. Kaur, G.; Guruprasad, K.; Temple BR, S.; Shirvanyants, D.G.; Dokholyan, N.V.; Pati, P.K. Structural complexity and functional diversity of plant NADPH oxidases. Amino Acids 2018, 50, 79–94. [Google Scholar] [CrossRef]
  53. Lightfoot, D.; Boettcher, A.; Little, A.; Shirley, N.; Able, A.J. Identification and characterisation of barley (Hordeum vulgare) respiratory burst oxidase homologue family members. Funct. Plant Biol. 2008, 35, 347–359. [Google Scholar] [CrossRef]
  54. Zhu, Y.; Du, B.; Qian, J.; Zou, B.; Hua, J. Disease Resistance gene-induced growth inhibition is enhanced by rcd1 independent of defense activation in Arabidopsis. Plant Physiol. 2013, 161, 2005–2013. [Google Scholar] [CrossRef]
  55. Kwak, J.M.; Mori, I.C.; Pei, Z.; Leonhardt, N.; Torres, M.A.; Dangl, J.L.; Bloom, R.E.; Bodde, S.; Jones JD, G.; Schroeder, J.I. NADPH oxidase AtrbohD and AtrbohF genes function in ROS-dependent ABA signaling in Arabidopsis. EMBO J. 2003, 22, 2623–2633. [Google Scholar] [CrossRef]
  56. Cepauskas, D.; Miliute, I.; Staniene, G.; Gelvonauskiene, D.; Stanys, V.; Jesaitis, A.J.; Baniulis, D. Characterization of apple NADPH oxidase genes and their expression associated with oxidative stress in shoot culture in vitro. Plant Cell. Tissue Organ Cult. 2016, 124, 621–633. [Google Scholar] [CrossRef]
  57. Kaur, G.; Pati, P.K. Analysis of cis-acting regulatory elements of Respiratory burst oxidase homolog (Rboh) gene families in Arabidopsis and rice provides clues for their diverse functions. Comput. Biol. Chem. 2016, 62, 104–118. [Google Scholar] [CrossRef] [PubMed]
  58. Shen, Q.; Ho, T.D. Functional dissection of an abscisic acid (ABA)-inducible gene reveals two independent ABA-responsive complexes each containing a G-Box and a nove1 cis-acting element. Plant Cell 1995, 7, 295–307. [Google Scholar] [CrossRef]
  59. Li, X.; Singh, J.; Qin, M.; Li, S.; Zhang, X.; Zhang, M.; Khan, A.; Zhang, S.; Wu, J. Development of an integrated 200K SNP genotyping array and application for genetic mapping, genome assembly improvement and genome wide association studies in pear (Pyrus). Plant Biotechnol. J. 2019. [Google Scholar] [CrossRef] [PubMed]
  60. Berens, M.L.; Berry, H.M.; Mine, A.; Argueso, C.T.; Tsuda, K. Evolution of hormone signaling networks in plant defense. Annu. Rev. Phytopathol. 2017, 55, 401–425. [Google Scholar] [CrossRef]
  61. Tao, S.; Wang, D.; Jin, C.; Sun, W.; Liu, X.; Zhang, S.; Gao, F.; Khanizadeh, S. Cinnamate-4-hydroxylase gene is involved in the step of lignin biosynthesis in Chinese white pear. J. Am. Soc. Hortic. Sci. 2015, 140, 573–579. [Google Scholar] [CrossRef]
  62. Wang, Y.; Zhang, X.; Wang, R.; Bai, Y.; Liu, C.; Yuan, Y.; Yang, Y.; Yang, S. Differential gene expression analysis of “Chili” (Pyrus bretschneideri) fruit pericarp with two types of bagging treatments. Hortic. Res. 2017, 4, 1–10. [Google Scholar] [CrossRef] [PubMed]
  63. Li, J.M.; Huang, X.S.; Li, L.T.; Zheng, D.M.; Xue, C.; Zhang, S.L.; Wu, J. Proteome analysis of pear reveals key genes associated with fruit development and quality. Planta 2015, 241, 1363–1379. [Google Scholar] [CrossRef] [PubMed]
Cells 08 00520 g001 550
Figure 1. A summary of the evolutionary relationships, genome size, and RBOH distribution of the five fruit-producing trees.
Figure 1. A summary of the evolutionary relationships, genome size, and RBOH distribution of the five fruit-producing trees.
Cells 08 00520 g001
Cells 08 00520 g002 550
Figure 2. The gene structures and conserved motifs of RBOH family members of five fruit-producing trees based on the evolutionary relationship. (A) Phylogenetic tree of RBOH family members of five fruit-producing trees. (B) Conserved motifs of RBOHs of five fruit-producing trees. (C) The exon-intron structure of PbRBOHs of five fruit-producing trees.
Figure 2. The gene structures and conserved motifs of RBOH family members of five fruit-producing trees based on the evolutionary relationship. (A) Phylogenetic tree of RBOH family members of five fruit-producing trees. (B) Conserved motifs of RBOHs of five fruit-producing trees. (C) The exon-intron structure of PbRBOHs of five fruit-producing trees.
Cells 08 00520 g002
Cells 08 00520 g003 550
Figure 3. Depiction of various domains in PbRBOHs.
Figure 3. Depiction of various domains in PbRBOHs.
Cells 08 00520 g003
Cells 08 00520 g004 550
Figure 4. Protein sequence multi-alignment and domain structure of RBOHs from pear and Arabidopsis. Conservative residues are highlighted by black asterisks. Histidine residues involved in heme binding and the gp91phox residues Pro-415 and Asp-500 are indicated with diamonds. EF-hands and conserved binding sites for FAD, NAD(P)H-adenine, and NAD(P)H-ribose are represented by straight lines or boxes. Transmembrane domains are indicated by straight lines.
Figure 4. Protein sequence multi-alignment and domain structure of RBOHs from pear and Arabidopsis. Conservative residues are highlighted by black asterisks. Histidine residues involved in heme binding and the gp91phox residues Pro-415 and Asp-500 are indicated with diamonds. EF-hands and conserved binding sites for FAD, NAD(P)H-adenine, and NAD(P)H-ribose are represented by straight lines or boxes. Transmembrane domains are indicated by straight lines.
Cells 08 00520 g004
Cells 08 00520 g005 550
Figure 5. Microsynteny maps of PbRBOHs and their flanking genes within syntenic chromosomal intervals. The relative positions of all flanking protein-coding genes were defined by the anchored RBOHs, highlighted in red. The chromosome segments are shown as horizontal lines, with arrows corresponding to individual genes and their transcriptional orientations. Conserved gene pairs among the segments are connected with lines.
Figure 5. Microsynteny maps of PbRBOHs and their flanking genes within syntenic chromosomal intervals. The relative positions of all flanking protein-coding genes were defined by the anchored RBOHs, highlighted in red. The chromosome segments are shown as horizontal lines, with arrows corresponding to individual genes and their transcriptional orientations. Conserved gene pairs among the segments are connected with lines.
Cells 08 00520 g005
Cells 08 00520 g006 550
Figure 6. Phylogenetic analysis of RBOHs from five fruit-producing trees (pear, sweet orange, mei, peach and grape) and other known RBOHs. Full-length polypeptide sequences were used to make the interspecific phylogenetic tree.
Figure 6. Phylogenetic analysis of RBOHs from five fruit-producing trees (pear, sweet orange, mei, peach and grape) and other known RBOHs. Full-length polypeptide sequences were used to make the interspecific phylogenetic tree.
Cells 08 00520 g006
Cells 08 00520 g007 550
Figure 7. Heatmap representation of PbRBOHs expression in response to salicylic acid (SA), methyl jasmonate (MeJA) and abscisic acid (ABA) treatments. The gene expression data of PbRBOHs at 0 h post-treatment (HPT), 1 HPT, 2 HPT and 3 HPT were obtained through reverse transcription-quantitative real-time polymerase chain reaction (qRT-PCR). The expression data of each gene at 0 HPT was used as a control sample to show the relative expression level. Purple and yellow boxes indicate lower or higher transcript levels, respectively. Each qRT-PCR analysis was performed in triplicate.
Figure 7. Heatmap representation of PbRBOHs expression in response to salicylic acid (SA), methyl jasmonate (MeJA) and abscisic acid (ABA) treatments. The gene expression data of PbRBOHs at 0 h post-treatment (HPT), 1 HPT, 2 HPT and 3 HPT were obtained through reverse transcription-quantitative real-time polymerase chain reaction (qRT-PCR). The expression data of each gene at 0 HPT was used as a control sample to show the relative expression level. Purple and yellow boxes indicate lower or higher transcript levels, respectively. Each qRT-PCR analysis was performed in triplicate.
Cells 08 00520 g007
Cells 08 00520 g008 550
Figure 8. Analysis of apoplastic superoxide content in pear fruit at different developmental stages. (A) The content of O2.− in pear fruits at different developmental stages. (B) Specific histochemical staining of apoplastic superoxide in pear pulp.
Figure 8. Analysis of apoplastic superoxide content in pear fruit at different developmental stages. (A) The content of O2.− in pear fruits at different developmental stages. (B) Specific histochemical staining of apoplastic superoxide in pear pulp.
Cells 08 00520 g008
Cells 08 00520 g009 550
Figure 9. Expression profiles of pear RBOHs across different developmental stages. The expression level of 10 PbRBOHs in pear fruits at eight developmental stages (15 DAF, 39 DAF, 47 DAF, 63 DAF, 79 DAF, and 145 DAF). The relative expression level was calculated by setting the expression value of PbRBOHs in fruit (15 DAF) of pear at 1.
Figure 9. Expression profiles of pear RBOHs across different developmental stages. The expression level of 10 PbRBOHs in pear fruits at eight developmental stages (15 DAF, 39 DAF, 47 DAF, 63 DAF, 79 DAF, and 145 DAF). The relative expression level was calculated by setting the expression value of PbRBOHs in fruit (15 DAF) of pear at 1.
Cells 08 00520 g009
Cells 08 00520 g010 550
Figure 10. Expression profiles of PbRBOHs in different pear varieties at different developmental stages. FPKM values were obtained by RNA-seq analysis and presented as a heatmap. The colour scale and circle size are provided with the heat map to indicate the differential pattern of expression. DS: Pyrus bretschneideri cv. Dangshan Su; LS: Pyrus bretschneideri cv. Lianglizaosu; 23: 23 DAF; 55: 55 DAF; MP: mature period. The heatmap is drawn using the TBTools v. 0.6644 software (http://www.tbtools.com/).
Figure 10. Expression profiles of PbRBOHs in different pear varieties at different developmental stages. FPKM values were obtained by RNA-seq analysis and presented as a heatmap. The colour scale and circle size are provided with the heat map to indicate the differential pattern of expression. DS: Pyrus bretschneideri cv. Dangshan Su; LS: Pyrus bretschneideri cv. Lianglizaosu; 23: 23 DAF; 55: 55 DAF; MP: mature period. The heatmap is drawn using the TBTools v. 0.6644 software (http://www.tbtools.com/).
Cells 08 00520 g010
Cells 08 00520 g011 550
Figure 11. Expression profiles of pear RBOHs across different tissues. The transcript level of the 10 PbRBOHs in different tissues (fruits, leaves, stems, buds, and flowers). The expression data of each gene in fruit (at 15 DAF) of pear was used as a control sample to show the relative expression level. Each qRT-PCR analysis was performed in triplicate. The heatmap is drawn using the TBTools v. 0.6644 software (https://github.com/CJ-Chen/TBtools/releases).
Figure 11. Expression profiles of pear RBOHs across different tissues. The transcript level of the 10 PbRBOHs in different tissues (fruits, leaves, stems, buds, and flowers). The expression data of each gene in fruit (at 15 DAF) of pear was used as a control sample to show the relative expression level. Each qRT-PCR analysis was performed in triplicate. The heatmap is drawn using the TBTools v. 0.6644 software (https://github.com/CJ-Chen/TBtools/releases).
Cells 08 00520 g011
Cells 08 00520 g012 550
Figure 12. Subcellular localization of PbRBOHA-GFP and PbRBOHD-GFP. (A) Subcellular localization of PbRBOHA-GFP. (B) Subcellular localization of PbRBOHD-GFP.
Figure 12. Subcellular localization of PbRBOHA-GFP and PbRBOHD-GFP. (A) Subcellular localization of PbRBOHA-GFP. (B) Subcellular localization of PbRBOHD-GFP.
Cells 08 00520 g012
Cells 08 00520 g013 550
Figure 13. A model of reactive oxygen species (ROS) metabolism and lignin monomer polymerization in pear fruit.
Figure 13. A model of reactive oxygen species (ROS) metabolism and lignin monomer polymerization in pear fruit.
Cells 08 00520 g013
Table
Table 1. Sequence information of the RBOH family genes of five fruit-producing trees.
Table 1. Sequence information of the RBOH family genes of five fruit-producing trees.
SpeciesSequence IDNameChrProtein
Length (aa)MW
(kDa)		pISubcellular Localization Predicted
PearPbr018609.1PbRBOHA13962108.079.02Plas
Pbr037815.1PbRBOHB16964108.429.11Plas
Pbr038667.1PbRBOHC6947106.188.07Plas
Pbr006277.1PbRBOHD6978111.349.37Cyto
Pbr036006.1PbRBOHE2915103.568.59Cyto
Pbr003403.1PbRBOHF/914103.118.59Cyto
Pbr023445.1PbRBOHG1185697.658.66Cyto
Pbr033955.1PbRBOHH386097.568.79Cyto
Pbr007212.1PbRBOHI1477989.079.06Plas
Pbr037399.1PbRBOHJ1570179.099.29Chlo
PeachPrupe.1G211000PpRBOHA1971108.879.18Chlo
Prupe.5G107400PpRBOHB5964109.479.32Chlo
Prupe.5G138300PpRBOHC5941105.228.97Plas
Prupe.5G204900PpRBOHD581092.488.85Plas
Prupe.6G088800PpRBOHE685998.488.88Nucl
Prupe.6G321500PpRBOHF6893101.948.87Nucl
Prupe.7G193000PpRBOHG7906102.848.93Cyto
MeiPm000665PmRBOHA111613.744.58Cyto
Pm000754PmRBOHB186098.628.83Cyto
Pm003305PmRBOHC1893102.019.00Nucl
Pm007134PmRBOHD2982110.119.02Nucl
Pm023955PmRBOHE71067121.099.29Chlo
Pm024297PmRBOHF7940105.508.73Chlo
Pm024959PmRBOHG781092.328.81Plas
Pm027130PmRBOHH8914103.698.74Plas
GrapeGSVIVT01001122001VvRBOHA182794.139.16Chlo
GSVIVT01001123001VvRBOHB1906102.809.13Nucl
GSVIVT01014350001VvRBOHC1984095.406.47Plas
GSVIVT01015025001VvRBOHD11917103.819.26Chlo
GSVIVT01019429001VvRBOHE2922104.759.26Cyto
GSVIVT01025074001VvRBOHF685297.149.12Cyto
GSVIVT01031128001VvRBOHG1487399.449.34Cyto
Sweet orangeCs3g14240.1CsRBOHA3889101.049.14Cyto
Cs4g06920.1CsRBOHB4910103.349.19Chlo
Cs5g02940.1CsRBOHC5946107.509.37Chlo
Cs5g11890.1CsRBOHD577788.708.57Plas
Cs7g19320.2CsRBOHE7970109.698.78Nucl
Cs7g19380.1CsRBOHF730835.236.72Plas
Cs8g12000.1CsRBOHG8915103.529.11Plas
Cs8g17640.1CsRBOHH884296.619.11Plas
“/” indicates that the specific chromosome is not known. Chloroplast (Chlo); Cytoplasmic (Cyto); Nucl (Nucleus); Plasmodesmata (Plas).

Share and Cite

MDPI and ACS Style

Cheng, X.; Li, G.; Manzoor, M.A.; Wang, H.; Abdullah, M.; Su, X.; Zhang, J.; Jiang, T.; Jin, Q.; Cai, Y.; et al. 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. Cells 2019, 8, 520. https://doi.org/10.3390/cells8060520

AMA Style

Cheng X, Li G, Manzoor MA, Wang H, Abdullah M, Su X, Zhang J, Jiang T, Jin Q, Cai Y, et al. 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. Cells. 2019; 8(6):520. https://doi.org/10.3390/cells8060520

Chicago/Turabian Style

Cheng, Xi, Guohui Li, Muhammad Aamir Manzoor, Han Wang, Muhammad Abdullah, Xueqiang Su, Jingyun Zhang, Taoshan Jiang, Qing Jin, Yongping Cai, and et al. 2019. "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" Cells 8, no. 6: 520. https://doi.org/10.3390/cells8060520

APA Style

Cheng, X., Li, G., Manzoor, M. A., Wang, H., Abdullah, M., Su, X., Zhang, J., Jiang, T., Jin, Q., Cai, Y., & Lin, Y. (2019). 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. Cells, 8(6), 520. https://doi.org/10.3390/cells8060520

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop