Skip to main content

Advertisement

Identifying the candidate genes involved in the calyx abscission process of 'Kuerlexiangli’ (Pyrus sinkiangensis Yu) by digital transcript abundance measurements

Article metrics

Abstract

Background

'Kuerlexiangli’ (Pyrus sinkiangensis Yu), a native pear of Xinjiang, China, is an important agricultural fruit and primary export to the international market. However, fruit with persistent calyxes affect fruit shape and quality. Although several studies have looked into the physiological aspects of the calyx abscission process, the underlying molecular mechanisms remain unknown. In order to better understand the molecular basis of the process of calyx abscission, materials at three critical stages of regulation, with 6000 × Flusilazole plus 300 × PBO treatment (calyx abscising treatment) and 50 mg.L-1GA3 treatment (calyx persisting treatment), were collected and cDNA fragments were sequenced using digital transcript abundance measurements to identify candidate genes.

Results

Digital transcript abundance measurements was performed using high-throughput Illumina GAII sequencing on seven samples that were collected at three important stages of the calyx abscission process with chemical agent treatments promoting calyx abscission and persistence. Altogether more than 251,123,845 high quality reads were obtained with approximately 8.0 M raw data for each library. The values of 69.85%-71.90% of clean data in the digital transcript abundance measurements could be mapped to the pear genome database. There were 12,054 differentially expressed genes having Gene Ontology (GO) terms and associating with 251 Kyoto Encyclopedia of Genes and Genomes (KEGG) defined pathways. The differentially expressed genes correlated with calyx abscission were mainly involved in photosynthesis, plant hormone signal transduction, cell wall modification, transcriptional regulation, and carbohydrate metabolism. Furthermore, candidate calyx abscission-specific genes, e.g. Inflorescence deficient in abscission gene, were identified. Quantitative real-time PCR was used to confirm the digital transcript abundance measurements results.

Conclusions

We identified candidate genes that showed highly dynamic changes in expression during the calyx abscission process. These genes are potential targets for future functional characterization and should be valuable for exploration of the mechanisms of calyx abscission, and eventually for developing methods based on small molecule application to induce calyx abscission in fruit production.

Background

The 'Kuerlexiangli’ pear is one of the characteristic and economically important fruit trees in Xinjiang Uygur Autonomous Region, China. Fruit of this cultivar is a rich source of juice and has few stone cells, good flavor, and tantalizing aroma. However, a persistent calyx seriously affects the fruit shape (Additional file1) as well as quality. Fruit with calyx abscission have higher titer of soluble sugars and vitamin C, lower titer of titratable acids compared with persistent calyx fruit[1]. Thus calyx persistence significantly affects the commercial importance of 'Kuerlexiangli’ and causes serious loss in economic value.

Previously, few studies have focused on calyx abscission of pear, which is clearly regulated by ethylene, whilst auxins apparently decrease the ratio of calyx abscission[2, 3]. On the other hand, the rate of calyx abscission in 'Kuerlexiangli’ fruit is varied among different pollinizer varieties, with higher calyx abscission rates obtained when flowers were pollinated with 'Xueqing’, 'Yali’, 'Zhongliyihao’ and 'Cuiguan’ pear[4]. Moreover, manually excising the calyx was helpful for enhancing the comprehensive quality of pear fruit[5].

Recent research has mainly focused on the physiological aspects of the calyx development process, with emphasis on anatomy, and the role of nutrition and hormone control. Microstructure data of calyx tubes at the calyx developing stage after spraying 300 mg kg-1 paclobutrazol (PP333) (calyx abscising treatment) or 50 mg kg-1 GA3 (calyx persisting treatment) showed that the average area of vascular bundle of calyx tube tissue was bigger and many sieve tube cells and idioblasts gradually appeared when the tree was treated with GA3. When the tree was treated with PP333, calyx tubes of young fruit only had vessels in vascular bundle. Abscission layer appeared at the late young fruit of calyx tube developing stage, and finally the calyx tube broke off and young fruit developed as calyx fruit[6]. The physiological evidence indicates a role of both nutrition and hormone regulation in calyx abscission. Sufficient Fe content is a guarantee for fruit calyx development[7]. High content of zeatin-riboside (ZR) and indoleacetic acid (IAA) in young fruit, high content of IAA and GA3 but low content of abscisic acid (ABA) in calyx, high ratio of IAA, ZR, GA3 between calyx and young fruit, and low ratio of ABA between calyx and young fruit all promoted calyx persistence[8]. However, the key genes that control or regulate fruit calyx abscission or persistence are still unknown. Thus the understanding of gene expression and regulation differences under calyx abscission and calyx persistence status and the responses to different chemical agents leading to different calyx status is necessary.

Next generation sequencing methods have emerged as a cost-effective high-throughput approach to the sequencing of a very large number of expressed genes even in small experiments. One of these, digital transcript abundance measurements method is a revolutionary approach for expression analysis competing to replace microarrays for analyzing transcriptome[9]. It is tag-based transcriptome sequencing for measuring relative gene expression levels, which can identify, quantify, and annotate expressed genes on the whole genome level with or even without prior sequence knowledge. It enables an entirely new scale of biological experimentation to reveal related pathways or identify target genes involved in different bio-processes. For example, digital transcript abundance measurements approach has been used to study gene expression in the poplar (Populus simonii × Populus nigra) under salt stress[10]. To identify the candidate genes for sex determination of papaya (Carica papaya), papaya male, female, and hermaphrodite plants were used for digital transcriptome analysis using high-throughput serial analysis of gene expression[11]. The molecular regulation mechanism of the physiological and biochemical response to potassium starvation in soybean roots and shoots was investigated by high-throughput tag-sequencing[12]. Differentially expressed genes in cucumber (Cucumis sativus L.) root under waterlogging stress have also been identified by digital transcript abundance measurements[13]. Overall, the digital transcript abundance measurements method has provided more valuable tools for qualitative and quantitative gene expression analysis than the earlier microarray based assays[14].

Most current molecular knowledge on the abscission process has been obtained from the model plant Arabidopsis thaliana[15], as well as tomato (Solanum lycopersicum)[16]. The apple cluster during immature fruit drop also represents an ideal system to study the shedding of actively growing organs[17]. Recently developed molecular approaches have been used in abscission process in horticulture crops. Previous studies have identified the transcriptomes associated with flower abscission in tomato[18], differential gene expression by the cDNA microarray technique during abscission of citrus leaves under ethylene treatment[19] and fruitlet abscission in apple[17]. The differentially expressed genes have been identified during shedding of immature apple fruits have been identified with a cDNA-AFLP approach[20], and the role of NAA and shading in apple fruit abscission have been analyzed by the transcriptome method[21]. However, a deep knowledge of the molecular events occurring during the early phases of calyx abscission induction is still lacking. The new release of the whole genome sequence of Asian pear 'Dangshansuli’[22], lays a good platform for genome-wide gene analysis. Here we report on the first use of genome-wide analysis to gain insight into the wide range of transcriptional responses associated with calyx abscission processes. Using Solexa/Illumina’s sequencing system, the transcriptomes were compared between chemical regulation of calyx abscission, 6000 × Flusilazole plus 300 × PBO inducing calyx abscission, and 50 mg.L-1 GA3 treatment to reduce calyx abscission. By investigating the expression of genes related to calyx abscission in 'Kuerlexiangli’, a number of pathways and candidate genes that are important in this process were identified.

