Revelation of enzyme/transporter-mediated metabolic regulatory model for high-quality terpene accumulation in developing fruits of Lindera glauca

Lindera glauca with rich resource and fruit terpene has emerged as potential material for utilization in China, but different germplasms show a variation for essential oil content and volatile profiling. This work aimed to determine key regulators (enzymes or transporters) and unravel mechanism of governing high production of essential oil of L. glauca fruit (EO-LGF). Temporal analysis of fruit growth and EO-LGF accumulation (yield, volatile compounds and contents) during development revealed a notable change in the contents of EO-LGF and its 45 compounds in developing fruits, and the major groups were monoterpene and sesquiterpene, showing good antioxidant and antimicrobial activities. To highlight molecular mechanism that govern such difference in terpene content and compound in developing fruits, Genome-wide assay was used to annotate 104 genes for terpene-synthesis pathway based on recent transcriptome data, and the comparative associations of terpene accumulative amount with gene transcriptional level were conducted on developing fruits to identify some crucial determinants (enzymes and transporters) with metabolic regulation model for high-quality terpene accumulation, involving in carbon allocation (sucrose cleavage, glycolysis and OPP pathway), metabolite transport, isoprene precursor production, C5-unit formation (MEP and MVA pathways), and mono-/sesqui-terpene synthesis. Our findings may present strategy for engineering terpene accumulation for utilization.

Characterization of 20 complete plastomes from the tribe Laureae (Lauraceae) and distribution of small inversions

Lindera Thunb. (Lauraceae) consists of approximately 100 species, mainly distributed in the temperate and tropical regions of East Asia. In this study, we report 20 new, complete plastome sequences including 17 Lindera species and three related species, Actinodaphne lancifolia, Litsea japonica and Sassafras tzumu. The complete plastomes of Lindera range from 152,502 bp (L. neesiana) to 154,314 bp (L. erythrocarpa) in length. Eleven small inversion (SI) sites are documented among the plastomes. Six of the 11 SI sites are newly reported and they locate in rpoB-trnC, psbC-trnS, petA-psbJ, rpoA and ycf2 regions. The distribution patterns of SIs are useful for species identification. An average of 83 simple sequence repeats (SSRs) were detected in each plastome. The mono-SSRs accounted for 72.7% of total SSRs, followed by di- (12.4%), tetra- (9.4%), tri- (4.2%), and penta-SSRs (1.3%). Of these SSRs, 64.6% were distributed in an intergenic spacer (IGS) region. In addition, 79.8% of the SSRs are located in a large single copy (LSC) region. In contrast, almost no SSRs are distributed in inverted repeat (IR) regions. The SSR loci are useful to identifying species but the phylogenetic value is low because the majority of them show autapomorphic status or highly homoplastic characteristics. The nucleotide diversity (Pi) values also indicated the conserved nature of the IR region compared to LSC and small single copy (SSC) regions. Five spacer regions with high Pi values, trnH-psbA, petA-psbJ and ndhF-rpl32, rpl32-trnL and Ψycf1-ndhF, have a potential use for the molecular identification study of Lindera and related species. Lindera species form a paraphyletic group in the plastome tree because of the inclusion of related genera such as Actinodaphne, Laurus, Litsea and Neolitsea. A former member of tribe Laureae, Sassafras, forms a clade with the tribe Cinnamomeae. The SIs do not affect the phylogenetic relationship of Laureae. This result indicated that ancient plastome captures may have contribute to the mixed intergeneric relationship of Laureae. Alternatively, the result may indicate that the morphological characters defined the genera of Lauraceae originated for several times.

De novo assembly and characterization of the floral transcriptome of an economically important tree species, Lindera glauca (Lauraceae), including the development of EST-SSR markers for population genetics

Lindera glauca (Lauraceae) is an economically important East Asian forest tree characterized by a dioecy in China and apomixis in Japan. However, patterns of population genetic diversity and structure of this species remain unknown for this species due to a lack of efficient molecular markers. In this study, we employed Illumina sequencing to analyze the transcriptomes of the female and male flower buds of L. glauca. We retrieved 59,753 and 75,075 unigenes for the female and male buds, respectively. Based on sequence similarity, 44,379 (74.27 %) unigenes for the female and 45,414 (60.49 %) unigenes for the male were matched to public databases. We identified 11,127 putative differentially expressed genes between the female and male buds and 20,048 expressed sequence tag-simple sequence repeats (EST-SSRs). From 3147 primer pairs designed successfully, 120 were selected for validation of polymorphism, and 13 could reliably amplify polymorphic bands and exhibited moderate levels of genetic diversity (e.g., N _A = 4.42; H _E = 0.56) when surveyed across 96 individuals of altogether six L. glauca populations from China and Japan. One of the three population genetic clusters identified in China was fixed in Japan, suggesting a historical population bottleneck following island immigration. The present study has generated a wealth of transcriptome data for future functional genomic research focused on the variable reproductive system of L. glauca (dioecy, apomixis) as well as EST-SSR markers for population genetics studies and its intriguing evolutionary shift from dioecy to apomixis in the wake of island colonization.

Integrated transcriptome sequencing and dynamic analysis reveal carbon source partitioning between terpenoid and oil accumulation in developing Lindera glauca fruits

Lindera glauca fruits (LGF) with the abundance of terpenoid and oil has emerged as a novel specific material for industrial and medicinal application in China, but the complex regulatory mechanisms of carbon source partitioning into terpenoid biosynthetic pathway (TBP) and oil biosynthetic pathway (OBP) in developing LGF is still unknown. Here we perform the analysis of contents and compositions of terpenoid and oil from 7 stages of developing LGF to characterize a dramatic difference in temporal accumulative patterns. The resulting 3 crucial samples at 50, 125 and 150 days after flowering (DAF) were selected for comparative deep transcriptome analysis. By Illumina sequencing, the obtained approximately 81 million reads are assembled into 69,160 unigenes, among which 174, 71, 81 and 155 unigenes are implicated in glycolysis, pentose phosphate pathway (PPP), TBP and OBP, respectively. Integrated differential expression profiling and qRT-PCR, we specifically characterize the key enzymes and transcription factors (TFs) involved in regulating carbon allocation ratios for terpenoid or oil accumulation in developing LGF. These results contribute to our understanding of the regulatory mechanisms of carbon source partitioning between terpenoid and oil in developing LGF, and to the improvement of resource utilization and molecular breeding for L. glauca.