Skip to main content
  • Research article
  • Open access
  • Published:

Transcriptomic analysis of early fruit development in Chinese white pear (Pyrus bretschneideri Rehd.) and functional identification of PbCCR1 in lignin biosynthesis

Abstract

Background

The content of stone cells and lignin is one of the key factors affecting the quality of pear fruit. In a previous study, we determined the developmental regularity of stone cells and lignin in ‘Dangshan Su’ pear fruit 15-145 days after pollination (DAP). However, the development of fruit stone cells and lignin before 15 DAP has not been heavily researched.

Results

In this study, we found that primordial stone cells began to appear at 7 DAP and that the fruit had formed a large number of stone cells at 15 DAP. Subsequently, transcriptome sequencing was performed on fruits at 0, 7, and 15 DAP and identified 3834 (0 vs. 7 DAP), 4049 (7 vs. 15 DAP) and 5763 (0 vs. 15 DAP) DEGs. During the 7-15 DAP period, a large number of key enzyme genes essential for lignin biosynthesis are gradually up-regulated, and their expression pattern is consistent with the accumulation of lignin in this period. Further analysis found that the biosynthesis of S-type lignin in ‘Dangshan Su’ pear does not depend on the catalytic activity of PbSAD but is primarily generated by the catalytic activity of caffeoyl-CoA through CCoAOMT, CCR, F5H, and CAD. We cloned PbCCR1, 2 and analysed their functions in Chinese white pear lignin biosynthesis. PbCCR1 and 2 have a degree of functional redundancy; both demonstrate the ability to participate in lignin biosynthesis. However, PbCCR1 may be the major gene for lignin biosynthesis, while PbCCR2 has little effect on lignin biosynthesis.

Conclusions

Our results revealed that ‘Dangshan Su’ pear began to form a large number of stone cells and produce lignin after 7 DAP and mainly accumulated materials from 0 to 7 DAP. PbCCR1 is mainly involved in the biosynthesis of lignin in ‘Dangshan Su’ pear and plays a positive role in lignin biosynthesis.

Background

Pear belongs to the family Rosaceae and is an important fruit crop. Because of its rich nutrition and good taste, pear fruit is popular with the public, and pear tree has a notably wide planting area around the world [1]. Pear has more than 2000 years of cultivation history in China [2]. ‘Dangshan Su’ pear (Pyrus bretschneideri cv. Dangshan Su) is a high-quality cultivar native to Daishan County (Dangshan County, Anhui Province, China). It has a long history of cultivation and is known for its unique flavour, rich nutrition and high medicinal value [3,4,5]. The ‘Dangshan Su’ pear has excellent fruit quality, high yield and storage resistance, and it is a late-maturing variety. However, ‘Dangshan Su’ pear has the drawback of having a relatively high number of stone cells, which strongly affect the quality of the fruit.

Stone cells are unique to pear fruit, and they influence the unique flavour and quality of pear fruits. Excessive stone cell content in pear fruit not only affects the accumulation of sugar and other nutrients but also impairs taste. Stone cells are thick-walled tissue cells composed mainly of lignin and cellulose, and those of pear are classified as short stone cells [6]. Numerous studies have shown that stone cells are derived from primordial cells. At 15 DAP, the cell wall of the parenchyma cells was unevenly thickened, while a large amount of lignin was transported to the thickened secondary wall and deposited to form stone cells [7]. After the stone cells are formed, they aggregate into a cluster to form a stone cell mass. The peak period of stone cell formation is 15-67 DAP, and then, their content gradually decreases [6, 8]. In the development of pear stone cells, the thickening of secondary walls and the accumulation of lignin are two critically important steps [9]. In pear stone cells, lignin content accounts for approximately 20-30% of the total dry weight [10, 11]. Therefore, a systematic study of the thickening of the secondary walls of parenchyma cells and the biosynthesis of lignin in pear fruit have important implications for the molecular mechanisms of pear stone cell formation.

Lignin is the most important component of stone cells, accounting for 18% of the stone cell content [12]. Lignin monomers can be divided into three different types: G-units (guaiacyl lignin), S-units (syringyl lignin) and H-units (hydroxyphenyl lignin). G-unit lignin is synthesized from coniferyl alcohol, S-unit lignin is polymerized from sinapyl alcohol, and H-unit lignin is synthesized from p-coumaryl alcohol [13]. Previous studies have found that the fruit lignin of ‘Dangshan Su’ pear consists of two components, G and S [14]. Research on lignin biosynthesis has been underway for more than a century, and the known lignin biosynthetic pathway has been continuously improved [15,16,17]. First, L-phenylalanine (L-Phe) enters the lignin biosynthetic pathway. Phenylalanine ammonia-lyase (PAL) catalyses the formation of phenylacrylic acid from L-Phe (non-oxidative deamination) [18]. P-coumaric acid is synthesized from phenylacrylic acid in a process catalysed by cinnamate 4-hydroxylase (C4H) and 4-hydroxycinnamate-CoA ligase (4CL) [19, 20]. Then, caffeic acid 3-O-methyltransferase (COMT), caffeoyl-CoA O-methyltransferase (CCoAOMT), cinnamoyl-CoA reductase (CCR), cinnamyl alcohol dehydrogenase (CAD) and other key enzymes are responsible for the formation of lignin monomers [21,22,23,24]. The lignin monomers are gradually polymerized under the action of peroxidase (POD) and laccase (LAC) to form lignin complexes [25, 26]. The lignin produced by these complex catalytic reactions is deposited on the thicker secondary wall to form stone cells.

Transcriptome analysis can comprehensively examine the function and structure of genes and reveal the molecular mechanisms of specific biological processes [27]. Transcriptome analysis has been widely used in the study of various metabolic pathways in plants. For example, Scutellaria viscidula was subjected to transcriptome analysis to analyse the biosynthetic mechanism of flavonoids [28], and the differential genetic pathways and pathways related to drought tolerance in Pennisetum glaucum were explored [29]. In pear, Zhang et al. revealed the molecular mechanism of the difference in stone cells between ‘Dangshan Su’ pear and ‘Lianglizaosu’ by comparative transcriptome analysis [5]. However, this previous study of lignin in pear fruit examined only the 23-145 DAP period, during which a large number of stone cells are produced. Therefore, we here attempt to advance the time point of study to find the critical period when stone cells first begin to form in large quantities. This initial development period will allow more effective exploration of the molecular mechanisms of pear stone cell formation and lignin biosynthesis.

In this study, we used fruits obtained at 0, 7, and 15 DAP as experimental materials to investigate the mechanism of stone cell formation in the early stage of pear fruit development. First, the stone cell and lignin contents of fruits in three periods were compared. Subsequently, transcriptome sequencing was performed to analyse differential gene expression in the fruits of these three periods. Finally, the functions of two key enzymes (PbCCR1 and 2) were characterized to clarify their roles in lignin biosynthesis. Our results reveal the molecular mechanism of stone cell formation and development in pear fruit, providing a theoretical basis for ‘Dangshan Su’ pear molecular genetic breeding and variety improvement.

Results

Analysis of stone cells and lignin content in fruits

For the fresh samples obtained (0, 7, 15 DAP fruits), we first stained the samples with phloroglucinol and then performed sectioning. After staining with phloroglucinol, we found that the results of tissue sectioning showed significant differences. As shown in the Fig. 1a, the 0 DAP fruit was not dyed red by the phloroglucinol dye but was red only in the vascular bundle area. The staining results of the 7 DAP fruit were similar to those of the 0 DAP fruit, but we found stone cell primordia stained red near the core of the fruit. This finding indicates that a few stone cell primordia are present in this period, and the development of stone cells begins to enter a period of relatively rapid growth. In the 15 DAP fruit, we found a large number of stone cell clusters dyed red by phloroglucinol near the core and pulp, indicating that a large number of stone cells were formed between 7 and 15 DAP and aggregated to form stone cell clusters. Then, we measured the stone cell content and lignin content of the fruits in the three periods, and the measured results were consistent with the staining results (Fig. 1b, c). Stone cells were scarce in the 0 DAP samples, and notably few stone cells were identified in the 7 DAP samples, while the content of stone cells increased greatly in the 15 DAP samples. The change trend in lignin content was consistent with the change in stone cell content, which was also a gradually increasing process.

Fig. 1
figure 1

Investigation of stone cell and lignin contents in 0, 7, 15 DAP fruits. a Fruit morphology and stone cell phloroglucinol staining of 0, 7, and 15 DAP fruits were performed. b Stone cell contents of 0, 7, 15 DAP fruits. c Lignin contents of 0, 7, and 15 DAP fruits. **Significant differences between the stone cell or lignin levels of the different developmental periods (P < 0.01)

Overview of RNA-Seq results

According to our research results, stone cells and lignin undergo considerable changes among the three time points 0, 7, 15 DAP, and the contents at 15 DAP exhibit a significant increase compared with the previous two periods. To clarify the transcriptional changes in ‘Dangshan Su’ fruit in the early stages of stone cell development, fruits from 0, 7, 15 DAP were selected for comparative analysis by Illumina sequencing, and three biological replicates were evaluated for each time point. In this study, we extracted total RNA from all nine samples (Additional file 9: Figure S1) and then constructed cDNA transcriptome sequencing library. Primary component analysis (PCA) and duplicate correlation checking scatter plot analysis showed high correlations among the three biologically replicated samples in each group (Additional file 10: Figure S2, Additional file 11: Figure S3). After trimming the raw data (deleting adaptors and low-quality sequences), 0 DAP fruit obtained the most reads (50313188, 53291526, 51885204), followed by 15 DAP fruit (45576288, 40812768, 45463438), 7 DAP fruit obtained the least reads (35573064, 31374310, 38135470) (Table 1). In addition, We found that the GC content of 9 libraries was in 48.27-49.11%, while the Q20 and Q30 values were greater than 96.37 and 91.75%, respectively. Because of the high quality of the sequencing results, most reads could be mapped to the reference genome of ‘Dangshan Su’. The percentages of mapped reads for the libraries were highly similar (about 80%), and the ratios of multiply mapped reads were 8.35-9.77%. The reads mapped to ‘+’ and reads mapped to ‘-’ showed greater than 35% coverage at all samples. The statistical results of the genome structure distribution showed that each sample had a very high similarity (Additional file 12: Figure S4). In single-site mapping, the intergenic region was the largest component. Of the total mapped reads, 56.36-58.80% of reads were mapped to intergenic regions, and 40.13-42.60% of reads were mapped to exon regions. The intron region had the smallest percentage, with only approximately 0.99-1.06% of reads being mapped to this region. Taken together, these results indicated that the RNA-Seq sequencing results were reliable and could be further analysed.