Results and discussion

Effects of different treatments on calyx abscission rate

Comparison of Flusilazole treatment and GA3 treatment as inducer/inhibitor of fruit abscission revealed significant differences in abscission rates (Figure 1). At 22 d after treatments, the rate of calyx abscission in the untreated control was 16.78%, and the Flusilazole treatment increased the calyx abscission rate to 91.25%, but the GA3 treatment decreased the calyx abscission ratio to 1.38%. Thus the application of Flusilazole treatment increased the calyx abscission rate by 4.4 times whereas GA3 treatment decreased the calyx abscission rate to one fourteenth.

Figure 1
figure1

Effect of Flusilazole and GA 3 treatments on calyx abscission of 'Kuerlexiangli’ calyx. 6000 × Flusilazole plus 300 × PBO increases the rate of calyx abscission, whereas GA3 decreases calyx abscission. Results represent the mean (±SD) of three replicates. Different letters indicate significant differences (P ≤ 0.05, Student t-test).

Digital transcript abundance measurements libraries sequencing

Illumina sequencing analysis was performed to obtain a global view of the calyx abscission of 'Kuerlexiangli’ through analysis of the transcriptomes of Flusilazole treatment and GA3 treatment. Seven samples from C1 to C7 were used as material for digital transcript abundance measurements analysis (Table 1). Calyx digital transcript abundance measurements libraries (Table 2) were deep sequenced altogether. The number of reads for each library ranged from 7.1 to 9.9 million. A total of 8,132,745, 7,613,791, 7,370,972, 9,924,450, 7,174,967, 8,179,209, and 8,050,331 reads of raw data were obtained for C1, C2, C3, C4, C5, C6, and C7 libraries, respectively. In this study, the tag sequences of the seven digital transcript abundance measurements libraries were mapped to the assembled pear genome of 'Dangshansuli’[22]. Finally, this generated 5,766,027, 5,381,267, 5,166,155, 6,932,208, 5,158,694, 5,653,121, and 5,511,540 reads corresponding to 70.90%, 70.68%, 70.09%, 69.85%, 71.90%, 69.12%, and 68.46% of all clean data in the seven libraries that could be mapped to the reference database. The Q30 percentages (A Q-score of 30 corresponds to an error rate of 1 per 1000) of all seven libraries were above 95%. All of these data showed that the throughput and sequencing quality was high enough for further analysis.

Table 1 Description of calyx tube abscission zone samples used in this study
Table 2 Statistics of digital transcript abundance measurements library sequencing and tag mapping

Analysis of differential gene expression

To compare differential expression patterns among seven libraries, we employed IDEG6 (a web tool for detection of differentially expressed genes in multiple tag sampling experiments) to identify mRNAs showing statistically significant differences based on their relative abundance, as reflected by total count of individual sequence reads, between all pairs of libraries. We compared the six test libraries with the control, two treatments during the same stage, and the same treatment during different stages, so that 15 pairs of comparisons were implemented. Among these comparisons, we found that 1,760 to 8,184 genes had significant changes in expression, and the average number was 5,994 (Figure 2). The differential expression pattern among libraries revealed that the largest differences occurred between C3 and C5. A total of 8,184 genes were significantly differentially expressed between C3 and C5. Of these genes, 2,353 genes were up-regulated and 5,831 were down-regulated in C3 compared with C5. The smallest difference was observed between C6 and C7, in which only 1,760 differentially expressed genes (DEGs) were identified.

Figure 2
figure2

Genes differentially expressed between different libraries. Up-regulated (red) and down-regulated (green) genes were quantified. The results of 15 comparisons are shown.

We also detected a large number of specifically expressed genes (SEGs) between each pair of libraries (Figure 3). Comparisons showed that there were 983 to 2,563 SEGs with an average number of 1,732 among 15 comparisons. These SEGs will help us to find genes correlated with calyx abscission process. Especially a number of stage-specific and treatment-specific expressed genes are likely to be key genes associated with calyx abscission.

Figure 3
figure3

Quantification of genes specifically expressed between different samples. The numbers of digital transcript abundance measurements genes of 15 comparisons between each two samples are shown in red and blue histograms (top vs. bottom).

Functional annotation of differentially expressed genes

Gene Ontology (GO) is an international standardized gene function classification system that describes properties of genes and their products in any organism. In this study, a total of 12,054 differentially expressed genes that could be categorized into 41 functional groups were found (Figure 4). The major subcategories were as follows: four subcategories for cellular component ('cell’, 'cell part’, 'organelle’, and 'macromolecular complex’); three subcategories for molecular function ('binding’, 'catalytic activity’, and 'transporter activity’); and seven subcategories for biological process ('metabolic process’, 'cellular process’, 'biological regulation’, 'pigmentation’, 'localization’, 'establishment of location’, and 'response to stimulus’). These results indicate that expressed genes functioning in 'binding’, 'catalytic activity’, 'metabolic process’, and 'cellular process’ are important during the calyx abscission process. Only a few genes were clustered in terms of 'synapse’, 'synapse part’, 'metallochaperone activity’, 'biological adhesion’ and 'immune system process’ and 'viral reproduction’.

Figure 4
figure4

GO categories of the genes identified. Y-axis (left) represents percentages of genes identified in this study, Y-axis (right) represent the actual gene number. The genes were annotated in three main categories: biological process, cellular component, and molecular function (X-axis).

To further investigate the function of different expressed genes during calyx abscission, significantly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were identified according to the P values and enrichment factor. Performing a BLAST search against the KEGG database indicated that expressed genes were involved in 251 pathways (Additional file2). As shown in Additional file3, nine KEGG pathways were observed to be significantly overrepresented in calyx abscission processes. Those genes correlated with calyx abscission mainly involved in photosynthesis, plant hormone signal transduction, cell wall modification, transcriptional regulation and carbohydrate metabolism were used for subsequent analysis. These trends were consistent with overall developmental activities during abscission processes[23]. In addition, several other biological processes that have not previously been reported to be associated with calyx abscission, such as flavonoid biosynthesis and flavone/flavonol biosynthesis, were found dramatically changed during calyx abscission processes. These might be novel genes that are relevant to the calyx abscission process in 'Kuerlexiangli’ fruit.

Cluster of calyx abscission-related genes

Impacts on photosynthesis

