De novo assembly and characterization of fruit transcriptome in Litchi chinensis Sonn and analysis of differentially regulated genes in fruit in response to shading
© Li et al.; licensee BioMed Central Ltd. 2013
Received: 7 January 2013
Accepted: 9 August 2013
Published: 14 August 2013
Litchi (Litchi chinensis Sonn.) is one of the most important fruit trees cultivated in tropical and subtropical areas. However, a lack of transcriptomic and genomic information hinders our understanding of the molecular mechanisms underlying fruit set and fruit development in litchi. Shading during early fruit development decreases fruit growth and induces fruit abscission. Here, high-throughput RNA sequencing (RNA-Seq) was employed for the de novo assembly and characterization of the fruit transcriptome in litchi, and differentially regulated genes, which are responsive to shading, were also investigated using digital transcript abundance(DTA)profiling.
More than 53 million paired-end reads were generated and assembled into 57,050 unigenes with an average length of 601 bp. These unigenes were annotated by querying against various public databases, with 34,029 unigenes found to be homologous to genes in the NCBI GenBank database and 22,945 unigenes annotated based on known proteins in the Swiss-Prot database. In further orthologous analyses, 5,885 unigenes were assigned with one or more Gene Ontology terms, 10,234 hits were aligned to the 24 Clusters of Orthologous Groups classifications and 15,330 unigenes were classified into 266 Kyoto Encyclopedia of Genes and Genomes pathways. Based on the newly assembled transcriptome, the DTA profiling approach was applied to investigate the differentially expressed genes related to shading stress. A total of 3.6 million and 3.5 million high-quality tags were generated from shaded and non-shaded libraries, respectively. As many as 1,039 unigenes were shown to be significantly differentially regulated. Eleven of the 14 differentially regulated unigenes, which were randomly selected for more detailed expression comparison during the course of shading treatment, were identified as being likely to be involved in the process of fruitlet abscission in litchi.
The assembled transcriptome of litchi fruit provides a global description of expressed genes in litchi fruit development, and could serve as an ideal repository for future functional characterization of specific genes. The DTA analysis revealed that more than 1000 differentially regulated unigenes respond to the shading signal, some of which might be involved in the fruitlet abscission process in litchi, shedding new light on the molecular mechanisms underlying organ abscission.
KeywordsLitchi chinensis Transcriptome Fruit RNA-Seq DTA Shade Abscission
Litchi (Litchi chinensis Sonn.), one of the most important tropical and subtropical plants in the Sapindaceae family was originated in China, where it has been cultivated for more than 2,300 years . The Sapindaceae is a relatively large family, including at least 125 genera and 1,000 species, which are widely distributed in the tropical and subtropical areas. Being the most widely cultivated fruit crop in this family, litchi becomes a significant contribution to the lives and economies of millions of people in Southeast Asia. In 2010, litchi production in China was 1,780,000 tons from 553,000 ha , providing major employment for the local people. The desirable characteristics of litchi fruit, such as a bright colour, exotic aroma, excellent flavour and rich nutrients, make it very attractive and popular in the international markets [3, 4]. Despite the economic and commercial importance of litchi, there have been few genomic studies on this crop. This is confirmed by the very few number of litchi gene sequences available in public databases, with only 354 nucleotide sequences available in the NCBI GenBank as of 6 Augest 2013 (http://www.ncbi.nlm.nih.gov/GenBank/). The lack of sequence data for the species of the Sapindaceae family has greatly hindered their research at the molecular level.
Next generation sequencing technology (NGS), such as high-throughput paired-end RNA sequencing (RNA-Seq) and digital transcript abundance(DTA)tag profiling, has greatly facilitated investigation of the functional complexity of transcriptomes for non-model organisms without a reference genome [5–7]. De novo assembly of a transcriptome from RNA-Seq produces a genome-scale transcription map that contains both the transcriptional structure and expression level for each gene.
Even though RNA-Seq and DTA have not yet been performed with the members of the Sapindaceae family, they have been widely applied to many other fruit crops for transcriptome analysis, gene discovery and molecular marker development. For examples: more than 59 million reads were generated and 6,695 unigenes were expressed in a stage-specific manner in three different developmental stages of berries in Vitis vinifera; 9,839 transcript assemblies were obtained and about 115 simple sequence repeat markers were identified in pomegranate (Punica granatum) ; 41,239 unigenes with a mean length of 531 bp were de novo assembled and more than 3,600 unigenes were differentially expressed in three ripening stages in Chinese bayberry (Myrica rubra) . In Siraitia grosvenorii, RNA-Seq was combined with DTA to understand mogroside biosynthesis. About 49 million high-quality reads were assembled into 43,891 unigenes, in which 85 cytochrome P450 and 90 UDP-glucosyltransferase unigenes were identified, and seven and five of these unigenes, respectively, were found to be involved in mogroside biosynthesis .
It is well known that litchi trees generally produce many more female flowers than necessary . Most of the inevitable massive abscission of flowers and fruitlet occurs after pollination, leading to a low fruit set . To date, rather than investigating at the molecular level, most biological studies on abscission in litchi have been carried out on the physiological level related to nutrients and hormones [14, 15].
During early fruit development, active sinks such as the growing shoots and fruit compete for limited carbohydrate and nutrient resources . Artificial shading over the whole canopy is a practical research approach used to induce fruit abscission and investigate its various underlying physiological and molecular mechanisms. Shading during this period not only rapidly decreases photosynthesis, reducing the availability of assimilates but further aggravates the competition among these sinks [17–19]. A decrease in the relative growth rates of fruit is apparent within 2 days after shading [18, 19]. Fruit abscission begins at 5 to 10 days after shading and peaks at 15 days [20–22]. Therefore, reduction in fruit growth is an earlier response to shading, while fruit abscission is the final effect of shading on fruit. Severe fruitlet abscission usually leads to a very low yield and causes significant economic losses for farmers.
During the last few years, several methods have been used to study the transcriptional regulation of fruit abscission, including suppression subtractive hybridization and cDNA-amplified fragment length polymorphism, as well as microarrays. With suppression subtractive hybridization technology, a total of 347 expressed sequence tags were obtained and only 112 unigenes were found in the cDNA library of fruitlet, in apples (Malus × domestica), after shading . cDNA-amplified fragment length polymorphism identified 227 differentially expressed clones isolated from the seeds, cortices and peduncles of abscising fruitlet and non-abscising fruitlet populations. Because of the limited throughput and high false-positive rate of the two methods mentioned above, microarray was applied to several species to obtain large-scale transcriptional regulation information during organ abscission in apples [22, 25], citrus (Citrus clementina)  and tomato (Solanum lycopersicum) . Recently, Gil-Amado and Gomez-Jimenez  compared the olive fruit abscission zone transcriptomes at two different stages (pre-abscission versus abscission) using the RNA-Seq technique; 148 Mb of sequences (443,811 good-quality sequence reads) were obtained and 4,728 differentially expressed genes were identified from these two samples. Of these 4,728 differentially expressed genes, 2,314 were up-regulated and 2,414 were down-regulated at the abscission stage in the fruit abscission zone .
RNA-Seq is a fast and comprehensive approach for direct sequencing at an extraordinary depth, while DTA is based on sequencing serial analysis of gene- expression libraries and generates a digital output proportional to the number of transcripts per mRNA . Over the last few years, studies in marine fish (Lateolabrax japonicus) , whitefly (Bemisia tabaci)  and Siraitia grosvenorii have demonstrated that the combination of these two technologies is very suitable for studying transcriptome profiles and provides a good understanding of the complexity of gene expression, regulation and networks. However, none of the techniques mentioned above have been applied in litchi to date, simply because of the lack of genomic sequence information. In this study, the litchi fruit transcriptome was first de novo assembled and annotated using RNA-Seq and various computation tools, following which large-scale differentially expressed transcripts in response to shading were explored using DTA. This information will be very helpful for improving the functional genomics studies in litchi and furthering our understanding of the molecular mechanisms behind shade-induced fruit abscission in fruit trees.
Paired-end sequencing and de novo assembly
Summary of the litchi fruit transcriptome
Mean length (bp)
Total nucleotides (bp)
Functional annotation and classification
Summary of functional annotation of assembled unigenes
Public protein database
Number of unigene hits
In addition to GO analysis, Clusters of Orthologous Groups (COG) analysis was used to further evaluate the function of our assembled unigenes. In total, out of 34,029 NR hits, 10,234 hits were aligned with the 24 COG classifications (Table 2). In the 24 COG categories, the cluster for ‘general function prediction only’ (24.40%) represented the largest group, followed by ‘transcription’ (12.54%) and ‘replication, recombination and repair’ (11.81%). Only a small portion of the unigenes were assigned to ‘nuclear structure’ (0.059%) or ‘cell motility’ (0.15%). It is worth noting that a large number of genes were assigned to ‘signal transduction mechanisms’ (9.55%), ‘carbohydrate transport and metabolism’ (8.76%) and ‘cell wall/membrane/envelope biogenesis’ (3.33%) (Figure 2B).
Evaluation of the DTA library and sequencing
To explore the differentially expressed genes induced by shading treatment, DTA profiling was conducted for both shaded and non-shaded samples. After data cleaning, a total of 3.6 million and 3.5 million high-quality tags, representing 126,156 and 128,505 tag entities with unique nucleotide sequences (distinct tags) were generated from shaded and non-shaded libraries respectively. Among these, 58,758 and 59,691 distinct tags identified in 24,181 and 23,881 unigenes, respectively, were unambiguously mapped to the litchi reference transcriptome database obtained in this study (Additional file 3). The distribution of total tags and distinct tags throughout different tag abundance categories showed significantly similar patterns in these two libraries (Additional file 4).
Differentially expressed unigenes between shaded and non-shaded libraries
Effects of shading on fruitlet abscission and expression of selected differentially regulated genes between two TDA libraries
First litchi fruit reference transcriptome generated by RNA-Seq
Computer-based de novo assembly tools (e.g., Trans-ABySS, Velvet-Oases, SOAPdenovo, Trinity) have been developed in conjunction with massively parallel sequencing and their usefulness in transcriptome assembly has been demonstrated in previous studies. Garg et al. found that transcriptome assembly with Velvet followed by Oases program was better than that with other programs based on various assessment parameters (e.g., N50 length, average contig length and sequence similarity with closely related species) . It was also noted that Oases was the most appropriate program for de novo assembly of the wheat (triticum aestivum) transcriptome, when compared with Trans-ABySS and Trinity .
Only the SOAPdenovo program was used for sequence assembly in the current study since many studies have demonstrated that SOAPdenovo is a powerful tool for annotating transcriptomes for plant species across taxa, including sweet potato (Ipomoea batatas) , Chinese bayberry (Myrica rubra) , Siraitia grosvenorii, Taxus mairei, Hevea brasiliensis, Salvia miltiorrhiza, citrus , Momordica cochinchinensis and peanut (Arachis hypogaea) . In this study, a total of 4,007,808,300 bp were sequenced and assembled into a reference transcriptome of litchi fruit, which is also the first library in the Sapindaceae family to be reported. Using the SOAPdenovo software, we generated a total of 57,050 unigenes (≥200 bp), whose average length (601 bp) was longer than that documented in other studies using similar technology and different research organisms, such as Momordica cochinchinensis (388 bp) , Hevea brasiliensis (436 bp) , sweet potato (581 bp) , Salvia miltiorrhiza (467 bp)  and Chinese bayberry (531 bp) . Our results showed that more than 59% of unigenes were matched with functional annotations in the NR database, which is comparable with results from other studies using the same approach (e.g., >61% in Salvia miltiorrhiza, 43% in Momordica cochinchinensis and 58% in Hevea brasiliensis). We found that more unigenes are involved in the starch and sucrose metabolism pathways than in other carbohydrate metabolism pathways that are active in the fruit. These results are consistent with the findings in Siraitia grosvenorii, Chinese bayberry  and date palm (Phoenix dactylifera) . Taken together, the litchi fruit reference transcriptome assembled in this study is comprehensive, accurate and useful for future genetic research of litchi fruit.
Discovery of more than 1,000 unigenes differentially regulated in fruit in response to shading using the DTA method
Fruit is an organ that depends mainly on nutrition imported from elsewhere in the plant, especially photosynthates from leaves . Shading during early fruit development decreases fruit growth and induces fruit abscission by reducing the availability of assimilates . Although many researchers have made great efforts to elucidate the mechanisms of shading effects on fruit growth and abscission, the symptomatic molecular events related to this treatment are still largely unknown [19, 22, 23, 44]. By employing the method of suppression subtractive hybridization, Zhou et al. isolated, from apple fruitlet, 112 shade-response unigenes belonging to eight functional categories . Through the use of apple oligonucleotide microarray, 1,057 selected genes belonging to 15 functional categories, from shade-treated fruit abscission zone, were discovered . In this study, we discovered that the expression of 1,039 unigenes, which fell into 14 groups, was impacted by shading.
A reduction in light intensity is thought to be the most direct effect of shading. It has been hypothesized that photosynthetic rates are significantly altered by different levels of shading [45, 46]. Zhu et al. found that about 6% of shade-response genes are related to photosynthesis and more than 90% of these genes are repressed in the apple fruit abscission zone . In the current study, about 2% of litchi unigenes induced by shading were found to be related to photosynthesis, and only 28% of the genes documented were repressed. The repressed unigenes were mainly involved in light harvesting and electron transport for photosystem I, while the remaining up-regulated unigenes included genes involved in chlorophyll biosynthesis, light harvesting and electron transport for photosystem II, and the chlorophyll DNA binding genes. Fruit setting and development are highly dependent on the carbohydrate supply. In response to shading, 36 differentially expressed unigenes encode enzymes involved in carbohydrate metabolisms, such as glycosidases, hydrolases, and transferases, most were up-regulated, which is in agreement with the findings in apples [22, 23]. The up-regulated expression of these genes is likely a direct or indirect early response of the fruitlet to the carbohydrate shortage.
Shading during early fruit development retards fruit growth by decreasing cell division and expansion . Reductions in cell division have been associated with coordinated alterations in the expression of some key cell-cycle genes [44, 47]. Even though the current study found that only CDK inhibitor was up-regulated, it also showed that six cell-cycle-related unigenes were repressed by shading, including cyclin-dependent kinase (CDK), cyclin-dependent kinase-activating kinase (CAK) and some cell-division control proteins. These results are also consistent with those from previous studies in apples [22, 44].
Some research groups have exploited the fact that dark-induced cell death is usually accompanied by chlorophyll degradation, protein hydrolysis, ROS generation, and cell-wall modification [48, 49]. In this study, eight unigenes related to cell death were up-regulated by shading, particularly cysteine protease and senescence-associated and autophagy-related genes. Shading also had a large impact on protein ubiquitination and degradation; for example, F-box proteins, 26S proteasome subunit proteins, ubiquitin E3 ligase complex including cullin and ubiquitin-conjugating enzymes were positively regulated in response to shading, which is in agreement with the findings in apples . Different levels of dark-induced ROS production have also been reported during the senescence process [49–51]. Only a core NADPH oxidase unigene was increased in the shaded litchi fruitlet. In contrast, some unigenes, involving in the cellular response to oxidative stress, were induced such as a glutathione S-transferase, two glutaredoxin and six thioredoxin genes. It has been previously pointed out that shading has an effect on cell-wall loosening and degradation [22, 23]. However, in our study, only two cell-wall degradation related genes (polygalacturonase and β-d-galactosidase) were induced after shading treatment. Three cellulose synthases, two cell-wall structure genes and a single callose synthase, involved in the secondary cell-wall formation and stress resistance [52–54] were induced. These results indicated that fruits active cellular oxidation resistance and cell-wall formation in response to low-light stress. More than 92% of the differentially expressed unigenes related to stress response were up-regulated, including a large group of resistant genes against diverse stress conditions in this study.
A complex combination of signalling pathways has been proposed to play important roles in the modulation of fruit development. Hormone signalling is key for the coordination of fruit growth, fruit abscission and plant defence [27, 55, 56]. As many as 18 differentially expressed unigenes were hormone-related and involved in the signal transduction pathways of auxins, brassinosteroids, gibberellins, jasmonic acid and polyamine. All of them were up-regulated, except for a single auxin efflux carrier gene, in response to the shading treatment. It is puzzling that no unigene, related to abscisic acid (ABA), ethylene biosynthesis and signalling, was found to be significantly changed by shading. In fact, ABA and ethylene have been implicated in the regulation of stress-induced senescence. Several genes involved in ABA and ethylene biosynthesis and signalling were induced in the shading-treated fruit abscission zone in apples . However, several groups of transcription factors associated with senescence and stress, such as AP2/ERF, NAC, WRKY and zinc finger (C3HC4 and C2H2 type) families, were all largely increased in the shaded fruits, which is consistent with the results on leaf senescence in Arabidopsis by shading . Another group of differentially expressed genes involved in ATP-binding cassette, calcium binding and membrane transport, as well as those involved in the transportation of sugars, amino acids, nucleosides, proteins, lipids, anions and cations. As previously mentioned above, shading during early fruit development reduces canopy photosynthesis and enhances competition between fruits. This results in the alteration of available carbohydrate and transportation of nutrients in favour of shoot growth and at the expense of fruit development [19, 58]. Therefore, 41 out of 55 transcription factors were positively induced by shading in our study, which might reflect the shifting function of the shaded fruitlet, as they become a source tissue for the mobilization of nutrients to other vegetative organs.
Differentially regulated genes involved in litchi fruitlet abscission in response to shading
The imposition of shading treatment during early fruit development promotes fruit abscission in several crops including apples [18, 20], peaches  and grapes . Shade is presumed to induce a transient reduction in the supply of carbohydrates to developing fruit during a period when the fruit is sensitive to such a stress. Shade-induced fruit abscission begins about 5 to 10 days after shading and peaks at 15 days [20–23]. Zhou et al. and Zhu et al. have exploited dark-induced fruit abscission to characterize many altered metabolic pathways and transcript abundance during this senescence process. In our study, shading treatment significantly stimulated fruit drop from 3 days after treatment and induced 88.5% of the fruit to abscise within 7 days, compared with a 62.2% fruit loss within the same time period in the control. Fourteen differentially regulated unigenes, between shaded and non-shaded DTA-tag profiling libraries, were randomly selected in order to compare their expression patterns with the fruit abscission dynamics during shading treatment. The qRT-PCR results of all the selected up-regulated unigenes obtained on day 3 agreed with the expression levels described by DTA. It is possible that transcriptional changes in response to shading occurred and exerted their effects at the early stage of the abscission process and thus became less apparent after 7 days.
Two factors for comparison between expression profiles were considered: whether the mRNA transcript abundance continued to be generally increased or decreased after days of treatment; and whether the relative expression ratio in samples undergoing shaded treatment was significantly higher than that in the controls at specific sampling dates, particularly at the early stage (1, 3 and 5 days) after treatment. According to these two criteria, 11 of the selected differentially regulated shading-responsive genes are regarded as being involved in the process of fruitlet abscission. Six of these genes are also documented as being related to the shedding of plant organs,including CHI[26, 27, 61], FFU, GH3, SAMDC, MYB and PMEI. The remaining five genes (AS, ACL, TNL1, UBCE2 and LL-DAP-AT) have not yet been reported in the literature to be involved in organ shedding. The expression of AS can be induced in a range of low-carbohydrate situations, such as dark periods and sugar-starved treatment [66, 67]. Moreover, an AS gene promoter has been identified as a potential carbohydrate-responsive element in Asparagus. After shading treatment, the expression of AS in litchi fruit was immediately up-regulated and was consistently increased during fruit abscission, indicating that AS expression might be a direct effect of the low-carbohydrate conditions. ACL catalyzes the ATP-dependant reaction of citrate and coenzyme A to form acetyl- coenzyme A and oxaloacetic acid . In our study, the expression of ACL displayed a down-regulated profile after shading treatment, showing that fruit need more energy to maintain the normal physiological functions under low-carbohydrate conditions. Some indirect evidence supported the possibility that ROS generation is related to the organ abscission . Thioredoxin is a ubiquitous small redox active protein in plants that plays a critical role in redox balance regulation through thiol-disulfide exchange reactions. The transcript level of TNL1 in this study was found to be significantly increased as shading-induced fruit abscission progressed, suggesting that it might be involved in regulating fruit drop by influencing stress responses. The transcript levels of ubiquitin-conjugating enzyme genes (GhUBC1 and GhUBC2) from allotetraploid cotton have been found to increase significantly in leaves and flowers at senescence, indicating that GhUBC1/2 may play a role in the degradation of target proteins that function in delaying the senescence program . The litchi UBCE2 unigene in our study was found to be significantly and consistently up-regulated in the shaded fruits, implying that shading-induced fruit abscission induced by shading treatment might be, to some extent, related to the organ senescence and triggered by UBCE2 up-regulated expression. LL-diaminopimelate aminotransferase, discovered in 2006, is a crucial enzyme in the plant lysine biosynthetic pathway, which converts tetrahydrodipicolinate to LL-diaminopimelate in a single step . Lysine is an amino acid essential for peptidoglycan synthesis. There is information in the literature concerning the relationship between LL-DAP-AT or peptidoglycan and organ shedding. However, our study showed that the LL-DAP-AT unigene transcript level in shaded fruit continued to be increased and significantly higher in shaded fruit compared with controls from day 1 to day 5 after shading treatment, indicating that LL-DAP-AT might have some role in regulating fruit abscission.
It is not surprising that some of the genes, differentially regulated in response to shading, are involved in fruitlet abscission since shading during early fruit development induces an increase in fruitlet drop. However, our study was carried out when fruitlet abscission was occurring, and the expression-pattern analysis of genes in both DTA and qRT-PCR experiments were compared between abscission and increased abscission. This would go some way to making accurate transcriptional changes difficult to measure. Thus, comparing between RNA samples from abscission induced by shading versus no abscission in the control should be suggested in the future experiments to find out the candidate genes specifically responsible for fruitlet abscission in litchi.
In this report, we present the sequencing, de novo assembly and analysis of the fruit transcriptome and provide a global description of transcript abundance profiles in fruits from shaded and unshaded trees using next-generation sequencing technology in Litchi chinensis. The transcriptome is described in details in the results section, with an emphasis on annotation using the NR, Swiss-Prot, KEGG and COG databases. We generated more than 53 million paired-end reads and assembled 57,050 unigenes. Using the DTA method, a total of 3.6 million and 3.5 million high-quality tags were obtained from the shaded and non-shaded libraries, respectively. Of these high quality tags, 3,209 regulated tag entities, which were mapped to 1,039 unigenes, showed significant differential expression, suggesting that these potential regulators might function in response to the shading signal. Furthermore, some of shading-induced differentially expressed unigenes were implicated in the fruitlet abscission process in litchi, shedding new light on the molecular mechanisms involved in organ abscission.
Plant materials, sample collection and determination of fruit abscission
Plant materials were collected in 2009 from an orchard located at the South China Agricultural University, Guangzhou, China. For transcriptome sequencing, 16 samples comprising fruitlet (six samples), pericarp (five samples), aril (three samples) and seed (two samples) were collected at different development stages (0, 14, 28, 42, 56, 70 and 84 days after female flower anthesis) from three different litchi cultivars including ‘Hehuadahongli’, ‘Nuomici’ and ‘Heiye’. For DTA library construction, six ‘Nuomici’ litchi trees with uniform vigour and moderate initial crop load were selected. Three were randomly selected as untreated controls while the other three were shaded to 18% of full sun with a neutral-density black-polypropylene shade cloth at 30 days after anthesis. This level of shading has been used to cause fruitlet abscission in many other studies, such as those of apples [23, 58]. Each tree was regarded as a biological replicate. Fruitlet samples were collected at 1, 3, 5 and 7 days after treatment from both control and treated trees. To determine the fruit abscission rate, 10 fruit-bearing shoots from each tree were tagged based on the uniformity of the shoots, and the fruit on tagged shoots were counted at 0, 1, 3, 5 and 7 days after treatment. All of the samples mentioned above were immediately frozen in liquid nitrogen and stored at −80°C until use.
RNA isolation and library construction for transcriptome analysis
Total RNA was isolated using a previously described method  and RNA integrity was confirmed using the 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Total RNA of each sample was then pooled with equivalent quantity. After RNA extraction, poly(A)-containing mRNA was purified using oligo(dT) magnetic beads. To avoid priming bias when synthesizing cDNA, the mRNA was first fragmented into small pieces using divalent cations under elevated temperature. The cleaved RNA fragments were then used to synthesize the first-strand cDNA using reverse transcriptase and random primers, followed by second-strand cDNA synthesis using DNA polymerase I and RNase H. Next, these cDNA fragments underwent processes including end repair, addition of a single ‘A’ base and ligation with adapters. Finally, these products were purified and amplified through PCR to create the library for sequencing.
Transcriptome sequencing, de novo assembly and functional annotation
The Solexa HiSeq™ 2000 was employed to sequence the library (BGI, Shenzhen, China). The cDNA fragments were about 200 bp in length and the fragments were sequenced with a paired-end strategy. The raw reads obtained were pre-processed by removing adaptor sequences, and discarding empty reads and low-quality sequences. All the reads were then used for transcriptome de novo assembly using the short read assembling program SOAPdenovo (version 1.04; http://soap.genomics.org.cn/soapdenovo.html)  at the parameters of “-K 29 -M 2 -L 50” by BGI. Short reads were first assembled into contigs with no gap, and the reads were mapped back to contigs. The resulting contigs were joined into scaffolds using the read mate pairs, with unknown sequences replaced with ‘N’s. Finally, paired-end reads were performed to fill the gap between different scaffolds in order to obtain unigenes with the least ‘N’s and longest sequences. The assembled unigenes were annotated by BLASTx alignment (E-value <0.00001) to protein databases such as the NCBI NR protein database (http://www.ncbi.nlm.nih.gov), the Swiss-Prot protein database (http://www.expasy.ch/sprot), the KEGG pathway database (http://www.genome.jp/kegg) and the COG database (http://www.ncbi.nlm.nih.gov/COG). The best-aligning results from the four databases were chosen to decide the sequence direction of unigenes. Concerning the contradictory results between different databases, a priority order of NR, Swiss-Prot, KEGG and COG was considered. As for other sequences that failed to be annotated to any one of the above databases, the ESTScan program (http://www.ch.embnet.org/software/ESTScan.html) was used to predict the ‘CDS’ (coding regions) and their orientations. Functional annotations of unigenes by GO were carried out using the Blast2GO software , followed by GO functional classifications using the WEGO software  to view the distribution of gene functions in Litchi chinensis at the macro level. For pathway-enrichment analysis, all unigenes were mapped to terms in the KEGG database. The datasets are available from the NCBI Short Read Archive (SRA) with the accession number SRX0255051. The assembled sequences have been stored in the NCBI’s Transcriptome Shotgun Assembly (TSA) database which can be searched for using the Gene-IDs listed in Additional file 9.
DTA library preparation and sequencing
The total RNA of fruitlet samples collected at 1, 3, 5 and 7 days after treatment from both shaded and non-shaded trees, was isolated and pooled with equivalent quantities from the same treatment. Each DTA library included 12 pooled samples from four harvesting times (1, 3, 5 and 7 days after treatment) of three biological replicate trees. The sequence tag library preparation for two samples (non-shaded and shaded) was performed in parallel using the Gene Expression Sample Prep Kit (Illumina).
First, total RNA was used for mRNA purification with magnetic oligo(dT) beads. After reverse transcription and second-strand synthesis, the bead-bound double-stranded cDNA was subsequently digested with endonuclease Nla III which recognizes and cuts at the CATG sites on cDNA. These cDNA fragments with 3′ ends were then purified by magnetic beads precipitation and Illumina adaptor 1 was added to the 5′ ends. The junction of Illumina adaptor 1 and the CATG site is the recognition site of Mme I, which cuts 17 bp downstream of the CATG site, producing tags with adapter 1. After removing 3′ fragments with magnetic beads precipitation, Illumina adaptor 2 was introduced at 3′ ends of tags, forming a tag library. After 15 cycles of linear PCR amplification, 85 base strips were purified by PAGE gel electrophoresis. Finally, these strips were digested, and the single-chain molecules were fixed onto the Solexa Sequencing Chip (flow cell) for sequencing. The datasets are deposited in the NCBI’s SRA database with the accession numbers SRX258094 and SRX258095 for the shaded and non-shaded libraries, respectively.
Analysis and mapping of DTA tags
Sequencing-received raw image data were transformed by base calling into raw reads. Clean tags were obtained by removing the adaptor sequence, empty tags (reads with only adaptor sequences but no tags), low-quality sequences (tags with unknown sequences ‘N’) and low-complexity sequences as well as tags with a copy number of 1 (probably due to sequencing error). All clean tags were mapped back to the litchi fruit transcriptome reference database by using SOAPaligner (Version 2.20) with a maximum of one nucleotide mismatch allowed at the parameters of “-m 0 -x 1000 -s 40 -l 35 -v 3 -r 2”. The meaning and selection principles of the parameters are available on the internet (http://soap.genomics.org.cn/soapaligner.html). Clean tags mapped to reference sequences, from multiple genes, were filtered and the remaining tags became unambiguous clean tags. For gene-expression analysis, the number of unambiguous clean tags was calculated and then normalized to TPM (Transcript Per Million clean tags) [76, 77].
Identification of differentially expressed unigenes
A rigorous algorithm was developed to identify differentially expressed unigenes between the two samples (non-shaded and shaded) based on a previously described method . The FDR was controlled to determine differentially expressed genes . In this study, FDR ≤0.001 and the absolute value of log2 ratio ≥1 (fold-change ≥2) was used as the threshold to judge the significance of unigene expression differences. For functional and pathway-enrichment analysis, all differentially expressed unigenes were mapped to terms in GO and the KEGG database. Unigenes with similar expression patterns usually mean functional correlation. We performed cluster analysis of unigene expression patterns with ‘cluster’  software and ‘Java TreeView’  software. The annotation of all significant unigenes was further supplemented with manual BLASTX, conserved domains, and literature searches. Using this combined information, a functionally driven classification was manually created.
Quantitative real-time PCR analysis
Fourteen genes were randomly selected, from significant differentially expressed DTA tags between the two libraries, to observe their expression patterns during shading treatment by qRT-PCR. The list of gene-specific primers, designed on the 3′ untranslated gene region, is shown in Additional file 8. Purified total RNA (2 μg) from each sample was used to synthesize cDNA in a 20 μL reaction using the M-MLV first-strand synthesis system (Promega). Each qRT-PCR reaction was performed using 40 ng of cDNA in a 20 μL reaction using a SYBR Green-based PCR assay (TaKaRa). The quantitative reactions were performed on the ABI 7500 Real-Time PCR System with 7500 System Software (Applied Biosystems, USA). The comparative C t (threshold cycle) method of quantification was used with litchi EF-1α as the reference which was amplified in parallel with the target unigene(s), allowing unigene expression normalization and providing quantification according to our previous study . The relative fold-change for each of the selected genes was determined for both the control and shaded trees. Three independent biological replicates of each sample and three technical replicates of each biological replicate were used for real-time PCR analysis. PCR amplifications included the following condition: 95°C for 1 min, followed by 40 cycles of 95°C for 5 s, 55°C for 30 s and 72°C for 30 s. Dissociation curves were run to determine the specificity of the amplification reactions. After validation tests, normalization to EF-1α was performed using the ΔΔC t method (Applied Biosystems).
Basic local alignment search tool
Cumulative fruit abscission rate
Clusters of orthologous groups
Digital transcript abundance
False discovery rate
Kyoto encyclopedia of genes and genomes
National center for biotechnology information
- Q20 percentage:
Percentage of bases whose quality was larger than 20 in clean reads
Expressed sequence tags
False discovery rate
Quantitative real-time reverse transcriotion polymerase chain reaction
Reactive oxygen species.
This work was supported by the National Natural Science Foundation of China (31171920) and China Agricultural Research System (CARS-33-11). We are grateful to anonymous reviewers for their comments and English language editing on the manuscript and Dr Yuhai Cui (Agriculture and Agri-Food Canada, Southern Crop Protection and Food Research Centre) for his critical English language editing. We thank Dr. Will Wadlington from the Department of Plant Biology, University of Illinois at Urbana-Champaign and Dr. Rui Xia from the Delaware Biotechnology Institute (USA) for helpful suggestions and critical reading.
- Li JG: The litchi. 2008, Beijing: China Agriculture PressGoogle Scholar
- Li JG: Chapter 2. Litchi. Southern China Fruits Cultivation. Edited by: Chen JZ. 2011, Beijing: Chinese Agriculture PressGoogle Scholar
- Wall MM: Ascorbic acid and mineral composition of longan (Dimocarpus longan), lychee (Litchi chinensis) and rambutan (Nephelium lappaceum) cultivars grown in Hawaii. J Food Compos Anal. 2006, 19: 655-663. 10.1016/j.jfca.2005.12.001.View ArticleGoogle Scholar
- Noh JS, Park CH, Yokozawa T: Treatment with oligonol, a low-molecular polyphenol derived from lychee fruit, attenuates diabetes-induced hepatic damage through regulation of oxidative stress and lipid metabolism. Br J Nutr. 2011, 106: 1013-1022. 10.1017/S0007114511001322.View ArticlePubMedGoogle Scholar
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10: 57-63. 10.1038/nrg2484.PubMed CentralView ArticlePubMedGoogle Scholar
- Anisimov SV: Serial Analysis of Gene Expression (SAGE): 13 years of application in research. Curr Pharm Biotechnol. 2008, 9: 338-350. 10.2174/138920108785915148.View ArticlePubMedGoogle Scholar
- Surget-Groba Y, Montoya-Burgos JI: Optimization of de novo transcriptome assembly from next-generation sequencing data. Genome Res. 2010, 20: 1432-1440. 10.1101/gr.103846.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Zenoni S, Ferrarini A, Giacomelli E, Xumerle L, Fasoli M, Malerba G, Bellin D, Pezzotti M, Delledonne M: Characterization of transcriptional complexity during berry development in Vitis vinifera using RNA-Seq. Plant Physiol. 2010, 152: 1787-1795. 10.1104/pp.109.149716.PubMed CentralView ArticlePubMedGoogle Scholar
- Ono NN, Britton MT, Fass JN, Nicolet CM, Lin DW, Tian L: Exploring the transcriptome landscape of pomegranate fruit peel for natural product biosynthetic gene and SSR marker discovery. J Integr Plant Biol. 2011, 53: 800-813. 10.1111/j.1744-7909.2011.01073.x.View ArticlePubMedGoogle Scholar
- Feng C, Chen M, Xu CJ, Bai L, Yin XR, Li X, Allan AC, Ferguson IB, Chen KS: Transcriptomic analysis of Chinese bayberry (Myrica rubra) fruit development and ripening using RNA-Seq. BMC Genomics. 2012, 13: 19-10.1186/1471-2164-13-19.PubMed CentralView ArticlePubMedGoogle Scholar
- 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
- Stern RA, Gazit S: The reproductive biology of the lychee. Hort Rev. 2002, 28: 393-453.Google Scholar
- Mitra SK, Pereira LS, Pathak PK, Majumdar D: Fruit abscission pattern of lychee cultivars. Acta Hort. 2005, 665: 215-218.View ArticleGoogle Scholar
- Yuan RC, Huang HB: Litchi fruit abscission: its patterns, effect of shading and relation to endogenous abscisic acid. Scientia Hort. 1988, 36: 281-292. 10.1016/0304-4238(88)90063-5.View ArticleGoogle Scholar
- Huang XM, Huang HB, Li JG, Wang HC, Gao FF: Why do detached litchi fruit senesce more rapidly?. Acta Hort. 2001, 558: 217-220.View ArticleGoogle Scholar
- Lakso AN, Wünsche JN, Palmer JW, Corelli Grappadelli L: Measurement and modeling of carbon balance of the apple tree. Hortscience. 1999, 34: 1040-1047.Google Scholar
- Polomski RF, Barden JA, Byers RE, Wolf DD: Apple fruit nonstructural carbohydrates and abscission as influenced by shade and terbacil. J Amer Soc Hort Sci. 1988, 113: 506-511.Google Scholar
- Zibordi M, Domingos S, Grappadelli LC: Thinning apples via shading: an appraisal under field conditions. J Hort Sci Biotech. 2009, 84: 138-144.Google Scholar
- Morandi B, Zibordi M, Losciale P, Manfrini L, Pierpaoli E, Grappadelli LC: Shading decreases the growth rate of young apple fruit by reducing their phloem import. Sci Hortic-Amsterdam. 2011, 127: 347-352. 10.1016/j.scienta.2010.11.002.View ArticleGoogle Scholar
- McArtney S, White M, Latter I, Campbell J: Individual and combined effects of shading and thinning chemicals on abscission and dry-matter accumulation of ‘Royal Gala’ apple fruit. J Hort Sci Biotechnol. 2004, 79: 441-448.Google Scholar
- Kolarič J, Pleško IM, Tojnko S, Stopar M: Apple fruitlet ethylene evolution and MdACO1, MdACS5A, and MdACS5B expression after application of naphthaleneacetic acid, 6-benzyladenine, ethephon, or shading. Hortscience. 2011, 46: 1381-1386.Google 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
- Zhou CJ, Lakso AN, Robinson TL, Gan SS: Isolation and characterization of genes associated with shade-induced apple abscission. Mol Genet Genomics. 2008, 280: 83-92. 10.1007/s00438-008-0348-z.View 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
- Botton A, Eccher G, Forcato C, Ferrarini A, Begheldo M, Zermiani M, Moscatello S, Battistelli A, Velasco R, Ruperti B: Signaling pathways mediating the induction of apple fruitlet abscission. Plant Physiol. 2011, 155: 185-208. 10.1104/pp.110.165779.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
- Meir S, Philosoph-Hadas S, Sundaresan S, Selvaraj KSV, 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
- Gil-Amado JA, Gomez-Jimenez MC: Transcriptome analysis of mature fruit abscission control in olive. Plant Cell Physiol. 2013, 54: 244-269. 10.1093/pcp/pcs179.View ArticlePubMedGoogle Scholar
- Veitch NJ, Johnson PC, Trivedi U, Terry S, Wildridge D, MacLeod A: Digital gene expression analysis of two life cycle stages of the human-infective parasite. Trypanosoma brucei gambiense reveals differentially expressed clusters of co-regulated genes. BMC Genomics. 2010, 11: 124-10.1186/1471-2164-11-124.PubMed CentralView ArticlePubMedGoogle Scholar
- Xiang LX, He D, Dong WR, Zhang YW, Shao JZ: Deep sequencing-based transcriptome profiling analysis of bacteria-challenged Lateolabrax japonicus reveals insight into the immune-relevant genes in marine fish. BMC Genomics. 2010, 11: 472-10.1186/1471-2164-11-472.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang XW, Luan JB, Li JM, Bao YY, Zhang CX, Liu SS: De novo characterization of a whitefly transcriptome and analysis of its gene expression during development. BMC Genomics. 2010, 11: 400-10.1186/1471-2164-11-400.PubMed CentralView ArticlePubMedGoogle Scholar
- Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7: 986-995.PubMedGoogle Scholar
- Garg R, Patel RK, Tyagi AK, Jain M: De novo assembly of chickpea transcriptome using short reads for gene discovery and marker identification. Dna Res. 2011, 18: 53-63. 10.1093/dnares/dsq028.PubMed CentralView ArticlePubMedGoogle Scholar
- Oono Y, Kobayashi F, Kawahara Y, Yazawa T, Handa H, Itoh T, Matsumoto T: Characterisation of the wheat (triticum aestivum L.) transcriptome by de novo assembly for the discovery of phosphate starvation-responsive genes: gene expression in Pi-stressed wheat. BMC genomics. 2013, 14: 77-10.1186/1471-2164-14-77.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang ZY, Fang BP, Chen JY, Zhang XJ, Luo ZX, Huang LF, Chen XL, Li YJ: De novo assembly and characterization of root transcriptome using Illumina paired-end sequencing and development of cSSR markers in sweet potato (Ipomoea batatas). BMC Genomics. 2010, 11: 726-10.1186/1471-2164-11-726.PubMed CentralView ArticlePubMedGoogle Scholar
- Hao DC, Ge GB, Xiao PG, Zhang YY, Yang L: The first insight into the tissue specific taxus transcriptome via Illumina second generation sequencing. PLoS One. 2011, 6: e21220-10.1371/journal.pone.0021220.PubMed CentralView ArticleGoogle Scholar
- Xia ZH, Xu HM, Zhai JL, Li DJ, Luo HL, He CZ, Huang X: RNA-Seq analysis and de novo transcriptome assembly of Hevea brasiliensis. Plant Mol Biol. 2011, 77: 299-308. 10.1007/s11103-011-9811-z.View ArticlePubMedGoogle Scholar
- Hua WP, Zhang Y, Song J, Zhao LJ, Wang ZZ: De novo transcriptome sequencing in Salvia miltiorrhiza to identify genes involved in the biosynthesis of active ingredients. Genomics. 2011, 98: 272-279. 10.1016/j.ygeno.2011.03.012.View ArticleGoogle Scholar
- Liu B, Jiang G, Zhang Y, Li J, Li X, Yue J, Chen F, Liu H, Li H, Zhu S: Analysis of transcriptome differences between resistant and susceptible strains of the citrus red mite Panonychus citri (Acari: Tetranychidae). PloS one. 2011, 6: e28516-10.1371/journal.pone.0028516.PubMed CentralView ArticlePubMedGoogle Scholar
- Hyun TK, Rim Y, Jang H, Kim CH, Park J, Kumar R, Lee S, Kim BC, Bhak J, Nguyen-Quoc B: De novo transcriptome sequencing of Momordica cochinchinensis to identify genes involved in the carotenoid biosynthesis. Plant Mol Biol. 2012, 79: 413-427. 10.1007/s11103-012-9919-9.View ArticlePubMedGoogle Scholar
- Zhang JA, Liang S, Duan JL, Wang J, Chen SL, Cheng ZS, Zhang Q, Liang XQ, Li YR: De novo assembly and characterisation of the transcriptome during seed development, and generation of genic-SSR markers in Peanut (Arachis hypogaea L.). BMC Genomics. 2012, 13: 90-10.1186/1471-2164-13-90.PubMed CentralView ArticlePubMedGoogle Scholar
- Yin YX, Zhang XW, Fang YJ, Pan LL, Sun GY, Xin CQ, Abdullah MMB, Yu XG, Hu SN, Al-Mssallem IS: High-throughput sequencing-based gene profiling on multi-staged fruit development of date palm (Phoenix dactylifera, L.). Plant Mol Biol. 2012, 78: 617-626. 10.1007/s11103-012-9890-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Hieke S, Menzel CM, Ludders P: Shoot development, chlorophyll, gas exchange and carbohydrates in lychee seedlings (Litchi chinensis). Tree Physiol. 2002, 22: 947-953. 10.1093/treephys/22.13.947.View ArticlePubMedGoogle Scholar
- Dash M, Johnson LK, Malladi A: Severe shading reduces early fruit growth in apple by decreasing cell production and expansion. J Am Soc Hortic Sci. 2012, 137: 275-282.Google Scholar
- Cai ZQ: Shade delayed flowering and decreased photosynthesis, growth and yield of Sacha Inchi (Plukenetia volubilis) plants. Ind Crop Prod. 2011, 34: 1235-1237. 10.1016/j.indcrop.2011.03.021.View ArticleGoogle Scholar
- Cartechini A, Palliotti A: Effect of shading on vine morphology and productivity and leaf gas exchange characteristics in grapevines in the field. Am J Enol Viticult. 1995, 46: 227-234.Google Scholar
- Malladi A, Johnson LK: Expression profiling of cell cycle genes reveals key facilitators of cell production during carpel development, fruit set, and fruit growth in apple (Malus× domestica Borkh.). J Exp Bot. 2011, 62: 205-219. 10.1093/jxb/erq258.PubMed CentralView ArticlePubMedGoogle Scholar
- Brouwer B, Ziolkowska A, Bagard M, Keech O, Gardeström P: The impact of light intensity on shade-induced leaf senescence. Plant Cell Environ. 2012, 35: 1084-1098. 10.1111/j.1365-3040.2011.02474.x.View ArticlePubMedGoogle Scholar
- Rosenwasser S, Rot I, Sollner E, Meyer AJ, Smith Y, Leviatan N, Fluhr R, Friedman H: Organelles contribute differentially to reactive oxygen species-related events during extended darkness. Plant Physiol. 2011, 156: 185-201. 10.1104/pp.110.169797.PubMed CentralView ArticlePubMedGoogle Scholar
- Pastori GM, Del Río LA: Natural senescence of pea leaves (an activated oxygen-mediated function for peroxisomes). Plant Physiol. 1997, 113: 411-418.PubMed CentralPubMedGoogle Scholar
- Causin HF, Roberts IN, Criado MV, Gallego SM, Pena LB, Ríos MDC, Barneix AJ: Changes in hydrogen peroxide homeostasis and cytokinin levels contribute to the regulation of shade-induced senescence in wheat leaves. Plant Sci. 2009, 177: 698-704. 10.1016/j.plantsci.2009.08.014.View ArticleGoogle Scholar
- Hernández-Blanco C, Feng DX, Hu J, Sánchez-Vallet A, Deslandes L, Llorente F, Berrocal-Lobo M, Keller H, Barlet X, Sánchez-Rodríguez C: Impairment of cellulose synthases required for Arabidopsis secondary cell wall formation enhances disease resistance. The Plant Cell Online. 2007, 19: 890-903. 10.1105/tpc.106.048058.View ArticleGoogle Scholar
- Priyanka B, Sekhar K, Reddy VD, Rao KV: Expression of pigeonpea hybrid-proline-rich protein encoding gene (CcHyPRP) in yeast and Arabidopsis affords multiple abiotic stress tolerance. Plant Biotechnol J. 2010, 8: 76-87. 10.1111/j.1467-7652.2009.00467.x.View ArticlePubMedGoogle Scholar
- Piršelová B, Matušíková I: Callose: the plant cell wall polysaccharide with multiple biological functions. Acta Physiol Plant. 2013, 35: 635-644. 10.1007/s11738-012-1103-y.View ArticleGoogle Scholar
- Robert-Seilaniantz A, Grant M, Jones JD: Hormone crosstalk in plant disease and defense: more than just jasmonate-salicylate antagonism. Annu Rev Phytopathol. 2011, 49: 317-343. 10.1146/annurev-phyto-073009-114447.View ArticlePubMedGoogle Scholar
- Davies PJ: The plant hormones: their nature, occurrence, and functions. Plant Hormones: Biosynthesis, Signal Transduction, Action. Edited by: Davies PJ. 2010, New York: Springer, 3View ArticleGoogle Scholar
- Parlitz S, Kunze R, Mueller-Roeber B, Balazadeh S: Regulation of photosynthesis and transcription factor expression by leaf shading and re-illumination in Arabidopsis thaliana leaves. J Plant Physiol. 2011, 168: 1311-1319. 10.1016/j.jplph.2011.02.001.View ArticlePubMedGoogle Scholar
- Bepete M, Lakso AN: Differential effects of shade on early-season fruit and shoot growth rates in ‘Empire’ apple. HortSci. 1998, 33: 823-825.Google Scholar
- Byers RE, Lyons CG, Yoder KS, Barden JA, Young RW: Peach and apple thinning by shading and photosynthetic inhibition. J Hortic Sci. 1985, 60: 465-472.Google Scholar
- Ferree DC, McArtney SJ, Scurlock DM: Influence of irradiance and period of exposure on fruit set of French-American hybrid grapes. J Am Soc Hortic Sci. 2001, 126: 283-290.Google Scholar
- Del Campillo E, Lewis LN: Occurrence of 9.5 cellulase and other hydrolases in flower reproductive organs undergoing major cell wall disruption. Plant Physiol. 1992, 99: 1015-1020. 10.1104/pp.99.3.1015.PubMed CentralView ArticlePubMedGoogle Scholar
- Ruiz R, Guardiola JL: Carbohydrate and mineral nutrition of orange fruitlets in relation to growth and abscission. Physiol Plant. 1994, 90: 27-36. 10.1111/j.1399-3054.1994.tb02188.x.View ArticleGoogle Scholar
- Kuang JF, Wu JY, Zhong HY, Li CQ, Chen JY, Lu WJ, Li JG: Carbohydrate stress affecting fruitlet abscission and expression of genes related to auxin signal transduction pathway in litchi. Int J Mol Sci. 2012, 13: 16084-16103. 10.3390/ijms131216084.PubMed CentralView ArticlePubMedGoogle Scholar
- Gil-Amado JA, Gomez-Jimenez MC: Regulation of polyamine metabolism and biosynthetic gene expression during olive mature-fruit abscission. Planta. 2012, 235: 1221-1237. 10.1007/s00425-011-1570-1.View ArticlePubMedGoogle Scholar
- Lashbrook CC, Cai S: Cell wall remodeling in Arabidopsis stamen abscission zones: temporal aspects of control inferred from transcriptional profiling. Plant Signal Behav. 2008, 3: 733-736. 10.4161/psb.3.9.6489.PubMed CentralView ArticlePubMedGoogle Scholar
- Davies KM, Seelye JF, Irving DE, Borst WM, Hurst PL, King GA: Sugar regulation of harvest-related genes in asparagus. Plant Physiol. 1996, 111: 877-883. 10.1104/pp.111.3.877.PubMed CentralView ArticlePubMedGoogle Scholar
- Buchanan-Wollaston V, Page T, Harrison E, Breeze E, Lim PO, Nam HG, Lin JF, Wu SH, Swidzinski J, Ishizaki K: Comparative transcriptome analysis reveals significant differences in gene expression and signalling pathways between developmental and dark/starvation-induced senescence in Arabidopsis. Plant J. 2005, 42: 567-585. 10.1111/j.1365-313X.2005.02399.x.View ArticlePubMedGoogle Scholar
- Winichayakul S, Moyle RL, Coupe SA, Davies KM, Farnden KJF: Analysis of the asparagus (Asparagus officinalis) asparagine synthetase gene promoter identifies evolutionarily conserved cis-regulatory elements that mediate Suc-repression. Func Plant Biol. 2004, 31: 63-72. 10.1071/FP03179.View ArticleGoogle Scholar
- Fatland BL, Ke J, Anderson MD, Mentzen WI, Cui LW, Allred CC, Johnston JL, Nikolau BJ, Wurtele ES: Molecular characterization of a heteromeric ATP-citrate lyase that generates cytosolic acetyl-coenzyme A in Arabidopsis. Plant Physiol. 2002, 130: 740-756. 10.1104/pp.008110.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang XD, Jenkins JN, Callahan FE, Creech RG, Si Y, McCarty JC, Saha S, Ma DP: Molecular cloning, differential expression, and functional characterization of a family of class I ubiquitin-conjugating enzyme (E2) genes in cotton (Gossypium). Biochim Biophys Acta. 2003, 1625: 269-279. 10.1016/S0167-4781(02)00623-1.View ArticlePubMedGoogle Scholar
- Hudson AO, Singh BK, Leustek T, Gilvarg C: An LL-diaminopimelate aminotransferase defines a novel variant of the lysine biosynthesis pathway in plants. Plant Physiol. 2006, 140: 292-301.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu JY, Peng G, Li CQ, Lu WJ, Wang ZH, Li JG: A new rapid and effective method for RNA isolation from litchi tissues of fruitlet and abscission zone. Acta Hortic Sin. 2011, 38: 1191-1196.Google Scholar
- Li RQ, Zhu HM, Ruan J, Qian WB, Fang XD, Shi ZB, Li YR, Li ST, Shan G, Kristiansen K: De novo assembly of human genomes with massively parallel short read sequencing. Genome Res. 2010, 20: 265-272. 10.1101/gr.097261.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMedGoogle Scholar
- Ye J, Fang L, Zheng HK, Zhang Y, Chen J, Zhang ZJ, Wang J, Li ST, Li RQ, Bolund L: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34: W293-W297. 10.1093/nar/gkl031.PubMed CentralView ArticlePubMedGoogle Scholar
- Morrissy AS, Morin RD, Delaney A, Zeng T, McDonald H, Jones S, Zhao Y, Hirst M, Marra MA: Next-generation tag sequencing for cancer gene expression profiling. Genome Res. 2009, 19: 1825-1835. 10.1101/gr.094482.109.PubMed CentralView ArticlePubMedGoogle Scholar
- ‘t Hoen PAC, Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RH, de Menezes RX, Boer JM, van Ommen GB, den Dunnen JT: Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008, 36: e141-10.1093/nar/gkn705.PubMed CentralView ArticlePubMedGoogle Scholar
- Benjamini Y, Yekutieli D: The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001, 29: 1165-1188. 10.1214/aos/1013699998.View ArticleGoogle Scholar
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.PubMed CentralView ArticlePubMedGoogle Scholar
- Saldanha AJ: Java Treeview-extensible visualization of microarray data. Bioinformatics. 2004, 20: 3246-3248. 10.1093/bioinformatics/bth349.View ArticlePubMedGoogle Scholar
- Zhong HY, Chen JW, Li CQ, Chen L, Wu JY, Chen JY, Lu WJ, Li JG: Selection of reliable reference genes for expression studies by reverse transcription quantitative real-time PCR in litchi under different experimental conditions. Plant Cell Rep. 2011, 30: 641-653. 10.1007/s00299-010-0992-8.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.