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.