Table 1 Summary of the RNA-Seq data in three stages of ‘Dangshan Su’ fruit

Gene expression analysis and cluster analysis of DEGs

According to the relationship between gene transcripts per million (TPM) deviation within 15% of the final value (based on 100% mapped reads) and percentage of mapped reads, when all samples were 3.5 < TPM < 15, the curve reached saturation state (Additional file 13: Figure S5). This result showed that our sequencing depth was good and could be used for gene expression analysis. By analysing the density distributions of gene expression, we gained a macroscopic understanding of the gene expression levels in these samples (Additional file 14: Figure S6). In the 0 DAP samples, the number of genes with low expression levels (log2(TPM) < − 1) was higher than those in the other two samples (7 DAP and 15 DAP), but the number of genes with high expression levels (log2(TPM) > 1) was significantly lower than those of the other two samples. However, the number of highly expressed genes at 7 DAP was less than that at 15 DAP. Therefore, low-expression genes were the most abundant at 0 DAP, while high-expression genes were significantly enriched at 15 DAP. Similar results were obtained in the hierarchical clustering heat map of transcript expression (Additional file 15: Figure S7). In the heat map, we can clearly observe that the transcripts of the three parallel samples in each group showed largely the same expression levels. However, the transcript expression levels changed with fruit development. Among the 0 DAP samples, a large number of transcripts (probably 60-70%) showed low expression levels, which represented the greatest proportion among the three groups. At 7 DAP compared with 0 DAP, the number of transcripts with low expression and high expression levels was not notably different, but more transcripts were observed at the medium expression level. Clearly, 15 DAP was the group with the highest expression of transcripts, and approximately 60% of its transcripts exhibited high expression levels.

After analysing the expression patterns of each group of transcripts, we summarized the numbers of DEGs in the three stages of fruit development. In all pairwise comparisons, a total of 5763 DEGs were identified in 0 DAP vs. 15 DAP, making this the most abundant group (3521 up-regulated and 2242 down-regulated). We identified 3834 DEGs in the 0 DAP vs. 7 DAP group (2289 up-regulated and 1545 down-regulated). A total of 4049 DEGs (2278 up-regulated and 1771 down-regulated) were discovered in the 7 DAP vs. 15 DAP group (Additional file 16: Figure S8). Venn diagrams were used to show the DEGs among all pairwise comparisons (Fig. 2a, b). In all pairwise comparisons, a total of 354 DEGs were up-regulated, and 125 DEGs were down-regulated. However, we found some genes that were specifically differentially expressed. For example, we identified 1593 distinctive DEGs in the group between the 0 DAP and 7 DAP stages (888 up-regulated and 705 down-regulated), while we found only 1397 (586 up-regulated and 811 down-regulated) and 1349 (782 up-regulated and 567 down-regulated) distinctive DEGs in the 7 DAP vs. 15 DAP and 0 DAP vs. 15 DAP groups, respectively. The greatest number of down-regulated genes was identified in the 0 DAP vs. 15 DAP group (2242), and the proportion of strongly down-regulated DEGs (with high expression fold changes) was also the highest among all the pairwise comparisons (Fig. 2c). Interestingly, in this group, we found that although it contained the greatest number of up-regulated genes, the proportion of strongly up-regulated DEGs was not the highest. Instead, the proportion of strongly up-regulated DEGs was the highest of the three groups in the 7 DAP vs. 15 DAP group.

Fig. 2
figure 2

Venn diagram of significantly up-regulated and down-regulated genes in 0, 7, 15 DAP fruits and cluster analysis of DEG expression based on log ratio q-value data. a, b Venn diagram analysis of stone cell development-related DEGs. c Clustering analysis of DEGs based on fold change

Functional annotation and pathway enrichment analysis of DEGs

To explore the molecular mechanism of stone cell development in ‘Dangshan Su’ pear, we used the EuKaryotic Orthologous Groups (KOG), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases to analyse DEGs. The results of the KOG enrichment analysis showed that the DEGs were mainly annotated in general function prediction only (R), posttranslational modification (O), signal transduction mechanisms (T), carbohydrate transport and metabolism (G), and amino acid transport and metabolism (E). Twenty-five KOG functional categories were represented (Additional file 17: Figure S9). However, we did not find any DEGs that were annotated in the cell motility (N) classes. In the three groups of pairwise comparisons, 157 (8.4%), 173 (8.4), and 211 (7.2%) DEGs were annotated to involve the function of secondary metabolite biosynthesis (Q class).

To further understand the functions of the DEGs, we also performed GO classification analysis (Additional file 18: Figure S10). Specifically, for 0 DAP vs. 7 DAP, 7 DAP vs. 15 DAP and 0 DAP vs. 15 DAP, the three groups had 3061 (79.8), 3290 (81.3%) and 4667 (80.1%) DEGs, respectively, annotated in the GO database. These DEGs were enriched in the three major functional categories of biological process (BP), cellular component (CC) and molecular function (MF). The above three functional categories could be divided into 67 more detailed subcategories, including 27 (BP1-27), 22 (CC1-22), and 18 (CC1-18) functional subcategories. Cellular process (BP-22), cell (CC-16), and binding (MF-15) were the most important functional subcategories of the three major functional categories in addition to these large numbers of functional subcategories. In all pairwise comparisons, significantly higher amounts of the up-regulated DEGs were annotated in single-organism process (BP-1), metabolic process (BP-8), membrane (CC-13), organelle (CC-17), catalytic activity (MF-4), and nucleic acid binding transcription factor activity (MF-1); the down-regulated DEGs were annotated in biological regulation (BP-19), response to stimulus (BP-4), membrane part (CC-11), organelle part (CC-18), transporter activity (MF-3), and signal transducer activity (MF-10).

Mapping DEGs to the KEGG database allows for the identification of significant metabolic or signalling pathways in which DEGs are enriched. The 0 DAP vs. 15 DAP group had 660 DEGs annotated to 177 KEGG pathways (Additional file 19: Figure S11, Additional file 2: Table S2). These KEGG pathways involved four major functional groups: cellular processes, environmental information processing, genetic information processing, and metabolism, and these functional groups were divided into four, three, four, and twelve functional subgroups, respectively (Additional file 3: Table S3). In the 0 DAP vs. 15 DAP group, the up-regulated DEGs were mainly concentrated in pathways such as carbon metabolism (ko01200), pyrimidine metabolism (ko00240), and biosynthesis of amino acids (ko01230); the down-regulated DEGs were concentrated in purine metabolism (ko00230), retinol metabolism (ko00830), glycerophospholipid metabolism (ko00564) and other KEGG pathways (Fig. 3). In addition, we found that some DEGs were enriched in the pathways of apoptosis (ko04210), phenylalanine metabolism (ko00360, ko00400) and phenylpropanoid biosynthesis (ko00940) (Fig. 3, Additional file 2: Table S2).

Fig. 3
figure 3

KEGG enriched scatter plots of up-regulated (left) and down-regulated (right) DEGs in the 0 DAP vs. 15 DAP comparison. Rich factor refers to the ratio of the number of DEGs enriched in the pathway to the number of annotated genes. The larger the Rich factor, the more significant the enrichment is. The q-value is the corrected P-value after multiple hypothesis testing, which ranges from 0 to 1. The closer to zero, the more significant the enrichment is

DEGs associated with stone cell lignin biosynthesis

The phenylpropanoid metabolism pathway is usually closely related to the biosynthesis of secondary walls and lignin [30]. RNA-seq data analysis showed that active transcription and expression were detected for some key enzyme genes in the phenylpropanoid metabolism pathway that control the entry of lignin biosynthesis carbon sources, including PAL (3), C4H (3), 4CL (2), HCT (5), CCoAOMT (1), F5H (2), CCR (3), CAD (3), POD (3) and other key enzyme genes in lignin biosynthesis (Additional file 20: Figure S12, Table 2). Coumarin and lignin are closely related in the phenylpropanoid metabolism pathway with a common precursor substance (coumaric acid). Interestingly, we did not identify the presence of DEGs in the coumarin biosynthetic pathway (Additional file 20: Figure S12).

Table 2 Expression level of key enzymes in lignin biosynthesis pathway of ‘Dangshan Su’ pear

Based on the results of transcriptome sequencing, we calculated the relative expression levels of these key lignin biosynthetic enzyme genes. In the 0 DAP vs. 7 DAP group, most of the key lignin biosynthetic enzyme genes had low expression levels. However, some key enzymes also showed up-regulated expression trends, such as PbPAL1, Pb4CL1, 2, PbCCR1, 2, PbCOMT1, 3, PbCAD3 and PbPOD2. It is worth mentioning that PbPOD2 is the key enzyme gene in this stage that is related to the polymerization of lignin monomers, and it showed up-regulated expression (fragments per kilobase of exon per million mapped reads (TPM) value as high as 1455.64 at 7 DAP). At the 7-15 DAP stage, at least one member of each key lignin biosynthetic enzyme gene showed an upregulation trend (Table 2).