Decrease in photosynthesis may be an important contributing factor for the abscission of flowers and fruitlets in the abscission processes[24], which was confirmed by our present experiment. Our results showed that 230 genes encoding photosynthesis-related genes were differential expression under different treatments (Additional file4). Altered expressions were found for numerous genes involved in carbon fixation, photosystem I, photosystem II, and photosynthetic electron transport. Flusilazole treatment decreased the expression of genes involved with PSI, PSII activity and photosynthetic electron transport in photosystem. The down-regulated genes included those that encode subunit V, PsaA/PsaB, F-type H+-transporting ATPase subunit gamma in photosystem I, and 13 kDa protein in photosystem II (Table 3). We identified that GA3 treatment led to strong induction in photosynthesis-related gene expression (Additional file4). These results indicate that photosynthesis might play an important role in the calyx abscission and samples with diminished photosynthesis function are more prone to abscission.

Table 3 Selected genes with altered expression during calyx abscission process

Plant hormone signal transduction

Many kinds of hormones regulate the process of abscission, such as IAA, ABA, GA, JA, and ethylene, among which IAA and ethylene play a pivotal role. IAA prevents, while ethylene accelerates the abscission processes[25]. Forty-three genes were found to be involved in ethylene synthesis, perception, and response in the present study. ACC synthase (ACS) and ACC oxidase (ACO) are rate-limiting enzymes in ethylene biosynthesis[26, 27]. Genes involved in ACC synthase (Pbr002233.1, Pbr007947.1, Pbr032688.1, Pbr001095.1, Pbr029891.1, Pbr030234.1) showed expressions increased by 2.109 to 6.236 times during the early calyx abscission process after Flusilazole treatment. Two ACC oxidase genes (Pbr005179.1 and Pbr031954.1) were up-regulated in later calyx abscission process after Flusilazole treatment. On the other hand, two ethylene receptor genes (Pbr002199.1 and Pbr023072.1) and three ethylene responsive genes (Pbr000435.1, Pbr025988.1, and Pbr030542.1) were up-regulated in later calyx abscission process after Flusilazole treatment (Additional file4). This is consistent with the correlation between abscission and increased expression of genes for ethylene synthesis and ethylene receptors in AZ (abscission zone), which has been reported in mature fruit of olive[28] and apple[21]. The general rule states that provided the flux of IAA to the abscission zone is maintained, cell separation is inhibited and abscission does not happen[29]. One hundred and two auxin-related genes were differentially expressed by the Flusilazole and GA3 treatments, including indole-3-acetic acid-amido synthetase, indole-3-acetic acid-induced protein, and auxin response factor. Nine genes encoding AUX/IAA protein were down-regulated at 6 d after Flusilazole treatment and 13 genes encoding auxin responsive factors were also down-regulated at 10 d after Flusilazole treatment. Genes related to polar auxin transporting were also affected by GA3 treatment, with auxin influx carriers induced and efflux carriers largely repressed. We found that an auxin transport-related gene (Pbr013424.1) was up-regulated at 6 d, but repressed at 10 d after Flusilazole treatment. In addition, a gene (Pbr011705.1) encoding spermidine synthase was up-regulated with Flusilazole treatment at onset of calyx abscission processes. As mentioned above, hormones seem to play a relatively important role during the early phases of abscission in the calyx, since a majority of the transcriptionally activated elements involved in hormone signaling seem to be downstream of the induction of abscission. It is especially interesting that expression of more genes associated with ethylene and auxin metabolism and signal transduction were altered during calyx expression process than previously reported.

Cell wall degradation-related genes are highly expressed in calyx process

Abscission involves progressive dissolution of the middle lamella. A key step in loss of adhesion between cells within a separation layer is the induction of cell wall-degrading enzymes such as polygalacturonases. The role of other wall-modifying proteins such as expansin, xyloglucan transglucosylase hydrolase, pectinesterase, etc. has also been studied during abscission processes[30]. Polygalacturonases have been studied in tomato abscission zones[31], in oilseed rape and Arabidopsis leaf and flower abscission zones[32], while celluloses have been studied in cotton[33], citrus fruit[34] and red raspberry abscission zones[35]. Reports have indicated that an increase in expansin[36], xyloglucan transglucosylase hydrolases[37] and pectate lyases[38] correlate with organ abscission. A large number of genes encoding cell wall hydrolases, expansin, and lyases were found to be over-represented with Flusilazole treatment (Additional file4). In this study, a gene (Pbr036382.1) encoding polygalacturonase and a gene (Pbr014900.1) encoding pectinesterase were up-regulated at 22 d after Flusilazole treatment. We also observed expansin expression in calyx abscission process. Three genes (Pbr000187.1, Pbr042842.1, and Pbr040203.1) encoding xyloglucan endotransglucosylase/hydrolase were up-regulated at 22 d after Flusilazole treatment (Table 3). Pectate lyases (PLs,EC:4.2.2.2) could facilitate cell wall disassembly by their action on demethylated pectin through beta elimination in the presence of Ca2+, leading to depolymerization of pectic polysaccharides[39]. Some candidate genes related with this reaction identified in this study may play roles in cell wall degradation, such as Pbr011950.1, Pbr021157.1, Pbr035305.1, and Pbr041275.2, all of which probably aid in later abscission processes.

Genes encoding transcription factors

Transcription factors (TFs) act as major switches of regulatory cascades during development, and alterations in the expression of such genes may affect various developmental processes. Our digital transcript abundance measurements results showed that 198 genes encoding transcription factors were differentially expressed, including ERF/APETALA2, MYB, and WRKY transcription factors, and MADS-box proteins (Additional file4). APETALA2 (AP2) is a floral homeotic factor that plays an important role in the control of Arabidopsis flower and seed development and encodes a putative TF that is distinguished by a novel DNA-binding motif referred to as the AP2 domain[40]. The expression of a gene (Pbr029902.1) encoding for a TF containing an AP2 domain was transiently up-regulated specifically at a later stage under GA3 treatment. Multiple up-regulated MYB factors at 22 d after Flusilazole treatment were identified, including MYB 105, 160, 220, 219, 12, and 190. Moreover, MYB 12 is a flavonol-specific activator of flavonoid biosynthesis, strongly triggering the promoters of genes encoding chalcone synthase (CHS), flavanone-3-hydroxylase (F3H), and flavonol synthase (FLS), which are all involved in the biosynthesis of flavonols[41, 42]. The up-regulation of MBY 12 (Pbr022669.1) consistently resulted in an enhanced transcript level of CHS at 22 d after Flusilazole treatment. Interestingly, auxin and ethylene induce flavonol accumulation through partly identical transcription networks[43]. We hypothesize that the calyx abscission zone contains higher levels of flavonols and flavonoids, which might have functions in stress defense after calyx abscission. The WRKY family is a superfamily of TFs, which hold central positions mediating fast positive and negative regulation of disease resistance[44]. Two WRKY genes were up-regulated transiently at 6 d after Flusilazole treatment (Table 3), which showed similar change pattern to that reported in tomato flower abscission zone[18]. On the whole, all above results suggest that many transcription factors play an important role in regulating calyx abscission processes in pear.

