De novo transcriptome assembly of the eight major organs of Sacha Inchi (Plukenetia volubilis) and the identification of genes involved in α-linolenic acid metabolism

Sacha Inchi (Plukenetia volubilis L.), which belongs to the Euphorbiaceae, has been considered a new potential oil crop because of its high content of polyunsaturated fatty acids in its seed oil. The seed oil especially contains high amounts of α-linolenic acid (ALA), which is useful for the prevention of various diseases. However, little is known about the genetic information and genome sequence of Sacha Inchi, which has largely hindered functional genomics and molecular breeding studies. In this study, a de novo transcriptome assembly based on transcripts sequenced in eight major organs, including roots, stems, shoot apexes, mature leaves, male flowers, female flowers, fruits, and seeds of Sacha Inchi was performed, resulting in a set of 124,750 non-redundant putative transcripts having an average length of 851 bp and an N50 value of 1909 bp. Organ-specific unigenes analysis revealed that the most organ-specific transcripts are found in female flowers (2244 unigenes), whereas a relatively small amount of unigenes are detected to be expressed specifically in other organs with the least in stems (24 unigenes). A total of 42,987 simple sequence repeats (SSRs) were detected, which will contribute to the marker assisted selection breeding of Sacha Inchi. We analyzed expression of genes related to the α-linolenic acid metabolism based on the de novo assembly and annotation transcriptome in Sacha Inchi. It appears that Sacha Inchi accumulates high level of ALA in seeds by strong expression of biosynthesis-related genes and weak expression of degradation-related genes. In particular, the up-regulation of FAD3 and FAD7 is consistent with high level of ALA in seeds of Sacha Inchi compared with in other organs. Meanwhile, several transcription factors (ABI3, LEC1 and FUS3) may regulate key genes involved in oil accumulation in seeds of Sacha Inchi. The transcriptome of major organs of Sacha Inchi has been sequenced and de novo assembled, which will expand the genetic information for functional genomic studies of Sacha Inchi. In addition, the identification of candidate genes involved in ALA metabolism will provide useful resources for the genetic improvement of Sacha Inchi and the metabolic engineering of ALA biosynthesis in other plants.