We then performed qRT-PCR analysis on these key enzyme genes. We found that the results were consistent with the RNA-seq sequencing results and mapped the results to the lignin biosynthetic metabolic pathway (Fig. 4). L-Phe, as a precursor of lignin biosynthesis, flows into the metabolic pathway. The expression of PbPAL1 was up-regulated at first, and PbPAL2 rapidly increased in expression level after 7 DAP. These genes, including PbC4H1, 2, PbC3H, PbCCoAOMT, PbF5H1, 2 and PbCAD2, exhibited similar expression patterns. These genes were expressed at lower levels before 7 DAP, but their expression after 7 DAP significantly and continuously increased, reaching the range of up-regulated differential expression. In addition, some genes began to increase from 0 DAP, such as Pb4CL1, 2, PbCCR1, 2, and PbCOMT1, 3. In contrast, genes such as PbCCR3 and PbCAD1 showed extremely low expression levels during the whole period from 0 to 15 DAP. Finally, of the genes responsible for the polymerization of lignin monomers, only PbLAC1 and PbPOD2 showed upregulation of transcript levels. Overall, the expression levels of most key enzyme genes in lignin biosynthesis at the 0-7 DAP stage were low, and their expression levels continued to increase significantly from 7 DAP.

Fig. 4
figure 4

Expression patterns of key enzymes in lignin biosynthesis by qRT-PCR

To further explore candidate genes related to stone cell development in pear fruits, 354 up-regulated DEGs were identified as common among the three groups. These DEGs include Dirigent-like protein (DIR, Pbr000712.1), KIP1-like protein (KIP1, Pbr001414.1), S-adenosylmethionine synthetase (S-AdoMet, Pbr001686.1), Glutamine amidotransferase (GATase, Pbr002524.1), GMC oxidoreductase (GMC oxred, Pbr011059.2), EamA-like transporter protein (EamA, Pbr024589.1), Terpene synthase (TPS, Pbr024829.1), Multicopper oxidase (Cu-oxidase, Pbr025074.1) and other genes. Both the transcriptome sequencing data and qRT-PCR results showed that these genes were up-regulated from 0 to 15 DAP (Additional file 21: Figure S13, Additional file 4: Table S4). This finding suggests that these genes may play an important role in the development and biosynthesis of ‘Dangshan Su’ pear stone cells and lignin.

Numerous studies have shown that transcription factors regulate the formation of plant cell walls and the biosynthesis of lignin [31,32,33]. Thus, we explored the expression patterns of three transcription factor families (MYB, AP2/ERF and NAC) (Additional file 22: Figure S14). In the MYB transcription factor family [34], we found 16 MYB gene family members that were gradually up-regulated and 8 that were gradually down-regulated. In addition, four MYB genes (PbMYB47, 73, 93, 125) showed a pattern of expression that decreased first and then increased. Fourteen members of the AP2/ERF gene family [35] showed up-regulated expression, whereas only four genes showed a progressively down-regulated expression pattern. The expression of the NAC gene family [36] seemed to be relatively silent from 0 to 15 DAP, while the expression of only 7 family members was gradually up-regulated.

Cinnamoyl-CoA reductase involved in lignin biosynthesis

The catalytic formation of cinnamoyl-CoA into cinnamaldehyde is a key step before the formation of lignin monomers. The CCR gene plays a major role in this critical step [37]. According to previous research results, 31 CCR genes were retained in the genome of Chinese white pear, and these CCR genes have conserved domains in this family [24]. We validated the accuracy of our transcriptome sequencing data by qRT-PCR experiments. After that step, we collected the FPKM values of all 31 CCR genes at 0, 7, 15 DAP and used the data to analyse the expression profiles of these genes (Fig. 5a). Approximately 74.2% of the CCR genes showed very low transcript levels at the 0-15 DAP stage. The expression patterns of PbCCR-like7, 11 were similar; both gene expression levels were high at 0 DAP and then decreased gradually. The expression levels of PbCCR-like12, 22 increased first and then decreased, although they differed in that PbCCR-like22 was highly expressed at 7 DAP. In contrast, the expression of PbCCR-like9 decreased first and then increased, reaching a low level of expression at 7 DAP. According to our research results, the lignin content of pear fruit increased gradually during the 0-15 DAP period. Previous studies have speculated that PbCCR1, 2, 3 may be related to lignin biosynthesis in pear fruit. However, the expression profiles and qRT-PCR results showed that only PbCCR1, 2 gradually increased in the 0-15 DAP period, while the expression of PbCCR3 was not detected at this stage.

Fig. 5
figure 5

Expression profiles of 31 CCR genes in ‘Dangshan Su’ pear and sequence alignment of PbCCR1, 2, 3 against genes from other plant species. a Analysis of the expression patterns of 31 CCR genes. b Using the crystal structure of MtCCR as a template, the blue underline indicates the NADP specificity site, and the green underline indicates the NAD(P) binding site. Black circles indicate catalytic triads. Orange triangles represent key catalytically active amino acid residues. The purple underline indicates the catalytic centre. The black wavy lines and arrows represent α-helices and β-sheets, respectively. ZmCCR1 (CAA74071); OsCCR19 and 20 (Os09g25150, Os08g34280); AtCCR1 and 2 (AAG46037, AAG53687); PtoCCR7 (KF145198)

To further understand the characteristics of PbCCR1, 2, 3, we collected and aligned some CCR genes that have been proven to be related to lignin biosynthesis (Fig. 5b). The obtained protein sequence contains conserved G-X-X-G-X-X-A and D-X-X-D NAD(P) binding sites. At the same time, the presence of the R-X-X-X-X-X-K NADP specificity motif was identified; this particular sequence distinguishes the CCR gene from other NAD(H)-dependent short-chain dehydrogenase superfamily (SDR) genes. More importantly, the NWYCYGK motif was also identified. This motif is a CCR signature motif that is highly conserved among these protein sequences. The Ser-Tyr-Lys catalytic triad, which represents active residues in SDRs, was found and marked with black circles. In addition, yellow triangles were used to identify key amino acid residues closely related to catalytic reactions. These key catalytic sites are also conserved, but there are exceptions. For example, the amino acid residues at position 43 associated with substrate and cofactor binding showed lower conservation, but the three genes in Chinese white pears showed strong species specificity (all are Leu at this location). Additionally, at the substrate binding site, the amino acid residue at position 155 was highly conserved, and only in OsCCR20 was the tyrosine (Y) at this site replaced by phenylalanine (F) [38].

Heterologous overexpression of PbCCR1 and 2 in Arabidopsis

We attempted to further explore the role of PbCCR1 and 2 in lignin biosynthesis in pear fruit. First, we obtained the full-length sequences of these two genes. We ligated these two genes into the pCAMBIA1304 vector and obtained complete plant expression vectors, and we subsequently transformed wild-type Arabidopsis (Fig. 6a). Using RT-PCR for detection (Fig. 6b), we finally obtained 4 (PbCCR1-OE1-4) and 4 (PbCCR2-OE1-4) active transgenic lines. We detected the expression levels of PbCCR1 and 2 in T3 generation transgenic plants. Based on the transcript levels (Fig. 6c, d), we found that the gusA gene had higher expression in PbCCR1-OE2, 4 and PbCCR2-OE2, 4.

Fig. 6
figure 6

Overexpression of PbCCR1 and PbCCR2 in Arabidopsis. a Construction strategy of the plant expression vectors pCAMBIA1304-PbCCR1 and 2. b RT-PCR analysis for detecting the GFP gene in transgenic Arabidopsis. M: DL2000 DNA Marker; 1: negative control; 2: positive control; 3-6: PbCCR1-OE1-4 transgenic lines; 7-10: PbCCR2-OE1-4 transgenic lines. c The expression levels of PbCCR1 and 2 in transgenic lines. Error bars represent the SD of three biological replicates. d β-Glucuronidase (GUS) histochemical staining of PbCCR1-OE2 and 4 and PbCCR2-OE2 and 4

To understand the distribution of lignin in the inflorescence stems of wild-type and transgenic Arabidopsis more intuitively, we chose PbCCR1-OE2 and PbCCR2-OE1, which had higher lignin contents, and performed histochemical (toluidine blue) staining (Fig. 7a–f). We compared the stained cross-sections of the inflorescence stems of the three groups of samples. The staining intensities of the xylem and the interfascicular fibres of PbCCR1-OE2 were the highest, followed by those of PbCCR2-OE1, and those of wild-type plants were the lowest. In addition, we found that compared with wild-type plants, PbCCR1-OE2 showed an increase in the number of cells in the xylem region, while there was no significant difference between wild-type plants and PbCCR2-OE1. Subsequently, we observed the cell walls of the xylem in the cross-section of the inflorescence stem of Arabidopsis (wild type, PbCCR1-OE2 and PbCCR2-OE1) by transmission electron microscope (TEM) (Fig. 8a–c). The cell wall thickness of the xylem in the inflorescence stem cross-sections of PbCCR1-OE2 increased significantly compared with that of wild-type plants. In PbCCR2-OE1, the cell wall thickness of these cells also seemed to increase, but the change range was not as obvious as that observed in PbCCR1-OE2 plants (Fig. 8d).

Fig. 7
figure 7

Observation of transverse section staining of inflorescence stems and determination of lignin contents in transgenic plants. a–f Toluidine blue staining of the inflorescence stems from wild-type (WT) and transgenic lines. a, c, e: bar = 100 μm. b, d, f: bar = 200 μm. F: interfascicular fibre; X: xylem. g–h Determination of lignin contents in inflorescence stems of Arabidopsis. Significant differences in lignin content levels between WT and transgenic overexpression lines of Arabidopsis were determined by a t test; *P < 0.05, **P < 0.01

Fig. 8
figure 8

