Abstract
Mountain pepper (Litseacubeba (Lour.) Pers.) (Lauraceae) is an important industrial crop as an ingredient in cosmetics, pesticides, food additives and potential biofuels. These properties are attributed to monoterpenes and sesquiterpenes. However, there is still no integrated model describing differentially expressed genes (DEGs) involved in terpenoid biosynthesis during the fruit development of L. cubeba. Here, we performed digital gene expression (DGE) using the Illumina NGS platform to evaluated changes in gene expression during fruit development in L. cubeba. DGE generated expression data for approximately 19354 genes. Fruit at 60 days after flowering (DAF) served as the control, and a total of 415, 1255, 449 and 811 up-regulated genes and 505, 1351, 1823 and 1850 down-regulated genes were identified at 75, 90, 105 and 135 DAF, respectively. Pathway analysis revealed 26 genes involved in terpenoid biosynthesis pathways. Three DEGs had continued increasing or declining trends during the fruit development. The quantitative real-time PCR (qRT-PCR) results of five differentially expressed genes were consistent with those obtained from Illumina sequencing. These results provide a comprehensive molecular biology background for research on fruit development, and information that should aid in metabolic engineering to increase the yields of L. cubeba essential oil.