TPIA2: an updated tea plant information archive for Camellia genomics.

The genus Camellia consists of about 200 species, which include many economically important species widely used for making tea, ornamental flowers and edible oil. Here, we present an updated tea plant information archive for Camellia genomics (TPIA2; http://tpia.teaplants.cn) by integrating more novel large-scale genomic, transcriptomic, metabolic and genetic variation datasets as well as a variety of useful tools. Specifically, TPIA2 hosts all currently available and well assembled 10 Camellia genomes and their comprehensive annotations from three major sections of Camellia. A collection of 15 million SNPs and 950 950 small indels from large-scale genome resequencing of 350 diverse tea accessions were newly incorporated, followed by the implementation of a novel 'Variation' module to facilitate data retrieval and analysis of the functionally annotated variome. Moreover, 116 Camellia transcriptomes were newly assembled and added, leading to a significant extension of expression profiles of Camellia genes to 13 developmental stages and eight abiotic/biotic treatments. An updated 'Expression' function has also been implemented to provide a comprehensive gene expression atlas for Camellia. Two novel analytic tools (e.g. Gene ID Convert and Population Genetic Analysis) were specifically designed to facilitate the data exchange and population genomics in Camellia. Collectively, TPIA2 provides diverse updated valuable genomic resources and powerful functions, and will continue to be an important gateway for functional genomics and population genetic studies in Camellia.

Integrated physiological, metabolite and proteomic analysis reveal the glyphosate stress response mechanism in tea plant (Camellia sinensis).

Glyphosate residues can tremendously impact the physiological mechanisms of tea plants, thus threatening tea security and human health. Herein, integrated physiological, metabolite, and proteomic analyses were performed to reveal the glyphosate stress response mechanism in tea plant. After exposure to glyphosate (≥1.25 kg ae/ha), the leaf ultrastructure was damaged, and chlorophyll content and relative fluorescence intensity decreased significantly. The characteristic metabolites catechins and theanine decreased significantly, and the 18 volatile compounds content varied significantly under glyphosate treatments. Subsequently, tandem mass tags (TMT)-based quantitative proteomics was employed to identify the differentially expressed proteins (DEPs) and to validate their biological functions at the proteome level. A total of 6287 proteins were identified and 326 DEPs were screened. These DEPs were mainly catalytic, binding, transporter and antioxidant active proteins, involved in photosynthesis and chlorophyll biosynthesis, phenylpropanoid and flavonoid biosynthesis, sugar and energy metabolism, amino acid metabolism, and stress/defense/detoxification pathway, etc. A total of 22 DEPs were validated by parallel reaction monitoring (PRM), demonstrating that the protein abundances were consistent between TMT and PRM data. These findings contribute to our understanding of the damage of glyphosate to tea leaves and molecular mechanism underlying the response of tea plants to glyphosate.

Evolutionary Landscape of Tea Circular RNAs and Its Contribution to Chilling Tolerance of Tea Plant.

Chilling stress threatens the yield and distribution pattern of global crops, including the tea plant (Camellia sinensis), one of the most important cash crops around the world. Circular RNA (circRNA) plays roles in regulating plant growth and biotic/abiotic stress responses. Understanding the evolutionary characteristics of circRNA and its feedbacks to chilling stress in the tea plant will help to elucidate the vital roles of circRNAs. In the current report, we systematically identified 2702 high-confidence circRNAs under chilling stress in the tea plant, and interestingly found that the generation of tea plant circRNAs was associated with the length of their flanking introns. Repetitive sequences annotation and DNA methylation analysis revealed that the longer flanking introns of circRNAs present more repetitive sequences and higher methylation levels, which suggested that repeat-elements-mediated DNA methylation might promote the circRNAs biogenesis in the tea plant. We further detected 250 differentially expressed circRNAs under chilling stress, which were functionally enriched in GO terms related to cold/stress responses. Constructing a circRNA-miRNA-mRNA interaction network discovered 139 differentially expressed circRNAs harboring potential miRNA binding sites, which further identified 14 circRNAs that might contribute to tea plant chilling responses. We further characterized a key circRNA, CSS-circFAB1, which was significantly induced under chilling stress. FISH and silencing experiments revealed that CSS-circFAB1 was potentially involved in chilling tolerance of the tea plant. Our study emphasizes the importance of circRNA and its preliminary role against low-temperature stress, providing new insights for tea plant cold tolerance breeding.

Effects of Different Chemicals on Sexual Regulation in Persimmon (Diospyros kaki Thunb.) Flowers.

Research on crop sexuality is important for establishing systems for germplasm innovation and cultivating improved varieties. In this study, androecious persimmon trees were treated with various concentrations of ethrel (100, 500, and 1,000 mg/L) and zeatin (1, 5, and 10 mg/L) to investigate the morphological, physiological, and molecular characteristics of persimmon. Ethrel at 1,000 mg/L and zeatin at 10 mg/L both significantly reduced the stamen length and pollen grain diameter in androecious trees. Ethrel treatment also led to reduced stamen development with degenerated cellular contents; zeatin treatment promoted the development of arrested pistils via maintaining relatively normal mitochondrial morphology. Both treatments altered carbohydrate, amino acid, and endogenous phytohormone contents, as well as genes associated with hormone production and floral organ development. Thereafter, we explored the combined effects of four chemicals, including ethrel and zeatin, as well as zebularine and 5-azacytidine, both of which are DNA methylation inhibitors, on androecious persimmon flower development. Morphological comparisons showed that stamen length, pollen viability, and pollen grain diameter were significantly inhibited after combined treatment. Large numbers of genes involving in carbohydrate metabolic, mitogen-activated protein kinase (MAPK) signaling, and ribosome pathways, and metabolites including uridine monophosphate (UMP) and cyclamic acid were identified in response to the treatment, indicating complex regulatory mechanisms. An association analysis of transcriptomic and metabolomic data indicated that ribosomal genes have distinct effects on UMP and cyclamic acid metabolites, explaining how male floral buds of androecious persimmon trees respond to these exogenous chemicals. These findings extend the knowledge concerning sexual differentiation in persimmon; they also provide a theoretical basis for molecular breeding, high-yield cultivation, and quality improvement in persimmon.

Comparative transcriptomic analysis unveils the deep phylogeny and secondary metabolite evolution of 116 Camellia plants.

Camellia plants include more than 200 species of great diversity and immense economic, ornamental, and cultural values. We sequenced the transcriptomes of 116 Camellia plants from almost all sections of the genus Camellia. We constructed a pan-transcriptome of Camellia plants with 89 394 gene families and then resolved the phylogeny of genus Camellia based on 405 high-quality low-copy core genes. Most of the inferred relationships are well supported by multiple nuclear gene trees and morphological traits. We provide strong evidence that Camellia plants shared a recent whole genome duplication event, followed by large expansions of transcription factor families associated with stress resistance and secondary metabolism. Secondary metabolites, particularly those associated with tea quality such as catechins and caffeine, were preferentially heavily accumulated in the Camellia plants from section Thea. We thoroughly examined the expression patterns of hundreds of genes associated with tea quality, and found that some of them exhibited significantly high expression and correlations with secondary metabolite accumulations in Thea species. We also released a web-accessible database for efficient retrieval of Camellia transcriptomes. The reported transcriptome sequences and obtained novel findings will facilitate the efficient conservation and utilization of Camellia germplasm towards a breeding program for cultivated tea, camellia, and oil-tea plants.

Metabolite Profiling and Transcriptome Analysis Revealed the Conserved Transcriptional Regulation Mechanism of Caffeine Biosynthesis in Tea and Coffee Plants.

Caffeine is a characteristic bioactive compound in tea and coffee plants, which is synthesized and accumulated extensively in leaves and seeds. However, little is known about the regulatory mechanism of caffeine synthesis in plants. This study compared the caffeine metabolite between tea and coffee plants. We found that tea leaves contained significantly higher caffeine than coffee leaves, which is perhaps due to more members of N-methyltransferase (NMT) genes as well as higher expression levels in tea plants. Substantial numbers of transcription factors were predicted to be involved in caffeine biosynthesis regulation, combining weighted gene co-expression network analysis and the cis-element of NMT promoter analysis in tea and coffee plants. Furthermore, analysis of the transcription factors from the caffeine-related modules suggested that the regulatory mechanism of caffeine biosynthesis was probably partly conservative in tea and coffee plants. This study provides an essential resource for the regulatory mechanism of caffeine biosynthesis in plants.

Divergent DNA methylation contributes to duplicated gene evolution and chilling response in tea plants.

The tea plant (Camellia sinensis) is a thermophilic cash crop and contains a highly duplicated and repeat-rich genome. It is still unclear how DNA methylation regulates the evolution of duplicated genes and chilling stress in tea plants. We therefore generated a single-base-resolution DNA methylation map of tea plants under chilling stress. We found that, compared with other plants, the tea plant genome is highly methylated in all three sequence contexts, including CG, CHG and CHH (where H = A, T, or C), which is further proven to be correlated with its repeat content and genome size. We show that DNA methylation in the gene body negatively regulates the gene expression of tea plants, whereas non-CG methylation in the flanking region enables a positive regulation of gene expression. We demonstrate that transposable element-mediated methylation dynamics significantly drives the expression divergence of duplicated genes in tea plants. The DNA methylation and expression divergence of duplicated genes in the tea plant increases with evolutionary age and selective pressure. Moreover, we detect thousands of differentially methylated genes, some of which are functionally associated with chilling stress. We also experimentally reveal that DNA methyltransferase genes of tea plants are significantly downregulated, whereas demethylase genes are upregulated at the initial stage of chilling stress, which is in line with the significant loss of DNA methylation of three well-known cold-responsive genes at their promoter and gene body regions. Overall, our findings underscore the importance of DNA methylation regulation and offer new insights into duplicated gene evolution and chilling tolerance in tea plants.

The Reference Genome of Tea Plant and Resequencing of 81 Diverse Accessions Provide Insights into Its Genome Evolution and Adaptation.

Tea plant is an important economic crop, which is used to produce the world's oldest and most widely consumed tea beverages. Here, we present a high-quality reference genome assembly of the tea plant (Camellia sinensis var. sinensis) consisting of 15 pseudo-chromosomes. LTR retrotransposons (LTR-RTs) account for 70.38% of the genome, and we present evidence that LTR-RTs play critical roles in genome size expansion and the transcriptional diversification of tea plant genes through preferential insertion in promoter regions and introns. Genes, particularly those coding for terpene biosynthesis proteins, associated with tea aroma and stress resistance were significantly amplified through recent tandem duplications and exist as gene clusters in tea plant genome. Phylogenetic analysis of the sequences of 81 tea plant accessions with diverse origins revealed three well-differentiated tea plant populations, supporting the proposition for the southwest origin of the Chinese cultivated tea plant and its later spread to western Asia through introduction. Domestication and modern breeding left significant signatures on hundreds of genes in the tea plant genome, particularly those associated with tea quality and stress resistance. The genomic sequences of the reported reference and resequenced tea plant accessions provide valuable resources for future functional genomics study and molecular breeding of improved cultivars of tea plants.

Transcriptomic analyses reveal a systemic defense role of the uninfested adjacent leaf in tea plant (Camellia sinensis) attacked by tea geometrids (Ectropis obliqua).

To get a more detailed understanding of the interaction between tea plant (Camellia sinensis) and tea geometrids (Ectropis obliqua), transcriptomic profile in undamaged adjacent leaf (TGL) of tea geometrids fed local leaves (LL) was investigated for the first time. Here, approximately 245 million clean reads contained 39.39 Gb of sequence data were obtained from TGL. Further analysis revealed that systemic response was induced in TGL after tea geometrids feeding on LL, although the defense response was weaker than that in LL. The differentially expressed genes (DEGs) identification analysis showed little overlap of DEGs between TGL and LL. Comparative transcriptome analysis suggested that JA signal regulated resistant pathway was induced in LL; whereas primary metabolism pathway was activated in TGL in response to tea geometrids feeding. This study reveals a novel resistance mechanism of TGL to tea geometrids feeding.

Theanine transporters identified in tea plants (Camellia sinensis L.).

Theanine, a unique non-proteinogenic amino acid, is an important component of tea, as it confers the umami taste and relaxation effect of tea as a beverage. Theanine is primarily synthesized in tea roots and is subsequently transported to young shoots, which are harvested for tea production. Currently, the mechanism for theanine transport in the tea plant remains unknown. Here, by screening a yeast mutant library, followed by functional analyses, we identified the glutamine permease, GNP1 as a specific transporter for theanine in yeast. Although there is no GNP1 homolog in the tea plant, we assessed the theanine transport ability of nine tea plant amino acid permease (AAP) family members, with six exhibiting transport activity. We further determined that CsAAP1, CsAAP2, CsAAP4, CsAAP5, CsAAP6, and CsAAP8 exhibited moderate theanine affinities and transport was H+ -dependent. The tissue-specific expression of these six CsAAPs in leaves, vascular tissues, and the root suggested their broad roles in theanine loading and unloading from the vascular system, and in targeting to sink tissues. Furthermore, expression of these CsAAPs was shown to be seasonally regulated, coincident with theanine transport within the tea plant. Finally, CsAAP1 expression in the root was highly correlated with root-to-bud transport of theanine, in seven tea plant cultivars. Taken together, these findings support the hypothesis that members of the CsAAP family transport theanine and participate in its root-to-shoot delivery in the tea plant.

Optimized sequencing depth and de novo assembler for deeply reconstructing the transcriptome of the tea plant, an economically important plant species.

BACKGROUND: Tea is the oldest and among the world's most popular non-alcoholic beverages, which has important economic, health and cultural values. Tea is commonly produced from the leaves of tea plants (Camellia sinensis), which belong to the genus Camellia of family Theaceae. In the last decade, many studies have generated the transcriptomes of tea plants at different developmental stages or under abiotic and/or biotic stresses to investigate the genetic basis of secondary metabolites that determine tea quality. However, these results exhibited large differences, particularly in the total number of reconstructed transcripts and the quality of the assembled transcriptomes. These differences largely result from limited knowledge regarding the optimized sequencing depth and assembler for transcriptome assembly of structurally complex plant species genomes. RESULTS: We employed different amounts of RNA-sequencing data, ranging from 4 to 84 Gb, to assemble the tea plant transcriptome using five well-known and representative transcript assemblers. Although the total number of assembled transcripts increased with increasing sequencing data, the proportion of unassembled transcripts became saturated as revealed by plant BUSCO datasets. Among the five representative assemblers, the Bridger package shows the best performance in both assembly completeness and accuracy as evaluated by the BUSCO datasets and genome alignment. In addition, we showed that Bridger and BinPacker harbored the shortest runtimes followed by SOAPdenovo and Trans-ABySS. CONCLUSIONS: The present study compares the performance of five representative transcript assemblers and investigates the key factors that affect the assembly quality of the transcriptome of the tea plants. This study will be of significance in helping the tea research community obtain better sequencing and assembly of tea plant transcriptomes under conditions of interest and may thus help to answer major biological questions currently facing the tea industry.

The tea plant reference genome and improved gene annotation using long-read and paired-end sequencing data.

Tea is a globally consumed non-alcohol beverage with great economic importance. However, lack of the reference genome has largely hampered the utilization of precious tea plant genetic resources towards breeding. To address this issue, we previously generated a high-quality reference genome of tea plant using Illumina and PacBio sequencing technology, which produced a total of 2,124 Gb short and 125 Gb long read data, respectively. A hybrid strategy was employed to assemble the tea genome that has been publicly released. We here described the data framework used to generate, annotate and validate the genome assembly. Besides, we re-predicted the protein-coding genes and annotated their putative functions using more comprehensive omics datasets with improved training models. We reassessed the assembly and annotation quality using the latest version of BUSCO. These data can be utilized to develop new methodologies/tools for better assembly of complex genomes, aid in finding of novel genes, variations and evolutionary clues associated with tea quality, thus help to breed new varieties with high yield and better quality in the future.

Transcriptomic and Phytochemical Analyses Reveal Root-Mediated Resource-Based Defense Response to Leaf Herbivory by Ectropis oblique in Tea Plant ( Camellia sinensis).

Leaf herbivory on tea plants ( Camellia sinensis) by tea geometrids ( Ectropis oblique) severely threaten the yield and quality of tea. In previous work, we found that local defense response was induced in damaged leaves by geometrids at the transcriptome level. Here, we investigated the systemic response triggered in undamaged roots and the potential role of roots in response to leaf herbivory. Comparative transcriptome analysis and carbohydrate dynamics indicated that leaf herbivory activated systemic carbon reallocation to enhance resource investment for local secondary metabolism. The crucial role of jasmonic acid and the involvement of other potential hormone signals for local and systemic signaling networks were supported by phytohormone quantification and dynamic expression analysis of phytohormone-related genes. This work represents a deep understanding of the interaction of tea plants and geometrids from the perspective of systems biology and reveals that tea plants have evolved an intricate root-mediated resource-based resistance strategy to cope with geometrid attack.

Tea Plant Information Archive: a comprehensive genomics and bioinformatics platform for tea plant.

Tea is the world's widely consumed nonalcohol beverage with essential economic and health benefits. Confronted with the increasing large-scale omics-data set particularly the genome sequence released in tea plant, the construction of a comprehensive knowledgebase is urgently needed to facilitate the utilization of these data sets towards molecular breeding. We hereby present the first integrative and specially designed web-accessible database, Tea Plant Information Archive (TPIA; http://tpia.teaplant.org). The current release of TPIA employs the comprehensively annotated tea plant genome as framework and incorporates with abundant well-organized transcriptomes, gene expressions (across species, tissues and stresses), orthologs and characteristic metabolites determining tea quality. It also hosts massive transcription factors, polymorphic simple sequence repeats, single nucleotide polymorphisms, correlations, manually curated functional genes and globally collected germplasm information. A variety of versatile analytic tools (e.g. JBrowse, blast, enrichment analysis, etc.) are established helping users to perform further comparative, evolutionary and functional analysis. We show a case application of TPIA that provides novel and interesting insights into the phytochemical content variation of section Thea of genus Camellia under a well-resolved phylogenetic framework. The constructed knowledgebase of tea plant will serve as a central gateway for global tea community to better understand the tea plant biology that largely benefits the whole tea industry.

Circular RNA architecture and differentiation during leaf bud to young leaf development in tea (Camellia sinensis).

MAIN CONCLUSION: Circular RNA (circRNA) discovery, expression patterns and experimental validation in developing tea leaves indicates its correlation with circRNA-parental genes and potential roles in ceRNA interaction network. Circular RNAs (circRNAs) have recently emerged as a novel class of abundant endogenous stable RNAs produced by circularization with regulatory potential. However, identification of circRNAs in plants, especially in non-model plants with large genomes, is challenging. In this study, we undertook a systematic identification of circRNAs from different stage tissues of tea plant (Camellia sinensis) leaf development using rRNA-depleted circular RNA-seq. By combining two state-of-the-art detecting tools, we characterized 3174 circRNAs, of which 342 were shared by each approach, and thus considered high-confidence circRNAs. A few predicted circRNAs were randomly chosen, and 20 out of 24 were experimental confirmed by PCR and Sanger sequencing. Similar in other plants, tissue-specific expression was also observed for many C. sinensis circRNAs. In addition, we found that circRNA abundances were positively correlated with the mRNA transcript abundances of their parental genes. qRT-PCR validated the differential expression patterns of circRNAs between leaf bud and young leaf, which also indicated the low expression abundance of circRNAs compared to the standard mRNAs from the parental genes. We predicted the circRNA-microRNA interaction networks, and 54 of the differentially expressed circRNAs were found to have potential tea plant miRNA binding sites. The gene sets encoding circRNAs were significantly enriched in chloroplasts related GO terms and photosynthesis/metabolites biosynthesis related KEGG pathways, suggesting the candidate roles of circRNAs in photosynthetic machinery and metabolites biosynthesis during leaf development.

Characterization and Alternative Splicing Profiles of the Lipoxygenase Gene Family in Tea Plant (Camellia sinensis).

Oxylipins, including jasmonic acid (JA) and volatiles, are important for signaling in plants, and these are formed by the lipoxygenase (LOX) enzyme family. There is a large gap in understanding of the underlying molecular basis of their roles in tea plants. Here, we identified 11 CsLOX genes from the tea plant (Camellia sinensis), and characterized their phylogeny, gene structure and protein features into three subclasses. We then examined their enzymatic activities, LOX expression and alternative splicing of transcripts during development and in response to abiotic or biotic stresses in tea plants. In vitro expressed protein assays showed that the CsLOX2, 3 and 9 enzymatically function to produce 9/13-HPOT, 13-HPOT and 9-HPOT, respectively. CsLOX2 and CsLOX9 green fluorescent protein (GFP) fusion proteins localized to chloroplasts and the cytoplasm, respectively. RNA sequencing, quantitative reverse transcription-PCR and Northern blot analysis suggested that CsLOX5, 6 and 9 were predominantly expressed in seeds, flowers and roots, respectively. CsLOX2, 3, 4, 6 and 7 were up-regulated after attack by the insect Ectropis oblique, while CsLOX1 was induced after infection with the pathogen Glomerella cingulata. CsLOX3, 7 and 10 were up-regulated by JA but not ABA or salicylic acid. Long-term cold stress down-regulated CsLOX expression while a short duration of cold induced the expression of CsLOX1, 6 and 7. Alternatively spliced transcripts of six CsLOX genes were dynamically regulated through time and varied in relative abundances under the investigated stresses; we propose a mechanism of competing or compensating regulation between isoforms. This study improves our understanding of evolution of LOXs and regulation of their diverse functions in plants.

Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality.

Tea, one of the world's most important beverage crops, provides numerous secondary metabolites that account for its rich taste and health benefits. Here we present a high-quality sequence of the genome of tea, Camellia sinensis var. sinensis (CSS), using both Illumina and PacBio sequencing technologies. At least 64% of the 3.1-Gb genome assembly consists of repetitive sequences, and the rest yields 33,932 high-confidence predictions of encoded proteins. Divergence between two major lineages, CSS and Camellia sinensis var. assamica (CSA), is calculated to ∼0.38 to 1.54 million years ago (Mya). Analysis of genic collinearity reveals that the tea genome is the product of two rounds of whole-genome duplications (WGDs) that occurred ∼30 to 40 and ∼90 to 100 Mya. We provide evidence that these WGD events, and subsequent paralogous duplications, had major impacts on the copy numbers of secondary metabolite genes, particularly genes critical to producing three key quality compounds: catechins, theanine, and caffeine. Analyses of transcriptome and phytochemistry data show that amplification and transcriptional divergence of genes encoding a large acyltransferase family and leucoanthocyanidin reductases are associated with the characteristic young leaf accumulation of monomeric galloylated catechins in tea, while functional divergence of a single member of the glutamine synthetase gene family yielded theanine synthetase. This genome sequence will facilitate understanding of tea genome evolution and tea metabolite pathways, and will promote germplasm utilization for breeding improved tea varieties.

Ectopic expression of EuWRI1, encoding a transcription factor in E. ulmoides, changes the seeds oil content in transgenic tobacco.

The WRINKLED1 (WRI1) gene is a well-established key transcriptional regulator involved in the regulation of fatty acid biosynthesis in developing seeds. In this study, a new WRI1 gene was isolated from seeds of Eucommia ulmoides and named EuWRI1. A close link between gibberellins signaling and EuWRI1 gene expression was suggested in this study. Functional characterization of EuWRI1 was elucidated through seed-specific expression in tobacco. In transgenic tobacco, the expression of EuWRI1 in eight independent transgenic lines was detected by semiquantitative RT-PCR. The relative mRNA accumulation of genes encoding enzymes involved in fatty acid biosynthesis (biotin carboxyl carrier protein and keto-ACP synthase 1) was also assayed in tobacco seeds. Analysis of the seeds oil content and starch content indicated that the transgenic lines showed a significant increase in seeds oil content, whereas starch content decreased significantly. Further analysis of the fatty acid composition revealed that palmitic acid (16:0), linoleic acid (18:2) and linolenic acid (18:3) increased significantly in seeds of transgenic tobacco lines, but stearic acid (18:0) levels significantly declined. © 2018 American Institute of Chemical Engineers Biotechnol. Prog., 34:337-346, 2018.

Characterization of Glycolytic Pathway Genes Using RNA-Seq in Developing Kernels of Eucommia ulmoides.

Eucommia ulmoides Oliver, the only member of the Eucommiaceae family, is a rare and valuable tree used to produce a highly valued traditional Chinese medicine and contains α-linolenic acid (ALA) up to 60% of the total fatty acids in the kernels (embryos). Glycolysis provides both cellular energy and the intermediates for other biosynthetic processes. However, nothing was known about the molecular basis of the glycolytic pathway in E. ulmoides kernels. The purposes of this study were to identify novel genes of E. ulmoides related to glycolytic metabolism and to analyze the expression patterns of selected genes in the kernels. Transcriptome sequencing based on the Illumina platform generated 96,469 unigenes in four cDNA libraries constructed using RNAs from 70 and 160 days after flowering kernels of both low- and high-ALA varieties. We identified and characterized the digital expression of 120 unigenes coding for 24 protein families involved in kernel glycolytic pathway. The expression levels of glycolytic genes were generally higher in younger kernels than in more mature kernels. Importantly, several unigenes from kernels of the high-ALA variety were expressed more than those from the low-ALA variety. The expression of 10 unigenes encoding key enzymes in the glycolytic pathway was validated by qPCR using RNAs from six kernel stages of each variety. The qPCR data were well consistent with their digital expression in transcriptomic analyses. This study identified a comprehensive set of genes for glycolytic metabolism and suggests that several glycolytic genes may play key roles in ALA accumulation in the kernels of E. ulmoides.