Micrographs of transmission electron microscopy on the cell wall of xylem of inflorescence stem of wild type and transgenic Arabidopsis (the scale = 2.0 μm). a Wild type; b PbCCR1-OE2; c PbCCR2-OE1; d Measurement of cell wall thickness (P < 0.05)

We measured the lignin content in the inflorescence stems of Arabidopsis, which supported our staining results (Fig. 7g, h). The lignin contents of PbCCR1-OE2 and PbCCR1-OE4 reached 14.20 and 13.92%, respectively, which were significantly higher than those of wild-type plants (10.02%). However, although the lignin content (11.78%) of PbCCR2-OE1 was significantly higher than that of wild-type plants, and the lignin contents of the other lines were higher than that of wild-type plants, the difference in the other lines was not significant. In transgenic plants, the expression levels of key enzyme genes involved in Arabidopsis lignin biosynthesis changed greatly (Fig. 7i). Overexpression of PbCCR1 in wild-type Arabidopsis increased the expression of most key lignin synthesis enzymes, especially At4CL, AtF5H, AtCOMT, AtCAD4 and other genes. In contrast, the expression of AtCCOMT and AtCCR1 was inhibited. However, in PbCCR2-OE1, the expression of many key lignin biosynthetic enzyme genes was also activated, but this activation was small. Similar to the results of PbCCR1-OE2, the expression of AtCCOMT and AtCCR1 was inhibited.

Discussion

The ‘Dangshan Su’ pear has a long history of cultivation and represents an important germplasm resources that occupies an important position in the world fruit market. However, in recent years, due to the deterioration of varieties and poor management, fruit quality has declined to some extent. The content of stone cells in pear fruit has increased, causing a less appealing taste. To solve this problem, many studies have been conducted on stone cells, and some progress has been made in revealing the biosynthetic and regulatory processes of pear stone cells and lignin [6, 39, 40]. However, the existing research has mainly concentrated on pear fruit after 15 DAP, and pear fruit in this period already contains a large number of stone cells [5, 41, 42]. The molecular mechanisms of stone cell formation and development in early fruit development (before 15 DAP) remain unclear. In this study, we selected fruit samples from the early developmental stages (0, 7, 15 DAP) of ‘Dangshan Su’ pear. We found that no stone cells were present near the fruit core at 0 DAP, indicating that the fruit did not form stone cell primordia during this period, and the metabolic pathways associated with stone cell formation and development may not be fully activated at this stage. At 7 DAP, sporadic stone cell primordia appeared near the core. At 15 DAP, the content of stone cells increased significantly (lignin content increased from 0.01 to 0.34%), and the metabolic network related to stone cell formation was fully activated (Fig. 1).

We then selected fruits obtained at 0, 7, 15 DAP for transcriptome sequencing to ascertain the molecular basis for the formation of large numbers of stone cells and lignin deposition at 0-15 DAP. We obtained a total of 35.2 Gb of data. For the 0, 7, 15 DAP fruit, we obtained at least 40506401, 25200858 and 32775750 mapped reads (Table 1), greater numbers than in the previous study of 23, 55 and 145 DAP fruit [5]. Because stone cells undergo a secondary wall thickening process accompanied by lignin deposition [43], we focused on cell wall development and lignin biosynthesis. The GO database contains three major categories, CC, MF, and BP, which can be applied to individual species to annotate the functions of genes [44]. For the 7-15 DAP comparison, in the BP classification (Additional file 18: Figure S10), a large number of DEGs were annotated to plant-type secondary cell wall biogenesis (44), phenylpropanoid biosynthetic process (44), lignin metabolic process (36), lignin biosynthetic process (32) and other functional categories. In the MF classification, many DEGs were annotated to GO nodes associated with plant cellulose biosynthesis, such as cellulose synthase activity (10) and cellulose synthase (UDP-forming) activity (10). All of the above are functional categories associated with cell wall development and lignin biosynthesis. However, in the 0-7 DAP period, where stone cells did not increase significantly, DEGs also fell into these functional categories, but the number of DEGs was significantly less than in the 7-15 DAP period. This finding indicated that the fruit began to produce stone cells on a large scale after 7 DAP.

The development of stone cells depends on the biosynthesis of lignin, which involves multiple metabolic pathways. The most important of these is the phenylpropanoid metabolic pathway (ko00360 and ko00940) [30]. Thirty-four key lignin biosynthetic enzyme genes were identified, and most of the genes showed an upward trend, such as PbPAL1, Pb4CL1, 2, PbHCT50, PbCCR1, 2, PbCOMT1, 3, PbPOD2 (Fig. 4, Table 2). These results showed that lignin biosynthesis entered an active stage after 7 DAP, and a large number of lignin monomers were synthesized. However, we found that there was no differential expression of PbSAD between 0 and 15 DAP, and this result was confirmed by qRT-PCR. SAD is involved in the late stage of lignin biosynthesis and can transform sinapaldehyde into sinapyl alcohol to produce syringyl monolignol [45]. In ‘Dangshan Su’ pear fruits at 23, 55 and 145 DAP, the expression of PbSAD was consistently very low. Perhaps the biosynthesis of syringyl monolignol in ‘Dangshan Su’ pear does not depend on the pathway of sinapyl alcohol catalysed by 4CL, CCR and SAD. The same results were also found in other species. For example, CAD may be involved in the synthesis of coniferyl alcohol and sinapyl alcohol in Arabidopsis and does not require any specific SAD to produce sinapyl alcohol [46]. A rice study also showed no evidence that a specific SAD is needed to produce sinapyl alcohol [47]. Therefore, we believe that the formation of syringyl monolignol in ‘Dangshan Su’ pear is mainly generated by caffeoyl-CoA catalysed by CCoAOMT, CCR, F5H and CAD. Among the genes involved in the polymerization of lignin monomers, the expression level of PbLAC1 gradually increased. This gene has been shown to be involved in the biosynthesis of lignin and thickening of the cell wall [26]. In PbPOD1, 2, 3, only the expression level of PbPOD2 (Pbr031894.1) increased as the lignin content increased. Therefore, we speculate that PbLAC1 and PbPOD2 are mainly responsible for the polymerization of lignin monomers in the fruit of ‘Dangshan Su’ pear. In addition to the phenylpropanoid pathway, we found that many metabolic pathways involved in lignin biosynthesis were up-regulated at 0-7 DAP. Many key enzyme genes (carbamoyl phosphate synthase, CPSase, Pbr026918.3; histidine phosphatase, His Phos, Pbr011138.1; D-isomer specific 2-hydroxyacid dehydrogenase, 2-Hacid dh, Pbr025526.1) in the carbon metabolic pathway (ko01200) are up-regulated during this period (Additional file 2: Table S2). The carbon metabolic pathway is activated to provide an important carbon source for downstream lignin biosynthesis [48]. At the same time, there are a number of up-regulated DEGs in the citrate cycle pathway (ko00020). Through the citrate cycle, sugars and other organic substances can be synthesized. These substances then enter glycolysis to form L-Phe [49]. This finding indicates that from 0 to 7 DAP, the fruit initiates upstream pathways to prepare for lignin biosynthesis, producing a large quantity of lignin biosynthetic precursors and accumulating materials for the large-scale synthesis of lignin after 7 DAP.

After 7 DAP, various pathways associated with lignin biosynthesis are activated, and the fruit begins to synthesize lignin in large quantities. At the same time, some metabolic pathways involved in the formation of plant cell walls also become active, for example, pentose and glucuronate interconversions (alcohol dehydrogenase GroES-like protein, ADH, Pbr013913.1; UDP-glucoronosyl and UDP-glucosyl transferase, UDPGT, Pbr009446.1; pentose phosphate pathway phosphofructokinase, PFK, Pbr008761.1) (Additional file 2: Table S2). The intermediates of these metabolic pathways are closely related to the biosynthesis of plant cell walls. In addition, we collected cellulose synthase genes from Arabidopsis thaliana and identified homologues of these genes in ‘Dangshan Su’ pear (Additional file 5: Table S5). According to previous studies, AtCESA1, 3, 6 are mainly involved in the biosynthesis of the plant primary wall [50,51,52], while AtCESA4, 7, 8 play an important role in the biosynthesis of plant secondary cell walls (SCWs) [53, 54]. We found that the homologous genes of AtCESA4, 7, 8 (PbCESA4-1, 7, 8) in ‘Dangshan Su’ pear were up-regulated (Additional file 5: Table S5). To further verify that these three genes are involved in lignin biosynthesis in fruit stone cells, we explored the expression patterns of PbCESA4-1, 7, 8 during the whole developmental period of ‘Dangshan Su’ pear (Additional file 23: Figure S15). The results showed that the expression patterns of these three genes were consistent with the developmental trends of lignin and stone cells in ‘Dangshan Su’ pear. This finding suggests that the three genes PbCESA4-1, 7, 8 may be involved in the development of stone cells in pear fruit and may be responsible for SCW formation.

