Identifying the candidate genes involved in the calyx abscission process of 'Kuerlexiangli’ (Pyrus sinkiangensis Yu) by digital transcript abundance measurements
© Qi et al.; licensee BioMed Central Ltd. 2013
Received: 8 April 2013
Accepted: 18 October 2013
Published: 23 October 2013
'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.
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.
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.
KeywordsPear Calyx abscission Digital transcript abundance measurements Differentially expressed genes
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. 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. Moreover, manually excising the calyx was helpful for enhancing the comprehensive quality of pear fruit.
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. 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. 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. 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. 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. 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. 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. Differentially expressed genes in cucumber (Cucumis sativus L.) root under waterlogging stress have also been identified by digital transcript abundance measurements. 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.
Most current molecular knowledge on the abscission process has been obtained from the model plant Arabidopsis thaliana, as well as tomato (Solanum lycopersicum). The apple cluster during immature fruit drop also represents an ideal system to study the shedding of actively growing organs. 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, differential gene expression by the cDNA microarray technique during abscission of citrus leaves under ethylene treatment and fruitlet abscission in apple. The differentially expressed genes have been identified during shedding of immature apple fruits have been identified with a cDNA-AFLP approach, and the role of NAA and shading in apple fruit abscission have been analyzed by the transcriptome method. 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’, 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
Digital transcript abundance measurements libraries sequencing
Description of calyx tube abscission zone samples used in this study
Calyx tubes at 6 days after 6000 × Flusilazole + 300 × PBO treatment
Calyx tubes at 6 days after 50 mg/L GA3 treatment
Yellow loop in abscission zone from calyx tube at 10 days after 6000 × Flusilazole + 300 × PBO treatment
Calyx tubes at 10 days after treated with 50 mg/L GA3
Abscission zone from calyx tubes at 22 days after 6000 × Flusilazole + 300 × PBO treatment with calyx removing
Calyx tubes were still present at 22 days after 6000 × Flusilazole + 300 × PBO treatment
Calyx tubes were still present with control at the same time as in C6
Statistics of digital transcript abundance measurements library sequencing and tag mapping
Total length (bp)
Single length (bp)
Data mapping to gene
Analysis of differential gene expression
Functional annotation of differentially expressed genes
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. 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
Selected genes with altered expression during calyx abscission process
TPM (transcripts per million clean tag)
Photosystem I subunit V
Photosystem I PsaA/PsaB
Photosystem I PsaA/PsaB
F-type H + -transporting ATPase subunit gamma
Photosystem II 13 kDa protein
ACC oxidase 7
Ethylene-responsive transcription factor
Ethylene-responsive transcription factor
Ethylene-responsive transcription factor 1B
PtrAUX6; auxin influx carrier component
PIN6; auxin efflux carrier component
Cell wall modification
AP2 domain-containing transcription factor
MYB 12 transcription factor
WRKY transcription factor
WRKY transcription factor
SCL domain class transcription factor
Carbohydrate biosynthesis and metabolism
Sucrose phosphate syntase
Sucrose phosphate syntase
Other differentially regulated genes
Probable E3 ubiquitin-protein ligase
Cysteine-rich receptor-like protein kinase
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. 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 and apple. 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. 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. Polygalacturonases have been studied in tomato abscission zones, in oilseed rape and Arabidopsis leaf and flower abscission zones, while celluloses have been studied in cotton, citrus fruit and red raspberry abscission zones. Reports have indicated that an increase in expansin, xyloglucan transglucosylase hydrolases and pectate lyases 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:184.108.40.206) 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. 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. 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. 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. 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. 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. This phenomenon has been also described for pistachio. 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: 220.127.116.11), beta-fructofuranosidase (EC: 18.104.22.168), alcohol dehydrogenase (EC: 22.214.171.124), and acetyl-CoA carboxylase (EC: 126.96.36.199) under Flusilazole treatment. A different gene set showed lower expression in calyx abscission. This set comprised genes coding for 2-isopropylmalate synthase (EC: 188.8.131.52), pyruvate decarboxylase (EC: 184.108.40.206), sorbitol dehydrogenase (EC: 220.127.116.11) and ribulose-1,5-bisphosphate carboxylase (EC: 18.104.22.168). 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 (EC22.214.171.124) was highly expressed under Flusilazole treatment in the abscission process (Table 3). In apple, UDP-glucosyltransferases were also expressed during fruit abscission. 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. 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. 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. 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. 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. 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.
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.
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.
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 using the software tool Bowtie. 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, we employed IDEG6 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 were annotated by InterProScan Release 36.0 and functional assignments were mapped onto Gene Ontology (GO). Furthermore, the GO classification and draw GO tree using WEGO. Genes were mapped to terms in the Kyoto Encyclopedia of Genes and Genomes database (KEGG, release) using BLASTX 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.
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. 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.
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).
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google 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
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Taylor JE, Whitelaw CA: Signals in abscission. New Phytol. 2001, 151: 323-340. 10.1046/j.0028-646x.2001.00194.x.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Wang KLC, Li H, Ecker JR: Ethylene biosynthesis and signaling networks. Plant Cell. 2002, 14 (suppl 1): S131-S151.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Y, Mao Z, Craker LE: Auxin transport and abscission in coleus. HortSci. 1991, 26: 686-686.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.Google Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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 issueView 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 (suppl 2): W293-W297.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 issuePubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.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.