Genes involved in carbohydrate metabolism

A strong connection between the carbohydrate amounts available for the fruitlet, especially soluble sugars, and their probability of abscission has been suggested[45]. This phenomenon has been also described for pistachio[46]. So abscission may be due to a lack of carbohydrates. In our digital transcript abundance measurements results, 486 carbohydrate metabolism genes had significantly altered expression profiles in response to Flusilazole treatment and GA3 treatment (Additional file4). Affected genes within this group include those associated with glycolysis/gluconeogenesis, fatty acid biosynthesis, and sucrose metabolic processes. For example, we found induction of alpha-amylase (EC: 3.2.1.1), beta-fructofuranosidase (EC: 3.2.1.26), alcohol dehydrogenase (EC: 1.1.1.1), and acetyl-CoA carboxylase (EC: 6.4.1.2) under Flusilazole treatment. A different gene set showed lower expression in calyx abscission. This set comprised genes coding for 2-isopropylmalate synthase (EC: 2.3.3.13), pyruvate decarboxylase (EC: 4.1.1.1), sorbitol dehydrogenase (EC: 1.1.1.14) and ribulose-1,5-bisphosphate carboxylase (EC: 4.1.1.39). Some of the genes are clearly connected to sugar mobilization, such as alpha-amylase and glycosyl hydrolase, which were up-regulated at 6 d after Flusilazole treatment. The gene (Pbr001279.1) encoding UDP-glucosyltransferases (EC2.4.1.115) was highly expressed under Flusilazole treatment in the abscission process (Table 3). In apple, UDP-glucosyltransferases were also expressed during fruit abscission[17]. Genes related to sucrose metabolism, e.g., sucrose-phosphate synthase (Pbr001745.1, Pbr008035.1), glycosyltransferase (Pbr014962.1 and Pbr033075.1), and glucose-1-phosphate adenyltransferase (Pbr036396.1) were up-regulated under GA3 treatment during the later calyx abscission process, while beta-fructofuranosidase was up-regulated under Flusilazole treatment at the onset of calyx abscission process. Levels of sugars like galactose, fructose, and glucose, as well as transcripts for their transporters, increase in senescing leaves analogously to what happens under environmental stresses, such as cold and dehydration[47]. A gene encoding a sugar transporter was up-regulated with Flusilazole treatment at the onset of the calyx abscission processes observed for Pbr008092.1. Sorbitol dehydrogenase (SDH) has been identified as a key enzyme in sorbitol metabolism, converting sorbitol into fructose[48, 49]. In this study, four genes encoding SDH (Pbr013916.1, Pbr032770.1, Pbr013912.1, and Pbr013913.1) were down-regulated with Flusilazole treatment in later calyx abscission processes, suggesting that sorbitol catabolism was largely inhibited, resulting in abscission.

Other differentially regulated genes

There were other genes that showed high-level differential expression related to calyx abscission (Table 3). E3 ubiquitin-protein ligase has been shown to play an important role in hormone regulation, photomorphogenesis, floral homeosis, senescence, and pathogen defense in plant[50]. Expression of Pbr009422.1, encoding E3 ubiquitin-protein ligase, was up-regulated under Flusilazole treatment during the calyx abscission processes. Expression of Pbr039793.1, encoding acyltransferase, was up-regulated in late calyx abscission processes. One gene (Pbr032411.1) encoding cysteine-rich receptor-like protein kinase 26-like was identified. Cysteine-rich receptor-like kinases is one of the largest groups of receptor-like kinases, which have been suggested to play important roles in the regulation of pathogen defense and programmed cell death[51]. It is notable that the gene Pbr013236.1 was up-regulated with Flusilazole treatment in both early and late calyx abscission processes; it showed 69% identity with IDA (Inflorescence deficient in abscission) protein of Arabidopsis, which is associated with regulation of floral organ abscission. IDA encodes a small protein with an N-terminal signal peptide. Analysis of ida mutant plants indicates that IDA regulates floral organ abscission through an ethylene insensitive pathway. Overexpression of IDA results in early abscission and production of arabinose and galactose in the floral AZs[15]. This suggests that the activity of IDA may be important to the onset and later stages of the calyx abscission process. However, further functional experiments are necessary to confirm this point. In addition, some differentially expressed genes (Pbr012802.1, Pbr012908.1 Pbr042050.1, and Pbr020727.1) without annotation were also found. We hypothesize that these genes are putative calyx abscission-related transcripts.

Confirmation of differentially expressed genes by qRT-PCR

In order to verify the genes that were actually differentially expressed during the calyx abscission processes, the expressive abundance of seven selected genes was analyzed by quantitative real-time PCR. The results showed that although the exact fold changes of six of the selected genes at several data points varied between digital transcript abundance measurements and qRT-PCR analysis, trends of gene expression change detected by the two different approaches were largely consistent. Only one gene (Pbr000187.1) did not show consistent expression between accurate quantification of expression and digital transcript abundance measurements (Additional file5). Pearson’s correlation coefficient (r) showed that both the digital transcript abundance measurements and qRT-PCR data were highly correlated, with the r-value ranging from 0.656 (Pbr036692.1) to 0.934 (Pbr001279.1), which was in agreement with previous report[52]. The qRT-PCR further demonstrated that genes related to photosystem reaction (Pbr041327.1), hormone-related transcripts (Pbr036692.1), carbohydrate metabolism (Pbr021608.1, Pbr008092.1, and Pbr001279.1) and other differentially regulated genes (Pbr027181.1) showed significant difference between treatments and participated in the process of calyx abscission or persistent processes.

Conclusions

The present results have demonstrated the usefulness of the digital transcript abundance measurements approach to identify differentially expressed genes between Flusilazole treatment and GA3 treatment. These differentially expressed genes may well be important for calyx abscission in fruit. In addition, a list of candidate target genes for functional studies involving calyx abscission process was generated. Among the isolated candidate genes, IDA appears to play a key role during calyx abscission processes. Further studies should be concentrated on functional characterization of these genes in the future. This study could lead to better understanding of the molecular mechanism of the phenotypic difference between calyx abscission and persistent fruits. Moreover, the findings of this study may facilitate the selection of new chemical agents and accelerate genetic strategies for the development of more effective pear calyx abscission for commercial purposes.

Methods

Plant materials and treatments