We investigated the expression patterns of MYB, NAC and AP2/ERF transcription factors in ‘Dangshan Su’ pear from 0 to 15 DAP to identify genes closely related to stone cell formation (Additional file 22: Figure S14). Arabidopsis NST1, 2, 3 have been shown to regulate SCW formation in plants [32, 55]. We identified the NAC gene most closely related to AtNST1, 2, 3 in ‘Dangshan Su’ pear and found that only the AtNST1 homologous gene (PbNAC46) showed up-regulated expression at 0-15 DAP (Additional file 6: Table S6). Thus, we believe that PbNAC46 is the most important NAC gene involved in SCW formation in ‘Dangshan Su’ pear stone cells. MYB transcription factors can induce, catalyse and hinder the biosynthesis of plant lignin [56]. In this study, we found that multiple MYB genes may be involved in the development of pear stone cells (Additional file 6: Table S6). We found that the homologous genes of AtMYB46 and 83 (functional redundancy, positive regulation of secondary wall, cellulose, lignin biosynthesis) in ‘Dangshan Su’ pear (PbMYB105, 104) were up-regulated, suggesting that there may be a regulatory network similar to that of Arabidopsis in ‘Dangshan Su’ pear. This regulatory mechanism is achieved through the regulation of downstream genes by PbMYB105 and 104 [57]. In addition, we identified several MYB genes that may be involved in the regulation of lignin biosynthesis (including the activation or inhibition of lignin biosynthesis), such as PbMYB103, 122, 127 (Additional file 6: Table S6). Compared with NAC and MYB transcription factors, AP2/ERF genes often regulate lignin biosynthesis by directly regulating the expression of key enzyme genes in lignin biosynthesis [58,59,60]. We found that the expression of PbSHN (Pbr027839.1) is consistent with the accumulation of lignin content, and we believe this gene may have similar functions to AtSHN in regulating the biosynthesis of lignin [33]. PbERF110 has very close affinities with some ERF genes that regulate lignin biosynthesis (PtERF34, DcERF2, OsERF71) [60,61,62]. We believe that PbERF110 is a potential AP2/ERF gene that positively regulates lignin biosynthesis in ‘Dangshan Su’ pear. Conversely, PbAP2-24 was discovered as a possible negative regulation of lignin biosynthetic genes. We found that the similarity between the amino acid sequences of PbAP2-24 and EjAP2-1 (a transcription factor that inhibits lignin biosynthesis) [63] was as high as 96.75% (Additional file 24: Figure S16). Both the transcriptome sequencing and qRT-PCR results showed that PbAP2-24 was consistently low from 0 to 15 DAP. This finding indicates that PbAP2-24 did not exert its biological function of inhibiting lignin synthesis during the activation period of stone cell development and lignin synthesis.

CCR, as a key enzyme gene regulating lignin monomer biosynthesis, plays an important role in the lignin biosynthetic pathway. There has been controversy regarding the types of substrates used by CCR genes in catalytic reactions, but numerous studies have shown that CCR is more likely to use feruloyl-CoA or caffeoyl-CoA rather than sinapoyl-CoA as a substrate in catalytic reactions [64]. In this study, we do not believe that the biosynthesis of S-type lignin in ‘Dangshan Su’ pear depends on the catalytic activity of PbSAD. The reaction of sinapoyl-CoA catalysed by PbCCR is located upstream of the reaction, which indicates that PbCCR in ‘Dangshan Su’ pear probably does not use sinapoyl-CoA as the substrate for catalytic reaction (Fig. 4). Thirty-one CCR genes were retained in the ‘Dangshan Su’ pear genome, and PbCCR1, 2, and 3 were likely to participate in lignin biosynthesis [24]. According to the results of the study of expression profiles, only PbCCR1 and 2 showed up-regulated expression patterns, and the expression level of PbCCR3 in this period was very low (Fig. 5a). We believe that PbCCR1 and 2 may be closely related to early stone cell development and lignin biosynthesis in pear fruit, while PbCCR3 has little effect on stone cell development and lignin biosynthesis in the early fruit development stage. To further clarify the roles of PbCCR1 and 2 in lignin biosynthesis, we overexpressed these two genes in wild Arabidopsis (Fig. 6). The functions of PbCCR1 and 2 were verified by overexpression analysis. We found that overexpression of PbCCR1 increased lignin content in transgenic plants and cell wall thickness in xylem tissues and interfascicular fibres (Fig. 8). Similar results have been reported in birch, where overexpression of BpCCR1 increased the lignin content of transgenic plants to 14.6% and increased the cell wall thickness of xylem vessel cells [65]. However, overexpression of PbCCR2 did not play a significant role in these indicators, perhaps because PbCCR2 showed low catalytic activity, so it could not cause significant changes in lignin content. This finding may be observed because overexpression of PbCCR1 activates the expression of most lignin biosynthetic genes in Arabidopsis, but the activation resulting from PbCCR2 overexpression is weak (Fig. 7i). Overall, PbCCR1 and 2 have a degree of functional redundancy, and both demonstrate the ability to participate in lignin biosynthesis. However, PbCCR1 may be the major gene for lignin biosynthesis, while PbCCR2 has little effect on lignin biosynthesis. This finding is consistent with the fact that there are multiple CCR genes with catalytic ability in a species but only one major gene [38, 66].

Conclusion

By investigating pear fruit at 0, 7 and 15 DAP, we found that primordial stone cells began to appear at 7 DAP, and the fruit had formed a large number of stone cells at 15 DAP. We performed transcriptome sequencing on the fruits at 0, 7 and 15 DAP and identified 3834 (0 vs. 7 DAP), 4049 (7 vs. 15 DAP) and 5763 (0 vs. 15 DAP) DEGs. The stone cells of ‘Dangshan Su’ pear fruit began to form in large quantities after 7 DAP, and the key enzyme genes in the lignin biosynthetic pathway were up-regulated widely. Before 7 DAP, the fruit primarily activates upstream material metabolic pathways (carbon metabolism and citrate cycle pathway) to provide important precursor substances for lignin biosynthesis. Further analysis found that the biosynthesis of S-type lignin in ‘Dangshan Su’ pear does not depend on the catalytic activity of PbSAD but is primarily generated by the catalytic activity of caffeoyl-CoA through CCoAOMT, CCR, F5H and CAD. In addition, some transcription factors regulating lignin biosynthesis were identified, such as PbNAC46, PbMYB105, 104, PbSHN and PbAP2-24. Subsequently, we cloned PbCCR1 and 2 and analysed their functions in the Chinese white pear lignin biosynthetic process. Overall, PbCCR1, 2 has a degree of functional redundancy, both of which demonstrate the ability to participate in lignin biosynthesis. However, PbCCR1 may be the major gene for lignin biosynthesis, while PbCCR2 has little effect on lignin biosynthesis.

Methods

Plant materials

The plant material used in this experiment was the ‘Dangshan Su’ pear, which grows in the Yuanyichang agricultural park (Dangshan County, Anhui Province, China). In the flowering period of ‘Dangshan Su’ pear (approximately April 1 of each year), ‘Cuiguan’ (Pyrus pyrifolia) was used as the pollen parent of the fruit. We first collected the pollen of ‘Cuiguan’ for pollinating ‘Dangshan Su’ pear. Selecting 40-year-old Pyrus bretschneideri cv. Dangshan Su trees with good growth, the stamens of the flowers were manually removed, and the pollen was applied to the stigma for artificial pollination. Immediately after artificial pollination, the stigma was wrapped in a bag and cultured for 7 days in the dark. After artificial pollination, we ensure that the fruits after pollination were kept under the same growth conditions, including water and fertilizer and management conditions. At 7 and 15 DAP, we collected the fruit for transcriptome sequencing along with the fruit from 0 DAP. For consistency with our previous study of the stone cells and lignin of ‘Dangshan Su’ pear, we collected the corresponding fruits for qRT-PCR analysis at 39, 47, 55, 63, 79, 102, and 145 DAP. Some of the 0, 7, and 15 DAP fresh fruit was used for sectioning, microscopic observation and determination of stone cell and lignin contents, and the remaining fruits were stored in a − 80 °C freezer for subsequent experiments. Three random biological repeats were set up for fruit collection in each period to eliminate the errors caused by differences between fruit samples.

Determination of stone cells and lignin in pear fruit

The fruits at 0, 7, and 15 DAP were selected for the determination of stone cell and lignin contents. Each group of samples was weighed to 5.0 g of fruit and frozen at − 20 °C for 1 day. The fruit samples were then homogenized at 20000 r·min− 1 for 5 min. Water was added to the homogenate, which was then allowed to stand. After all of the contents sank to the bottom of the beaker, the suspended matter in the upper layer was poured out (this operation was repeated five times). Finally, the stone cells were dried and weighed, and the stone cell content of each fruit sample was calculated as follows: formula weight of stone cells (g DW)/weight of pulp (g FW) × 100 = stone cell content (%) [5].

The lignin content of ‘Dangshan Su’ pear fruit was also determined at 0, 7, and 15 DAP. The fruit epidermis was removed and dried in an oven at 37 °C, and the dried fruit was ground finely. The powder was extracted with methanol after grinding, and then, the extracted residue was dried. Next, 0.2 g dried residue was extracted for 1 h in 15 ml 70% H2SO4 at 30 °C. After extraction, 115 mL of distilled water was added, and the solution was boiled for 1 h; the total volume of the sample remained unchanged during boiling. The boiled mixture was filtered with filter paper and rinsed with 70 °C water until the rinse water was clear. The lignin residue was dried and weighed, and each sample was repeated three times.

Extraction of total RNA from plant material and qRT-PCR

The total RNA of all the used ‘Dangshan Su’ pear fruits in this study was extracted using a plant RNA extraction kit from Tiangen (Beijing, China) (Additional file 9: Figure S1). Reverse transcription after total RNA extraction from the fruit was performed using a PrimeScriptâ„¢ RT reagent kit with gDNA Eraser (Takara, Kusatsu, Japan). All qRT-PCR primers in this study were designed using Beacon Designer 7 software, and liquid primers were ordered from Sangon Biotech (Shanghai, China) (Additional file 7: Table S7). We tested the amplification efficiency of all primers by reference to the method of Zhu et al. (Additional file 7: Table S7) [67]. At the same time, the standard curves of each primer in the detection of amplification efficiency were listed (Additional file 25: Figure S17). The qRT-PCR system consisted of 10 μL SYBR® Premix Ex TaqTM II (2X), 6.4 μL water, 0.8 μL upstream and downstream primers and 2 μL cDNA (20 μL total). For the qRT-PCR experiment, we chose the tubulin gene (No. AB239680.1) as the internal reference gene [68]. Each gene was evaluated with three biological replicates during qRT-PCR, and the 2-△△CT method was used to process the data and calculate the relative expression levels of the genes [69].

