Transcriptome profiling of fruit development and maturation in Chinese white pear (Pyrus bretschneideri Rehd)
- Min Xie†1,
- Ying Huang†1,
- Yanping Zhang1,
- Xin Wang1,
- Hua Yang2,
- Oliver Yu2,
- Wenhao Dai3 and
- Congbing Fang1, 2Email author
© Xie et al.; licensee BioMed Central Ltd. 2013
Received: 9 August 2013
Accepted: 20 November 2013
Published: 23 November 2013
Pear (Pyrus spp) is an important fruit species worldwide; however, its genetics and genomic information is limited. Combining the Solexa/Illumina RNA-seq high-throughput sequencing approach (RNA-seq) with Digital Gene Expression (DGE) analysis would be a powerful tool for transcriptomic study. This paper reports the transcriptome profiling analysis of Chinese white pear (P. bretschneideri) using RNA-seq and DGE to better understand the molecular mechanisms in fruit development and maturation of Chinese white pear.
De novo transcriptome assembly and gene expression analysis of Chinese white pear were performed in an unprecedented depth (5.47 gigabase pairs) using high-throughput Illumina RNA-seq combined with a tag-based Digital Gene Expression (DGE) system. Approximately, 60.77 million reads were sequenced, trimmed, and assembled into 90,227 unigenes. These unigenes comprised 17,619 contigs and 72,608 singletons with an average length of 508 bp and had an N50 of 635 bp. Sequence similarity analyses against six public databases (Uniprot, NR, and COGs at NCBI, Pfam, InterPro, and KEGG) found that 61,636 unigenes can be annotated with gene descriptions, conserved protein domains, or gene ontology terms. By BLASTing all 61,636 unigenes in KEGG, a total of 31,215 unigenes were annotated into 121 known metabolic or signaling pathways in which a few primary, intermediate, and secondary metabolic pathways are directly related to pear fruit quality. DGE libraries were constructed for each of the five fruit developmental stages. Variations in gene expression among all developmental stages of pear fruit were significantly different in a large amount of unigenes.
Extensive transcriptome and DGE profiling data at five fruit developmental stages of Chinese white pear have been obtained from a deep sequencing, which provides comprehensive gene expression information at the transcriptional level. This could facilitate understanding of the molecular mechanisms in fruit development and maturation. Such a database can also be used as a public information platform for research on molecular biology and functional genomics in pear and other related species.
The genus Pyrus is one of the most important genera in Rosaceae family for fruit production. Pyrus species are widely used for commercial fruit production in 76 countries or regions around the world (http://faostat.fao.org) and their economic importance has been well recognized . There are four major edible species in Pyrus: P. communis L. is mainly cultivated in Europe, North America, South America, Africa, and Australia. The other three species, P. bretschneideri Rehd., P. ussuriensis Maxim., and P. pyrifolia (Burm.) Nakai., are grown in East Asian [1, 2]. The world production of pear has doubled in the past 16 years. China is the largest pear producer. In 2010, China produced 15.23 million tons (Mt) of pear fruits that accounted of 67.26% of world pear production (22.64 Mt) (FAOSTAT, 2012).
As of July 2013, 26,696 nucleotide sequences, 4,413 expressed sequence tags (ESTs), 52 genome survey sequences (GSS), and 2,636 proteins from the Pyrus genus have been deposited in GenBank. These sequences mostly derived from cDNA cloning and EST sequencing [3–7] provide useful information for transcriptional analysis, candidate gene discovery, and gene functional analysis; however, a comprehensive description of genes that expressed in pear fruit during the fruit development and maturation period remains unavailable.
Recent years, use of RNA-seq, the next generation sequencing technology, has generated over one billion bp of high-quality DNA sequence per analysis/experiment  and has dramatically improved the efficiency of gene discovery and functional analysis [9, 10], which largely facilitates the investigation of the functional complexity of transcriptomes [11, 12]. Illumina sequencing of transcriptomes for yeast [13, 14], Arabidopsis, mouse [16, 17], and human cells [18, 19] has confirmed that next generation sequencing is well-suited for surveying transcriptome profiles in eukaryotes. Recently Illumina RNA-seq system has been used to identify genes related to bud dormancy in pear (P. pyrifolia) . RNA-seq is not limited to detecting transcripts for organisms with existing genomic sequences; it also can be used to sequence non-model organisms that lack of genomic information [21–25]. DGE is a tag-based transcriptome sequencing approach in which the expression level of all genes in a sample is measured by counting the number of individual mRNA molecules produced from each gene, which enables the DGE protocol more suitable and affordable for comparative gene expression studies [25–31]. Despite its obvious potential, the next generation sequencing has not yet been applied to pear research.
In this study, 5.47 gigabase pairs of high quality DNA sequence were generated using Illumina technology to survey the poly (A)+ transcriptome of P. bretschneideri. A total of 90,227 unigenes were assembled. All known homologous genes in major metabolic pathways related to fruit development and maturation were identified. Furthermore, five DGE libraries were constructed and gene expression profiles at different fruit developmental stages were analyzed. These annotated transcriptome sequences and gene expression profiles provide useful information for identification of genes involved in quality trait development during fruit development and maturation in pear species.
Plant material and RNA extraction
cDNA library construction, Illumina sequencing and De novoassembly
To obtain a complete gene expression profile, RNA from each tissue sample was pooled. The poly (A)+ RNA was isolated from 20 μg of the pooled RNA using Dynal oligo(dT) 25 beads (Invitrogen) according to the Illumina manufacturer’s instructions. The mRNA was then cleaved into short sequences using divalent cations at 70°C for 5 min. The cleaved RNA fragments were used for first strand cDNA synthesis using SuperScript III reverse transcriptase (Invitrogen) and N6 random hexamers (Takara). Second strand cDNA synthesis was performed using RNase H (Invitrogen) and DNA polymerase (Invitrogen). Subsequently, cDNA fragments were ligated to adapters after an end repair process. These products were purified and enriched using the QIAquick PCR Purification Kit to develop a cDNA library.
The cDNA library was sequenced from both 5′ and 3′ ends on the Illumina sequencing platform (HiSeqTM 2000). Image deconvolution and Q-value (quality value assigned to each base) calculation were carried out in the Illumina data processing pipeline (version 1.6). Before assembly, adaptor sequences, empty reads, low quality sequences with N percentage (i.e., the percentage of nucleotides in read which could not be sequenced) over 5%, and those containing more than 20% bases with Q-value ≤ 10 were removed. De novo assembly of the cleaned reads was performed using SOAPdenovo software (version 1.03, http://soap.genomics.org.cn). The de Bruijn graph was firstly used to generate contigs . The reads were then mapped back to the contigs, and the paired-end reads and contigs from the same transcript were further joined into scaffolds. The complete scaffolds were subsequently obtained after the paired-end reads were again used for gap fillings . To reduce any sequence redundancy, the scaffolds were further assembled into unigenes using TGICL (copyright(c), http://www.tigr.org/tdb/tgi) . Among them, the scaffolds with more than 70% similarity after multiple alignments were grouped into clusters, and others that cannot reach the threshold set and fall into any assembly should remain as a list of singletons. The assembled unigenes (larger than 200 bp) were deposited in the Transcriptome Shotgun Assembly Sequence Database (TSA) at NCBI with the accession number from JR595427 to JR673747.
Functional annotation and classification
Unigenes that were larger than 150 bp were used for BLAST search and annotation against the NCBI nonredundant (nr) database with the BLASTX algorithm (http://www.ncbi.nlm.nih.gov/) using an E-value cut-off of 10-5. Functional annotation by gene ontology terms (GO; http://www.geneontology.org) was carried out based on the best hits from nr annotation using Blast2go software (http://www.blast2go.de/) , and the resultant GO id were further used for GO classification by WEGO  (http://wego.genomics.org.cn/cgi-bin/wego/index.pl). Annotation of COG and KEGG pathways was performed by sequence comparisons using BLASTX algorithm against Cluster of Orthologous Groups database and Kyoto Encyclopedia of Genes and Genomes database with an E-value threshold of 10-5. Annotations of all assembled sequences were deposited in the Transcriptome Shotgun Assembly Sequence Database (TSA) at NCBI and can be searched using the Gene-ID listed in Additional file 1.
Digital gene expression (DGE) library construction and sequencing
Total RNA was extracted from fruit samples collected at five fruit developmental stages. The tag library was prepared using the Illumina gene expression sample prep kit. Briefly, the poly(A)+ RNA was purified from 6 μg of RNA using oligo(dT) magnetic beads. Double-stranded cDNA was directly synthesized on the beads and subsequently digested with NlaIII. All cDNA fragments with 3′ ends were purified with magnetic beads before their 5′ ends were ligated to Illumina adapter 1. A 21 bp tag with the adaptor 1 was produced after digested with MmeI (an enzyme that recognizes the junction of the Illumina adaptor 1 and the CATG site). The Illumina adapter 2 was then ligated to the 3′ ends of the tag. The library was amplified by PCR for 15 cycles and 90 bp strips were purified by PAGE gel electrophoresis. These strips were then digested and single-chain molecules were attached to the Illumina sequencing chip for sequencing. All raw tag data are available in Gene Expression Omnibus (GEO) at NCBI with the accession number: GSE33388.
Analysis and mapping of DGE tags
To map DGE tags, sequencing-received raw image data were filtered to remove low quality tags (tags with unknown sequences ‘N’), empty tags (sequence with only the adaptors but no tags), and tags with only one copy number (probable sequencing error). For annotation, cleaned tags with CATG and the 21 bp tag sequence were mapped to our transcriptome reference database with no more than 1 nucleotide mismatch. All tags that mapped to reference sequences of multiple genes were filtered out and the remaining tags were designated as unambiguous tags for gene expression analysis. The number of unambiguous tags of each gene was calculated and then normalized to TPM (number of transcripts per million clean tags). The differentially expressed tags were used for mapping and annotation. A complete list of all differentially expressed genes is shown in Additional file 2, Additional file 3, Additional file 4 and Additional file 5.
Evaluation of DGE libraries
To compare the gene expression in different fruit developmental stages, the frequency of each tag in different DGE libraries was statistically analyzed using the method of Audic et al.. The false discovery rate (FDR) was used to determine the threshold P-value in multiple tests. The threshold determining the significance of differentially expressed genes was set at FDR ≤ 0.001 and log2 foldchange ≥ 2.
Pathway enrichment analysis
Here N is the number of all genes that with KEGG annotation, n is the number of DEGs in N, M is the number of all genes that were annotated to specific pathways, and m is number of DEGs in M.
Illumina sequencing and de novoassembly
Summary of sequence assembly of P. bretschneideri mRNA pooled from ten different tissues samples
Base pairs (bp)
Mean length (bp)
Raw sequencing reads
Contigs (≥100 bp)
Clusters (≥150 bp)
Singletons (≥150 bp)
Total unigenes (≥150 bp)
Annotation of predicted proteins
Summary of annotations of the P. bretschneideri unigenes
All assembled unigenes
Gene annotations against NR plant proteins
Gene annotations against NR Arabidopsis proteins
Unique gene annotations against NR
Gene annotations against UniProt
Gene annotations against InterPro
Gene annotations against Pfam
Gene annotations against COG
Gene annotations against KEGG
GO annotations for NR protein hits
3 main categories, 44 sub-categories
All annotated Unigenes
Conserved domain annotation
Thirty most frequently occurring InterPro domains/families in the P. bretschneideri unigenes
Protein kinase, catalytic domain
Serine/threonine-protein kinase-like domain
Serine/threonine-protein kinase domain
Serine/threonine-protein kinase, active site
Tyrosine-protein kinase, catalytic domain
Protein kinase, ATP binding site
RNA recognition motif domain
WD40 repeat, subgroup
Zinc finger, RING-type
WD40 repeat 2
Zinc finger, C3HC4 RING-type
Leucine-rich repeat-containing N-terminal, type 2
HTH transcriptional regulator, Myb-type, DNA-binding
SANT domain, DNA binding
F-box domain, cyclin-like
Cytochrome P450, E-class, group I
EF-Hand 1, calcium-binding site
WD40 repeat, conserved site
Ankyrin repeat-containing domain
Helicase, superfamily 1/2, ATP-binding domain
Gene ontology (GO), clusters of orthologous groups (COG) and Kyoto Encyclopedia of genes and genomes (KEGG) ontology (KO) classification
A total of 90,227 annotated sequences were mapped to canonical pathways in Kyoto Encyclopedia of Genes and Genomes (KEGG)  to identify active pathways in white pear. Of which, 24,169 sequences were assigned to 121 KEGG pathways (Additional file 9). The most represented pathways by the unique sequences were metabolic pathways (5,488 members), biosynthesis of secondary metabolites (2,904 members), plant-pathogen interaction (2,389 members), and spliceosome (1,121 members). These annotations provide a valuable resource for investigating specific processes, functions, and pathways in pear growth and development.
Analysis of metabolic pathway genes using P. bretschneideri unigenes
The overall quality of pear fruits is closely related to flavor, grit (stone cell) content, flesh texture, skin russet, and appearance [1, 46]. During fruit development and maturation, pear fruits undergo a series of physiological and biochemical changes including expansion of size, accumulation of soluble solids, change of pigments, and formation of aromatic volatiles [1, 2]. Most of these traits are inherited in a complex polygenic manner and controlled by multiple QTLs, which poses a significant challenge to traditional breeding [2, 47]. In order to better understand the genetic and molecular basis of these changes related to quality formation, a few pathways including several primary, intermediate, and secondary metabolic pathways that are related to pear quality development were analyzed. A total of 61,636 annotated unigenes were used to analyze metabolic pathway genes with simple keyword searches. Each search result was confirmed with a BLAST search.
Digital gene expression (DGE) library sequencing
Statistics of P. bretschneideri DGE sequencing
Total number of raw tags
Total number of distinct raw tags
Total number of clean tags
Total number of distinct clean tags
Clean tag/Raw tag
Total number of tags mapping to genes
% of tags mapping to genes
Total number of distinct tags mapping to gene
% of distinct tags mapping to gene
Number of all tag-mapped genes
% of all tag-mapped genes
Number of unambiguous tags mapping to gene
% of unambiguous tags mapping to gene
Number of distinct unambiguous tags mapping to gene
% of distinct unambiguous tags mapping to gene
Number of unambiguous tag-mapped genes
% of unambiguous tag-mapped genes
Number of unknown tags
% of unknown tags
Number of distinct unknown tags
% of distinct unknown tags
Tag sequences in the five DGE libraries were mapped to the transcriptome reference database generated in the above-mentioned Illumina sequencing. The reference database contains 162,456 distinct sequences with 142,331 unambiguous reference tags. Among the distinct tags (68,916 to 77,039) generated from the Illumina sequencing in the five libraries of fruit developmental stages, 32,705 to 36,335 distinct tags were mapped to individual genes in the reference database (Table 4). Tags mapped to a unique sequence are the most critical subset of the DGE libraries as they can explicitly identify a transcript. Up to 31.2% (28,175) of the sequences in the transcriptome reference tag database could be unequivocally identified by a unique tag (Table 4).
Gene expression profile changes between developmental stages
To profile the gene expression pattern during different developmental stages of pear fruit, the number of clean tags for each gene was calculated and its differentially expressed tags between two samples of two adjacent stages were identified using an algorithm developed by Audic et al..
The comparative analysis between fruit stage 3 (FS3) and FS2 libraries revealed that 2,135 genes showed significant expression changes. Among these genes, 923 and 1,212 genes were up-regulated and down-regulated, respectively, in FS3 compared to FS2 (Figure 10, Additional file 3). Among the ten most up-regulated and ten most down-regulated genes, seven annotated genes were up-regulated in FS3: two prunin 2 precursor [Prunus dulcis], one unnamed protein product [Vitis vinifera], one dehydration-responsive protein RD22 [Prunus persica], one mitochondrial carrier protein [Ricinus communis], one conserved hypothetical protein [Ricinus communis], and one tonoplast intrinsic protein (putative) [Ricinus communis]. Seven functionally defined genes were down-regulated in FS3, including a flavonoid 3-hydroxylase gene [Ricinus communis], a RNA binding protein (putative) [Ricinus communis], an APETALA2 (AP2) domain class transcription factor [Malus x domestica], a CBF/DREB1 transcription factor [Malus baccata], an anthocyanidin synthase gene [Pyrus communis], an Auxin-binding protein ABP19a, and an NADH dehydrogenase subunit 6 [Nicotiana tabacum]. Six genes were also blasted without annotation (Additional file 12). Based on the GO functional classification, almost all genes were involved in cellular components, i.e., extracellular region, membrane part, anchored to membrane and external encapsulating structure. Among genes enriched in KO, significant changes were observed in genes of metabolic pathways, such as biosynthesis of phenylpropanoids, biosynthesis of plant hormones, flavonoid biosynthesis, and starch and sucrose metabolism (Additional file 14).
Between the FS4 and FS3 libraries, 1,824 genes demonstrated significant changes. In the FS4 library, 797 and 1,027 genes were up-regulated and down-regulated, respectively, in comparison with the FS3 library (Figure 10, Additional file 4). Nine genes among the ten most differentially up- and down-regulated genes showed no similarity (Additional file 12). Among the up-regulated genes, five were hypothetical protein from Vitis vinifera and Populus trichocarpa, in which only one matched a flavonoid 3-hydroxylase gene from Ricinus communis, two down-regulated genes were predicted to encode proteins found in Ricinus communis and Populus trichocarpa, an unnamed protein product [Vitis vinifera], a laccase 1a [Populus trichocarpa] and an Nonspecific lipid-transfer protein precursor [Ricinus communis] (Additional file 12). Almost all gene sets enriched in GO were correlated to plastid, plastid part, thylakoid, and chloroplast, while the gene sets enriched in KO were related to metabolic pathways (Additional file 15).
The comparison between the FS5 and FS4 libraries also revealed significant variations in gene expression. A total of 1,411 genes, including 484 up-regulated and 927 down-regulated, were identified in FS5 compared to FS4 (Figure 10, Additional file 5). Eight genes among the ten most differentially up- and down-regulated genes showed no similarity (Additional file 12). Among the ten up-regulated genes, one aligned with a beta-D-xylosidase gene [Pyrus pyrifolia], one with an AP2 domain class transcription factor [Malus x domestica], and two genes were unnamed protein product [Vitis vinifera] and conserved hypothetical protein [Ricinus communis]. Among the ten down-regulated genes, three aligned with dehydration-responsive protein RD22, two with prunin precursors, one with a delta-12 oleate desaturase [Gossypium hirsutum], and two were unnamed protein product or hypothetical protein. Most genes enriched in GO were related (or functional in) chloroplast stroma and photosystem. Genes in enriched KO showing significant changes between the FS5 and FS4 libraries were related to cysteine and methionine metabolism, photosynthesis-antenna proteins and isoquinoline alkaloid biosynthesis (Additional file 16).
Expression profiling during ripening
In this study, differences in the number and expression profile of DEGs were observed at five developmental stages of pear fruit. 28,743 unigenes were expressed during fruit ripening, with 810, 2,135, 1,824, and 1,411 showing differential expression between 25 and 55 DPA, 55 and 85 DPA, 85 and 115 DPA, and 115 and 145 DPA, respectively. It showed that fewer genes expressed in the cell division stage (25 to 55 DPA) than in the period of cell expansion (55 to 145 DPA) and fruit ripening (115 to 145 DPA) (Figures 1 and 10).
To characterize the levels of gene expression, the number of unambiguous tags of each gene was calculated and normalized to the number of TPM. Based on this analysis, the gene expression levels in the five fruit developmental stages were categorized to rare (TPM < 5), low (TPM > 5–50), moderate (TPM > 50–100), and high (TPM >100) (Additional file 17). The largest portion of transcripts (20,431 out of 28,744, 71.08%) exhibited rare expression, and followed by low expression (6,483 out of 28,744, 22.55%). Only a small fraction of transcripts was expressed at moderate (3.21%) and high (3.15%) levels.
Selected genes with different expression trends in five fruit developmental stages
TPM (transcript per million clean tag)
Sucrose phosphate synthase
Auxin influx transport protein
COL domain class transcription factor
UDPG flavonoid 7-O-glucosyltransferase
NAD-dependent sorbitol dehydrogenase
NAD-dependent sorbitol dehydrogenase
S6 ribosomal protein
Structural molecule activity
DEAD-box RNA helicase-like protein
Large ribosomal subunit
Plasma membrane ATPase
Plasma membrane intrinsic protein
Plasma intrinsic protein
BHLH domain class transcription factor
60S ribosomal protein L35a
Ribosomal protein S21e
Illumina sequencing and function annotation
In this study, one ‘Dangshansuli’ pear cDNA library and five DGE libraries (from fruit samples collected at 25, 55, 85, 115 and 145 DPA) were constructed using RNA-Seq technology. These constructed libraries were subjected to comparative gene expression studies. As a result, we obtained 90,227 unigenes including 17,619 clusters and 72,608 singletons. 62,077 out of 90,227 assembled unigenes were annotated. Of which, 61,624 were annotated to the non-redundant (nr) NCBI database and 28,114, 33,205, and 24,196 were annotated to GO, COG and KEGG databases, respectively. For the remaining 28,150 unigenes (31.2% of all 90,227 assembled unigenes), the absence of significant homology to existing genes could be caused by several factors. Obviously, the length and completeness of the assembled unigenes were one main factor, which can be seen from the increasing proportion of sequences with matches in the NR database as the length of uingenes increased (Figure 3A). However, the assembled unigenes shorter than 500 bp added up to 64,812 (71.8% of all 90,227 assembled unigenes) (Figure 2B), and some of them were too short to allow statistically meaningful matches. In addition, for some unigenes, lack of homologous sequences in the public databases may indicate specific roles for them in P. bretschneideri.
Unigenes identified through conserved domain annotation showed that most abundant InterPro domains/families were Protein kinase and its subcategories that are known to regulate the majority of cellular pathways. Cytochrome P450, Myb-DNA-binding families (MYB), and “WD40-repeat” domains were also very abundant. MYB and WD40-repeat proteins contributed to the extensive modification and regulation of various secondary compounds in many pathways have been validated to be key roles in anthocyanin biosynthesis by regulating anthocyanin structural genes . MYB proteins also appear to be versatile in regulating secondary metabolism, cellular morphogenesis, meristem formation, and the cell cycle . Recent research showed that Cytochrome P450 encoded enzymes that are responsible for the conversion of the mogroside backbone to various mogrosides .
ESTs sequences had been intensively used for transcriptional analysis, candidate gene discovery, and gene functional analysis. However, since RNA-Seq technology was developed, only 4,413 expressed sequence tags (ESTs) of Pyrus plants have been deposited in GenBank due to the fact that EST analyses can only identify limited number of candidate genes that are involved in complex biosynthetic pathways. Using RNA-seq, over one billion nucleotides of high-quality DNA sequence per analysis/experiment can be generated , which has dramatically improved the efficiency of gene discovery and functional analysis [9, 10]. Furthermore, RNA-Seq was less expensive, more efficient, and allowed faster gene discovery than traditional EST analysis. RNA-seq recently has been widely used for transcriptome profiling analysis in many plant species [15, 26, 41, 48–52].
Genes related to metabolic pathway
In the present study, a total of 31,215 unigenes were assigned to 121 KEGG pathways, among which 2,417 unigenes were identified to be associated with carbohydrate metabolisms. A total of 34.8% and 29.2% of all 2,417 unigenes are involved in flavonoid biosynthesis and starch and sucrose metabolism, respectively. Invertase (25 unigenes), SS (24 unigenes), ADP glucose pyrophosphorylase (24 unigenes), SPS (23 unigenes), and starch synthase (23 unigenes) were ranked top five enzymes in 18 key enzymes (see Additional file 10) involved in starch and sucrose metabolism pathways.
Sugar is one of the most important biochemical components that determine fruit quality. A series of enzymes control sucrose metabolism during fruit development and maturation. Invertase, SS, and SPS are three key enzymes that are deeply involved in fruit sucrose matabolism. For example, invertase (β-D-fructofuranosidase) cleaves sucrose to glucose and fructose irreversible, while soluble acid invertase (AIV) presumably hydrolyzes sucrose to hexose for cell growth and development [53, 54]. A decrease of soluble AIV activity was correlated with a rapid increase of sucrose during fruit maturation in Japanese pear [55, 56]. SS is involved in the movement of sucrose to diverse pathways important for metabolic structure and storage functions in plant cells . It has been reported that SS plays a role in sucrose cleavage rather than sucrose synthesis [58, 59]; however, one report suggested that SS was also involved in sucrose synthesis in mature peach fruit where a high level of sucrose was accumulated . In pear fruit, SS also appears to be involved in sucrose synthesis since SS activity increased along with a sucrose accumulation in pear fruit . Recent research in potato also showed that SS strongly determines the intracellular levels of UDP-glucose, ADP-glucose, and starch, and total yield in potato tubers . SPS also plays a major role in sucrose biosynthesis. Both SPS and SS are two important determinants of sucrose accumulation in Asian pear fruit. In 23 pear cultivars, the activity of SS was closely correlated with sucrose content, while SPS showed a weak correlation .
In the family Rosaceae, sorbitol is a major carbohydrate of translocated photosynthates. Sorbitol is converted into other sugars in the fruit by sorbitol-metabolizing enzymes, in which sorbitol oxidase and NAD-dependent sorbitol dehydrogenase are two major players [62, 63]. Research also revealed that three gene families, sorbitol transport (SOT), sorbitol dehydrogenase (SDH), and sorbitol-6-phosphate dehydrogenase (S6PDH), showed more impacts on sugar metabolism in pear than in non-rosaceous species . Wu et al. (2013)  identified four S6PDH genes through genome sequencing; however, no unigenes coding S6PDH were detected in this study. This may be caused by the difference in gene assembly or annotation. Yamaki and Moriguchi (1989) reported that NAD-dependent sorbitol dehydrogenase that converts sorbitol to fructose showed a high activity throughout the entire fruit developing period in Japanese pear fruit. However, the activity of sorbitol oxidase activity that is about one-tenth of NAD-dependent sorbitol dehydrogenase, was high in immature fruit, but decreased during the fruit expansion period and increased again during the fruit maturation stage . In peach fruit, sorbitol oxidase activity was relatively high, but other sorbitol-related enzymes were barely detectable .
In this study, two important enzymes, ADP glucose pyrophosphorylase (24 unigenes) and starch synthase (23 unigenes) that are involved in starch metabolic pathways were identified. In plants, regulation of starch metabolism is complex. In synthesis of starch, ADP-glucose was the glucosyl donor for the elongation of α-1, 4-glucosidic chains . The first committed step is ADP-glucose synthesis that is catalyzed by ADP-glucose pyrophosphorylase (ADPGlc PPase) . As a glucose donor, ADPG molecule is transferred to non-reducing end of the (1–4) glucan primer by starch synthase catalyses, thus formation of a long-chain amylase. Ghosh and Preiss (1966) showed that the reaction catalyzed by ADPGlc PPase is a regulation step in starch synthesis in higher plants as well as in the cyanobacteria [68–70]. Most of the plant ADPGlc PPases are allosterically regulated by 3-PGA and inorganic orthophosphate (Pi) [69, 70]. Some reports also suggested that in higher plants the enzyme activity can also be regulated by its reductive state [71–73].
In this study, multiple unigenes in the main flavonoid biosynthesis pathways were annotated according to KEGG (Figure 6B). Flavonoids including phenylpropanoids, flavones, flavonols, and anthocyanins are important secondary metabolites that are directly involved in the development of fruit quality, such as color, flavor, and health beneficial ingredients [1, 46, 74]. These compounds are also involved in the formation of undesirable brown pigments in fruits following bruising or cutting and/or storage . Information on phenylpropanoid metabolism in pear is limited. Nishitani et al. (2010) have studied the importance of phenylpropanoid metabolism in pear fruit ripening using oligoarray analysis . In present study, a large number (497) of unigenes involved in phenylpropanoid biosynthesis were detected. Lignin biosynthesis in the phenylpropanoid pathway is one of the important factors for pear fruit quality. Lignin is the primary component of stone cells in pear fruit  and its synthesis has direct influence on formation and content of stone cells, ultimately influencing pear fruit quality . Anthocyanin biosynthesis is essential for fruit coloration. Red coloration is determined by the content and composition of anthocyanins in the fruit skin . Most enzymes in the anthocyanin biosynthetic pathway in pears are well studied . For example, Feng et al. (2010) reported that anthocyanin biosynthesis in pears is regulated by a R2R3-MYB transcription factor PyMYB10 .
Digital gene expression (DGE) at different stages of pear fruits
Analyses of KEGG pathways showed that DEGs were observed in several different pathways including some metabolic pathways, such as biosynthesis of phenylpropanoids, plant hormones, and flavonoids; starch and sucrose metabolism; cysteine and methionine metabolism; biosynthesis of photosynthesis-antenna proteins and isoquinoline alkaloid. These are closely related to fruit development in pear and other species [60, 61, 66, 74, 79].
Comparative analysis showed that dehydration-responsive protein RD22 was significantly up-regulated in pear fruit between FS3 and FS2 and significantly down-regulated between the FS3 and FS2 libraries. This showed that this gene involved in the whole expansion stage of pear fruit, so we infer that it may play important roles in the process of pear fruit growth and development. Although some research reported that dehydration-responsive protein RD22 was related to stress/defense response [80–82]. The role of this gene in pear fruit development remains unknown. AP2 was well known for its roles in floral organ identity and develop. In this study, we found that AP2 domain class transcription factor was significantly down-regulated and up-regulated in the fruit developmental stage of FS3/FS2 and FS5/FS4, respectively, which indicated that the AP2 domain class transcription factor was involved in fruit growth and development in pear. Recent research showed that AP2 was involved in many aspects of fruit development in other species. For example, Chung et al. (2010) proved that a tomato (Solanum lycopersicum) APETALA2/ERF gene (SlAP2a) is a negative regulator of fruit ripening . AP2 genes are also involved in seed mass and yield development via regulation of embryo cell number and cell size . Rohrmann et al. (2011) reported that some genes in the AP2-EREBP family responsive to ethylene also showed the altered expression from the green fruit developmental stage onwards . Ripoll et al. (2011) found that AP2 acts to prevent overgrowth of replum by negatively regulating BP and RPL, two genes that normally act to promote replum formation in Arabidopsis. AP2 also acts to prevent overgrowth of the valve margin by repressing the expression of the valve margin identity gene . These studies indicated that AP2 domain class transcription factor is an important regulatory factor in fruit development. Other genes identified in this study, such as metallothionein-like protein type 3, nonspecific lipid-transfer protein precursor, cytochrome P450, zinc finger protein, prunin 2 precursor, flavonoid 3-hydroxylase gene, CBF/DREB1 transcription factor, anthocyanidin synthase gene, etc. may have their own roles in the process of pear fruit growth and development. These roles need further verification.
In this study, the transcriptome profile of Chinese white pear (P. bretschneideri) was investigated using Solexa/Illumina RNA-seq and DGE deep sequencing technologies. A total of 90,227 unigenes were assembled and 62,077 unigenes were annotated. The results demonstrated the feasibility of using Illumina sequencing-based DGE system for gene expression profiling and provided new directions for functional analysis of genes involved in pear fruit development. These findings provide a substantial contribution to existing sequence resources of pear species and will certainly valuable for elucidation of molecular mechanisms in fruit development and maturation in pear or related species. Therefore, this study clearly evidenced that Illumina sequencing technology could be applied as a rapid and cost-effective method for de novo transcriptome analysis of non-model plant species that lack of prior genome annotation.
Availability of supporting data
The assembled unigenes (larger than 200 bp) were deposited in the Transcriptome Shotgun Assembly Sequence Database (TSA) at NCBI, with accession numbers from JR595427 to JR673747.
The project was funded by the National Natural Science Foundation of China (grants no. 31071756 and 31201583), and the Program for Innovative Research Team in Anhui Agricultural University.
- Bell RL, Quamme HA, Layne REC, Skirvin RM: Pears. Fruit breeding, volume I: tree and tropical fruits. Edited by: Janick J, Moore JN. 1996, UK: John Wiley and Sons Inc, 441-514.Google Scholar
- Bell RL: Pears (Pyrus). Genetic resources of temperate fruit and nut crops I. Edited by: Moore JN, Ballington JRJr. 1990, Wageningen, The Netherlands: International Society for Horticultural Science, 655-697.Google Scholar
- Zhang XD, Andrew CA, Yi Q, Chen LM, Li KZ, Shu Q, Su J: Differential gene expression analysis of Yunnan red pear, Pyrus pyrifolia, during fruit skin coloration. Plant Mol Biol Rep. 2011, 29: 305-314. 10.1007/s11105-010-0231-z.View ArticleGoogle Scholar
- Nishitani C, Terakami S, Sawamura Y, Takada N, Yamamoto T: Development of novel EST-SSR markers derived from Japanese pear (Pyrus pyrifolia). Breed Sci. 2009, 59: 391-400. 10.1270/jsbbs.59.391.View ArticleGoogle Scholar
- Zhang HL, Wang YJ, Zhang CH, Wang XP, Li HE, Xu WR: Isolation, characterization and expression analysis of resistance gene candidates in pear (Pyrus spp). Sci Hortic. 2011, 127: 282-289. 10.1016/j.scienta.2010.10.016.View ArticleGoogle Scholar
- Nahla B, Joseph DP: Identification of European and Asian pears using EST-SSRs from Pyrus. Genet Resour Crop Evol. 2010, 57: 357-370. 10.1007/s10722-009-9474-7.View ArticleGoogle Scholar
- Javier S: Genomic characterization of self-incompatibility ribonucleases (S-RNases) in European pear cultivars and development of PCR detection for 20 alleles. Tree Genetics & Genomes. 2009, 5: 393-405. 10.1007/s11295-008-0194-5.View ArticleGoogle Scholar
- Huang W, Marth G: EagleView: a genome assembly viewer for next generation sequencing technologies. Genome Res. 2008, 18: 1538-1543. 10.1101/gr.076067.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Schuster SC: Next-generation sequencing transforms today’s biology. Nat Methods. 2008, 5: 16-18.View ArticlePubMedGoogle Scholar
- Ansorge WJ: Next-generation DNA sequencing techniques. N Biotechnol. 2009, 25: 195-203. 10.1016/j.nbt.2008.12.009.View ArticlePubMedGoogle Scholar
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10: 57-63. 10.1038/nrg2484.PubMed CentralView ArticlePubMedGoogle Scholar
- Morozova O, Marra MA: Applications of next-generation sequencing technologies in functional genomics. Genomics. 2008, 92: 255-264. 10.1016/j.ygeno.2008.07.001.View ArticlePubMedGoogle Scholar
- Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, Gerstein M, Snyder M: The transcriptional landscape of the yeast genome defined by RNA sequencing. Science. 2008, 320: 1344-1349. 10.1126/science.1158441.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilhelm BT, Marguerat S, Watt S, Schubert F, Wood V, Goodhead I, Penkett CJ, Rogers J, Bähler J: Dynamic repertoire of a eukaryotic transcriptome surveyed at single-nucleotide resolution. Nature. 2008, 453: 1239-1243. 10.1038/nature07002.View ArticlePubMedGoogle Scholar
- Lister R, O’Malley RC, Tonti-Filippini J, Gregory BD, Berry CC, Millar AH, Ecker JA: Highly integrated single-base resolution maps of the epigenome in Arabidopsis. Cell. 2008, 133: 523-536. 10.1016/j.cell.2008.03.029.PubMed CentralView ArticlePubMedGoogle Scholar
- Cloonan N, Forrest ARR, Kolle G, Gardiner BBA, Faulkner GJ, Brown MK, Taylor DF, Steptoe AL, Wani S, Bethel G: Stem cell transcriptome profiling via massive-scale mRNA sequencing. Nat Methods. 2008, 5: 613-619. 10.1038/nmeth.1223.View ArticlePubMedGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5: 621-628. 10.1038/nmeth.1226.View ArticlePubMedGoogle Scholar
- Marioni J, Mason C, Mane S, Stephens M, Gilad Y: RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008, 18: 1509-1517. 10.1101/gr.079558.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Sultan M, Schulz MH, Richard H, Magen A, Klingenhoff A, Scherf M, Seifert M, Borodina T, Soldatov A, Parkhomchuk D, Schmidt D, O’Keeffe S, Haas S, Vingron M, Lehrach H, Yaspo ML: A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science. 2008, 321: 956-960. 10.1126/science.1160342.View ArticlePubMedGoogle Scholar
- Liu GQ, Li WS, Zheng PH, Xu T, Chen LJ, Liu DF, Hussain S, Teng YW: Transcriptomic analysis of ‘Suli’ pear (Pyrus pyrifolia white pear group) buds during the dormancy by RNA-Seq. BMC Genomics. 2012, 13: 700-10.1186/1471-2164-13-700.PubMed CentralView ArticlePubMedGoogle Scholar
- Morin R, Bainbridge M, Fejes A, Hirst M, Krzywinski M, Pugh T, McDonald H, Varhol R, Jones S, Marra M: Profiling the HeLa S3 transcriptome using randomly primed cDNA and massively parallel short-read sequencing. Biotechniques. 2008, 45: 81-94. 10.2144/000112900.View ArticlePubMedGoogle Scholar
- Vera JC, Wheat CW, Fescemyer HW, Frilander MJ, Crawford DL, Hanski I, Marden JH: Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing. Mol Ecol. 2008, 17: 1636-1647. 10.1111/j.1365-294X.2008.03666.x.View ArticlePubMedGoogle Scholar
- Kristiansson E, Asker N, Forlin L, Larsson DGJ: Characterization of the Zoarces viviparus liver transcriptome using massively parallel pyrosequencing. BMC Genomics. 2009, 10: 345-10.1186/1471-2164-10-345.PubMed CentralView ArticlePubMedGoogle Scholar
- Meyer E, Aglyamova GV, Wang S, Buchanan-Carter J, Abrego D, Colbourne JK, Willis BL, Matz MV: Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx. BMC Genomics. 2009, 10: 219-10.1186/1471-2164-10-219.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang XW, Luan JB, Li JM, Bao YY, Zhang CX, Liu SS: De novo characterization of a whitefly transcriptome and analysis its gene expression during development. BMC Genomics. 2010, 11: 400-10.1186/1471-2164-11-400.PubMed CentralView ArticlePubMedGoogle Scholar
- Tang Q, Ma XJ, Mo CM, Wilson IW, Song C, Zhao H, Yang YF, Fu W, Qiu DY: An efficient approach to finding Siraitia grosvenorii triterpene biosynthetic genes by RNA-seq and digital gene expression analysis. BMC Genomics. 2011, 12: 343-10.1186/1471-2164-12-343.PubMed CentralView ArticlePubMedGoogle Scholar
- Xue J, Bao YY, Li BL, Cheng YB, Peng ZY, Liu H, Xu HJ, Zhu ZR, Lou YG, Cheng JA, Zhang CX: Transcriptome analysis of the brown planthopper Nilaparvata lugens. PloS ONE. 2010, 5: e14233-10.1371/journal.pone.0014233.PubMed CentralView ArticlePubMedGoogle Scholar
- Xiang LX, He D, Dong WR, Zhang YW, Shao JZ: Deep sequencing-based transcriptome profiling analysis of bacteria-challenged Lateolabrax japonicus reveals insight into the immune-relevant genes in marine fish. BMC Genomics. 2010, 11: 472-10.1186/1471-2164-11-472.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu J, Zhang YL, Zhang HQ, Huang H, Kevin MF, Lu J: Whole genome wide expression profiles of Vitis amurensis grape responding to downy mildew by using Solexa sequencing technology. BMC Plant Biol. 2010, 10: 234-10.1186/1471-2229-10-234.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu T, Qin ZW, Zhou XY, Feng Z, Du YL: Transcriptome profile analysis of floral sex determination in cucumber. J Plant Physiol. 2010, 167: 905-913. 10.1016/j.jplph.2010.02.004.View ArticlePubMedGoogle Scholar
- Wang QQ, Liu F, Chen XS, Ma XJ, Zeng HQ, Yang ZM: Transcriptome profiling of early developing cotton fiber by deep-sequencing reveals significantly differential expression of genes in a fuzzless/ lintless mutant. Genomics. 2010, 96: 369-376. 10.1016/j.ygeno.2010.08.009.View ArticlePubMedGoogle Scholar
- Fátima RC, Noé DF, Jean PVC, Andrés CH, Nayelli MM, Stefan F: A simple and efficient method for isolating small RNAs from different plant species. Plant Methods. 2011, 7: 4-10. 10.1186/1746-4811-7-4.View ArticleGoogle Scholar
- Carra A, Gambino G, Schubert A: A cetyltrimethylammonium bromide-based method to extract low-molecular-weight RNA from polysaccharide-rich plant tissues. Anal Biochem. 2007, 360: 318-320. 10.1016/j.ab.2006.09.022.View ArticlePubMedGoogle Scholar
- Li R, Zhu H, Ruan J, Qian W, Fang X, Shi Z, Li Y, Li S, Shan G, Kristiansen K, Li S, Yang H, Wang J, Wang J: De novo assembly of human genomes with massively parallel short read sequencing. Genome Res. 2010, 20: 265-272. 10.1101/gr.097261.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, Tsai J, Quackenbush J: TIGR gene indices clustering tools (TGICL):a software system for fast clustering of large EST datasets. Bioinformatics. 2003, 19: 651-652. 10.1093/bioinformatics/btg034.View ArticlePubMedGoogle Scholar
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO:a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMedGoogle Scholar
- Ye J, Fang L, Zheng HK, Zhang Y, Chen J, Zhang ZJ, Wang J, Li ST, Li RQ, Bolund L, Wang J: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34: W293-W297. 10.1093/nar/gkl031.PubMed CentralView ArticlePubMedGoogle Scholar
- Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7: 986-995.PubMedGoogle Scholar
- Shi CY, Yang H, Wei CL, Yu O, Zhang ZZ, Jiang CJ, Sun J, Li YY, Chen Q, Xia T, Wan XC: Deep sequencing of the Camellia sinensis transcriptome revealed candidate genes for major metabolic pathways of tea-specific compounds. BMC Genomics. 2011, 12: 131-139. 10.1186/1471-2164-12-131.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang X, Sun Q, McGrath SD, Mardis ER, Soloway PD, Clark AG: Transcriptome-wide identification of novel imprinted genes in neonatal mouse brain. PLoS ONE. 2008, 3: e3839-10.1371/journal.pone.0003839.PubMed CentralView ArticlePubMedGoogle Scholar
- Feng C, Chen M, Xu CJ, Bai L, Yin XR, Li X, Allan AC, Ferguson IB, Chen KS: Transcriptomic analysis of Chinese bayberry (Myrica rubra) fruit development and ripening using RNA-Seq. BMC Genomics. 2012, 13: 19-10.1186/1471-2164-13-19.PubMed CentralView ArticlePubMedGoogle Scholar
- Hunter S, Apweiler R, Attwood TK, Bairoch A, Bateman A, Binns D, Bork P, Das U, Daugherty L, Duquenne L, Finn RD, Gough J, Haft D, Hulo N, Kahn D, Kelly E, Laugraud A, Letunic I, Lonsdale D, Lopez R, Madera M, Maslen J, McAnulla C, McDowall J, Mistry J, Mitchell A, Mulder N, Natale D, Orengo C, Quinn AF, Selengut JD, Sigrist CJ, Thimma M, Thomas PD, Valentin F, Wilson D, Wu CH, Yeats C: InterPro: the integrative protein signature database. Nucleic Acids Res. 2009, 37: D211-D215. 10.1093/nar/gkn785.PubMed CentralView ArticlePubMedGoogle Scholar
- Sompompailin K, Makita Y, Saito K: A WD-repeat-containing putative regulatory protein in anthocyanin biosynthesis in Perilla frutescens. Plant Mol Biol. 2002, 50: 485-495. 10.1023/A:1019850921627.View ArticleGoogle Scholar
- Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K, Holm L, Sonnhammer ELL, Eddy SR, Bateman A: The Pfam protein families database. Nucleic Acids Res. 2010, 38: D211-D222. 10.1093/nar/gkp985.PubMed CentralView ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M: The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004, 32: D277-D280. 10.1093/nar/gkh063.PubMed CentralView ArticlePubMedGoogle Scholar
- Bell RL, Janick J: Quantitative genetic analysis of fruit quality in pear. J Amer Soc Hort Sci. 1990, 115: 829-834.Google Scholar
- Yamamoto T, Chevreau E: Pear genomics. Plant genetics/ genomics, volume 6: genetics and genomics of rosaceae. Edited by: Kevin MF, Susan EG. 2009, LLC, USA: Springer Science + Business Media, 163-186.View ArticleGoogle Scholar
- Xie FL, Burklew CE, Yang YF, Liu M, Xiao P, Zhang BH, Qiu DY: De novo sequencing and a comprehensive analysis of purple sweet potato (Impomoea batatas L.) transcriptome. Planta. 2012, 236: 101-113. 10.1007/s00425-012-1591-4.View ArticlePubMedGoogle Scholar
- Jin H, Martin C: Multifunctionality and diversity within the plant MYB-gene family. Plant Mol Biol. 1999, 41: 577-585. 10.1023/A:1006319732410.View ArticlePubMedGoogle Scholar
- Xu H, Gao Y, Wang JB: Transcriptomic analysis of rice (Oryza sativa) developing embryos using the RNA-Seq technique. PLoS ONE. 2012, 7: e30646-10.1371/journal.pone.0030646.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang YG, Zhu J, Dai HY: Characterization of transcriptional differences between columnar and standard apple trees using RNA-Seq. Plant Mol Biol Rep. 2012, 30: 957-965. 10.1007/s11105-011-0396-0.View ArticleGoogle Scholar
- Mizrachi E, Hefer CA, Ranik M, Joubert F, Myburg AA: De novo assembled expressed gene catalog of a fast-growing Eucalyptus tree produced by Illumina mRNA-Seq. BMC Genomics. 2010, 11: 681-10.1186/1471-2164-11-681.PubMed CentralView ArticlePubMedGoogle Scholar
- Pfeiffer I, Kutschera U: Sucrose metabolism and cell elongation in developing sunflower hypocotyls. J Exp Bot. 1995, 46: 631-638. 10.1093/jxb/46.6.631.View ArticleGoogle Scholar
- Tang X, Ruffner HP, Scholes JD, Rolfe SA: Purification and characterisation of soluble invertases from leaves of Arabidopsis thaliana. Planta. 1996, 198: 17-23.View ArticlePubMedGoogle Scholar
- Moriguchi T, Abe K, Sanada T, Yamaki S: Levels and role of suscrose synthase, sucrose-phosphate synthase, and acid invertase in sucrose accumulation in fruit of Asian pear. J Am Soc Hort Sci. 1992, 117: 274-278.Google Scholar
- Tanase K, Yamaki S: Sucrose synthase isozymes related to sucrose accumulation during fruit development of Japanese pear (Pyrus pyrifolia Nakai). J Japan Soc Hort Sci. 2000, 69: 671-676. 10.2503/jjshs.69.671.View ArticleGoogle Scholar
- Komatsu A, Moriguchi T, Koyama K, Omura M, Akihama T: Analysis of sucrose synthase genes in citrus suggests different roles and phylogenetic relationships. J Exp Bot. 2002, 53: 61-71. 10.1093/jexbot/53.366.61.View ArticlePubMedGoogle Scholar
- Hawker JS: Sucrose. Biochemistry of storage carbohydrates in green plants. Edited by: Dey PM, Dixon RA. 1985, New York, USA: Academic, 1-51.Google Scholar
- Turner JF, Turner DH: The regulation of carbohydrate metabolism. Annu Rev Plant Physiol. 1975, 26: 159-186. 10.1146/annurev.pp.26.060175.001111.View ArticleGoogle Scholar
- Moriguchi T, Sanada T, Yamakil S: Seasonal fluctuations of some enzymes relating to sucrose and sorbitol metabolism in peach fruit. J Amer Soc Hort Sci. 1990, 115: 278-281.Google Scholar
- Baroja-Fernández E, Munoz FJ, Montero M, Etxeberria E, Sesma MT, Ovecka M, Bahaji A, Ezquer I, Li J, Prat S, Pozueta-Romero J: Enhancing sucrose synthase activity in transgenic potato (Solanum tuberosum L.) tubers results in increased levels of starch, ADPglucose and UDPglucose and total yield. Plant Cell Physiol. 2009, 50: 1651-1662. 10.1093/pcp/pcp108.View ArticlePubMedGoogle Scholar
- Yamaki S, Ishikawa K: Roles of four sorbitol related enzymes and invertase in the seasonal alteration of sugar metabolism in apple tissue. J Amer Soc Hort Sci. 1986, 111: 134-137.Google Scholar
- Yamaki S, Moriguchi T: Seasonal fluctuation of sorbitol related enzymes and invertase activities accompanying maturation of Japanese pear (Pyrus serotina Rehder var. cults Rehder) fruit. J Jpn Soc Hort Sci. 1989, 57: 602-607.View ArticleGoogle Scholar
- Wu J, Wang Z, Shi Z, Zhang S, Ming R, Zhu S, Khan MA, Tao S, Korban SS, Wang H, Chen NJ, Nishio T, et al: The genome of pear (Pyrus bretschneideri Rehd). Genome Res. 2012, doi:10.1101/gr.144311.112Google Scholar
- Moriguchi T, Sanada T, Yamaki S: Seasonal fluctuations of some enzymes relating to sucrose and sorbitol metabolism in peach fruit. J Amer Hort Sci. 1990, 115: 278-281.Google Scholar
- Ballicora MA, Iglesias AA, Preiss J: ADP-glucose pyrophosphorylase: a regulatory enzyme for plant starch synthesis. Photosynth Res. 2004, 79: 1-24.View ArticlePubMedGoogle Scholar
- Martin C, Smith AM: Starch biosynthesis. Plant Cell. 1995, 7: 971-985.PubMed CentralView ArticlePubMedGoogle Scholar
- Ghosh HP, Preiss J: Adenosine diphosphate glucose pyrophosphorylase: a regulatory enzyme in the biosynthesis of starch in spinach leaf chloroplasts. J Biol Chem. 1966, 241: 4491-4504.PubMedGoogle Scholar
- Preiss J, Sivak MN: Biochemistry, molecular biology and regulation of starch synthesis. Genet Eng. 1998, 20: 177-223.View ArticleGoogle Scholar
- Preiss J: Biosynthesis of bacterial and mammalian glycogen and plant starch synthesis and their regulation. Biorganic chemistry: carbohydrates. Edited by: Hecht SM. 1999, Oxford: Oxford University Press, 489-554.Google Scholar
- Fu Y, Ballicora MA, Leykam JF, Preiss J: Mechanism of reductive activation of potato tuber ADP-glucose pyrophosphorylase. J Biol Chem. 1998, 273: 25045-25052. 10.1074/jbc.273.39.25045.View ArticlePubMedGoogle Scholar
- Ballicora MA, Frueauf JB, Fu Y, Schurmann P, Preiss J: Activation of the potato tuber ADP-glucose pyrophosphorylase by thioredoxin. J Biol Chem. 2000, 275: 1315-1320. 10.1074/jbc.275.2.1315.View ArticlePubMedGoogle Scholar
- Tiessen A, Hendriks JH, Stitt M, Branscheid A, Gibon Y, Farre EM, Geigenberger P: Starch synthesis in potato tubers is regulated by post-translational redox modification of ADP-glucose pyrophosphorylase: a novel regulatory mechanism linking starch synthesis to the sucrose supply. Plant Cell. 2002, 14: 2191-2213. 10.1105/tpc.003640.PubMed CentralView ArticlePubMedGoogle Scholar
- Fischer TC, Gosch C, Pfeiffer J, Halbwirth H, Halle C, Stich K, Forkmann G: Flavonoid genes of pear (Pyrus communis). Trees. 2007, 21: 521-529. 10.1007/s00468-007-0145-z.View ArticleGoogle Scholar
- Nishitani C, Shimizu T, Fujii H, Hosaka F, Terakami S, Nakamura Y, Itai A, Yamaguchi-Nakamura A, Yamamoto T: Oligoarray analysis of gene expression in ripening Japanese pear fruit. Sci Hortic. 2010, 124: 195-203. 10.1016/j.scienta.2009.12.032.View ArticleGoogle Scholar
- 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. 10.1016/j.plantsci.2008.12.011.View ArticleGoogle Scholar
- Steyn WJ, Wand SJE, Holcroft DM, Jacobs G: Red colour development and loss in pears. Acta Hortic. 2005, 671: 79-85.View ArticleGoogle Scholar
- Feng SQ, Wang YL, Yang S, Xu YT, Chen XS: Anthocyanin biosynthesis in pears is regulated by a R2R3-MYB transcription factor PyMYB10. Planta. 2010, 232: 245-255. 10.1007/s00425-010-1170-5.View ArticlePubMedGoogle Scholar
- Baur AH, Yang SF: Methionine metabolism in apple tissue in relation to ethylene biosynthesis. Phytochemistry. 1972, 11: 3207-3214. 10.1016/S0031-9422(00)86375-X.View ArticleGoogle Scholar
- Fernandez L, Torregrosa L, Terrier N, Sreekantan L, Grimplet J, Davies C, Thomas MR, Romieu C, Ageorges A: Identification of genes associated with flesh morphogenesis during grapevine fruit development. Plant Mol Biol. 2007, 63: 307-323. 10.1007/s11103-006-9090-2.View ArticlePubMedGoogle Scholar
- Ziosi V, Bonghi C, Bregoli AM, Trainotti L, Biondi S, Sutthiwal S, Kondo S, Costa G, Torrigiani P: Jasmonate-induced transcriptional changes suggest a negative interference with the ripening syndrome in peach fruit. J Exp Bot. 2008, 59: 563-573. 10.1093/jxb/erm331.View ArticlePubMedGoogle Scholar
- Dhanaraj AL, Alkharouf NW, Beard HS, Chouikha IB, Matthews BF, Wei H, Arora R, Rowland LJ: Major diVerences observed in transcript proWles of blueberry during cold acclimation under weld and cold room conditions. Planta. 2007, 225: 735-751. 10.1007/s00425-006-0382-1.View ArticlePubMedGoogle Scholar
- Chung MY, Vrebalov J, Alba R, Lee J, McQuinn R, Chung JD, Klein P, Giovannoni J: A tomato (Solanum lycopersicum) APETALA2/ERF gene, SlAP2a, is a negative regulator of fruit ripening. Plant J. 2010, 64: 936-947. 10.1111/j.1365-313X.2010.04384.x.View ArticlePubMedGoogle Scholar
- Ohto M, Fischer RL, Goldberg RB, Nakamura K, Harada JJ: Control of seed mass by APETALA2. Proc Natl Acad Sci. 2005, 102: 3123-3128. 10.1073/pnas.0409858102.PubMed CentralView ArticlePubMedGoogle Scholar
- Rohrmann J, Tohge T, Alba R, Osorio S, Caldana C, McQuinn R, Arvidsson S, Merwe MJ, Rianõ-Pachón DM, Mueller-Roeber B, Fei Z, Nesi AN, Giovannoni JJ, Fernie AR: Combined transcription factor profiling, microarray analysis and metabolite profiling reveals the transcriptional control of metabolic shifts occurring during tomato fruit development. Plant J. 2011, 68: 999-1013. 10.1111/j.1365-313X.2011.04750.x.View ArticlePubMedGoogle Scholar
- Ripoll JJ, Roeder AHK, Ditta GS, Yanofsky MF: A novel role for the floral homeotic gene APETALA2 during Arabidopsis fruit development. Development. 2011, 138: 5167-5176. 10.1242/dev.073031.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.