The plant materials used in this study were obtained from the Research Institute of Pomology, Chinese Academy of Agricultural Sciences (CAAS), Xingcheng, Liaoning province during the 2012 growing season. Five uniform fifty-year-old 'Kuerlexiangli’ trees were selected and divided into three blocks of six branches each. Two branches from each block were treated with: 6000 × Flusilazole + 300 × PBO (calyx abscising treatment) or GA3 50 mg.L-1 (calyx persisting treatment) sprayed at 0 d after full bloom, with plants with no treatment as the control. Under these treatments, calyx tube abscission symptoms (yellow loop in abscission zone) were observed within 10 d and only a few abscission symptoms were found with calyx persisting treatment. Therefore, we selected two early time points after treatment (6 and 10 d) in order to detect genes responsible for the early abscission events. An additional sample was collected 22 d after the treatment time point because at that stage, no further calyx abscission occurs. Samples of the two treatments at defined time points were collected for digital transcript abundance measurements (Table 1). A pear fruit with calyx tube at 22 days after full bloom is shown in Additional file6. A total of seven independent libraries were sequenced. At each time point, about 100 fruits were collected from each branch. The calyx abscission zone (AZ) tissues, containing a few layers of AZ cells on the proximal side of the separation line and adjacent cells, were manually dissected from the calyx tube samples, using a razor blade of 1 mm3. The AZ tissues were collected and frozen in liquid nitrogen and kept at -80°C until RNA isolation.

Field test of calyx abscission rate induced by different chemical agents

To determine the calyx abscission rate, all flowers on branches marked at the beginning of experiment, were counted and recorded at 22 d after full bloom in different treatments. Calyx abscission rate = number of fruits with calyx abscission/ number of all fruits tested.

RNA isolation and Solexa/Illumina sequencing

Solexa/Illumina sequencing was carried out by CapitalBio Corporation, Beijing, China. The total RNA was extracted from the samples using Plant RNA Isolation Kit (AutoLab), followed by RNA purification with RNeasy MiniElute Cleanup Kit (Qiagen), according to the manufacturer’s instruction. Total RNA content, purity and degradation were assessed by Nanodrop2000 spectrophotometer (Thermo scientific, USA) and quality of RNA was confirmed by agarose gel electrophoresis before proceeding. For mRNA library construction and deep sequencing, RNA samples were prepared using the TruSeq RNA Sample Preparation Kit according to the manufacturer’s protocol. Briefly, the poly-A-containing mRNA molecules were purified from 3 μg of total RNA using poly-T oligo-attached magnetic beads with two rounds of purification. For the second round elution of the poly-A RNA, the RNA was fragmented using divalent cations under 95°C. For Solexa/Illumina sequencing, cDNA synthesis was carried out with the broken RNA fragments and these RNA fragments reversely transcribed into first strand cDNA using random hexamers. Second-strand cDNA synthesis using DNA Polymerase I and RNase H. The cDNA fragments were put through an end repair process to convert the overhangs into blunt ends using an End Repair (ERP) mix. The 3′ to 5′ exonuclease activity of this mix removes the 3′ overhangs and the polymerase activity fills in the 5′ overhangs. A single 'A’ nucleotide was then added to the 3′ ends of the blunt fragments to prevent them from ligating to one another during the adapter ligation reaction. A corresponding single 'T’ nucleotide on the 3′ end of the adapter provides a complementary overhang for ligating the adapter to the fragment. This strategy ensures a low rate of chimera (concatenated template) formation. The multiple indexing adapters were ligated to the ends of the double-stranded cDNA, making them for hybridization onto the Illumina Sequencing Chip (flow cell). PCR was used to selectively enrich those DNA fragments that have adapter molecules on both ends and to amplify the amount of DNA in the library, and was minimized to twelve cycles to avoid skewing the representation of the library. A gel purification procedure was carried out to select the fragments sized from 300 to 400 bp to produce the library for cluster generation and sequencing. The libraries were checked for quality by Agilent 2100 bioanalyzer and quantified by Qubit and qPCR. Cluster formation and sequencing on the GAIIx platform were performed following the manufacturer’s standard cBot and sequencing protocols. For the multiplexing sequencing, 35 cycles of single read were used to sequence the RNA, followed by 7 cycles of index identification.

Data analysis