Transcriptome sequencing

Referring to the methods of total RNA extraction and reverse transcription in plants, a DNA library of fruit samples from 0, 7, 15 DAP stages was constructed. Three sets of biological replicates were set for each period of fruit (a total of 9 cDNA libraries were constructed for transcriptome sequencing) to eliminate errors caused by differences between fruits. The constructed library was analysed using the Illumina HiSeqâ„¢ 2500 sequencing platform (Sangon Biotech, Shanghai, China). From the raw data, we first filtered out the adaptor sequences and low-quality reads. The remaining sequences were considered clean reads and were mapped to the reference genome of Chinese white pear using HISAT2 software. Finally, unigene expression was calculated as FPKM with the software package Cufflinks [70]. The raw data has been uploaded to the SRA database of NCBI (Additional file 1: Table S1).

Identification, classification and metabolic pathway analysis of differentially expressed genes

In this experiment, either 0 DAP was used as the control and 7, 15 DAP were used as the experimental groups or 7 DAP was used as the control and 15 DAP was used as the experimental group. Differential expression analysis was performed to find genes with different expression levels between the two samples. A direct expression of a gene’s expression level is the abundance of its transcript. The higher the transcript abundance, the higher the gene expression level is. We used the TPM value to calculate the expression level of each gene. The BlastX and BlastN programs were used to compare the assembled unigenes with protein and nucleic acid databases (e-value< 10− 5). These databases include the nonredundant protein database (http://ncbi.nlm.nih.gov/), the nonredundant nucleic acid sequence database (http://ncbi.nlm.nih.gov/), the KEGG database (http://www.kegg.jp), the Cluster of Orthologous Groups (COG) database (https://www.ncbi.nlm.nih.gov/COG/), and the GO annotation database (http://www.geneontology.org). GO annotations, including MF, BP and CC, were also analysed. The best alignment results were selected for functional annotation of the gene and to determine the metabolic pathway in which the gene is involved.

Cloning of genes and construction of plant expression vectors

Full-length sequence-specific primers were designed using Primer Premier 6.0 software based on the full-length sequences of PbCCR1, 2. RT-PCR was performed using ‘Dangshan Su’ pear fruit cDNA as template. Primers with restriction sites were also designed using Primer Premier 6.0 software (PbCCR1 and 2 both used NcoI and SpeI restriction endonuclease sites) (Additional file 8: Table S8). Each gene fragment was ligated into the pCAMBIA1304 vector using T4 ligase at 16 °C for 3 h to obtain complete pCAMBIA1304-PbCCR1, 2 recombinant plasmids.

Screening of transgenic Arabidopsis and determination of lignin content

Enrichment culture of Agrobacterium liquid carrying pCAMBIA1304-PbCCR1 and 2 plasmids was performed. After the enrichment culture was completed, the supernatant was discarded after centrifugation, and the cell pellet was retained. The cells were resuspended in invasion buffer (0.02% Silwet L-77, 1/2 MS, 5% sucrose) such that the OD600 value of the solution was between 0.7 and 0.8 for subsequent infection. Wild-type Arabidopsis plants with a large number of flower buds and good growth condition were selected for infection (Arabidopsis seeds were purchased from the Arabidopsis Biological Resource Center). The whole flower buds were soaked in the Agrobacterium solution for approximately 45 s. After infection, the Arabidopsis plants were cultured in a dark environment for 24 h and then moved out of the dark environment and grown under normal light. The seeds were collected after maturation, sterilised, and germinated on MS medium containing hygromycin. The plants that could grow normally were potential positive plants, and the positive transgenic plants were obtained after a GUS staining test.

According to the results of Anderson et al., the acetyl bromide method was used to detect lignin content [71]. The inflorescence stems of the 60-day-old T3 generation of transgenic Arabidopsis were sliced by hand, stained with 1% toluidine blue on slides and observed directly under a microscope. TEM analysis was performed with reference to the method of Wu et al. [72].

Availability of data and materials

Data supporting the results of this study can be obtained from [NCBI], but the availability of these data is limited (accession number: SRR9998933, SRR9998932, SRR9998940, SRR9998936, SRR9998935, SRR9998934, SRR9998939, SRR9998938, SRR9998937). These data are used under the permission of the current study and are therefore not publicly available. We chose to publish our data on June 20, 2020. Until such a time, the data is available from the corresponding author upon reasonable request.

Abbreviations

4CL:

4-hydroxycinnamate-CoA ligase

BP:

Biological process

C4H:

Cinnamate 4-hydroxylase

CAD:

Cinnamyl alcohol dehydrogenase

CC:

Cellular component

CCoAOMT:

Caffeoyl-CoA O-methyltransferase

CCR:

Cinnamoyl-CoA reductase

CESA :

Cellulose synthase

COMT:

Caffeic acid 3-O-methyltransferase

DAP:

Days after pollination

DEGs:

Differentially expressed genes

F5H:

Ferulate-5-Hydroxylase

GO:

Gene ontology

KEGG:

Kyoto encyclopedia of genes and genomes

KOG:

Karyotic ortholog groups

LAC:

Laccase

L-Phe:

L-phenylalanine

MF:

Molecular function

PAL:

Phenylalanine ammonia-lyase

POD:

Peroxidase

SAD:

Sinapyl alcohol dehydrogenase

SCW:

Secondary cell wall

SDRs:

Short-chain dehydrogenase superfamilies

TPM:

Transcripts per million

References

  1. Lace B, Lacis G. Evaluation of pear (Pyrus communis L.) cultivars in Latvia. Hortic Sci. 2015;42:107–13.

    Article  Google Scholar 

  2. Xue B, Yao JL, Qin MF, Zhang MY, Allan AC, Wang DF, Wu J. PbrmiR397a regulates lignification during stone cell development in pear fruit. Plant Biotechnol J. 2019;17:103–17.

    Article  CAS  PubMed  Google Scholar 

  3. Hamauzu Y, Forest F, Hiramatsu K, Sugimoto M. Effect of pear (Pyrus communis L.) procyanidins on gastric lesions induced by HCl/ethanol in rats. Food Chem. 2007;100:255–63.

    Article  CAS  Google Scholar 

  4. Huang LJ, Gao WY, Li X, Zhao WS, Huang LQ, Liu CX. Valuation of the in vivo anti-inflammatory effects of extracts from Pyrus bretschneideri Rehd. J Agr Food Chem. 2010;58:8983–7.

    Article  CAS  Google Scholar 

  5. Zhang JY, Cheng X, Jin Q, Su XQ, Li ML, Yan CC, Jiao XY, Li DH, Lin Y, Cai YP. Comparison of the transcriptomic analysis between two Chinese white pear (Pyrus bretschneideri Rehd.) genotypes of different stone cells contents. Plos One. 2017;12:e0187114.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Cai YP, Li GQ, Nie JQ, Lin Y, Nie F, Zhang JY, Xu YL. Study of the structure and biosynthetic pathway of lignin in stone cells of pear. Sci Hortic. 2010;125:374–9.

    Article  CAS  Google Scholar 

  7. Jin Q, Yan CC, Qiu JX, Zhang N, Lin Y, Cai YP. Structural characterization and deposition of stone cell lignin in Dangshan Su pear. Sci Hortic. 2013;155:123–30.

    Article  CAS  Google Scholar 

  8. Yan CC, Yin M, Zhang N, Jin Q, Fang Z, Lin Y, Cai YP. Stone cell distribution and lignin structure in various pear varieties. Sci Hortic. 2014;174:142–50.

    Article  CAS  Google Scholar 

  9. Nii N, Kawahara T, Nakao Y. The development of stone cells in Japanese pear fruit. J Hortic Sci Biotechnol. 2008;83:148–53.

    Article  Google Scholar 

  10. Tao ST, Khanizadeh S, Zhang H, Zhang SL. Anatomy, ultrastructure and lignin distribution of stone cells in two Pyrus species. Plant Sci. 2009;176:413–9.

    Article  CAS  Google Scholar 

  11. Brahem M, Renard CMGC, Gouble B, Bureau S, Bourvellec CL. Characterization of tissue specific differences in cell wall polysaccharides of ripe and overripe pear fruit. Carbohyd Polym. 2017;156:152–64.

    Article  CAS  Google Scholar 

  12. Zhao SG, Zhang JG, Zhao YP, Zhang YX. New discoveries of stone cell differentiation in fruitlets of ‘Yali’ pears (Pyrus bretschneideri Rehd.). J Food Agric Enviro. 2017;11:937–42.

    Google Scholar 

  13. Ruben V, Barbara DM, John R, Wout B. Lignin biosynthesis and its integration into metabolism. Curr Opin Biotech. 2019;56:230–9.

    Article  CAS  Google Scholar 

  14. Baucher M, Halpin C, Petit-Conil M, Boerjan W. Lignin: genetic engineering and impact on pulping. Crit Rev Biochem Mol Biol. 2003;38:305–50.

    Article  CAS  PubMed  Google Scholar 

  15. Reddy M, Chen F, Shadle G, Jackson L, Aljoe H, Dixon RA. Targeted down-regulation of cytochrome P450 enzymes for forage quality improvement in alfalfa (Medicago sativa L.). P Natl Acad Sci USA. 2005;102:16573–8.

    Article  CAS  Google Scholar 

  16. Chen F, Reddy M, Temple S, Jackson L, Shadle G, Dixon RA. Multi-site genetic modulation of monolignol biosynthesis suggests new routes for formation of syringyl lignin and wall-bound ferulic acid in alfalfa (Medicago sativa L.). Plant J. 2006;48:113–24.

    Article  CAS  PubMed  Google Scholar 

  17. Zhao J. Lignification: flexibility, biosynthesis and regulation. Trends Plant Sci. 2016;21:713–21.

    Article  CAS  PubMed  Google Scholar 

  18. Liang XW, Dron M, Cramer CL, Dixon RA, Lamb CJ. Differential regulation of phenylalanine ammonia-lyase genes during plant development and by environmental cues. J Biol Chem. 1989;264:14486–92.

    CAS  PubMed  Google Scholar 

  19. Gou MY, Ran XZ, Martin DW, Liu CJ. The scaffold proteins of lignin biosynthetic cytochrome P450 enzymes. Nat Plants. 2018;4:299–310.

    Article  CAS  PubMed  Google Scholar 

  20. Wei H, Xu C, Movahedi A, Sun W, Qiang ZG. Characterization, expression profiling, and functional analyses of a 4CL-like gene of Populus trichocarpa. Processes. 2019;7:45.

    Article  Google Scholar 

  21. Wu YY, Wang NF, Hisano H, Cao YP, Wu FY, Liu WW, Bao Y, Wang ZY, Fu CX. Simultaneous regulation of F5H in COMT-RNAi transgenic switchgrass alters effects of COMT suppression on syringyl lignin biosynthesis. Plant Biotechnol J. 2019;17:836–45.

    Article  CAS  PubMed  Google Scholar 

  22. Tetreault HM, Scully ED, Gries T, Palmer NA, Funnell-Harris DL, Baird L, Seravalli J, Dien BS, Sarath G, Clemente TE, Sattler SE. Overexpression of the Sorghum bicolor SbCCoAOMT alters cell wall associated hydroxycinnamoyl groups. PLoS One. 2018;13:e0204153.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Yan XJ, Liu J, Kim H, Liu BG, Huang X, Yang ZC, Lin YCJ, Chen H, Yang CM, Wang JP, Muddiman DC, Ralph J, Sederoff RR, Li QZ, Chiang VL. CAD1 and CCR2 protein complex formation in monolignol biosynthesis in Populus trichocarpa. New Phytol. 2019;222:244–60.

    Article  CAS  PubMed  Google Scholar 

  24. Cheng X, Li ML, Li DH, Zhang JY, Jin Q, Sheng LL, Cai YP, 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–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Liu CJ. Deciphering the enigma of lignification: precursor transport, oxidation, and the topochemistry of lignin assembly. Mol Plant. 2012;5:304–17.

    Article  CAS  PubMed  Google Scholar 

  26. Cheng X, Li GH, Ma CH, Abdullah M, Zhang JY, Zhao H, Jin Q, Cai YP, 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng QD, Chen ZH, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–U130.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Bai CK, Xu J, Cao B, Li X, Li GS. Transcriptomic analysis and dynamic expression of genes reveal flavonoid synthesis in Scutellaria viscidula. Acta Physiol Plant. 2018;40:161.

    Article  CAS  Google Scholar 

  29. Jaiswal S, Antala TJ, Mandavia MK, Chopra M, Jasrotia RS, Tomar RS, Kheni J, Angadi UB, Iquebal MA, Golakia BA, Rai A, Kumar D. Transcriptomic signature of drought response in pearl millet (Pennisetum glaucum L.) and development of web-genomic resources. Sci Rep. 2018;8:3382.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Pascual MB, El-Azaz J, de la Torre FN, Canas RA, Avila C, Canovas FM. Biosynthesis and metabolic fate of phenylalanine in conifers. Front Plant Sci. 2016;7:1030.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Ko JH, Jeon HW, Kim WC, Kim JY, Han KH. The MYB46/MYB83-mediated transcriptional regulatory programme is a gatekeeper of secondary wall biosynthesis. Ann Bot. 2014;114:1099–107.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Mitsuda N, Seki M, Shinozaki K, Ohme-Takagi M. The NAC transcription factors NST1 and NST2 of Arabidopsis regulate secondary wall thickenings and are required for anther dehiscence. Plant Cell. 2005;7:1030.

    Google Scholar 

  33. Ambavaram MMR, Krishnan A, Trijatmiko KR, Pereira A. Coordinated activation of cellulose and repression of lignin biosynthesis pathways in rice. Plant Physiol. 2011;155:916–31.

    Article  CAS  PubMed  Google Scholar 

  34. Cao YP, Han YH, Li DH, Lin Y, Cai YP. MYB transcription factors in chinese pear (Pyrus bretschneideri Rehd.): genome-wide identification, classification, and expression profiling during fruit development. Front Plant Sci. 2016;7:577.

    PubMed  PubMed Central  Google Scholar 

  35. Li XL, Tao S, Wei SW, Ming ML, Huang XS, Zhang SL, Wu J. The mining and evolutionary investigation of AP2/ERF genes in pear (Pyrus). BMC Plant Biol. 2018;18:46.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Ahmad M, Yan XH, Li JZ, Yang QS, Jamil W, Teng YW, Bai SL. Genome wide identification and predicted functional analyses of NAC transcription factors in Asian pears. BMC Plant Biol. 2018;18:214.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Lacombe E, Hawkins S, Van Doorsselaere J, Piquemal J, Goffner D, Poeydomenge O, Boudet AM, Grima-Pettenati J. Cinnamoyl CoA reductase, the first committed enzyme of the lignin branch biosynthetic pathway: cloning, expression and phylogenetic relationships. Plant J. 1997;11:429–41.

    Article  CAS  PubMed  Google Scholar 

  38. Chao N, Li N, Qi Q, Li S, Lv T, Jiang XN, Gai Y. Characterization of the cinnamoyl-CoA reductase (CCR) gene family in Populus tomentosa reveals the enzymatic active sites and evolution of CCR. Planta. 2017;245:61–75.

    Article  CAS  PubMed  Google Scholar 

  39. Li SM, Su XQ, Jin Q, Li GH, Sun YM, Abdullah M, Cai YP, Lin Y. iTRAQ-based identification of proteins related to lignin synthesis in the pear pollinated with pollen from different varieties. Molecules. 2018;23:548.

    Article  PubMed Central  CAS  Google Scholar 

  40. Xue C, Yao JL, Xue YS, Su GQ, Wang L, Lin LK, Allan AC, Zhang SL, Wu J. PbrMYB169 positively regulates lignification of stone cells in pear fruit. J Exp Bot. 2019;70:1801–14.

    Article  PubMed  Google Scholar 

  41. Li XL, Liu L, Ming ML, Hu HJ, Zhang MY, Fan J, Song BB, Zhang SL, Wu J. Comparative transcriptomic analysis provides insight into the domestication and improvement of pear (P. pyrifolia) fruit. Plant Physiol. 2019. https://doi.org/10.1104/pp.18.01322.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Li SM, Su XQ, Abdullah M, Sun YM, Li GH, Cheng X, Lin Y, Cai YP. Effects of different pollens on primary metabolism and lignin biosynthesis in pear. Int J Mol Sci. 2018;19:2273.

    Article  PubMed Central  CAS  Google Scholar 

  43. Choi JH, Lee SH. Distribution of stone cell in Asian, Chinese, and European pear fruit and its morphological changes. J Appl Bot Food Qual. 2013;86:185–9.

    Google Scholar 

  44. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Li LG, Cheng XF, Leshkevich J, Umezawa T, Harding SA, Chiang VL. The last step of syringyl monolignol biosynthesis in angiosperms is regulated by a novel gene encoding sinapyl alcohol dehydrogenase. Plant Cell. 2001;13:1567–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Sibout R, Eudes A, Mouille G, Pollet B, Lapierre C, Jouanin L, Seguin A. CINNAMYL ALCOHOL DEHYDROGENASE-C and -D are the primary genes involved in lignin biosynthesis in the floral stem of Arabidopsis. Plant Cell. 2005;17:2059–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Zhang KW, Qian Q, Huang ZJ, Wang YQ, Li M, Hong LL, Zeng DL, Gu MH, Chu CC, Cheng ZK. GOLD HULL AND INTERNODE2 encodes a primarily multifunctional cinnamyl-alcohol dehydrogenase in rice1. Plant Physiol. 2006;140:972–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Rueda-Lopez M, Canas RA, Canales J, Canovas FM, Avila C. The overexpression of the pine transcription factor PpDof5 in Arabidopsis leads to increased lignin content and affects carbon and nitrogen metabolism. Physiol Plantarum. 2015;155:369–83.

    Article  CAS  Google Scholar 

  49. Rohde A, Morreel K, Ralph J, Goeminne G, Hostyn V, De Rycke R, Kushnir S, Van Doorsselaere J, Joseleau JP, Vuylsteke M, Van Driessche G, Van Beeumen J, Messens E, Boerjan W. Molecular phenotyping of the pal1 and pal2 mutants of Arabidopsis thaliana reveals far-reaching consequences on phenylpropanoid, amino acid, and carbohydrate metabolism. Plant Cell. 2004;16:2749–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Arioli T, Peng L, Betzner AS, Burn J, Wittke W, Herth W, Camilleri C, Hofte H, Plazinski J, Birch R, Cork A, Glover J, Redmond J, Williamson RE. Molecular analysis of cellulose biosynthesis in Arabidopsis. Science. 1998;279:717–20.

    Article  CAS  PubMed  Google Scholar 

  51. Fagard M, Desnos T, Desprez T, Goubet F, Refregier G, Mouille G, McCann M, Rayon C, Vernhettes S, Hofte H. PROCUSTE1 encodes a cellulose synthase required for normal cell elongation specifically in roots and dark-grown hypocotyls of Arabidopsis. Plant Cell. 2000;12:2409–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Scheible WR, Eshed R, Richmond T, Delmer D, Somerville C. Modifications of cellulose synthase confer resistance to isoxaben and thiazolidinone herbicides in Arabidopsis Ixr1 mutants. P Natl Acad Sci USA. 2001;98:10079–84.

    Article  CAS  Google Scholar 

  53. Taylor NG, Howells RM, Huttly AK, Vickers K, Turner SR. Interactions among three distinct CesA proteins essential for cellulose synthesis. Proc Natl Acad Sci U S A. 2003;100:1450–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Taylor NG, Laurie S, Turner SR. Multiple cellulose synthase catalytic subunits are required for cellulose synthesis in Arabidopsis. Plant Cell. 2000;12:2529–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Yang C, Xu Z, Song J, Conner K, Vizcay Barrena G, Wilson ZA. Arabidopsis MYB26/MALE STERILE35 regulates secondary thickening in the endothecium and is essential for anther dehiscence. Plant Cell. 2007;19:534–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Gray J, Caparros-Ruiz D, Grotewold E. Grass phenylpropanoids: regulate before using! Plant Sci. 2012;184:112–20.

    Article  CAS  PubMed  Google Scholar 

  57. Zhong RQ, Ye ZH. MYB46 and MYB83 bind to the SMRE sites and directly activate a suite of transcription factors and secondary wall biosynthetic genes. Plant Cell Physiol. 2012;53:368–80.

    Article  CAS  PubMed  Google Scholar 

  58. Guo WF, Jin L, Miao YH, He X, Hu Q, Guo K, Zhu LF, Zhang XL. An ethylene response-related factor, GbERF1-like, from Gossypium barbadense improves resistance to Verticillium dahliae via activating lignin synthesis. Plant Mol Biol. 2016;91:305–18.

    Article  CAS  PubMed  Google Scholar 

  59. Ma RF, Xiao Y, Lv ZY, Tan HX, Chen RB, Li Q, Chen JF, Wang Y, Yin J, Zhang L, Chen WS. AP2/ERF transcription factor, Ii049, positively regulates Lignan biosynthesis in Isatis indigotica through activating salicylic acid signaling and Lignan/lignin pathway genes. Front Plant Sci. 2017;8:1361.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Lee DK, Jung H, Jang G, Jeong JS, Kim YS, Ha SH, Choi DY, Kim JK. Overexpression of the OsERF71 transcription factor alters Rice root structure and drought resistance. Plant Physiol. 2016;172:575–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Jorma V, Judith F, Jonathan L, Andras G, Lorenz G, Airi L, Jaakko K, Bjorn S. A genome-wide screen for ethylene-induced ethylene response factors (ERFs) in hybrid aspen stem identifies ERF genes that modify stem growth and wood properties. New Phytol. 2013;200:511–22.

    Article  CAS  Google Scholar 

  62. Soichi K, Yukie C, Masayuki K, Kazuhiro M, Yoshihiro O. Upregulation of the promoter activity of the carrot (Daucus carota) phenylalanine ammonia-lyase gene (DcPAL3) is caused by new members of the transcriptional regulatory proteins, DcERF1 and DcERF2, which bind to the GCC-box homolog and act as an activator to the DcPAL3 promoter. J Plant Res. 2008;121:499–508.

    Article  CAS  Google Scholar 

  63. Zeng JK, Xian L, Qian X, Chen JY, Yin XR, Ferguson IB, Chen KS. EjAP2-1, an AP2/ERF gene, is a novel regulator of fruit lignification induced by chilling injury, via interaction with EjMYB transcription factors. Plant Biotechnol J. 2015;13:1325–34.

    Article  CAS  PubMed  Google Scholar 

  64. Sun QN, Foston M, Meng XZ, Sawada D, Pingali SV, O'Neill HM, Li HJ, Wyman CE, Langan P, Ragauskas AJ, Kumar R. Effect of lignin content on changes occurring in poplar cellulose ultrastructure during dilute acid pretreatment. Biotechnol Biofuels. 2014;7:150.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  65. Zhang WB, Wei R, Chena S, Jiang J, Lia HY, Huang HJ, Yang G, Wang S, Wei HR, Liu GF. Functional characterization of CCR in birch (Betula platyphylla × Betula pendula) through overexpression and suppression analysis. Physiol Plantarum. 2015;154:283–96.

    Article  CAS  Google Scholar 

  66. Park HL, Bhoo SH, Kwon M, Lee SW, Cho MH. Biochemical and expression analyses of the Rice Cinnamoyl-CoA reductase gene family. Front Plant Sci. 2017;8:2099.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Zhu LF, Yang CQ, You YH, Liang W, Wang NN, Ma FW, Li CY. Validation of reference genes for qRT-PCR analysis in peel and flesh of six apple cultivars (Malus domestica) at diverse stages of fruit development. Sci Hortic. 2019;244:165–71.

    Article  CAS  Google Scholar 

  68. Wu T, Zhang RP, Gu C, Wu JY, Wan HJ, Zhang SJ, Zhang SL. Evaluation of candidate reference genes for real time quantitative PCR normalization in pear fruit. Afr J Agr Res. 2012;7:3701–4.

    Google Scholar 

  69. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-△△CT method. Methods. 2001;25:402–8.

    Article  CAS  PubMed  Google Scholar 

  70. Li DH, Wu D, Li SZ, Guo N, Gao JS, Sun X, Cai YP. Transcriptomic profiling identifies differentially expressed genes associated with programmed cell death of nucellar cells in Ginkgo biloba L. BMC Plant Biol. 2019;19:91.

    Article  PubMed  PubMed Central  Google Scholar 

  71. Anderson NA, Tobimatsu Y, Ciesielski PN, Ximenes E, Ralph J, Donohoe BS, Michael L, Clint C. Manipulation of Guaiacyl and Syringyl monomer biosynthesis in an Arabidopsis Cinnamyl alcohol dehydrogenase mutant results in atypical lignin biosynthesis and modified cell wall structure. Plant Cell. 2015;27:2195–209.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Wu XW, Lu XP, Riaz M, Yan L, Jiang CC. Boron toxicity induced specific changes of cell ultrastructure and architecture of components in leaf center and tip of trifoliate orange (Poncirus trifoliata (L.) Raf.). J Environ Manag. 2019;246:426–33.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We would like to thank Dr. Yunpeng Cao and Dr. Yongzhi Hua of Anhui Agricultural University for their help.

Funding

The research was supported by the National Natural Science Foundation of China (Grant #31640068), the Anhui Provincial Natural Science Foundation (Grant #1808085QC79) and the Graduate innovation fund of Anhui Agricultural University (Grant #2019yjs-54). The Funding bodies were not involved in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

XQS and YPC designed the research project; XQS, YZ, GHL, HW, XC and QJ carried out the experimental work; XQS and YZ analyzed the data; XQS wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Yongping Cai.

Ethics declarations

Ethics approval and consent to participate

The experiments did not involve endangered or protected species. No specific permits were required for these locations/activities because the pears used in this study were obtained from Yuanyichang agricultural park in Dangshan County, which were demonstration orchards at Auhui Agricultural University.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Accession number of RNA-Seq raw data in NCBI SRA database.

Additional file 2: Table S2.

Summary of KEGG pathways.

Additional file 3: Table S3.

KEGG pathway classification.

Additional file 4: Table S4.

The expression level of putative novel genes related to stone cell development.

Additional file 5: Table S5.

Cellulose synthase gene related to secondary wall biosynthesis in Chinese white pear.

Additional file 6: Table S6.

Putative functions of pear MYB, NAC, AP2/ERF transcription factors.

Additional file 7: Table S7.

Sequence of primers for qRT-PCR.

Additional file 8: Table S8.

Primer sequences used in transgenic Arabidopsis.

Additional file 9: Figure S1.

RNA quality testing of all nine samples.

Additional file 10: Figure S2.

PCA analysis of biological repeats of 0, 7, 15 DAP samples.

Additional file 11: Figure S3.

Duplicate correlation checking scatter plot analysis of biological repeats of 0, 7, 15 DAP samples.

Additional file 12: Figure S4.

Read distribution in different regions blasted against the reference genome.

Additional file 13: Figure S5.

Relationship between gene TPM deviation within 15% of the final value (based on 100% mapped reads) and percentage of mapped reads.

Additional file 14: Figure S6.

Density distribution of gene expression in 0, 7, 15 DAP fruit. Gene expression levels are represented as log2(TPM).

Additional file 15: Figure S7.

Cluster analysis of transcript expression.

Additional file 16: Figure S8.

Volcano plot of differentially expressed genes.

Additional file 17: Figure S9.

Unigene KOG function categories.

Additional file 18: Figure S10.

Unigene GO function categories.

Additional file 19: Figure S11.

Visual KEGG path classification.

Additional file 20: Figure S12.

DEGs related to phenylpropanoid metabolism in the 7 DAP vs. 15 DAP comparison. Red indicates upregulation of expression, while yellow indicates no change in expression.

Additional file 21: Figure S13.

qRT-PCR validation of the expression levels of putative novel genes related to stone cell development.

Additional file 22: Figure S14.

Expression profiles of the MYB, NAC, and AP2/ERF genes in Chinese white pear fruit at 0, 7, and 15 DAP. The fragments per kilobase of exon per million mapped reads (TPM) values were obtained by RNA sequencing analysis and are presented as a heat map. The colour scale shows different levels of gene expression. Red indicates a high level of expression, white signifies a medium level of expression, and blue denotes a low level of expression.

Additional file 23: Figure S15.

Expression patterns of PbCESA4-1, 7, 8 in fruit at different developmental stages.

Additional file 24: Figure S16.

PbAP2-24 alignment with EjAP2-1 proteins.

Additional file 25: Figure S17.

Standard curve of qRT-PCR primers.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Su, X., Zhao, Y., Wang, H. et al. Transcriptomic analysis of early fruit development in Chinese white pear (Pyrus bretschneideri Rehd.) and functional identification of PbCCR1 in lignin biosynthesis. BMC Plant Biol 19, 417 (2019). https://doi.org/10.1186/s12870-019-2046-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-019-2046-x

Keywords