Background
Sacha Inchi (Plukenetia volubilis L.), a member of the Euphorbiaceae family [1], is native to the rainforests of South America [2] . Sacha Inchi is also known as Inca peanut, wild peanut, Sacha peanut or mountain peanut, having been cultivated for centuries by the indigenous population [3]. And thus, the magnitude of several evolutionary forces like selection, genetic drift and gene flow on population structure of Sacha Inchi was decided by strong anthropogenic influence [4]. Based on cytological study, the most common chromosome number is 2n = 58 [5]. Sacha Inchi seeds contain 41-54% oil [2,3,6], which is characterized predominantly by high levels of polyunsaturated fatty acids (PUFAs), especially α-linolenic acid (ALA, C18:3 cis Δ9, 12, 15, ω-3) and linoleic acid (LA, C18:2 cis Δ9, 12, ω-6), which represent approximately 50 and 35% of the total oil, respectively [2,7,8]. In addition, Sacha Inchi seeds contain substantial amounts of total tocopherols (137 mg/100 g) [7], which are a class of chemical compounds consisting of various methylated phenols and display strong antioxidant activity [7,9,10]. Considerable amounts of phytosterols (75.7-86.2 mg/100 g) and 15 polyphenolic compounds, belonging to phenyl alcohol, flavonoid, seicoridoid, and lignan classes, were positively identified in seeds [7,11]. Particularly, condensed tannins are the main family of phenolic compounds which might be indicative of potential high antioxidant properties [12]. Leaf extracts were characterized by phenolic compounds, steroids, and/or terpenoids, which resulted in the antioxidant activities in leaves of Sacha Inchi [13]. These results indicate that Sacha Inchi could be considered as an important material for production of antioxidant phenolic compounds and phytosterols. ALA and LA are essential fatty acids [14], which are useful in the prevention of coronary heart disease, hypertension, diabetes, arthritis, high cholesterol, cancer, and inflammatory and autoimmune disorders [15][16][17][18][19]. Serum parameters in rats treated with Sacha Inchi oil indicated lower levels of cholesterol and triglycerides, and higher levels of high density lipoprotein in comparison with the control group [20]. The lipid profiles of patients with hypercholesterolemia who intake seed oil of Sacha Inchi for four months indicated a decreased in the values of total cholesterol and non-esterified fatty acids, and a rise in high density lipoprotein and the insulin levels [21]. Moreover, ALA and LA are also critical for the development of infants during pregnancy and breastfeeding periods [8]. However, the amount of ALA in most human diets is insufficient. The ALA/LA ratio in the Western diet is approximately 1:15, which much lower than the ratio recommended by the WHO (1:2 to 1:6) [15,22]. The Western diet is very high in ω-6 FAs and low in ω-3 FAs because of the unwise recommendation to substitute ω-6 FAs for saturated fats to lower serum cholesterol concentrations, resulting in the production of products rich in ω-6 and poor in ω-3 FAs [23]. Thus, it is necessary to increase the intake of ALA, which can be extracted from the seeds of some oil plants, such as Sacha Inchi. Triacylglycerol (TAG) is the primary unit of energy storage in eukaryotic cells. Normally, TAGs are derived either from the glycerol-3-phosphate (G3P) pathway (also known as the Kennedy pathway) or the acyl-dihydroxyacetone phosphate (acyl-DHAP) pathway [24,25]. TAG degradation through oxidation (primarily beta-oxidation) releases FAs [26]. One aim of this study is to analyze the expression of genes involved in the G3P pathway and beta-oxidation, which are major pathways for TAG biosynthesis and degradation in most tissues or organisms.
Thus far, most studies on Sacha Inchi have dealt with plant development and physiology [27][28][29], the characterization of seed oil [3,6,8], in vitro regeneration systems [30], and potential applications in biofuel production [31] and in cosmetic, pharmaceutical, and food industries [7,[32][33][34]. However, the genetic information and molecular mechanisms underlying ALA metabolism in Sacha Inchi have rarely been studied, especially absence of a genome. Only one transcriptome analysis and the expression of oil biosynthesis genes in Sacha Inchi seeds have been published [35,36]. In the transcriptome analysis, the developing seeds of two stages from two-year-old Sacha Inchi were sequenced and unigenes that may be involved in de novo FA and triacylglycerol biosynthesis were identified [35]. In another study, expression profiles of genes controlling unsaturated fatty acids biosynthesis and oil deposition in developing seeds of Sacha Inchi were investigated by quantitative real-time PCR [35,36]. However, only a few genes contributing to the high level of ALA have been characterized in Sacha Inchi. Our results provide priority candidates for future research.
In this study, we sequenced and de novo assembled the transcriptome of 8 major organs of Sacha Inchi using next-generation sequencing technology. The assembled transcriptome sequences will expand the genetic information for functional genomic studies of Sacha Inchi. In addition, the identification of candidate genes involved in ALA metabolism will provide useful resources for the genetic improvement of Sacha Inchi and the metabolic engineering of ALA metabolism in other plants.

Transcriptome sequencing and de novo assembly of Sacha Inchi
To comprehensively construct the transcriptome of Sacha Inchi, eight major organs, including roots, stems, shoot apexes, mature leaves, male flowers, female flowers, fruits, and seeds, were sampled for RNA isolation. Distinct cDNA libraries of those organs were constructed and sequenced, resulting in a total of 164 G 150-bp paired-end raw reads. After the removal of adapters, poly-N-containing reads and low-quality sequences from the raw data, approximately 162 G clean reads were retained and used for transcriptome assembly and analysis (Additional file 1: Table S1). The Trinity [37] assembly program, which has been shown to be the best single k-mer assembler for RNA-Seq short reads, was used for de novo assembly [38]. As a result, an assembly of 349,951 contigs was established for the Sacha Inchi transcriptome. To reduce redundancy and potential assembly errors, the candidate unigenes that most likely had the longest ORFs (Open Reading Frames) were chosen from the assembly result, and then those transcripts were filtered by their fragments per kilobase per million mapped base pairs of sequenced (FPKM) values less than 0.1. Finally, a set of 124,750 unigenes with an average length of 851 bp and an N50 value of 1909 bp was obtained (Table 1). We have compared our data with the seed transcriptome reported by Wang et al. [35] using BLASTn (E-value≤1e-10). The alignment rate of the reported seed transcriptome to our data is 93.8%, indicating the 6.2% of the seed transcriptome was not detected in our dataset, which probably resulted from the seed transcriptome reported by Wang et al. [35] containing transcripts of seeds at two developmental stages, whereas our dataset containing transcripts of seeds at a single developmental stage. The size distribution of unigenes is shown in Additional file 2: Figure S1a. Principal component analysis (PCA) was conducted using R package, the distance assessment reveals that all three independent biological replicates of each sample have good reproducibility, and the seed showed the most distinctive expression patterns in all tested organs (Additional file 3: Figure S2), which is in accord with the result in Arabidopsis [39].
To evaluate and summarize the reliability and quality of the assembly, the clean reads were mapped back to the Trinity-assembled transcriptome. The overall alignment rate was 83.39%, indicating that a high-quality de novo assembled transcriptome was obtained.

Gene ontology and KEGG pathway analysis of the Sacha Inchi transcriptome
We utilized the best BLASTx hit from the NR database to functionally classify the Sacha Inchi unigenes. The best hits were subjected to Blast2GO [41,42] to analyze the gene ontology (GO) terms and enzyme commission (EC) numbers. The result showed that 45,319 unigenes were assigned to 240,186 GO term annotations, with an average of 5.3 GO terms for each unigene. A total of 45,319 unigene gene functions were described under three main divisions (biological processes, cellular components and molecular functions). The predominant group in each of the biological processes, cellular components and molecular functions was metabolic process (GO: 0008152), cell (GO: 0005623) and catalytic activity (GO: 0003824), respectively (Additional file 6: Figure S5). To further understand the biological functions and interactions of transcripts, the unigenes of assembled sequences were assigned by the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The result showed that a total of 24,678 unigenes were involved in 323 different pathways (see Additional file 7).

Analysis of organ-specific unigenes
Examination of the unigenes expressed in each organ type revealed that female flowers expressed the most unigenes (16%), while stems and seeds expressed the least number of unigenes (11% each) (Additional file 8: Figure S6a). We designated a gene as the organ-specific with the criteria that the expression value (FPKM) was at least 10 in one organ while less than 1 in other organs. Using the criteria, 495, 24,249,390,195,2244,287, and 388 unigenes were specifically found in roots, stems, shoot apexes, mature leaves, male flowers, female flowers, fruits, and seeds, respectively (Additional file 8: Figure S6b). Detailed annotation and FPKM information of organ-specific unigenes are listed in Additional file 9. To evaluate the functional properties of these organ-specific unigenes, GO and KEGG annotation enrichment analyses were carried out. As a result, transcripts were significantly enriched in 7, 18, 10, 18 and 22 GO terms in male flowers, female flowers, roots, seeds and stems, respectively (Additional file 10). However, no significantly enriched GO term was found in fruits or mature leaves. Through KEGG pathway enrichment analysis of female flower, most of those pathways are related to metabolism pathways, including "Amino acid metabolism," "Biosynthesis of other secondary metabolites," "Carbohydrate metabolism," "Energy metabolism," "Glycan biosynthesis and metabolism," "Lipid metabolism," "Metabolism of cofactors and vitamins," "Metabolism of terpenoids and polyketides" and "Nucleotide metabolism," followed by pathways related to genetic information processing (Additional file 11: Figure S7), and the genes specially expressed in female flowers were significantly enriched in 17 pathways (Additional file 12). Genes specifically expressed in female flowers were enriched not only in pathways related to sugar, lipid, and amino acid metabolism but also in plant hormone signal transduction and phagosome pathways, which is consistent with the GO analysis(Additional file 11: Figure S7).

Identification and characterization of genes involved in ALA metabolism
Based on the de novo assembly and annotation of the Sacha Inchi transcriptome, a total of 211 transcripts were identified as candidate unigenes for 13 enzymes of ALA biosynthesis (Table 2 and Additional file 13). The ALA biosynthesis pathway ( Fig. 1) is composed of two sub-pathways: fatty acid (FA) elongation and desaturation [43,44]. The initiation and acyl chain elongation steps of de novo FA biosynthesis use acetyl-CoA. Acetyl-CoA is initially catalyzed by acetyl-CoA carboxylase (ACC) to form malonyl-CoA. Then, malonyl-ACP, which is the primary substrate for subsequent elongation, is generated from malonyl-CoA by the malonyl-CoA: ACP malonyl-transferase (MCMT). Four enzymes played a significant role during the addition of two carbons: ketoacyl-ACP synthase III (KAS III), ketoacyl-ACP reductase (KAR), 3-hydroxyacyl-ACP dehydratase (HAD, EC: 4.2.1.-) and enoyl-ACP reductase (EAR) [45,46]. After four reactions, C4:0-ACP, which is the substrate for further elongation, is produced. Next, ketoacyl-ACP synthase I (KAS I) is used for the elongation from C4 to C16. However, the reaction from C16 to C18 in the ALA biosynthesis pathway is catalyzed by ketoacyl-ACP synthase II (KAS II), and the enzyme stearoyl-ACP desaturase (SAD) removes two hydrogen atoms from stearic acid (18C:0) to form oleic acid (18C:1) during the process of unsaturated FA formation. The enzyme fatty acid desaturase 2 (FAD2) and fatty acid desaturase 6 (FAD6) desaturate oleic acid (C18:1Δ 9 ) to generate LA [47]. Fatty acid desaturase 3 (FAD3), FAD7, and FAD8 enzyme genes act on the ω6 fatty acid LA [48,49] to catalyze the biosynthesis of ALA from LA in Sacha Inchi. In the ALA biosynthesis pathway, the genes encoding FATA, KAS II, FAD2, FAD3, and FAD7 were identified and showed significant up-regulation in seeds compared with other organs (Fig. 2), indicating that the unsaturated FA biosynthesis pathway might be blocked in those organs. However, the genes encoding HAD, MCMT (Malonyl-CoA ACP transferase), and KAS I were not found, which might result from the single development stage of seeds sampled in this study.
It is worthy to note that phosphatidylcholine diacylglycerol cholinephosphotransferase (PDCT), which was not found in the previous report of transcriptome of Sacha Inchi seeds [35], but has been shown to play an important role in PUFA accumulation by catalyzing the interconversion between phosphatidylcholine (PC) and diacylglycerol (DAG) [50][51][52], was found to have high level of expression in female flowers and seeds in this study (Fig. 3).

Analysis of differentially expressed genes (DEGs) involved in TAG synthesis
Differentially expressed genes (DEGs) were analyzed in the eight major organs, especially the genes involved in TAG synthesis, glycolysis/gluconeogenesis pathway (ko00010) and pentose phosphate pathway (ko00030). Genes having a false discovery rate (FDR) value < 0.001 and |log 2 (fold change)| ≥2 found by edgeR were regarded as DEGs [53]. As a result, a total of 24,196 DEGs were detected in our dataset. Among these DEGs, thirty-three genes were involved in TAG biosynthesis (Table 2 and Additional file 13). Glycerol-3-phosphate acyltransferase (GPAT) catalyze the acylation of glycerol-3-phosphate to produce 1-acyl-sn-glycerol-3-phosphate (lysophosphatididic acid, LPA) (Fig. 4) [54,55]. Phosphatidic acid (PA) is synthesized de novo from the acylation of LPA in a reaction catalyzed by acyl-CoA: lysophosphatidic acid   . Diacylglycerol (DAG) is produced from PA catalyzed by phosphatidate phosphatase (PP). The DAG is then available for two reactions: diacylglycerol O-acyltransferase (DGAT) transfer acyl-CoAs to the sn-3 position of DAG to produce TAG; alternatively phospholipid: diacylglycerol acyltransferase (PDAT) transfer the sn-2 acyl group from PA to DAG, forming TAG [56,57]. In addition to the DAG derived from Kennedy pathway, a second PC-derived DAG pool has been identified [58]. The GPAT and DGAT were up-regulated in female flowers, while LPAAT, PP, and PDAT were up-regulated in shoot apexes, seeds and fruits, respectively (Fig. 2). The glycolysis/gluconeogenesis pathway produces glycerol-3-phosphate which is the source of triacylglycerol (TAG) biosynthesis, and the pentose phosphate pathway is a process of glucose turnover that produces precursor materials and NADPH for FA biosynthesis [54][55][56][57][58]. One hundred and thirty-seven and 64 unigenes were identified as involved in glycolysis/gluconeogenesis pathway and pentose phosphate pathway, respectively, most of which were up-regulated in female flowers and seeds (Additional file 14: Figure S8). In addition to the key enzymes, we found that some transcription factors related to TAG biosynthesis were highly up-regulated in the seed. ABSCISIC ACID-INSENSITIVE 3 (ABI3), LEAFY COTYLEDON1 (LEC1) and B3-domain transcription factor FUSCA3 (FUS3) showed a 100-fold or more expression difference in seeds compared to other organs (Table 3). These transcription factors probably regulate the seed oil by directly or indirectly regulating FA biosynthesis or TAG accumulation [59].

Identification and characterization of genes involved in FA catabolism
The major pathway of FA catabolism, beta-oxidation pathway, is composed of acyl-CoA oxidase (ACX), ketoacyl-CoA thiolase (KAT) and multifunctional protein (MFP) (Fig. 5). Based on the de novo assembly and annotation of the Sacha Inchi transcriptome, a total of 89 transcripts were identified as candidate unigenes for 4 enzymes of fatty acid catabolism pathway, including 27 unigenes encoding long-chain acyl-CoA synthetase (LACS) that catalyzes initial reactions of FA catabolism, 16 unigenes encoding ACX, 15 unigenes encoding MFP, and 31 unigenes encoding KAT (Table 2 and Additional file 13). Most of the FA catabolism-related genes exhibited weak expression in seeds, for example, ACX and MFP that catalyze the first and second steps of the β-oxidation of fatty acids, generating acetyl-CoA and energy [60,61] were down-regulated in seeds (Fig. 2).

Validation of gene expression profiles using qRT-PCR
To experimentally confirm the expression patterns of genes identified by transcriptome sequencing, 9 key unigenes involved in ALA metabolism (Fig. 6a), 7 unigenes detected in our data and 8 unigenes showing organ-specific expression (Fig. 6b) in Sacha Inchi were chosen for qRT-PCR analysis. The detailed information of the selected unigenes is presented in Additional file 16. The results showed that the expression patterns of the most genes tested by qPCR and RNA-Seq were consistent (Fig. 6a, b). Overall, a highly significant correlation (Pearson correlation coefficient r = 0.835) existed between qRT-PCR and RNA-Seq results regarding the ratios of gene expression levels (Additional file 17: Figure S9), suggesting that our transcriptome data reflect the expression patterns of most genes in Sacha Inchi.

Discussion
In this study, we provide a large number of organ-specific unigenes by analyzing the gene expression profile of individual organs. As shown in Additional file 8: Figure S6b, it is apparent that different organs expressed distinct sets of genes. The vast majority of the organ-specific transcripts are found in female flowers, whereas relative small amount of unigenes are predicted to be expressed specifically in other organs. These findings indicate that the metabolisms, such as lipid metabolism, amino acid metabolism in female flowers are active and involve many more specifically expressed genes. This result also provides a large number of candidate specific promoters of female flowers. Considering ALA has potential applications in food and pharmaceutical industries [7], the main objective of our research was to comprehensively study the ALA metabolism and investigate the molecular basis of this pathway in Sacha Inchi. The omega-3 fatty acid desaturases, especially FAD3, are key enzymes for the formation of ALA from the desaturation of LA [62,63]. Overexpression of FAD3 in roots and seeds led to the increase of ALA in a fad3-2 mutant of A. thaliana [64]. Meanwhile, overexpression of FAD7 increased levels of ALA that led to series of physiological alterations, such as electrolyte leakage and malondialdehyde contents in tomato [65]. Taken together, we suggested that high transcript levels of FAD3 and FAD7 is consistent with high level of ALA content in seeds of Sacha Inchi. Our results showed that the predominant genes related to these core pathways are sequence conserved. Furthermore, previous research showed that when the flax PDCTs were co-expressed with FAD2 and FAD3, PUFA levels increased in Saccharomyces cerevisiae [50]. The PDCT from flax are capable of increasing C18-PUFA levels substantially in metabolically engineered yeast and transgenic A. thaliana seeds [50]. These data strongly indicate that those genes appear to play an important role in the determination of PUFA content in TAG synthesis. DEGs that directly and indirectly regulate TAG biosynthesis were identified in our data. In addition to a number of genes related to glycolysis/ gluconeogenesis pathway and pentose phosphate pathway were up-regulated expression in female flower and seed, the key gene phosphatidate phosphatase (PP) which catalyzes DAG from PA had a high expression level in seed. The transcription factors also regulate the TAG biosynthesis directly or indirectly in plants. Overexpression of maize LEC1 can increase seed oil as much as 48% [66]. ABI3 and FUS3 are key regulators in phase transition and seed development, maturation and dormancy, and thereby indirectly regulate oil biosynthesis [26]. The genes involved in ALA biosynthesis are regulated by miRNAs as shown in previous study: KASII and KASIII are regulated by the miR159; KAR is regulated by the miR156b, miR156c, miR156g and miR6029; FATA is regulated by miR801, miR298, miR1430 and miR828; FATB is regulated by miR555 and miR113; and SAD is regulated by miR2163 [67].
Based on the gene expression pattern analysis, genes coding enzymes related to the ALA biosynthesis strongly express in the seeds, whereas FA catabolism-related

Conclusions
In the present study, the complete transcriptome of Sacha Inchi was de novo-assembled and annotated for the first time, generating a total of 124,750 non-redundant transcripts, of which 70,142 could be functionally annotated. Among the eight organs analyzed in this study, the largest number of specifically expressed genes was found in female flowers while the least was found in stems. We identified 211 unigenes and 89 unigenes potentially involved in the ALA biosynthesis and FA catabolism pathways, respectively. Compared with other organs, most of the unigenes related to ALA biosynthesis metabolism were up-regulated, whereas most of those enzymes related to FA catabolism were down-regulated in seeds of Sacha Inchi. In particular, the up-regulation of FAD3 and FAD7 may play an important role in high level accumulation of ALA in seeds of Sacha Inchi. Some transcription factors are highly up-regulated in seeds, which are potentially related to TAG accumulation. The transcriptome data reported here provide the foundation for the functional genomics research and genetic improvement of Sacha Inchi.
In conclusion, we present a high-quality transcriptome sequence for the Sacha Inchi. The sequences of genes related to organ-specific and ALA metabolism are obtained based on large-scale transcriptomic data, which Values are means ± standard deviations (n = 3). b Organ-specific genes. HPG, hypothetical protein GLOINDRAFT_348631; GDS, glutamine-dependent NAD(+) synthetase; CBL, CBL-interacting serine/threonine-protein kinase; PP, peroxidase 64 precursor; HPA, hypothetical protein AALP_AAs53585U000100; HSP90, HSP90 co-chaperone CPR7; LP, phospholipid transfer protein; RABC2a, ras-relatedprotein RABC2a-like. FF, female flowers; Fr, fruits; MF, male flowers; ML, mature leaves; Ro, roots; SA, shoot apexes; Se, seeds; St, stems. The levels of the detected amplicons were normalized using the amplified product of ACTIN of Sacha Inchi (n = 3) Fig. 7 The ratio of expression of genes encoding enzymes for ALA biosynthesis relative to fatty acid catabolism. FF, female flowers; Fr, fruits; MF, male flowers; ML, mature leaves; Ro, roots; SA, shoot apexes; Se, seeds; St, stems will enable further metabolomic and gene functional study. This ALA-rich species was studied to form a more diversified set of ALA to eventually increase storage ALA production and satisfy more human need worldwide.  [68]. The samples were collected at 60 days after pollination (DAP) when fruits and seeds have reached full size. Three independent biological replicates of each sample were collected from three individual plants. All samples were frozen immediately in liquid nitrogen and then stored at − 80°C for RNA extraction. A total amount of 3 μg of RNA per sample was used as input material for RNA sample preparation. Using poly-T oligo-attached magnetic beads, mRNA was purified from total RNA by following the manufacturer's recommendations (NEB, USA). Then, fragmentation was carried out using divalent cations under elevated temperature. Using these short fragments (about 200 bp) as templates, first-strand cDNA was synthesized using random hexamer primers and MMuLV Reverse Transcriptase (RNase H). Second-strand cDNA synthesis was subsequently performed using DNA polymerase I and RNase H. Then, these cDNA fragments were processed by an end-repair and the ligation of adapters, according to the manufacturer's protocol (Beckman Coulter, Beverly, USA). The products were purified and enriched with PCR for preparing the final sequencing library. Finally, the library quality was assessed using an Agilent Bioanalyzer 2100 system. A TruSeq SR Cluster Kit v3-cBot-HS (Illumina) was used to perform the clustering of index-coded samples. After cluster generation, the library preparations were sequenced on an Illumina Hi-Seq 4000 platform, and 150-bp paired-end reads were generated at the Novogene Bioinformatics Institute, Beijing, P. R. China.

Methods
Sequence data processing and de novo assembly First, raw reads in fastq format were processed through in-house Perl scripts. Then, the clean data were de novo assembled using Trinity software [37], with parameter settings of "-trimmomatic," which removes reads containing adapters, reads containing poly-N and reads of low quality [69].

Functional annotation and pathway assignments
Unigenes were aligned by BLASTx to the NCBI nonredundant (NR), Swiss-Prot, TAIR10, UniRef90 and KOG databases with a threshold E-value of 10 − 5 . For each unigene, the best BLASTx hit from the NR database was submitted to the BLAST2GO platform [42,70], and GO terms were obtained based on annotations between gene names and GO terms. To determine metabolic pathways, the unigenes were also submitted to the online KEGG (Kyoto Encyclopedia of Genes and Genomes) Automatic Annotation Server (KAAS), with single-directional best hit method [71].

Differential expressed genes (DEGs) analysis
The expression level of each unigene was measured using the fragments per kilobase of transcript sequence per millions base pairs sequenced (FPKM) method. The clean reads were aligned to the assembly using Bowtie [72], and the resulting alignments were used to estimate expression abundances in FPKM [37]. All read counts were normalized to FPKM. Differential expression analysis of the samples with three biological replications was performed using the edgeR [53]. Genes with FDR value < 0.001 and |log2(fold change)| ≥2 calculated by edgeR were regarded differentially expressed.

Validation by quantitative real-time PCR (qRT-PCR)
Gene-specific primer pairs were designed using Primer Premier 5.0 software, and the amplified PCR products varied from 100 to 250 bp (Additional file 15). qRT-PCR assays were performed as previously described [74]. The reference gene was chosen based on the methods of Niu et al. [75].

Additional files
Additional file 1: Table S1. Overview of the sequencing data of the Sacha Inchi transcriptomes. (XLS 25 kb) Additional file 2: Figure S1. Overview of Sacha Inchi transcriptome assembly and the characteristics of the homology search of unigenes against the NR database by BLAST (cut-off E-value of 1.0E-5). (a) Size distribution of the assembled unigenes. (b) Similarity distribution of the best BLAST hits for each unigene. (c) E-value distributions of the best BLAST hits for each unigene against the NR database. (d) Species distribution of the best BLAST hit for each unigene. (TIF 935 kb) Additional file 3: Figure S2. The Pearson correlation coefficient (r) was used to estimate the difference between the replicates of each tissue.