Primary data analysis and base calling were performed using the Illumina instrument software. The following sequencing data were excluded from the analysis: low quality sequences such as the 3’ adaptor sequence; tags which were too long or too short; tags with unknown sequence; single copy tags. The remaining high quality sequences (clean data) were mapped to the pear gene set[22] using the software tool Bowtie[53]. A Perl script was written to process the mapping results and generate the gene expression profile. Similar to credibility interval approaches reported for the analysis of SAGE data[54], we employed IDEG6[55] to identify mRNAs showing statistically significant differences based on their relative abundance (as reflected by total count of individual sequence reads) between the two libraries. The general chi-square test was performed, as it has been proven to be one of the most efficient tests (http://bioservices.capitalbio.com/xzzq/rj/3885.shtml). Finally, genes with a P value < = 0.01 and Fold Change > =2 were marked significantly different between the two libraries. InterPro domains[56] were annotated by InterProScan[57] Release 36.0 and functional assignments were mapped onto Gene Ontology (GO)[58]. Furthermore, the GO classification and draw GO tree using WEGO[59]. Genes were mapped to terms in the Kyoto Encyclopedia of Genes and Genomes database (KEGG, release[60]) using BLASTX[61] at E values < = 1e-10. A Perl script was used to retrieve KO (KEGG Ontology) information from the blast search result and then establish pathway associations between unigene and database.

Real-time PCR

To validate the expression patterns revealed by digital transcript abundance measurements results, seven genes identified through digital transcript abundance measurements were analyzed using quantitative real-time PCR. The RNA samples for digital transcript abundance measurements were also used for quantitative real-time PCR. Primer sequences for the real-timer PCR assay were designed using eprimer 3 (from the program EMBOSS Explorer) and listed in Additional file7. Total RNA was treated with DNase I to remove genomic DNA contamination. Approximately 1 μg of total RNA was used as a template for reverse transcription using ReverTra Ace-αFirst Strand cDNA Synthesis Kit (TOYOBO,TOYOBO Biotech Co. Ltd, Japan ) according to the manufacturers’ instructions. The qRT-PCRs were performed using SYBR® Green Premix kit (TaKaRa Biotechnology. Dalian, China). The composition of PCR mix was as follows: 10 μl 2 × SYBR Premix ExTaq™, 0.4 μl each primer and 1μl of cDNA template in a final volume of 20 μl. All reactions were run as duplicates in 96-well plates. The qRT-PCR was performed on the Lightcycle-480 (Roche). Each cDNA was analyzed in triplicate, after which the average threshold cycle (Ct) was calculated per sample. The relative expression levels were calculated with the 2-ΔΔCt method[62]. Tubulin (AB239681) was used as the internal control. The primer pairs were: forward: 5′ TGGGCTTTGCTCCTCTTAC 3′; reverse: 5′ CCTTCGTGCTCATCTTACC 3′. The protocol of real-time PCR was as follows: initiation with a 10 min denaturation at 95°C followed by 40 cycles of amplification with 15 s of denaturation at 95°C, 15 s of annealing at 60°C, 20 s of extension at 72°C and reading the plate for fluorescence data collection at 60°C. A melting curve was performed from 60 to 95°C to check the specificity to the amplified product.

References

  1. 1.

    Yusan K, Qiman Y, Li J, Qin WM: Comparison of fruit quality of normal and roughbark pear in persistent calyx and dropping calyx of Korla fragrant pear. Xingjiang Agric Sci. 2012, 49: 1028-1034.

  2. 2.

    Gang MH, Qiman Y, Qin WM, Li J, Li J: Effect of different plant growth regulated on Korla fragrant pear quality and morphological characteristics. Xingjiang Agric Sci. 2010, 47: 736-740.

  3. 3.

    Gang MH, Qiman Y, Qin WM, Li YH, Li H, Li J: Effect of different plant growth regualtion on set fruit rate, calyx dropping fruit rate, fruit shape and chlorophyll content of Korla fragrant pear. J Xinjiang Agric Univ. 2009, 32: 26-30.

  4. 4.

    Liu N, Tao ST, Zhang SL, Cao YF, Wu HQ, Wang JZ, Wu T: Effect of different pollinizer varieties on calyx retention and quality for 'Kuerlexiangli’ fruit. J Nanjing Agric Univ. 2011, 34: 43-47.

  5. 5.

    Jiang YC, Cao YF, Zhang SL, Wu T, Tian LM, Dong XG: Effects of excession of fruit calyx on fruit quality, sclereid and the activities of related enzymes in Pyrus bretschneideri cv. Laiyang Cili. J Fruit Sci. 2010, 27: 877-881.

  6. 6.

    Ma HC, Wang YL, Wen X, Qi M, Li J: Effects of different regent treatments on the microstructure of calyx tube with or without calyx of Korla fragrant pear. J Fruit Sci. 2011, 28: 518-520.

  7. 7.

    Jia B, Zhu LW, Zhang SL: Preliminary study on the action of the mineral element Fe in fruit calyx development of 'DangshanSuli’ pear. Acta Hortic Sin. 2012, 39: 159-167.

  8. 8.

    Niu JX, He ZS: Dynamic changes of phytohormone content in pear calyx and young fruit during calyx growth and development. J Fruit Sci. 2009, 26: 431-434.

  9. 9.

    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.

  10. 10.

    Chen S, Jiang J, Li HY, Liu GF: The salt-responsive transcriptome of Populus simonii × Populus nigra via DGE. Gene. 2012, 504: 203-212. 10.1016/j.gene.2012.05.023.

  11. 11.

    Urasaki N, Tarora K, Shudo A, Ueno H, Tamaki M, Miyagi N, Adaniya S, Matsumura H: Digital transcriptome analysis of putative sex-determination genes in papaya (Carica papaya). PLoS One. 2012, 7: e40904-10.1371/journal.pone.0040904.

  12. 12.

    Wang C, Chen HF, Hao QN, Sha AH, Shan ZH, Chen LM, Zhou R, Zhi HJ, Zhou XA: Transcript profile of the response of two soybean genotypes to potassium deficiency. PLoS One. 2012, 7: e39856-10.1371/journal.pone.0039856.

  13. 13.

    Qi XH, Xu XW, Lin XJ, Zhang WJ, Chen XH: Identification of differentially expressed genes in cucumber (Cucumis sativus L.) root under waterlogging stress by digital gene expression profile. Genomics. 2012, 99: 160-168. 10.1016/j.ygeno.2011.12.008.

  14. 14.

    Hao QN, Zhou XA, Sha AH, Wang C, Zhou R, Chen SL: Identification of genes associated with nitrogen-use efficiency by genome-wide transcriptional analysis of two soybean genotypes. BMC Genomics. 2011, 12: 525-10.1186/1471-2164-12-525.

  15. 15.

    Cho SK, Larue CT, Chevalier D, Wang HC, Jinn TL, Zhang SQ, Walker JC: Regulation of floral organ abscission in Arabidopsis thaliana. Proc Natl Acad Sci USA. 2008, 105: 15629-15634. 10.1073/pnas.0805539105.

  16. 16.

    Mao L, Begum D, Chuang HW, Budiman MA, Szymkowiak EJ, Irish EE, Wing RA: JOINTLESS is a MADS-box gene controlling tomato flower abscission zone development. Nature. 2000, 406: 910-913. 10.1038/35022611.

  17. 17.

    Botton A, Eccher G, Forcato C, Ferrarini A, Begheldo M, Zermiani M, Moscatello S, Battistelli A, Velasco R, Ruperti B, Ramina A: Signaling pathways mediating the induction of apple fruitlet abscission. Plant Physiol. 2011, 155: 185-208. 10.1104/pp.110.165779.

  18. 18.

    Meir S, Philosoph-Hadas S, Sundaresan S, Selvaraj KS, Burd S, Ophir R, Kochanek B, Reid MS, Jiang CZ, Lers A: Microarray analysis of the abscission-related transcriptome in the tomato flower abscission zone in response to auxin depletion. Plant Physiol. 2010, 154: 1929-1956. 10.1104/pp.110.160697.

  19. 19.

    Agustí J, Merelo P, Cercós M, Tadeo FR, Talón M: Ethylene-induced differential gene expression during abscission of citrus leaves. J Exp Bot. 2008, 59: 2717-2733. 10.1093/jxb/ern138.

  20. 20.

    Dal Cin V, Barbaro E, Danesin M, Murayama H, Velasco R, Ramina A: Fruitlet abscission: a cDNA-AFLP approach to study genes differentially expressed during shedding of immature fruits reveals the involvement of a putative auxin hydrogen symporter in apple (malus domestica L. Borkh). Gene. 2009, 442: 26-36. 10.1016/j.gene.2009.04.009.

  21. 21.

    Zhu H, Dardick CD, Beers EP, Callanhan AM, Xia R, Yuan RC: Transcriptomics of shading-induced and NAA-induced abscission in apple (Malus domestica) reveals a shared pathway involving reduced photosynthesis, alterations in carbohydrate transport and signaling and hormone crosstalk. BMC Plant Biol. 2011, 11: 138-10.1186/1471-2229-11-138.

  22. 22.

    Wu J, Wang ZW, Shi ZB, Zhang S, Ming R, Zhu SL, Khan MA, Tao ST, Korban SS, Wang H, Chen NJ, Nishio T, Xu X, Cong L, Qi KJ, Huang XS, Wang YT, Zhao X, Wu JY, Deng C, Gou CY, Zhou WL, Yin H, Qin GH, Sha YH, Tao Y, Chen H, Yang YN, Song Y, Zhan DL: The genome of pear (Pyrus bretschneideri Rehd.). Genome Res. 2012, 23: 396-408.

  23. 23.

    Taylor JE, Whitelaw CA: Signals in abscission. New Phytol. 2001, 151: 323-340. 10.1046/j.0028-646x.2001.00194.x.

  24. 24.

    Mesejo C, Rosito S, Reig C, Martínez-Fuentes A, Agustí M: Synthetic auxin 3, 5, 6-TPA provokes citrus clementina (Hort. ex Tan) fruitlet abscission by reducing photosynthate availability. J Plant Growth Regul. 2012, 31: 186-194. 10.1007/s00344-011-9230-z.

  25. 25.

    Dal Cin V, Velasco R, Ramina A: Dominance induction of fruitlet shedding in Malus x domestica (L. Borkh): molecular changes associated with polar auxin transport. BMC Plant Biol. 2009, 9: 139-10.1186/1471-2229-9-139.

  26. 26.

    Alexander L, Grierson D: Ethylene biosynthesis and action in tomato: a model for climacteric fruit ripening. J Exp Bot. 2002, 53: 2039-2055. 10.1093/jxb/erf072.

  27. 27.

    Wang KLC, Li H, Ecker JR: Ethylene biosynthesis and signaling networks. Plant Cell. 2002, 14 (suppl 1): S131-S151.

  28. 28.

    Parra-Lobato MC, Gomez-Jimenez MC: Polyamine-induced modulation of genes involved in ethylene biosynthesis and signalling pathways and nitric oxide production during olive mature fruit abscission. J Exp Bot. 2011, 62: 4447-4465. 10.1093/jxb/err124.

  29. 29.

    Wang Y, Mao Z, Craker LE: Auxin transport and abscission in coleus. HortSci. 1991, 26: 686-686.

  30. 30.

    Tucker ML, Burke A, Murphy CA, Thai VK, Ehrenfried ML: Gene expression profiles for cell wall-modifying proteins associated with soybean cyst nematode infection, petiole abscission, root tips, flowers, apical buds, and leaves. J Exp Bot. 2007, 58: 3395-3406. 10.1093/jxb/erm188.

  31. 31.

    Hong SB, Sexton R, Tucker ML: Analysis of gene promoters for two tomato polygalacturonases expressed in abscission zones and the stigma. Plant Physiol. 2000, 123: 869-881. 10.1104/pp.123.3.869.

  32. 32.

    González-Carranza ZH, Whitelaw CA, Swarup R, Roberts JA: Temporal and spatial expression of a polygalacturonase during leaf and flower abscission in oilseed rape and Arabidopsis. Plant Physiol. 2002, 128: 534-543. 10.1104/pp.010610.

  33. 33.

    Mishra A, Khare S, Trivedi PK, Nath P: Ethylene induced cotton leaf abscission is associated with higher expression of cellulase (GhCel1) and increased activities of ethylene biosynthesis enzymes in abscission zone. Plant Physiol Biochem. 2008, 46: 54-63. 10.1016/j.plaphy.2007.09.002.

  34. 34.

    Kazokas WC, Burns JK: Cellulase activity and gene expression in citrus fruit abscission zones during and after ethylene treatment. J Amer Soc Hort Sci. 1998, 123: 781-786.

  35. 35.

    Sexton R, Palmer JM, Whyte NA, Littlejohns S: Cellulase, fruit softening and abscission in red raspberry Rubus idaeus L. cv glen clova. Ann Bot (london). 1997, 80: 371-376. 10.1006/anbo.1997.0465.

  36. 36.

    Sane AP, Tripathi SK, Nath P: Petal abscission in rose (Rosa bourboniana var Gruss an Teplitz) is associated with the enhanced expression of an alpha expansin gene, RbEXPA 1. Plant Sci. 2007, 172: 481-487. 10.1016/j.plantsci.2006.10.005.

  37. 37.

    Singh AP, Tripathi SK, Nath P, Sane AP: Petal abscission in rose is associated with the differential expression of two ethylene-responsive xyloglucan endotransglucosylase/hydrolase genes, RbXTH1 and RbXTH2. J Exp Bot. 2011, 62: 5091-5103. 10.1093/jxb/err209.

  38. 38.

    Sun LX, Nocker SV: Analysis of promoter activity of members of the PECTATE LYASE-LIKE (PLL) gene family in cell separation in Arabidopsis. BMC Plant Biol. 2010, 10: 152-10.1186/1471-2229-10-152.

  39. 39.

    Singh AP, Pandey SP, Rajluxmi Pandey S, Nath P, Sane AP: Transcriptional activation of a pectate lyase gene, RbPel1, during petal abscission in rose. Postharvest Biol Tec. 2011, 60: 143-148. 10.1016/j.postharvbio.2010.12.014.

  40. 40.

    Okamuro JK, Caster B, Villarroel R, Van Montagu M, Jofuku KD: The AP2 domain of APETALA2 defines a large new family of DNA binding proteins in Arabidopsis. Proc Natl Acad Sci USA. 1997, 94: 7076-7081. 10.1073/pnas.94.13.7076.

  41. 41.

    Brenner WG, Schmülling T: Transcript profiling of cytokinin action in Arabidopsis roots and shoots discovers largely similar but also organ-specific responses. BMC Plant Biol. 2012, 12: 112-10.1186/1471-2229-12-112.

  42. 42.

    Mehrtens F, Kranz H, Bednarek P, Weisshaar B: The Arabidopsis transcription factor MYB12 is a flavonol-specific regulator of phenylpropanoid biosynthesis. Plant Physiol. 2005, 138: 1083-1096. 10.1104/pp.104.058032.

  43. 43.

    Lewis DR, Ramirez MV, Miller ND, Vallabhaneni P, Ray WK, Helm RF, Winkel BSJ, Muday GK: Auxin and ethylene induce flavonol accumulation through distinct transcriptional networks. Plant Physiol. 2011, 156: 144-164. 10.1104/pp.111.172502.

  44. 44.

    Robatzek S, Somssich IE: A new member of the Arabidopsis WRKY transcription factor family, At WRKY6, is associated with both senescence-and defence-related processes. Plant J. 2001, 28: 123-133. 10.1046/j.1365-313X.2001.01131.x.

  45. 45.

    Iglesias DJ, Tadeo FR, Primo-Millo E, Talon M: Fruit set dependence on carbohydrate availability in citrus trees. Tree Physiol. 2003, 23: 199-204. 10.1093/treephys/23.3.199.

  46. 46.

    Nzima MDS, Martin GC, Nishijima C: Effect of fall defoliation and spring shading on shoot carbohydrate and growth parameters among individual branches of alternate bearing 'Kerman’ pistachio trees. J Amer Soc Hortic Sci. 1999, 124: 52-60.

  47. 47.

    Kiyosue T, Abe H, Yamaguchi-Shinozaki K, Shinozaki K: ERD6, a cDNA clone for an early dehydration-induced gene of Arabidopsis, encodes a putative sugar transporter. BBA-Biomembranes. 1998, 1370: 187-191. 10.1016/S0005-2736(98)00007-8.

  48. 48.

    Berüter J, Feusi MES, Rüedi P: Sorbitol and sucrose partitioning in the growing apple fruit. J Plant Physiol. 1997, 151: 269-276. 10.1016/S0176-1617(97)80252-0.

  49. 49.

    Yamaguchi H, Kanayama Y, Soejima J, Yamaki S: Changes in the amounts of the NAD-dependent sorbitol dehydrogenase and its involvement in the development of apple fruit. J Amer Soc Hort Sci. 1996, 121: 848-852.

  50. 50.

    Zhang YY, Yang CW, Li Y, Zheng NY, Chen H, Zhao QZ, Gao T, Guo HS, Xie Q: SDIR1 is a ring finger E3 ligase that positively regulates stress-responsive abscisic acid signaling in Arabidopsis. Plant Cell. 2007, 19: 1912-1929. 10.1105/tpc.106.048488.

  51. 51.

    Wrzaczek M, Brosché M, Salojärvi J, Kangasjärvi S, Idänheimo N, Mersmann S, Robatzek S, Karpiński S, Karpińska B, Kangasjärvi J: Transcriptional regulation of the CRK/DUF26 group of receptor-like protein kinases by ozone and plant hormones in Arabidopsis. BMC Plant Biol. 2010, 10: 95-10.1186/1471-2229-10-95.

  52. 52.

    Pan FG, Zhao YY, Zhu S, Sun CJ, Lei LC, Feng X, Han WY: Different transcriptional profiles of RAW264. 7 infected with Mycobacterium tuberculosis H37Rv and BCG identified via deep sequencing. PLoS One. 2012, 7: e51988-10.1371/journal.pone.0051988.

  53. 53.

    Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.

  54. 54.

    Vêncio RZN, Brentani H, Pereira CAB: Using credibility intervals instead of hypothesis tests in SAGE analysis. Bioinformatics. 2003, 19: 2461-2464. 10.1093/bioinformatics/btg357.

  55. 55.

    Romualdi C, Bortoluzzi S, D’Alessi F, Danieli GA: IDEG6: a web tool for detection of differentially expressed genes in multiple tag sampling experiments. Physiol Genomics. 2003, 12: 159-162.

  56. 56.

    Mulder NJ, Apweiler R, Attwood TK, Bairoch A, Barrell D, Bateman A, Binns D, Biswas M, Bradley P, Bork P, Bucher P, Copley RR, Courcelle E, Das U, Durbin R, Falquet L, Fleischmann W, Griffiths-Jonses S, Haft D, Harte N, Hulo N, Kahn D, Kanapin A, Krestyaninova M, Lopez R, Letunic I, Lonsdale D, Silventoinen V, Orchard SE, Pagni M: The InterPro Database, 2003 brings increased coverage and new features. Nucleic Acids Res. 2003, 31: 315-318. 10.1093/nar/gkg046.

  57. 57.

    Zdobnov EM, Apwelier R: InterProScan an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001, 17: 847-848. 10.1093/bioinformatics/17.9.847.

  58. 58.

    Harris M, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C, Richter J, Rubin GM, Blake JA, Bult C, Dolan M, Drabkin H, Eppig JT, Hill DP, Ni L, Ringwald M, Balakrishnan R, Cherry JM, Christie KR, Costanzo MC, Dwight SS, Engel S, Fisk DG, Hirschman JE, Hong EL, Nash RS: The gene ontology (GO) database and informatics resource. Nucleic Acids Res. 2004, 32: D258-10.1093/nar/gkh036. Database issue

  59. 59.

    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 (suppl 2): W293-W297.

  60. 60.

    Kanehisa M, Goto S, Hattori M, Aoki-Kinoshita KF, Itoh M, Kawashima S, Katayama T, Araki M, Hirakawa M: From genomics to chemical genomics: new developments in KEGG. Nucleic Acids Res. 2006, 34: D354-D357. 10.1093/nar/gkj102. Database issue

  61. 61.

    Altschul SF, Madden TL, Schäffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.

  62. 62.

    Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2- ΔΔCT method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.

Download references

Acknowledgements

This work was financially supported by National Natural Science Foundation of China (31230063), the Earmarked Fund for China Agriculture Research System (CARS-29), and Independent Innovation of Agricultural Sciences in Jiangsu Province (CX (11) 1013).

Author information

Correspondence to Jun Wu or Shaoling Zhang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

XXQ contributed to the experimental design and management, data analysis, and manuscript preparation. LFW contributed to quantitative qRT-PCR. LTL contributed to data analysis. YFC, LMT and XGD prepared the plant materials and tissue collection. JW and SLZ designed and managed the experiments, organize and reviewed the manuscript. All authors have read and approved the final manuscript.

Electronic supplementary material

Additional file 1:Pear fruit without calyx (A) and pear fruit with calyx (white arrow) (B) during the fruit mature period.(DOC 758 KB)

Additional file 2:Pathway enrichment analysis of differentially expressed genes. A BLAST search against the KEGG database indicated that expressed genes from C1 to C7 were involved in 251 pathways. (DOC 220 KB)

Additional file 3:List of significantly altered pathways with enriched differentially expressed genes after Flusilazole treatment and GA 3 treatment. The table provides the list of significantly alter biochemical pathways during calyx abscission process: Pathway ID, pathway name, the total number of genes with pathway annotation, the number of differentially expressed genes with pathway annotation of each comparison (Fifteen pairs comparison in total in this table, e.g. C1 vs C7 indicates the difference expressed genes in the comparison between C1and C7). (DOC 44 KB)

Additional file 4:Categorization of significant genes encoding enzymes and proteins with a variety of biological functions. In this table, five functional categories of genes showed differential expression patterns after Flusilazole treatment and GA3 treatment. The details of gene expression analysis of 15 pairs of comparisons: Gene ID, gene description and TPM (transcripts copies per million tags) of genes. (DOC 2 MB)

Additional file 5:Real-time quantitative PCR analysis of selected differential genes detected via digital transcript abundance measurements. Representative genes selected for the analysis were those involved in photosynthesis, plant hormone signal transduction, carbohydrate metabolism, cell wall degradation, and other processes. (DOC 88 KB)

Additional file 6:A pear fruit with calyx tube (white arrow) at 22 d after full bloom.(DOC 764 KB)

Additional file 7:Primer details for genes selected for quantitative real-time PCR analysis from results of digital transcript abundance measurements. This is the primer list of seven genes selected for quantitative real-time PCR assay to confirm the reliability of digital transcript abundance measurements. Gene ID, forward and reverse primers are shown. (DOC 30 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Qi, X., Wu, J., Wang, L. et al. Identifying the candidate genes involved in the calyx abscission process of 'Kuerlexiangli’ (Pyrus sinkiangensis Yu) by digital transcript abundance measurements. BMC Genomics 14, 727 (2013) doi:10.1186/1471-2164-14-727

Download citation

Keywords

  • Pear
  • Calyx abscission
  • Digital transcript abundance measurements
  • Differentially expressed genes