- Research article
- Open Access
De novo transcriptome assembly of the eight major organs of Sacha Inchi (Plukenetia volubilis) and the identification of genes involved in α-linolenic acid metabolism
BMC Genomics volume 19, Article number: 380 (2018)
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.
Sacha Inchi (Plukenetia volubilis L.), a member of the Euphorbiaceae family , is native to the rainforests of South America  . Sacha Inchi is also known as Inca peanut, wild peanut, Sacha peanut or mountain peanut, having been cultivated for centuries by the indigenous population . 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 . Based on cytological study, the most common chromosome number is 2n = 58 . 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) , 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 . Leaf extracts were characterized by phenolic compounds, steroids, and/or terpenoids, which resulted in the antioxidant activities in leaves of Sacha Inchi . 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 , 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 . 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 . Moreover, ALA and LA are also critical for the development of infants during pregnancy and breastfeeding periods . 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 . 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 . 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 , and potential applications in biofuel production  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 . 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  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 . 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.  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.  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 .
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.
Functional annotation of the Sacha Inchi transcriptome
After assembly, the 124,750 non-redundant transcripts were subjected to a BLAST search to predict the gene function against five public databases, NR, TAIR10, UniRef90, KOG and Swiss-Prot, and a 10− 5 e-value cut-off value was used . We annotated 67,832 (54.37%) unigenes against the NR database, 68,063 (54.56%) unigenes against the UniRef90 database, 42,109 (33.75%) unigenes against the TAIR10 database, 55,403 (44.41%) unigenes against the KOG database, and 43,501 (34.87%) unigenes against the Swiss-Prot database (Table 1). In total, 70,142 (56.23%) unigenes had at least one homologous match from these databases, whereas 35,381 (28.36%) unigenes had significant BLAST matches to proteins in all of the five databases, as shown in a Venn diagram (Additional file 4: Figure S3). The similarity distribution of the top hits showed that 33.99% of the mapped sequences had similarities higher than 80%, while 62.87% of the hits had similarities ranging from 40 to 80% (Additional file 2: Figure S1b). The E-value distribution had a comparable pattern with 38.47% of the mapped sequences with high homologies (< 1e-50), whereas 61.53% of the homologous sequences ranged between 1e-5 and 1e-50 (Additional file 2: Figure S1c). The species distribution of NR BLAST matches is shown in Additional file 2: Figure S1d. The top-scoring BLASTx hits against the NR protein database showed that the top three species were Ricinus communis (29.27%), Jatropha curcas (9.46%) and Malus domestica (3.20%). We further analyzed the unigenes that showed high similarity with those of the KOG database. Among the 25 categories, the largest group was “General function prediction” (11,069, 19.98%), followed by “Posttranslational modification, protein turnover, chaperones” (5827, 10.52%) and “Signal transduction mechanisms” (5752, 10.38%) (Additional file 5: Figure S4). The categories “Cell motility” (66, 0.12%), “Nuclear structure” (160. 0.29%), “Extracellular structures” (208, 0.38%) and “Coenzyme transport and metabolism” (511, 0.92%) accounted for relatively low proportions (Additional file 5: Figure S4).
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 . 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 , 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 |log2(fold change)| ≥2 found by edgeR were regarded as DEGs . 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 acyltransferase (LPAAT). 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 . 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 .
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).
Identification of simple sequence repeats (SSRs)
To develop molecular markers for genetic analysis and marker-assisted selection breeding of Sacha Inchi, simple sequence repeats (SSRs) were identified in our transcriptome. Here, 42,987 SSRs were detected in all of the 124,750 assembled unigenes using MISA software (Additional file 15). Of the SSRs unigenes, 9902 sequences contained more than one SSR and 4788 SSRs were presented in compound formation (Table 4). Of the 42,987 detected SSRs, monomer nucleotide repeats were the most abundant type (30,088; 69.99%), followed by dimer nucleotides (7345; 17.09%), trimer nucleotides (5061; 11.77%), tetramer nucleotides (356; 0.83%), pentamer nucleotides (46; 0.11%) and hexamer nucleotides (91; 0.21%). Among the 42,987 nucleotide repeats, A/T (69.40%) was the most abundant motifs, followed by AT/AT (9.02%) and AG/CT (5.64%).
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.
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 , 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 . Meanwhile, overexpression of FAD7 increased levels of ALA that led to series of physiological alterations, such as electrolyte leakage and malondialdehyde contents in tomato . 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 . The PDCT from flax are capable of increasing C18-PUFA levels substantially in metabolically engineered yeast and transgenic A. thaliana seeds . 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% . ABI3 and FUS3 are key regulators in phase transition and seed development, maturation and dormancy, and thereby indirectly regulate oil biosynthesis . 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 .
Based on the gene expression pattern analysis, genes coding enzymes related to the ALA biosynthesis strongly express in the seeds, whereas FA catabolism-related genes exhibited weak expression in seeds compared to those in other organs (Figs. 2 and 7),which might related to high ALA content.
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 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.
Plant materials, cDNA library construction and sequencing
Eight organs, including roots (Ro), stems (St), shoot apexes (SA), mature leaves (ML), male flowers (MF), female flowers (FF), fruits (Fr), and seeds (Se), were harvested in August from the 1-year-old plants of Sacha Inchi grown at the Xishuangbanna Tropical Botanical Garden (21o54’_N, 101o46’_E, 580 m above sea level), Chinese Academy of Sciences, Mengla, Yunnan, China under natural climate conditions. August has a mean temperature of 25.1 °C and mean monthly precipitation greater than 200 mm in Xishuangbanna Tropical Botanical Garden . 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.
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 , with parameter settings of “-trimmomatic,” which removes reads containing adapters, reads containing poly-N and reads of low quality .
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 .
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 , and the resulting alignments were used to estimate expression abundances in FPKM . All read counts were normalized to FPKM. Differential expression analysis of the samples with three biological replications was performed using the edgeR . Genes with FDR value < 0.001 and |log2(fold change)| ≥2 calculated by edgeR were regarded differentially expressed.
Simple sequence repeat (SSR) analysis
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 . The reference gene was chosen based on the methods of Niu et al. .
Differentially expressed genes
Fatty acid desaturase 2
Fatty acid desaturase 3
Fatty acid desaturase 6
Fatty acid desaturase 7
Fatty acid desaturase 8
Acyl-ACP thioesterase A
Acyl-ACP thioesterase B
Fragments per kilobase per million mapped base pairs of sequenced
Ketoacyl-ACP synthase II
Ketoacyl-ACP synthase III
Kyoto Encyclopedia of Genes and Genomes
Lysophosphatidic acid acyltransferase
National Center for Biotechnology Information
Phospholipid: diacylglycerol acyltransferase
Phosphatidylcholine diacylglycerol cholinephosphotransferase
Polyunsaturated fatty acids
Gillespie LJ. A revision of paleotropical Plukenetia (Euphorbiaceae) including two new species from Madagascar. Syst Bot. 2007;32(4):780–802.
Hamaker BR, Valles C, Gilman R, Hardmeier RM, Clark D, Garcia HH, Gonzales AE, Kohlstad I, Castro M, Valdivia R. Amino acid and fatty acid profiles of the Inca peanut (Plukenetia volubilis). Cereal Chem. 1992;69(4):461–3.
Gutierrez LF, Rosada LM, Jimenez A. Chemical composition of Sacha Inchi (Plukenetia volubilis L.) seeds and characteristics of their lipid fraction. Grasas Aceites. 2011;62(1):76–83.
Vasek J, Hlasna Cepkova P, Viehmannova I, Ocelak M, Cachique Huansi D, Vejl P. Dealing with AFLP genotyping errors to reveal genetic structure in Plukenetia volubilis (Euphorbiaceae) in the Peruvian Amazon. PLoS One. 2017;12(9):e0184259.
Cai ZQ, Zhang T, Jian HY. Chromosome number variation in a promising oilseed woody crop,Plukenetia volubilisL. (Euphorbiaceae). Caryologia. 2013;66(1):54–8.
Niu LJ, Li JL, Chen MS, Xu ZF. Determination of oil contents in Sacha Inchi (Plukenetia volubilis) seeds at different developmental stages by two methods: Soxhlet extraction and time-domain nuclear magnetic resonance. Ind Crop Prod. 2014;56:187–90.
Chirinos R, Zuloeta G, Pedreschi R, Mignolet E, Larondelle Y, Campos D. Sacha Inchi (Plukenetia volubilis): a seed source of polyunsaturated fatty acids, tocopherols, phytosterols, phenolic compounds and antioxidant capacity. Food Chem. 2013;141(3):1732–9.
Follegatti-Romero LA, Piantino CR, Grimaldi R, Cabral FA. Supercritical CO2 extraction of omega-3 rich oil from Sacha Inchi (Plukenetia volubilis L.) seeds. J Supercrit Fluid. 2009;49(3):323–9.
Hounsome N, Hounsome B, Tomos D, Edwards-Jones G. Plant metabolites and nutritional quality of vegetables. J Food Sci. 2008;73(4):R48–65.
Kornsteiner M, Wagner KH, Elmadfa I. Tocopherols and total phenolics in 10 different nut types. Food Chem. 2006;98(2):381–7.
Fanali C, Dugo L, Cacciola F, Beccaria M, Grasso S, Dachà M, Dugo P, Mondello L. Chemical characterization of Sacha Inchi (Plukenetia volubilis L.) oil. J Agric Food Chem. 2011;59(24):13043–9.
Chirinos R, Necochea O, Pedreschi R, Campos D. Sacha Inchi (Plukenetia volubilis L.) shell: an alternative source of phenolic compounds and antioxidants. Int J Food Sci Tech. 2016;51(4):986–93.
Nascimento AK, Melo-Silveira RF, Dantas-Santos N, Fernandes JM, Zucolotto SM, Rocha HA, Scortecci KC. Antioxidant and Antiproliferative activities of leaf extracts from Plukenetia volubilis Linneo (Euphorbiaceae). Evid Based Compl Alt. 2013;2013(6):1–10.
Burdge GC. Metabolism of alpha-linolenic acid in humans. Prostag Leukotr Ess. 2006;75(3):161–8.
Simopoulos AP. Human requirement for n-3 polyunsaturated fatty acids. Poultry Sci. 2000;79(7):961–70.
Kremer JM. Effects of modulation of inflammatory and immune parameters in patients with rheumatic and inflammatory disease receiving dietary supplementation of n-3 and n-6 fatty acids. Lipids. 1996;31(Suppl):S243–7.
Belluzzi A, Brignola C, Campieri M, Pera A, Boschi S, Miglioli M. Effect of an enteric-coated fish-oil preparation on relapses in Crohn's disease. New Engl J Med. 1996;334(24):1557–60.
Ramaprasad TR, Srinivasan K, Baskaran V, Sambaiah K, Lokesh BR. Spray-dried milk supplemented with α-linolenic acid or eicosapentaenoic acid and docosahexaenoic acid decreases HMG CoA reductase activity and increases biliary secretion of lipids in rats. Steroids. 2006;71(5):409–15.
Simopoulos AP. Evolutionary aspects of diet: the Omega-6/Omega-3 ratio and the brain. Mol Neurobiol. 2011;44(2):203–15.
Gorriti A, Arroyo J, Quispe F, Cisneros B, Condorhuamán M, Almora Y, Chumpitaz V. Oral toxicity at 60-days of Sacha Inchi oil (Plukenetia volubilis L.) and linseed (Linum usitatissimum L.), and determination of lethal dose 50 in rodents. Rev peru med exp salud publica. 2010;27(3):352. Article in Spanish
Garmendia F, Pando R, Ronceros G. Effect of Sacha Inchi oil (Plukenetia volubilis L.) on the lipid profile of patients with hyperlipoproteinemia. Rev peru med exp salud publica. 2011;28(4):628–32. Article in Spanish
Wijendran V, Hayes KC. Dietary n-6 and n-3 fatty acid balance and cardiovascular health. Annu Rev Nutr. 2004;24(24):597–615.
Cleeman JI. Report of the national cholesterol education-program expert panel on detection, evaluation, and treatment of high blood cholesterol in adults. Arch Intern Med. 1988;148(1):36–69.
Hajra AK, Larkins LK, Das AK, Hemati N, Erickson RL, MacDougald OA. Induction of the peroxisomal glycerolipid-synthesizing enzymes during differentiation of 3T3-L1 adipocytes - Role in triacylglycerol synthesis. J Biol Chem. 2000;275(13):9441–6.
Storch J, Zhou YX, Lagakos WS. Metabolism of apical versus basolateral sn-2-monoacylglycerol and fatty acids in rodent small intestine. J Lipid Res. 2008;49(8):1762–9.
Manan S, Chen BB, She GB, Wan XC, Zhao J. Transport and transcriptional regulation of oil production in plants. Crit Rev Biotechnol. 2017;37(5):641–55.
Cai ZQ. Shade delayed flowering and decreased photosynthesis, growth and yield of Sacha Inchi (Plukenetia volubilis) plants. Ind Crop Prod. 2011;34(1):1235–7.
Luo YL, Su ZL, Bi TJ, Cui XL, Lan QY. Salicylic acid improves chilling tolerance by affecting antioxidant enzymes and osmoregulators in sacha inchi (Plukenetia volubilis). Braz J Bot. 2014;37(3):357–63.
Fu QT, Niu LJ, Zhang QF, Pan BZ, He HY, Xu ZF. Benzyladenine treatment promotes floral feminization and fruiting in a promising oilseed crop Plukenetia volubilis. Ind Crop Prod. 2014;59(59):295–8.
Dong Y, Chen M, Wang X, Niu L, Fu Q, Xu Z. Establishment of in vitro regeneration system of woody oil crop Plukenetia volubilis. Molecular Plant Breeding. 2016;14(2):462–70. Article in Chinese
Zuleta EC, Rios LA, Benjumea PN. Oxidative stability and cold flow behavior of palm, sacha-inchi, jatropha and castor oil biodiesel blends. Fuel Process Technol. 2012;102:96–101.
Hanssen H-P, Schmitz-Huebsch M. Sacha Inchi (Plukenetia volubilis L.) nut oil and its therapeutic and nutritional uses. In: Preedy VR, Watson RR, Patel VB, editors. Nuts and Seeds in Health and Disease Prevention. San Diego, USA: Elsevier academic press Inc; 2011. p. 991–4.
Gonzalez-Aspajo G, Belkhelfa H, Haddioui-Hbabi L, Bourdy G, Deharo E. Sacha Inchi oil (Plukenetia volubilis L.), effect on adherence of Staphylococus aureus to human skin explant and keratinocytes in vitro. J Ethnopharmacol. 2015;171:330–4.
Gonzales GF, Gonzales C. A randomized, double-blind placebo-controlled study on acceptability, safety and efficacy of oral administration of sacha inchi oil (Plukenetia volubilis L.) in adult human subjects. Food Chem Toxicol. 2014;65(0):168–76.
Wang XJ, Xu RH, Wang RL, Liu AZ. Transcriptome analysis of Sacha Inchi (Plukenetia volubilis L.) seeds at two developmental stages. BMC Genomics. 2012;13(1):716.
Wang X, Liu A. Expression of genes controlling unsaturated fatty acids biosynthesis and oil deposition in developing seeds of Sacha Inchi (Plukenetia volubilis L.). Lipids. 2014;49(10):1019–31.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng QD, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–U130.
Zhao QY, Wang Y, Kong YM, Luo D, Li X, Hao P. Optimizing de novo transcriptome assembly from short-read RNA-Seq data: a comparative study. BMC bioinformatics. 2011;12(Suppl 14(14)):S2.
Ma L, Sun N, Liu X, Jiao Y, Zhao H, Deng XW. Organ-specific expression of Arabidopsis genome during development. Plant Physiol. 2005;138(1):80–91.
Altschul SF, Madden TL, Schaffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.
Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talon M, Dopazo J, Conesa A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008;36(10):3420–35.
Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.
Bates PD. Understanding the control of acyl flux through the lipid metabolic network of plant oil biosynthesis. BBA-Mol Cell Biol L. 2016;1861(9):1214–25.
Li-Beisson Y, Shorrosh B, Beisson F, Andersson MX, Arondel V, Bates PD, Baud S, Bird D, DeBono A, Durrett TP. Acyl-lipid metabolism. The Arabidopsis Book. 2013;11:e0161.
Ullman AH, Iii HBW.  assay of fatty acid synthase using a bicyclic dione as substrate. Methods Enzymol. 1981;72:303–6.
Kunst L, Samuels AL. Biosynthesis and secretion of plant cuticular wax. Prog Lipid Res. 2003;42(1):51–80.
Zhou XR, Green AG, Singh SP. Caenorhabditis elegans delta 12-desaturase fat-2 is a bifunctional desaturase able to desaturate a diverse range of fatty acid substrates at the delta 12 and delta 15 positions. J Biol Chem. 2011;286(51):43644–50.
Watts JL, Browse J. Genetic dissection of polyunsaturated fatty acid synthesis in Caenorhabditis elegans. P Natl Acad Sci USA. 2002;99(9):5854–9.
Meesapyodsuk D, Reed DW, Savile CK, Buist PH, Schafer UA, Ambrose SJ, Covello PS. Substrate specificity, regioselectivity and cryptoregiochemistry of plant and animal omega-3 fatty acid desaturases. Biochem Soc T. 2000;28(6):632–5.
Wickramarathna AD, Siloto RMP, Mietkiewska E, Singer SD, Pan X, Weselake RJ. Heterologous expression of flax PHOSPHOLIPID:DIACYLGLYCEROL CHOLINEPHOSPHOTRANSFERASE (PDCT) increases polyunsaturated fatty acid content in yeast and Arabidopsis seeds. BMC Biotechnol. 2015;15(1):63.
Hu ZH, Ren ZH, Lu CF. The phosphatidylcholine diacylglycerol cholinephosphotransferase is required for efficient hydroxy fatty acid accumulation in transgenic Arabidopsis. Plant Physiol. 2012;158(4):1944–54.
Chen G, Woodfield HK, Pan X, Harwood JL, Weselake RJ. Acyl-trafficking during plant oil accumulation. Lipids. 2015;50(11):1057–68.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Kornberg A, Jr PW. enzymatic synthesis of the coenzyme a derivatives of long chain fatty acids. J Biol Chem. 1953;204(1):329–43.
Kent C. Eukaryotic phospholipid biosynthesis. Annu Rev Biochem. 1995;64(64):315–43.
Stahl U, Carlsson AS, Lenman M, Dahlqvist A, Huang BQ, Banas W, Banas A, Stymne S. Cloning and functional characterization of a phospholipid : diacylglycerol acyltransferase from Arabidopsis. Plant Physiol. 2004;135(3):1324–35.
Dahlqvist A, Stahl U, Lenman M, Banas A, Lee M, Sandager L, Ronne H, Stymne H. Phospholipid : diacylglycerol acyltransferase: an enzyme that catalyzes the acyl-CoA-independent formation of triacylglycerol in yeast and plants. P Natl Acad Sci USA. 2000;97(12):6487–92.
Bates PD, Browse J. The pathway of triacylglycerol synthesis through phosphatidylcholine in Arabidopsis produces a bottleneck for the accumulation of unusual fatty acids in transgenic seeds. Plant J. 2011;68(3):387–99.
Chiu RS, Nahal H, Provart NJ, Gazzarrini S. The role of the Arabidopsis FUSCA3 transcription factor during inhibition of seed germination at high temperature. BMC Plant Biol. 2012;12(1):15.
Bahnson BJ, Anderson VE, Petsko GA. Structural mechanism of enoyl-CoA hydratase: three atoms from a single water are added in either an E1cb stepwise or concerted fashion. Biochemistry-Us. 2002;41(8):2621–9.
Carrozzo R, Bellini C, Lucio S, Deodato F, Cassandrii D, Cassanello M, Caruso U, Rizzo C, Rizza T, Napolitano ML, et al. Peroxisomal acyl-CoA-oxidase deficiency: two new cases. Am J Med Genet A. 2008;146A(13):1676–81.
Browse J, Mcconn M, James D, Miquel M. Mutants of Arabidopsis deficient in the synthesis of alpha-linolenate. Biochemical and genetic characterization of the endoplasmic reticulum linoleoyl desaturase. J Biol Chem. 1993;268(22):16345–51.
Lemieux B, Miquel M, Somerville C, Browse J. Mutants of Arabidopsis with alterations in seed lipid fatty acid composition. Theor Appl Genet. 1990;80(2):234–40.
Shah S, Xin ZG, Browse J. Overexpression of the FAD3 desaturase gene in a mutant of Arabidopsis. Plant Physiol. 1997;114(4):1533–9.
Liu XY, Teng YB, Li B, Meng QW. Enhancement of low-temperature tolerance in transgenic tomato plants overexpressing Lefad7 through regulation of trienoic fatty acids. Photosynthetica. 2013;51(2):238–44.
Shen B, Allen WB, Zheng PZ, Li CJ, Glassman K, Ranch J, Nubel D, Tarczynski MC. Expression of ZmLEC1 and ZmWRI1 increases seed oil production in maize. Plant Physiol. 2010;153(3):980–7.
Wang J, Jian H, Wang T, Wei L, Li J, Li C, Liu L. Identification of microRNAs actively involved in fatty acid biosynthesis in developing Brassica napus seeds using high-throughput sequencing. Front Plant Sci. 2016;7
Cao M, Zou X, Warren M, Zhu H. Tropical forests of Xishuangbanna, China. Biotropica. 2010;38(3):306–9.
Hansen KD, Brenner SE, Dudoit S. Biases in Illumina transcriptome sequencing caused by random hexamer priming. Nucleic Acids Res. 2010;38(12):e131.
Conesa A, Gotz S: Blast2GO: a comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics 2008, 2008:619832.
Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007;35:W182–5.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9.
Thiel T, Michalek W, Varshney R, Graner A. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor Appl Genet. 2003;106(3):411–22.
Pan BZ, Chen MS, Ni J, Xu ZF. Transcriptome of the inflorescence meristems of the biofuel plant Jatropha curcas treated with cytokinin. BMC Genomics. 2014;15(1):974.
Niu LJ, Tao YB, Chen MS, Fu QT, Li CQ, Dong YL, Wang XL, He HY, Xu ZF. Selection of reliable reference genes for gene expression studies of a promising oilseed crop, Plukenetia volubilis, by real-time quantitative PCR. Int J Mol Sci. 2015;16(6):12513–30.
The authors gratefully acknowledge the Central Laboratory of the Xishuangbanna Tropical Botanical Garden for providing the research facilities.
X-DH and B-ZP designed the experiments, performed the experiments, analyzed the results, and wrote the manuscript. Z-FX conceived and designed the study, analyzed the results, and wrote the manuscript. QF and LN performed some of the experiments and revised the manuscript. M-SC revised the manuscript. All authors reviewed and approved the manuscript.
This work was supported by the Natural Science Foundation of Yunnan Province (2014FB186 and 2016FB051) and the CAS 135 program (2017XTBG-T02). Funding sources played no role in the design of this study or the collection, analysis, and the interpretation of data and in writing the manuscript.
Availability of data and materials
The data sets of sequencing reads have been deposited in the National Center for Biotechnology Information database (NCBI) and can be retrieved under the SRA accession number SRP101395.
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Overview of the sequencing data of the Sacha Inchi transcriptomes. (XLS 25 kb)
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)
Figure S2. The Pearson correlation coefficient (r) was used to estimate the difference between the replicates of each tissue. The number between these two samples is given in the plot. The color represents r value, which shows high correlation in red between two samples, while low correlation in blue. (TIF 4328 kb)
Figure S3. Venn diagram showing the BLAST searches of the Sacha Inchi transcriptome against the five public databases. De novo unigene sequences were used to search against the following public databases: NR, UniRef90, TAIR10, KOG and Swiss-Prot. The numbers of unigenes that have significant hits against the five databases are shown in each intersection in the Venn diagram. (TIF 574 kb)
Figure S4. Histogram presentation of clusters of orthologous group classification of assembled unigenes. A total of 124,750 unigenes were classified into 25 functional categories. (TIF 816 kb)
Figure S5. Distribution of gene ontology (GO) categories of unigenes for Sacha Inchi. GO functional annotations are summarized into three main categories: biological process, cellular component, and molecular function. The number of unigenes in each category is shown on the y-axis. (TIF 1684 kb)
Pathways identified in the Sacha Inchi transcriptome. Three hundred and fifty-three KEGG pathways identified in Sacha Inchi and the corresponding unigene numbers of each pathway are shown. (XLS 60 kb)
Figure S6. The statistics of unigenes in eight organs. a) The percentage of unigenes expressed in each organ. b) Number of organ-specific unigenes. FF, female flowers; Fr, fruits; MF, male flowers; ML, mature leaves; Ro, roots; SA, shoot apexes; Se, seeds; St, stems. (TIF 2536 kb)
Detailed FPKM and Annotation information of organ-specific unigenes. (XLS 646 kb)
The list of GO terms that were significantly enriched in male flowers, female flowers, roots, seeds and stems. Gene ontology (GO) terms were assigned to organ-specific unigenes based on the top hits against the NR database. (XLS 33 kb)
Figure S7. Histogram presentation of the KEGG pathway annotation of female flower (FF)-specific genes. (TIF 1210 kb)
The list of KEGG pathways were enriched in each organ. (XLS 70 kb)
List of FA and TAG metabolism-related genes detected in Sacha Inchi transcriptome. A) Unigenes related to ALA biosynthesis. B) Unigenes related to the fatty acid catabolism pathway. C) Unigenes related to TAG biosynthesis. (XLS 69 kb)
Figure S8. Heat map representation and hierarchical clustering of putative genes involved in glycolysis/gluconeogenesis pathway and pentose phosphate pathway. A: glycolysis/gluconeogenesis pathway (ko00010). B: pentose phosphate pathway (ko00030). (TIF 2539 kb)
Overview of the SSRs detected in the assembled unigenes of Sacha Inchi. (XLS 4872 kb)
Primers for qRT-PCR in this study. (XLS 31 kb)
Figure S9. Pearson correlation analysis of the gene expression ratios obtained from RNA-Seq and qPCR data. The qPCR log10 values (expression ratios; y-axis) were plotted against the RNA-Seq log10 values (x-axis). The Pearson correlation coefficient (r) is given in the plot, and the circle indicates the extremely significant difference at p < 0.01. (TIF 337 kb)