Transcriptome analysis of Sacha Inchi (Plukenetia volubilis L.) seeds at two developmental stages

Sacha Inchi (Plukenetia volubilis L.) Protein Hydrolysate as a New Ingredient of Functional Foods

Sacha inchi (Plukenetia volubilis L.) is an under-exploited crop with great potential due to its nutritional and medicinal characteristics. A Sacha inchi protein isolate (SII), obtained from defatted Sacha inchi flour (SIF), was hydrolyzed by Bioprotease LA 660 under specific conditions. The hydrolysates were characterized chemically, and their digestibility and antioxidant capacity were evaluated by in vitro cell-free experiments to select the hydrolysate with major antioxidant activity. Sacha inchi protein hydrolysate at 20 min (SIH20B) was selected, and the anti-inflammatory capacity was evaluated by RT-qPCR and ELISA techniques, using two different doses in monocytes THP-1 stimulated with lipopolysaccharide (LPS). The results obtained showed that the in vitro administration of SIH20B down-regulated the TNF-α gene and reduced the release of this cytokine, whereas the anti-inflammatory cytokines IL-10 and IL-4 were up-regulated in LPS-stimulated monocytes and co-administrated with SIH20B. The peptides contained in SIH20B were identified, and the 20 more relatively abundant peptides with a mass by 1 kDa were subjected to in silico analysis to hypothesize those that could be responsible for the bioactivity reported in the hydrolysate. From the identified peptides, the peptides AAGALKKFL and LGVKFKGGL, among others, are proposed as the most biologically actives. In conclusion, SIH20B is a novel, natural source of high-value-added biopeptides that could be used as an ingredient in formulations of food or nutraceutical compounds.

The transcription factor MYB1 activates DGAT2 transcription to promote triacylglycerol accumulation in sacha inchi (Plukenetia volubilis L.) leaves under heat stress

Triacylglycerol (TAG) accumulation is frequently triggered in vegetative tissues experiencing heat stress, which may increases plant basal plant thermo-tolerance by sequestering the toxic lipid intermediates that contribute to membrane damage or cell death under stress conditions. However, stress-responsive TAG biosynthesis and the underlying regulatory mechanisms are not fully understood. Here, we investigated the lipidomic and transcriptomic landscape under heat stress in the leaves of sacha inchi (Plukenetia volubilis L.), an important oilseed crop in tropical regions. Under heat stress (45 °C), the content of polyunsaturated TAGs (e.g., TAG18:2 and TAG18:3) and total TAGs were significantly higher, while those of unsaturated sterol esters, including ZyE 28:4, SiE 18:2 and SiE 18:3, were dramatically lower. Transcriptome analysis showed that the expression of PvDGAT2-2, encoding a type II diacylglycerol acyltransferase (DGAT) that is critical for TAG biosynthesis, was substantially induced under heat stress. We confirmed the function of PvDGAT2-2 in TAG production by complementing a yeast mutant defective in TAG biosynthesis. Importantly, we also identified the heat-induced transcription factor PvMYB1 as an upstream activator of PvDGAT2-2 transcription. Our findings on the molecular mechanism leading to TAG biosynthesis in leaves exposed to heat stress have implications for improving the biotechnological production of TAGs in vegetative tissues, offering an alternative to seeds.

Full-length transcriptome revealed the accumulation of polyunsaturated fatty acids in developing seeds of Plukenetia volubilis

Background

Plukenetia volubilis is cultivated as a valuable oilseed crop, and its mature seeds are rich in polyunsaturated fatty acids (FAs), which are widely used in food and pharmaceutical industries. Recently, next-generation sequencing (NGS) transcriptome studies in P. volubilis indicated that some candidate genes were involved in oil biosynthesis. The NGS were inaccuracies in assembly of some candidate genes, leading to unknown errors in date analyses. However, single molecular real-time (SMRT) sequencing can overcome these assembled errors. Unfortunately, this technique has not been reported in P. volubilis.

Methods

The total oil content of P. volubilis seed (PVS) was determined using Soxhlet extraction system. The FA composition were analyzed by gas chromatography. Combining PacBio SMRT and Illumina technologies, the transcriptome analysis of developing PVS was performed. Functional annotation and differential expression were performed by BLAST software (version 2.2.26) and RSEM software (version 1.2.31), respectively. The lncRNA-targeted transcripts were predicted in developing PVS using LncTar tool.

Results

By Soxhlet extraction system, the oil content of superior plant-type (SPT) was 13.47% higher than that of inferior plant-type (IPT) at mature PVS. The most abundant FAs were C18:2 and C18:3, among which C18:3 content of SPT was 1.11-fold higher than that of IPT. Combined with PacBio and Illumina platform, 68,971 non-redundant genes were obtained, among which 7,823 long non-coding RNAs (lncRNAs) and 7,798 lncRNA-targeted genes were predicted. In developing seed, the expressions of 57 TFs showed a significantly positive correlation with oil contents, including WRI1-like1, LEC1-like1, and MYB44-like. Comparative analysis of expression profiles between SPT and IPT implied that orthologs of FAD3, PDCT, PDAT, and DAGT2 were possibly important for the accumulation of polyunsaturated FAs. Together, these results provide a reference for oil biosynthesis of P. volubilis and genetic improvement of oil plants.

Integrated Transcriptome and Proteome Analysis Provides Insight into the Ribosome Inactivating Proteins in Plukenetia volubilis Seeds

Plukenetia volubilis is a highly promising plant with high nutritional and economic values. In our previous studies, the expression levels of ricin encoded transcripts were the highest in the maturation stage of P. volubilis seeds. The present study investigated the transcriptome and proteome profiles of seeds at two developmental stages (Pv-1 and Pv-2) using RNA-Seq and iTRAQ technologies. A total of 53,224 unigenes and 6026 proteins were identified, with functional enrichment analyses, including GO, KEGG, and KOG annotations. At two development stages of P. volubilis seeds, 8815 unique differentially expressed genes (DEGs) and 4983 unique differentially abundant proteins (DAPs) were identified. Omics-based association analysis showed that ribosome-inactivating protein (RIP) transcripts had the highest expression and abundance levels in Pv-2, and those DEGs/DAPs of RIPs in the GO category were involved in hydrolase activity. Furthermore, 21 RIP genes and their corresponding amino acid sequences were obtained from libraries produced with transcriptome analysis. The analysis of physicochemical properties showed that 21 RIPs of P. volubilis contained ricin, the ricin_B_lectin domain, or RIP domains and could be divided into three subfamilies, with the largest number for type II RIPs. The expression patterns of 10 RIP genes indicated that they were mostly highly expressed in Pv-2 and 4 transcripts encoding ricin_B_like lectins had very low expression levels during the seed development of P. volubilis. This finding would represent valuable evidence for the safety of oil production from P. volubilis for human consumption. It is also notable that the expression level of the Unigene0030485 encoding type I RIP was the highest in roots, which would be related to the antiviral activity of RIPs. This study provides a comprehensive analysis of the physicochemical properties and expression patterns of RIPs in different organs of P. volubilis and lays a theoretical foundation for further research and utilization of RIPs in P. volubilis.

Lipid composition of the Amazonian 'Mountain Sacha Inchis' including Plukenetia carolis-vegae Bussmann, Paniagua & C.Téllez

Several Amazonian species of Plukenetia are remarkably rich sources of polyunsaturated fatty acids, in particular α-linolenic acid. The lipid composition of the large-seeded, recently described ‘Mountain Sacha Inchi’ Plukenetia carolis-vegae is reported here for the first time, and compared with Plukenetia huayllabambana, two cultivars of Plukenetia volubilis, and a newly developed hybrid between P. volubilis and P. carolis-vegae. All species and cultivars had a very high content of polyunsaturated fatty acids, 82.6–86.7% of total fatty acids, and at least 46.6% α-linolenic acid of total fatty acids. The highest content was found in P. carolis-vegae which had 57.4%. The exceptionally high α-linolenic acid content suggests that P. carolis-vegae may be an important plant-derived dietary source of this essential fatty acid and that the species has considerable potential for further domestication and commercialisation of its seeds and seed oil. A TAG analysis was carried out for the two P. volubilis cultivars, in which LLnLn and LnLL were most prevalent, and for P. huayllabambana, in which LLnLn constituted the largest fraction, followed by LnLnLn, indicating that this large-seeded species also has interesting dietary properties.

Overexpression of PvFAD3 Gene from Plukenetia volubilis Promotes the Biosynthesis of α-Linolenic Acid in Transgenic Tobacco Seeds

The ω-3 fatty acid desaturase (FAD3) gene encodes a rate-limiting enzyme in the synthesis of α-linolenic acid. In this study, homologous cloning was used to obtain the full-length sequence of the PvFAD3 gene of Plukenetia volubilis. The full-length DNA sequence was 1871 bp long, with 8 exons and 7 introns. The structural analysis of the amino acid sequence revealed that the PvFAD3 protein contained three histidine-conserved regions and an endoplasmic reticulum retention signal. The real-time reverse transcription-polymerase chain reaction performed for determining the expression patterns of the PvFAD3 gene in different tissues of P. volubilis showed that PvFAD3 expression was highly expressed in the fast oil accumulation stage of seed. The analysis of subcellular localization assay in epidermal cells of tobacco (Nicotiana benthamiana) leaves showed that the PvFAD3 protein was mainly localized in the endoplasmic reticulum. Seed-specific overexpression vectors were constructed, and Agrobacterium-mediated genetic transformation was performed to obtain transgenic tobacco plants overexpressing PvFAD3. The results of fatty acid assays performed using harvested seeds showed a significant increase in α-linolenic acid content, a dramatic decrease in linoleic acid content, and an obvious increase in oil content in transgenic tobacco seeds. Collectively, the PvFAD3 gene of P. volubilis was confirmed as a key enzyme gene for α-linolenic acid synthesis; thus, indicating that the PvFAD3 gene can be used for fatty acid fraction improvement in oilseed plants.

Genome-Wide Analysis of Alternative Splicing and Non-Coding RNAs Reveal Complicated Transcriptional Regulation in Cannabis sativa L

It is of significance to mine the structural genes related to the biosynthetic pathway of fatty acid (FA) and cellulose as well as explore the regulatory mechanism of alternative splicing (AS), microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) in the biosynthesis of cannabinoids, FA and cellulose, which would enhance the knowledge of gene expression and regulation at post-transcriptional level in Cannabis sativa L. In this study, transcriptome, small RNA and degradome libraries of hemp 'Yunma No.1' were established, and comprehensive analysis was performed. As a result, a total of 154, 32 and 331 transcripts encoding key enzymes involved in the biosynthesis of cannabinoids, FA and cellulose were predicted, respectively, among which AS occurred in 368 transcripts. Moreover, 183 conserved miRNAs, 380 C. sativa-specific miRNAs and 7783 lncRNAs were predicted. Among them, 70 miRNAs and 17 lncRNAs potentially targeted 13 and 17 transcripts, respectively, encoding key enzymes or transporters involved in the biosynthesis of cannabinoids, cellulose or FA. Finally, the crosstalk between AS and miRNAs or lncRNAs involved in cannabinoids and cellulose was also predicted. In summary, all these results provided insights into the complicated network of gene expression and regulation in C. sativa.

Sacha inchi (Plukenetia volubilis L.): An emerging source of nutrients, omega-3 fatty acid and phytochemicals

Sacha inchi (Plukenetia volubilis) (SI) is an oleaginous plant producing oil and protein-rich seeds. It has been cultivated for centuries and is native to the tropical rainforest of the Amazon region of South America including parts of Peru and northwestern Brazil. At present, SI seeds are emerging as a potential source of macro- and micronutrients, α-linolenic acid and phytochemicals. This review attempts to elucidate the nutrients, phytonutrients, safety, toxicity, health benefits and food applications of SI seed. Recent scientific studies have associated the consumption of SI seed/oil with reduced risk of chronic inflammatory diseases. However, lack of awareness and in-depth understanding has resulted in it being neglected both at the consumer and industrial level. In all, SI is an underutilized and undervalued oleaginous crop which not only has the potential to mitigate food and nutritional insecurity but also offers humongous opportunities for the development of novel value-added food products.

Transcriptome analyses reveals the dynamic nature of oil accumulation during seed development of Plukenetia volubilis L

Sacha inchi (Plukenetia volubilis L.) is a shrub native to Amazon rainforests that’s of commercial interest as its seeds contain 35–60% edible oil (dry weight). This oil is one of the healthiest vegetable oils due to its high polyunsaturated fatty acid content and favourable ratio of omega-6 to omega-3 fatty acids. De novo transcriptome assembly and comparative analyses were performed on sacha inchi seeds from five stages of seed development in order to identifying genes associated with oil accumulation and fatty acid production. Of 30,189 unigenes that could be annotated in public databases, 20,446 were differentially expressed unigenes. A total of 14 KEGG pathways related to lipid metabolism were found, and 86 unigenes encoding enzymes involved in α-linolenic acid (ALA) biosynthesis were obtained including five unigenes encoding FATA (Unigene0008403), SAD (Unigene0012943), DHLAT (Unigene0014324), α-CT (Unigene0022151) and KAS II (Unigene0024371) that were significantly up-regulated in the final stage of seed development. A total of 66 unigenes encoding key enzymes involved in the synthesis of triacylglycerols (TAGs) were found, along with seven unigenes encoding PDCT (Unigene0000909), LPCAT (Unigene0007846), Oleosin3 (Unigene0010027), PDAT1 (Unigene0016056), GPDH (Unigene0022660), FAD2 (Unigene0037808) and FAD3 (Unigene0044238); these also proved to be up-regulated in the final stage of seed development.

Isolation and functional analyses of PvFAD2 and PvFAD3 involved in the biosynthesis of polyunsaturated fatty acids from Sacha Inchi (Plukenetia volubilis)

The development of ω-3 fatty acid-rich vegetable oils is essential to enrich the production of functional foods. Sacha Inchi (Plukenetia volubilis L.) is a unique oilseed crop with much potential. Its seeds contain rich polyunsaturated fatty acids (PUFAs), especially linoleic acid (LA, C18:2) and α-linolenic acid (ALA, C18:3). Endoplasmic reticulum -located ω-6 and ω-3 fatty acid desaturases (FAD) are responsible for the biosynthesis of LA and ALA, respectively, in plant seeds. Here, we isolated two full-length FAD genes from Sacha Inchi, named PvFAD2 and PvFAD3, which encoded predicted amino acid residues of 384 and 379 in protein, respectively. Protein sequence and subcellular localization analysis revealed that they were located in the endoplasmic reticulum (ER). Heterologous expression in Saccharomyces cerevisiae confirmed that PvFAD2 and PvFAD3 could catalyze LA and ALA synthesis, respectively. The stability and catalytic efficiency of the PvFAD3 protein may be closely related to temperature. In transgenic tobacco, using seed-specific expression promoters, PvFAD2 and PvFAD3 significantly promotes the production of LA (from 68% to 70.5%) and ALA (from 0.7% to 3.1%) in seed oil. These results show that PvFAD2 and PvFAD3 do, indeed, function as crucial enzymes for PUFAs biosynthesis, and provide a key gene source for the sustainable production of lipids with tailored fatty acid compositions via genetic engineering in other oil crops.

Sacha inchi (Plukenetia volubilis L.)-from lost crop of the Incas to part of the solution to global challenges?

The underutilized, oleaginous crop Plukenetia volubilis L. has a remarkable lipid composition and a large potential for further domestication, alleviation of malnutrition, and integration into sustainable food production systems. Current global challenges include climate change, increasing population size, lack of food security, malnutrition, and degradation of arable lands. In this context, a reformation of our food production systems is imperative. Underutilized crops, or orphan crops, can provide valuable traits for this purpose, e.g., climate change resilience, nutritional benefits, cultivability on marginal lands, and improvement of income opportunities for smallholders. Plukenetia volubilis L. (Euphorbiaceae)-sacha inchi-is a 'lost crop' of the Incas native to the Amazon basin. Its oleaginous seeds are large, with a high content of ω-3, and -6 fatty acids (ca. 50.5, and 34.1%, of the lipid fraction, respectively), protein, and antioxidants. Culinarily, the seeds are nut-like and the crop has been associated with humans since Incan times. Research has particularly been undertaken in seed biochemistry, and to some extent in phylogeny, genetics, and cultivation ecology, but attention has been unevenly distributed, causing knowledge gaps in areas such as ethnobotany, allergenicity, and sustainable cultivation practices. Recently, seed size evolution and molecular drivers of the fatty acid synthesis and composition have been studied, however, further research into the lipid biosynthesis is desirable. Targeted breeding has not been undertaken but might be especially relevant for yield, sensory qualities, and cultivation with low environmental impact. Similarly, studies of integration into sustainable management systems are of highest importance. Here, present knowledge on P. volubilis is reviewed and a general framework for conducting research on underutilized crops with the aim of integration into sustainable food production systems is presented.

Genome-wide analysis of ethylene-response factor family in adzuki bean and functional determination of VaERF3 under saline-alkaline stress

Ethylene-response factor (ERF) proteins are members of a transcription factor family involved in plant growth and environmental stress responses, but the biological functions of ERF members in adzuki bean (Vigna angularis var. angularis) remain unknown. In addition, it is unclear whether these proteins have a role in regulating responses to abiotic stressors. Here, we identified 47 ERF genes by analyzing the adzuki bean genome. Whole-transcriptome analyses of plants under saline-alkaline stress suggested that the expression of 13 ERF genes was induced in response to saline-alkaline stress. Analysis of the cis-acting elements showed that the promoters of these saline-alkaline stress-inducible ERF genes contained LTRs, DREs, MYBs, ABREs, MYCs, CGTCA-, and TGACG-motifs, which are involved in abiotic stress responses. The expression of VaERF3 was induced by NaHCO₃, polyethylene glycol 6000, NaCl, and ABA (abscisic acid), as determined by qRT-PCR. Overexpression of VaERF3 in transgenic Arabidopsis resulted in higher levels of proline accumulation and lower malondialdehyde and reactive oxygen species contents in plants grown under saline-alkaline stress conditions. Moreover, VaERF3 encoded a nuclear-localized transcriptional activator that promoted the expression of stress-responsive genes. Collectively, these results are of great significance in elucidating the mechanisms of saline-alkaline stress responses in adzuki bean.

Melatonin-Induced Transcriptome Variation of Rapeseed Seedlings under Salt Stress

Salt stress inhibits the production of all crop species, including rapeseed (Brassica napus L.), the second most widely planted oil crop species. Although melatonin was confirmed to alleviate salt stress in rapeseed seedlings recently, the mechanism governing the expression levels remains unknown. Therefore, the melatonin-induced transcriptome variation of salt-stressed seedlings was explored. In this study, the transcriptomes of leaves and roots under control (CK), salt (125 mM NaCl, ST) and melatonin (125 mM NaCl plus 50 µM melatonin, MS) treatments were evaluated by using next-generation sequencing techniques. After conducting comparisons of gene expression in the roots and leaves between MS and ST, the differentially expressed gene (DEG) pools were screened. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses highlighted the significant pathways, which were mainly related to plant hormone synthesis and signal transduction, lignin and fatty acid metabolism. The functional genes in the objective KEGG pathways were identified. Furthermore, members of several transcription factor (TF) families participated in the response process. Combined with the hormone (campesterol (CS), jasmonic acid (JA), and gibberellic acid 3 (GA3)) contents measured in the seedlings, it could be concluded that melatonin induced changes in the intrinsic hormone metabolic network, which promoted seedling growth. Thus, this study identified new candidate genes and pathways active during the interactions between melatonin and salt stress, which provide clues for disclosing melatonin’s function in resistance to salt injury. Our results contribute to developing a practical method for sustainable agriculture on saline lands.

Seed size evolution and biogeography of Plukenetia (Euphorbiaceae), a pantropical genus with traditionally cultivated oilseed species

BACKGROUND

Plukenetia is a small pantropical genus of lianas and vines with variably sized edible oil-rich seeds that presents an ideal system to investigate neotropical and pantropical diversification patterns and seed size evolution. We assessed the biogeography and seed evolution of Plukenetia through phylogenetic analyses of a 5069 character molecular dataset comprising five nuclear and two plastid markers for 86 terminals in subtribe Plukenetiinae (representing 20 of ~ 23 Plukenetia species). Two nuclear genes, KEA1 and TEB, were used for phylogenetic reconstruction for the first time. Our goals were: (1) produce a robust, time-dependent evolutionary framework for Plukenetia using BEAST; (2) reconstruct its biogeographical history with ancestral range estimation in BIOGEOBEARS; (3) define seed size categories; (4) identify patterns of seed size evolution using ancestral state estimation; and (5) conduct regression analyses with putative drivers of seed size using the threshold model.

RESULTS

Plukenetia was resolved into two major groups, which we refer to as the pinnately- and palmately-veined clades. Our analyses suggest Plukenetia originated in the Amazon or Atlantic Forest of Brazil during the Oligocene (28.7 Mya) and migrated/dispersed between those regions and Central America/Mexico throughout the Miocene. Trans-oceanic dispersals explain the pantropical distribution of Plukenetia, including from the Amazon to Africa in the Early Miocene (17.4 Mya), followed by Africa to Madagascar and Africa to Southeast Asia in the Late Miocene (9.4 Mya) and Pliocene (4.5 Mya), respectively. We infer a single origin of large seeds in the ancestor of Plukenetia. Seed size fits a Brownian motion model of trait evolution and is moderately to strongly associated with plant size, fruit type/dispersal syndrome, and seedling ecology. Biome shifts were not drivers of seed size, although there was a weak association with a transition to fire prone semi-arid savannas.

CONCLUSIONS

The major relationships among the species of Plukenetia are now well-resolved. Our biogeographical analyses support growing evidence that many pantropical distributions developed by periodic trans-oceanic dispersals throughout the Miocene and Pliocene. Selection on a combination of traits contributed to seed size variation, while movement between forest edge/light gap and canopy niches likely contributed to the seed size extremes in Plukenetia.

Comparative Transcriptome Analysis Reveals an Efficient Mechanism of α-Linolenic Acid in Tree Peony Seeds

Tree peony (Paeonia section Moutan DC.) species are woody oil crops with high unsaturated fatty acid content, including α-linolenic acid (ALA/18:3; >40% of the total fatty acid). Comparative transcriptome analyses were carried out to uncover the underlying mechanisms responsible for high and low ALA content in the developing seeds of P. rockii and P. lutea, respectively. Expression analysis of acyl lipid metabolism genes revealed upregulation of select genes involved in plastidial fatty acid synthesis, acyl editing, desaturation, and triacylglycerol assembly in seeds of P. rockii relative to P. lutea. Also, in association with ALA content in seeds, transcript levels for fatty acid desaturases (SAD, FAD2, and FAD3), which encode enzymes necessary for polyunsaturated fatty acid synthesis, were higher in P. rockii compared to P. lutea. Furthermore, the overexpression of PrFAD2 and PrFAD3 in Arabidopsis increased linoleic and ALA content, respectively, and modulated the final ratio 18:2/18:3 in the seed oil. In conclusion, we identified the key steps and validated the necessary desaturases that contribute to efficient ALA synthesis in a woody oil crop. Together, these results will aid to increase essential fatty acid content in seeds of tree peonies and other crops of agronomic interest.

Transcriptome analysis reveals a composite molecular map linked to unique seed oil profile of Neocinnamomum caudatum (Nees) Merr

BACKGROUND

Neocinnamomum caudatum (Nees) Merr., a biodiesel tree species in the subtropical areas of South China, India and Burma, is distinctive from other species in Lauraceae family and its seed oil is rich in linoleic acid (18:2) and stearic acid (18:0). However, there is little genetic information about this species so far. In this study, a transcriptomic analysis on developing seeds of N. caudatum was conducted in an attempt to discern the molecular mechanisms involving the control of the fatty acid (FA) and triacylglycerol (TAG) biosynthesis.

RESULTS

Transcriptome analysis revealed 239,703 unigenes with an average length of 436 bp and 137 putative biomarkers that are related to FA formation and TAG biosynthesis. The expression patterns of genes encoding β-ketoacyl-acyl carrier protein synthase I (KASI), β- ketoacyl-acyl carrier protein synthase II (KASII), stearoyl-ACP desaturase (SAD), fatty acid desaturase 2 (FAD2), fatty acid desaturase 8 (FAD8) and acyl-ACP thioesterase A/B (FATA/B) were further validated by qRT-PCR. These genes displayed a very similar expression pattern in two distinct assays. Moreover, sequence analysis of different FATBs from diverse plant species revealed that NcFATB is structurally different from its counterpart in other species in producing medium-chain saturated FAs. Concertedly, heterologous expression of NcFATB in E. coli BL21 (DE3) strain showed that this corresponding expressed protein, NcFATB, prefers long-chain saturated fatty acyl-ACP over medium-chain fatty acyl-ACP as substrate.

CONCLUSIONS

Transcriptome analysis of developing N. caudatum seeds revealed a composite molecular map linked to the FA formation and oil biosynthesis in this biodiesel tree species. The substrate preference of NcFATB for long-chain saturated FAs is likely to contribute to its unique seed oil profile rich in stearic acid. Our findings demonstrate that in the tree species of Lauraceae family, the FATB enzymes producing long-chain FAs are structurally distinct from those producing medium-chain FAs, thereby suggesting that the FATB genes may serve as a biomarker for the classification of tree species of Lauraceae family.

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

Background

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.

Results

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.

Conclusions

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.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4774-y) contains supplementary material, which is available to authorized users.

Transcriptomic analysis of topping-induced axillary shoot outgrowth in Nicotiana tabacum

Topping is an important agronomic practice that significantly impacts the yield of various crop plants. Topping and the regulation of axillary shoot outgrowth are common agronomic practices in tobacco. However, the effects of topping on gene expression in tobacco remain unknown. We applied the Illumina HiSeq™ 2000 platform and analyzed differentially expressed genes (DEGs) from untopped and topped plants to study the global changes in gene expression in response to topping. We found that the number of DEGs varied from 7609 to 18,770 based on the reads per kilobase per million mapped reads (RPKM) values. The Gene Ontology (GO) enrichment analysis revealed that the cellular carbohydrate metabolic process and the disaccharide metabolic process, which may contribute to starch accumulation and stress/defense, were overrepresented terms for the DEGs. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that many DEGs were involved in starch and sucrose metabolism, glycolysis/gluconeogenesis, pyruvate metabolism, and plant hormone signal transduction, among other processes. The knowledge gained will improve our understanding of the processes of axillary shoot formation and enlargement at the transcriptional level. This study lays a solid foundation for future studies on molecular mechanisms underlying the growth of axillary shoots.

De novo transcriptome assembly and comparative analysis between male and benzyladenine-induced female inflorescence buds of Plukenetia volubilis

Plukenetia volubilis is a promising oilseed crop due to its seeds being rich in unsaturated fatty acids, especially alpha-linolenic acid. P. volubilis is monoecious, with separate male and female flowers on the same inflorescence. We previously reported that male flowers were converted to female flowers by exogenous cytokinin (6-benzyladenine, 6-BA) treatment in P. volubilis. To identify candidate genes associated with floral sex differentiation of P. volubilis, we performed de novo transcriptome assembly and comparative analysis on control male inflorescence buds (MIB) and female inflorescence buds (FIB) induced by 6-BA using Illumina sequencing technology. A total of 57,664 unigenes with an average length of 979 bp were assembled from 104.1 million clean reads, and 45,235 (78.45%) unigenes were successfully annotated in the public databases. Notably, Gene Ontology analyses revealed that 4193 and 3880 unigenes were enriched in the categories of reproduction and reproductive processes, respectively. Differential expression analysis identified 1385 differentially expressed unigenes between MIB and FIB, of which six unigenes related to cytokinin and auxin signaling pathways and 16 important transcription factor (TF) genes including MADS-box family members were identified. In particular, several unigenes encoding important TFs, such as homologs of CRABS CLAW, RADIALIS-like 1, RADIALIS-like 2, HECATE 2, WUSCHEL-related homeobox 9, and SUPERMAN, were expressed at higher levels in FIB than in MIB. The expression patterns of the 36 selected unigenes revealed by transcriptome analysis were successfully validated by quantitative real-time PCR. This study not only provides comprehensive gene expression profiles of P. volubilis inflorescence buds, but also lays the foundation for research on the molecular mechanism of floral sex determination in P. volubilis and other monoecious plants.

Using Transcriptome Analysis to Identify Genes Involved in Switchgrass Flower Reversion

Floral reversion is a process in which differentiated floral organs revert back to vegetative organs. Although this phenomenon has been described for decades, the underlying molecular mechanisms remain unclear. In this study, we found that immature switchgrass (Panicum virgatum) inflorescences can revert to neonatal shoots when incubated on a basal medium with benzylaminopurine. We used anatomical and histological methods to verify that these shoots were formed from floret primordia through flower reversion. To further explore the gene regulation of floral reversion in switchgrass, the transcriptome of reversed, unreversed, and uncultured immature inflorescences were analyzed and 517 genes were identified as participating in flower reversion. Annotation using non-redundant databases revealed that these genes are involved in plant hormone biosynthesis and signal transduction, starch and sucrose metabolism, DNA replication and modification, and other processes crucial for switchgrass flower reversion. When four of the genes were overexpressed in Arabidopsis thaliana, vegetative growth was facilitated and reproductive growth was inhibited in transgenic plants. This study provides a basic understanding of genes regulating the floral transition in switchgrass and will promote the research of floral reversion and flower maintenance.

Fatty Acid and Associated Gene Expression Analyses of Three Tree Peony Species Reveal Key Genes for α-Linolenic Acid Synthesis in Seeds

The increasing demand for healthy edible oil has generated the need to identify promising oil crops. Tree peony (Paeonia section Moutan DC.) is a woody oil crop with α-linolenic acid (ALA) contributing for 45% of the total fatty acid (FA) content in seeds. Molecular and genetic differences that contribute to varied FA content and composition among the wild peony species are, however, poorly understood. Analyses of FA content and composition during seed development in three tree peony species (Paeonia rockii, P. potaninii, and P. lutea) showed varied FA content among them with highest in P. rockii, followed by P. potaninii, and P. lutea. Total FA content among these species increased with seed development and reached its maximum in its final stage. Seed FA composition analysis of the three species also revealed that ALA (C18:3) was the most abundant, followed by oleic (C18:1) and linoleic (C18:2) acids. Additionally, quantitative real-time RT-PCR analyses of 10 key seed oil synthesis genes in the three tree peony species revealed that FAD3, FAD2, β-PDHC, LPAAT, and Oleosin gene expression levels positively correlate with total FA content and rate of accumulation. Specifically, the abundance of FAD3 transcripts in P. rockii compared with P. potaninii, and P. lutea suggests that FAD3 might play an important role in synthesis of ALA via phosphatidylcholine-derived pathway. Overall, comparative analyses of FA content and composition in three different peony species revealed a correlation between efficient lipid accumulation and lipid gene expression during seed development. Further characterization and metabolic engineering of these key genes from peonies will allow for subsequent improvement of tree peony oil quality and production.

RNA Sequencing and Coexpression Analysis Reveal Key Genes Involved in α-Linolenic Acid Biosynthesis in Perilla frutescens Seed

Perilla frutescen is used as traditional food and medicine in East Asia. Its seeds contain high levels of α-linolenic acid (ALA), which is important for health, but is scarce in our daily meals. Previous reports on RNA-seq of perilla seed had identified fatty acid (FA) and triacylglycerol (TAG) synthesis genes, but the underlying mechanism of ALA biosynthesis and its regulation still need to be further explored. So we conducted Illumina RNA-sequencing in seven temporal developmental stages of perilla seeds. Sequencing generated a total of 127 million clean reads, containing 15.88 Gb of valid data. The de novo assembly of sequence reads yielded 64,156 unigenes with an average length of 777 bp. A total of 39,760 unigenes were annotated and 11,693 unigenes were found to be differentially expressed in all samples. According to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, 486 unigenes were annotated in the “lipid metabolism” pathway. Of these, 150 unigenes were found to be involved in fatty acid (FA) biosynthesis and triacylglycerol (TAG) assembly in perilla seeds. A coexpression analysis showed that a total of 104 genes were highly coexpressed (r > 0.95). The coexpression network could be divided into two main subnetworks showing over expression in the medium or earlier and late phases, respectively. In order to identify the putative regulatory genes, a transcription factor (TF) analysis was performed. This led to the identification of 45 gene families, mainly including the AP2-EREBP, bHLH, MYB, and NAC families, etc. After coexpression analysis of TFs with highly expression of FAD2 and FAD3 genes, 162 TFs were found to be significantly associated with two FAD genes (r > 0.95). Those TFs were predicted to be the key regulatory factors in ALA biosynthesis in perilla seed. The qRT-PCR analysis also verified the relevance of expression pattern between two FAD genes and partial candidate TFs. Although it has been reported that some TFs are involved in seed development, more direct evidence is still needed to verify their function. However, these findings can provide clues to reveal the possible molecular mechanisms of ALA biosynthesis and its regulation in perilla seed.

Unravelling a stearidonic acid-rich triacylglycerol biosynthetic pathway in the developing seeds of Buglossoides arvensis: A transcriptomic landscape

Buglossoides arvensis is an emerging oilseed crop that is rich in stearidonic acid (SDA) and has several potential applications in human health and nutrition. The molecular basis of SDA biosynthesis in this plant remains unknown due to lack of genomic information. To unravel key genes involved in SDA-rich triacylglycerol (TAG) biosynthesis, we performed transcriptome sequencing of pooled mRNA from five different developmental stages of B. arvensis seeds using Illumina NextSeq platform. De novo transcriptome assembly generated 102,888 clustered transcripts from 39.83 million high-quality reads. Of these, 62.1% and 55.54% of transcripts were functionally annotated using Uniprot-Viridiplantae and KOG databases, respectively. A total of 10,021 SSR-containing sequences were identified using the MISA tool. Deep mining of transcriptome assembly using in silico tools led to the identification of genes involved in fatty acid and TAG biosynthesis. Expression profiling of 17 key transcripts involved in fatty acid desaturation and TAG biosynthesis showed expression patterns specific to the development stage that positively correlated with polyunsaturated fatty acid accumulation in the developing seeds. This first comprehensive transcriptome analysis provides the basis for future research on understanding molecular mechanisms of SDA-rich TAG accumulation in B. arvensis and aids in biotechnological production of SDA in other oilseed crops.

Transcriptome analysis and identification of genes associated with ω-3 fatty acid biosynthesis in Perilla frutescens (L.) var. frutescens

Background

Perilla (Perilla frutescens (L.) var frutescens) produces high levels of α-linolenic acid (ALA), a ω-3 fatty acid important to health and development. To uncover key genes involved in fatty acid (FA) and triacylglycerol (TAG) synthesis in perilla, we conducted deep sequencing of cDNAs from developing seeds and leaves for understanding the mechanism underlying ALA and seed TAG biosynthesis.

Results

Perilla cultivar Dayudeulkkae contains 66.0 and 56.2 % ALA in seeds and leaves, respectively. Using Illumina HiSeq 2000, we have generated a total of 392 megabases of raw sequences from four mRNA samples of seeds at different developmental stages and one mature leaf sample of Dayudeulkkae. De novo assembly of these sequences revealed 54,079 unique transcripts, of which 32,237 belong to previously annotated genes. Among the annotated genes, 66.5 % (21,429 out of 32,237) showed highest sequences homology with the genes from Mimulus guttatus, a species placed under the same Lamiales order as perilla. Using Arabidopsis acyl-lipid genes as queries, we searched the transcriptome and identified 540 unique perilla genes involved in all known pathways of acyl-lipid metabolism. We characterized the expression profiles of 43 genes involved in FA and TAG synthesis using quantitative PCR. Key genes were identified through sequence and gene expression analyses.

Conclusions

This work is the first report on building transcriptomes from perilla seeds. The work also provides the first comprehensive expression profiles for genes involved in seed oil biosynthesis. Bioinformatic analysis indicated that our sequence collection represented a major transcriptomic resource for perilla that added valuable genetic information in order Lamiales. Our results provide critical information not only for studies of the mechanisms involved in ALA synthesis, but also for biotechnological production of ALA in other oilseeds.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2805-0) contains supplementary material, which is available to authorized users.

Transcriptome analysis and identification of genes associated with ω-3 fatty acid biosynthesis in Perilla frutescens (L.) var. frutescens

BACKGROUND

Perilla (Perilla frutescens (L.) var frutescens) produces high levels of α-linolenic acid (ALA), a ω-3 fatty acid important to health and development. To uncover key genes involved in fatty acid (FA) and triacylglycerol (TAG) synthesis in perilla, we conducted deep sequencing of cDNAs from developing seeds and leaves for understanding the mechanism underlying ALA and seed TAG biosynthesis.

RESULTS

Perilla cultivar Dayudeulkkae contains 66.0 and 56.2 % ALA in seeds and leaves, respectively. Using Illumina HiSeq 2000, we have generated a total of 392 megabases of raw sequences from four mRNA samples of seeds at different developmental stages and one mature leaf sample of Dayudeulkkae. De novo assembly of these sequences revealed 54,079 unique transcripts, of which 32,237 belong to previously annotated genes. Among the annotated genes, 66.5 % (21,429 out of 32,237) showed highest sequences homology with the genes from Mimulus guttatus, a species placed under the same Lamiales order as perilla. Using Arabidopsis acyl-lipid genes as queries, we searched the transcriptome and identified 540 unique perilla genes involved in all known pathways of acyl-lipid metabolism. We characterized the expression profiles of 43 genes involved in FA and TAG synthesis using quantitative PCR. Key genes were identified through sequence and gene expression analyses.

CONCLUSIONS

This work is the first report on building transcriptomes from perilla seeds. The work also provides the first comprehensive expression profiles for genes involved in seed oil biosynthesis. Bioinformatic analysis indicated that our sequence collection represented a major transcriptomic resource for perilla that added valuable genetic information in order Lamiales. Our results provide critical information not only for studies of the mechanisms involved in ALA synthesis, but also for biotechnological production of ALA in other oilseeds.

High light exposure on seed coat increases lipid accumulation in seeds of castor bean (Ricinus communis L.), a nongreen oilseed crop

Little was known on how sunlight affects the seed metabolism in nongreen seeds. Castor bean (Ricinus communis L.) is a typical nongreen oilseed crop and its seed oil is an important feedstock in industry. In this study, photosynthetic activity of seed coat tissues of castor bean in natural conditions was evaluated in comparison to shaded conditions. Our results indicate that exposure to high light enhances photosynthetic activity in seed coats and consequently increases oil accumulation. Consistent results were also reached using cultured seeds. High-throughput RNA-Seq analyses further revealed that genes involved in photosynthesis and carbon conversion in both the Calvin-Benson cycle and malate transport were differentially expressed between seeds cultured under light and dark conditions, implying several venues potentially contributing to light-enhanced lipid accumulation such as increased reducing power and CO2 refixation which underlie the overall lipid biosynthesis. This study demonstrated the effects of light exposure on oil accumulation in nongreen oilseeds and greatly expands our understanding of the physiological roles that light may play during seed development in nongreen oilseeds. Essentially, our studies suggest that potential exists to enhance castor oil yield through increasing exposure of the inflorescences to sunlight either by genetically changing the plant architecture (smart canopy) or its growing environment.

Selecting Superior De Novo Transcriptome Assemblies: Lessons Learned by Leveraging the Best Plant Genome

Whereas de novo assemblies of RNA-Seq data are being published for a growing number of species across the tree of life, there are currently no broadly accepted methods for evaluating such assemblies. Here we present a detailed comparison of 99 transcriptome assemblies, generated with 6 de novo assemblers including CLC, Trinity, SOAP, Oases, ABySS and NextGENe. Controlled analyses of de novo assemblies for Arabidopsis thaliana and Oryza sativa transcriptomes provide new insights into the strengths and limitations of transcriptome assembly strategies. We find that the leading assemblers generate reassuringly accurate assemblies for the majority of transcripts. At the same time, we find a propensity for assemblers to fail to fully assemble highly expressed genes. Surprisingly, the instance of true chimeric assemblies is very low for all assemblers. Normalized libraries are reduced in highly abundant transcripts, but they also lack 1000s of low abundance transcripts. We conclude that the quality of de novo transcriptome assemblies is best assessed through consideration of a combination of metrics: 1) proportion of reads mapping to an assembly 2) recovery of conserved, widely expressed genes, 3) N50 length statistics, and 4) the total number of unigenes. We provide benchmark Illumina transcriptome data and introduce SCERNA, a broadly applicable modular protocol for de novo assembly improvement. Finally, our de novo assembly of the Arabidopsis leaf transcriptome revealed ~20 putative Arabidopsis genes lacking in the current annotation.

Transcriptome Characterization for Non-Model Endangered Lycaenids, Protantigius superans and Spindasis takanosis, Using Illumina HiSeq 2500 Sequencing

The Lycaenidae butterflies, Protantigius superans and Spindasis takanosis, are endangered insects in Korea known for their symbiotic association with ants. However, necessary genomic and transcriptomics data are lacking in these species, limiting conservation efforts. In this study, the P. superans and S. takanosis transcriptomes were deciphered using Illumina HiSeq 2500 sequencing. The P. superans and S. takanosis transcriptome data included a total of 254,340,693 and 245,110,582 clean reads assembled into 159,074 and 170,449 contigs and 107,950 and 121,140 unigenes, respectively. BLASTX hits (E-value of 1.0 × 10⁻⁵) against the known protein databases annotated a total of 46,754 and 51,908 transcripts for P. superans and S. takanosis. Approximately 41.25% and 38.68% of the unigenes for P. superans and S. takanosis found homologous sequences in Protostome DB (PANM-DB). BLAST2GO analysis confirmed 18,611 unigenes representing Gene Ontology (GO) terms and a total of 5259 unigenes assigned to 116 pathways for P. superans. For S. takanosis, a total of 6697 unigenes were assigned to 119 pathways using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. Additionally, 382,164 and 390,516 Simple Sequence Repeats (SSRs) were compiled from the unigenes of P. superans and S. takanosis, respectively. This is the first report to record new genes and their utilization for conservation of lycaenid species population and as a reference information for closely related species.

De novo sequencing and assembly analysis of transcriptome in the Sodom apple (Calotropis gigantea)

Background

The Sodom apple (Calotropis gigantea), a member of the Asclepiadaceae family, is a large evergreen shrub native to continental Asia and northern Africa. As an important medicinal shrub and a fiber resource plant, there is an urgent need for developing molecular markers to facilitate breeding and genetic improvement of varieties.

Results

In this study, using the Illumina high throughput sequencing technique we obtained about 45 million paired end sequencing reads, De novo assembled and generated a total of 133,634 transcripts with a mean of 1837.47 bp in length. Based on protein homology searches against available databases, a total of 21,851 unigenes were functionally annotated. In particular, many transcripts that encode for putative proteins involved in fiber and secondary metabolite biosynthesis were identified and analyzed. Key fiber genes identified were validated experimentally through Real-Time PCR technique. Various transcription factors involved in regulating plant response to abiotic stress were also identified. In addition, based on the unigene sequences assembled, 11,623 microsatellites loci were detected, which provide very useful resources for developing microsatellite molecular markers.

Conclusion

This study is the first report on transcriptome information in the Calotropis species and provides rich gene transcript resources for conducting further studies on understanding the molecular basis of fiber and secondary metabolite biosynthesis, serving the genetic improvement and resource utilization in Calotropis plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1908-3) contains supplementary material, which is available to authorized users.

De Novo Transcriptome Analysis of Warburgia ugandensis to Identify Genes Involved in Terpenoids and Unsaturated Fatty Acids Biosynthesis

The bark of Warburgia ugandensis (Canellaceae family) has been used as a medicinal source for a long history in many African countries. The presence of diverse terpenoids and abundant polyunsaturated fatty acids (PUFAs) in this organ contributes to its broad range of pharmacological properties. Despite its medicinal and economic importance, the knowledge on the biosynthesis of terpenoid and unsaturated fatty acid in W. ugandensis bark remains largely unknown. Therefore, it is necessary to construct a genomic and/or transcriptomic database for the functional genomics study on W. ugandensis. The chemical profiles of terpenoids and fatty acids between the bark and leaves of W. ugandensis were compared by gas chromatography-mass spectrometry (GC-MS) analysis. Meanwhile, the transcriptome database derived from both tissues was created using Illumina sequencing technology. In total, about 17.1 G clean nucleotides were obtained, and de novo assembled into 72,591 unigenes, of which about 38.06% can be aligned to the NCBI non-redundant protein database. Many candidate genes in the biosynthetic pathways of terpenoids and unsaturated fatty acids were identified, including 14 unigenes for terpene synthases. Furthermore, 2,324 unigenes were discovered to be differentially expressed between both tissues; the functions of those differentially expressed genes (DEGs) were predicted by gene ontology enrichment and metabolic pathway enrichment analyses. In addition, the expression of 12 DEGs with putative roles in terpenoid and unsaturated fatty acid metabolic pathways was confirmed by qRT-PCRs, which was consistent with the data of the RNA-sequencing. In conclusion, we constructed a comprehensive transcriptome dataset derived from the bark and leaf of W. ugandensis, which forms the basis for functional genomics studies on this plant species. Particularly, the comparative analysis of the transcriptome data between the bark and leaf will provide critical clues to reveal the regulatory mechanisms underlying the biosynthesis of terpenoids and PUFAs in W. ugandensis.

Transcriptomic Analysis of the Regulation of Rhizome Formation in Temperate and Tropical Lotus (Nelumbo nucifera)

Rhizome is the storage organ of lotus derived from modified stems. The development of rhizome is a complex process and depends on the balanced expression of the genes that is controlled by environmental and endogenous factors. However, little is known about the mechanism that regulates rhizome girth enlargement. In this study, using RNA-seq, transcriptomic analyses were performed at three rhizome developmental stages-the stolon, middle swelling and later swelling stage -in the cultivars 'ZO' (temperate lotus with enlarged rhizome) and 'RL' (tropical lotus with stolon). About 348 million high-quality reads were generated, and 88.5% of the data were mapped to the reference genome. Of 26783 genes identified, 24069 genes were previously predicted in the reference, and 2714 genes were novel transcripts. Moreover, 8821 genes were differentially expressed between the cultivars at the three stages. Functional analysis identified that these genes were significantly enriched in pathways carbohydrate metabolism and plant hormone signal transduction. Twenty-two genes involved in photoperiod pathway, starch metabolism and hormone signal transduction were candidate genes inducing rhizome girth enlargement. Comparative transcriptomic analysis detected several differentially expressed genes and potential candidate genes required for rhizome girth enlargement, which lay a foundation for future studies on molecular mechanisms underlying rhizome formation.

Selection of Reliable Reference Genes for Gene Expression Studies of a Promising Oilseed Crop, Plukenetia volubilis, by Real-Time Quantitative PCR

Real-time quantitative PCR (RT-qPCR) is a reliable and widely used method for gene expression analysis. The accuracy of the determination of a target gene expression level by RT-qPCR demands the use of appropriate reference genes to normalize the mRNA levels among different samples. However, suitable reference genes for RT-qPCR have not been identified in Sacha inchi (Plukenetia volubilis), a promising oilseed crop known for its polyunsaturated fatty acid (PUFA)-rich seeds. In this study, using RT-qPCR, twelve candidate reference genes were examined in seedlings and adult plants, during flower and seed development and for the entire growth cycle of Sacha inchi. Four statistical algorithms (delta cycle threshold (ΔC_t), BestKeeper, geNorm, and NormFinder) were used to assess the expression stabilities of the candidate genes. The results showed that ubiquitin-conjugating enzyme (UCE), actin (ACT) and phospholipase A22 (PLA) were the most stable genes in Sacha inchi seedlings. For roots, stems, leaves, flowers, and seeds from adult plants, 30S ribosomal protein S13 (RPS13), cyclophilin (CYC) and elongation factor-1alpha (EF1α) were recommended as reference genes for RT-qPCR. During the development of reproductive organs, PLA, ACT and UCE were the optimal reference genes for flower development, whereas UCE, RPS13 and RNA polymerase II subunit (RPII) were optimal for seed development. Considering the entire growth cycle of Sacha inchi, UCE, ACT and EF1α were sufficient for the purpose of normalization. Our results provide useful guidelines for the selection of reliable reference genes for the normalization of RT-qPCR data for seedlings and adult plants, for reproductive organs, and for the entire growth cycle of Sacha inchi.

Fatty acid composition of developing tree peony (Paeonia section Moutan DC.) seeds and transcriptome analysis during seed development

BACKGROUND

Tree peony (Paeonia section Moutan DC.) is known for its excellent ornamental and medicinal values. In 2011, seeds from P. ostii have been identified as novel resource of α-linolenic acid (ALA) for seed oil production and development in China. However, the molecular mechanism on biosynthesis of unsaturated fatty acids in tree peony seeds remains unknown. Therefore, transcriptome data is needed to better understand the underlying mechanisms.

RESULTS

In this study, lipid accumulation contents were measured using GC-MS methods across developing tree peony seeds, which exhibited an extraordinary ALA content (49.3%) in P. ostii mature seeds. Transcriptome analysis was performed using Illumina sequencing platform. A total of 144 million 100-bp paired-end reads were generated from six libraries, which identified 175,874 contigs. In the KEGG Orthology enrichment of differentially expressed genes, lipid metabolism pathways were highly represented categories. Using this data we identified 388 unigenes that may be involved in de novo fatty acid and triacylglycerol biosynthesis. In particular, three unigenes (SAD, FAD2 and FAD8) encoding fatty acid desaturase with high expression levels in the fast oil accumulation stage compared with the initial stage of seed development were identified.

CONCLUSIONS

This study provides the first comprehensive genomic resources characterizing tree peony seeds gene expression at the transcriptional level. These data lay the foundation for further understanding of molecular mechanism responsible for lipid biosynthesis and the high unsaturated fatty acids (especially ALA) accumulation. Meanwhile, it provides theoretical base for potential oilseed application in the respect of n-6 to n-3 ratio for human diets and future regulation of target healthy components of oils.

Fatty acid composition of developing tree peony (Paeonia section Moutan DC.) seeds and transcriptome analysis during seed development

BACKGROUND

Tree peony (Paeonia section Moutan DC.) is known for its excellent ornamental and medicinal values. In 2011, seeds from P. ostii have been identified as novel resource of α-linolenic acid (ALA) for seed oil production and development in China. However, the molecular mechanism on biosynthesis of unsaturated fatty acids in tree peony seeds remains unknown. Therefore, transcriptome data is needed to better understand the underlying mechanisms.

RESULTS

In this study, lipid accumulation contents were measured using GC-MS methods across developing tree peony seeds, which exhibited an extraordinary ALA content (49.3%) in P. ostii mature seeds. Transcriptome analysis was performed using Illumina sequencing platform. A total of 144 million 100-bp paired-end reads were generated from six libraries, which identified 175,874 contigs. In the KEGG Orthology enrichment of differentially expressed genes, lipid metabolism pathways were highly represented categories. Using this data we identified 388 unigenes that may be involved in de novo fatty acid and triacylglycerol biosynthesis. In particular, three unigenes (SAD, FAD2 and FAD8) encoding fatty acid desaturase with high expression levels in the fast oil accumulation stage compared with the initial stage of seed development were identified.

CONCLUSIONS

This study provides the first comprehensive genomic resources characterizing tree peony seeds gene expression at the transcriptional level. These data lay the foundation for further understanding of molecular mechanism responsible for lipid biosynthesis and the high unsaturated fatty acids (especially ALA) accumulation. Meanwhile, it provides theoretical base for potential oilseed application in the respect of n-6 to n-3 ratio for human diets and future regulation of target healthy components of oils.

Transcriptome analysis of canola (Brassica napus) under salt stress at the germination stage

Canola (Brassica napus) is one of the most important oil crops in the world. However, its yield has been constrained by salt stress. In this study, transcriptome profiles were explored using Digital Gene Expression (DGE) at 0, 3, 12 and 24 hours after H2O (control) and NaCl treatments on B. napus roots at the germination stage. Comparisons of gene-expression between the control and the treatment were conducted after tag-mapping to the sequenced Brassica rapa genome. The differentially expressed genes during the time course of salt stress were focused on, and 163 genes were identified to be differentially expressed at all the time points. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed that some of the genes were involved in proline metabolism, inositol metabolism, carbohydrate metabolic processes and oxidation-reduction processes and may play vital roles in the salt-stress response at the germination stage. Thus, this study provides new candidate salt stress responding genes, which may function in novel putative nodes in the molecular pathways of salt stress resistance.

Characterization of the rainbow trout spleen transcriptome and identification of immune-related genes

Resistance against diseases affects profitability of rainbow trout. Limited information is available about functions and mechanisms of teleost immune pathways. Immunogenomics provides powerful tools to determine disease resistance genes/gene pathways and develop genetic markers for genomic selection. RNA-Seq sequencing of the rainbow trout spleen yielded 93,532,200 reads (100 bp). High quality reads were assembled into 43,047 contigs. 26,333 (61.17%) of the contigs had hits to the NR protein database and 7024 (16.32%) had hits to the KEGG database. Gene ontology showed significant percentages of transcripts assigned to binding (51%), signaling (7%), response to stimuli (9%) and receptor activity (4%) suggesting existence of many immune-related genes. KEGG annotation revealed 2825 sequences belonging to "organismal systems" with the highest number of sequences, 842 (29.81%), assigned to immune system. A number of sequences were identified for the first time in rainbow trout belonging to Toll-like receptor signaling (35), B cell receptor signaling pathway (44), T cell receptor signaling pathway (56), chemokine signaling pathway (73), Fc gamma R-mediated phagocytosis (52), leukocyte transendothelial migration (60) and NK cell mediated cytotoxicity (42). In addition, 51 transcripts were identified as spleen-specific genes. The list includes 277 full-length cDNAs. The presence of a large number of immune-related genes and pathways similar to other vertebrates suggests that innate and adaptive immunity in fish are conserved. This study provides deep-sequence data of rainbow trout spleen transcriptome and identifies many new immune-related genes and full-length cDNAs. This data will help identify allelic variations suitable for genomic selection and genetic manipulation in aquaculture.

Characterization of the miiuy croaker (Miichthys miiuy) transcriptome and development of immune-relevant genes and molecular markers

Background

The miiuy croaker (Miichthys miiuy) is an important species of marine fish that supports capture fisheries and aquaculture. At present commercial scale aquaculture of this species is limited due to diseases caused by pathogens and parasites which restrict production and limit commercial value. The lack of transcriptomic and genomic information for the miiuy croaker limits the ability of researchers to study the pathogenesis and immune system of this species. In this study we constructed a cDNA library from liver, spleen and kidney which was sequenced using Illumina paired-end sequencing to enable gene discovery and molecular marker development.

Principal Findings

In our study, a total of 69,071 unigenes with an average length of 572 bp were obtained. Of these, 45,676 (66.13%) were successfully annotated in public databases. The unigenes were also annotated with Gene Ontology, Clusters of Orthologous Groups and KEGG pathways. Additionally, 498 immune-relevant genes were identified and classified. Furthermore, 14,885 putative simple sequence repeats (cSSRs) and 8,510 putative single nucleotide polymorphisms (SNPs) were identified from the 69,071 unigenes.

Conclusion

The miiuy croaker (Miichthys miiuy) transcriptome data provides a large resource to identify new genes involved in many processes including those involved in the response to pathogens and diseases. Furthermore, the thousands of potential cSSR and SNP markers found in this study are important resources with respect to future development of molecular marker assisted breeding programs for the miiuy croaker.

High-throughput sequencing and de novo assembly of Brassica oleracea var. Capitata L. for transcriptome analysis

Background

The cabbage, Brassica oleracea var. capitata L., has a distinguishable phenotype within the genus Brassica. Despite the economic and genetic importance of cabbage, there is little genomic data for cabbage, and most studies of Brassica are focused on other species or other B. oleracea subspecies. The lack of genomic data for cabbage, a non-model organism, hinders research on its molecular biology. Hence, the construction of reliable transcriptomic data based on high-throughput sequencing technologies is needed to enhance our understanding of cabbage and provide genomic information for future work.

Methodology/Principal Findings

We constructed cDNAs from total RNA isolated from the roots, leaves, flowers, seedlings, and calcium-limited seedling tissues of two cabbage genotypes: 102043 and 107140. We sequenced a total of six different samples using the Illumina HiSeq platform, producing 40.5 Gbp of sequence data comprising 401,454,986 short reads. We assembled 205,046 transcripts (≥ 200 bp) using the Velvet and Oases assembler and predicted 53,562 loci from the transcripts. We annotated 35,274 of the loci with 55,916 plant peptides in the Phytozome database. The average length of the annotated loci was 1,419 bp. We confirmed the reliability of the sequencing assembly using reverse-transcriptase PCR to identify tissue-specific gene candidates among the annotated loci.

Conclusion

Our study provides valuable transcriptome sequence data for B. oleracea var. capitata L., offering a new resource for studying B. oleracea and closely related species. Our transcriptomic sequences will enhance the quality of gene annotation and functional analysis of the cabbage genome and serve as a material basis for future genomic research on cabbage. The sequencing data from this study can be used to develop molecular markers and to identify the extreme differences among the phenotypes of different species in the genus Brassica.

Transcriptomic analysis of endangered Chinese salamander: identification of immune, sex and reproduction-related genes and genetic markers

BACKGROUND

The Chinese salamander (Hynobius chinensis), an endangered amphibian species of salamander endemic to China, has attracted much attention because of its value of studying paleontology evolutionary history and decreasing population size. Despite increasing interest in the Hynobius chinensis genome, genomic resources for the species are still very limited. A comprehensive transcriptome of Hynobius chinensis, which will provide a resource for genome annotation, candidate genes identification and molecular marker development should be generated to supplement it.

PRINCIPAL FINDINGS

We performed a de novo assembly of Hynobius chinensis transcriptome by Illumina sequencing. A total of 148,510 nonredundant unigenes with an average length of approximately 580 bp were obtained. In all, 60,388 (40.66%) unigenes showed homologous matches in at least one database and 33,537 (22.58%) unigenes were annotated by all four databases. In total, 41,553 unigenes were categorized into 62 sub-categories by BLAST2GO search, and 19,468 transcripts were assigned to 140 KEGG pathways. A large number of unigenes involved in immune system, local adaptation, reproduction and sex determination were identified, as well as 31,982 simple sequence repeats (SSRs) and 460,923 putative single nucleotide polymorphisms (SNPs).

CONCLUSION

This dataset represents the first transcriptome analysis of the Chinese salamander (Hynobius chinensis), an endangered species, to be also the first time of hynobiidae. The transcriptome will provide valuable resource for further research in discovery of new genes, protection of population, adaptive evolution and survey of various pathways, as well as development of molecule markers in Chinese salamander; and reference information for closely related species.

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

The study was designed to assess acceptability and side-effects of consumption of sacha inchi oil, rich in α-linolenic acid and sunflower oil, rich in linoleic acid, in adult human subjects. Thirty subjects received 10 or 15ml daily of sacha inchi or sunflower oil for 4months. Acceptability was assessed with daily self-report and with a Likert test at the end of the study. Safety was assessed with self- recording of side-effects and with hepatic and renal markers. Primary efficacy variables were the change in lipid profile. Subjects reported low acceptability of sacha inchi oil at week-1 (37.5%). However, since week-6, acceptability was significantly increased to 81.25-93.75%. No differences were observed in acceptability with respect to sex or oil volume (P>0.05). Most frequent adverse effects during first weeks of consuming sacha inchi oil or sunflower oil were nauseas. The side-effects were reduced with time. Biochemical markers of hepatic and kidney function were maintained unchanged. Serum total cholesterol and LDL cholesterol levels and arterial blood pressure were lowered with both oils (P<0.05). Higher HDL-cholesterol was observed with sacha inchi oil at month-4. In conclusion, sacha inchi oil consumed has good acceptability after week-1 of consumption and it is safety.

Mining genes involved in the stratification of Paris polyphylla seeds using high-throughput embryo transcriptome sequencing

Background

Paris polyphylla var. yunnanensis is an important medicinal plant. Seed dormancy is one of the main factors restricting artificial cultivation. The molecular mechanisms of seed dormancy remain unclear, and little genomic or transcriptome data are available for this plant.

Results

In this study, massive parallel pyrosequencing on the Roche 454-GS FLX Titanium platform was used to generate a substantial sequence dataset for the P. polyphylla embryo. 369,496 high quality reads were obtained, ranging from 50 to 1146 bp, with a mean of 219 bp. These reads were assembled into 47,768 unigenes, which included 16,069 contigs and 31,699 singletons. Using BLASTX searches of public databases, 15,757 (32.3%) unique transcripts were identified. Gene Ontology and Cluster of Orthologous Groups of proteins annotations revealed that these transcripts were broadly representative of the P. polyphylla embryo transcriptome. The Kyoto Encyclopedia of Genes and Genomes assigned 5961 of the unique sequences to specific metabolic pathways. Relative expression levels analysis showed that eleven phytohormone-related genes and five other genes have different expression patterns in the embryo and endosperm in the seed stratification process.

Conclusions

Gene annotation and quantitative RT-PCR expression analysis identified 464 transcripts that may be involved in phytohormone catabolism and biosynthesis, hormone signal, seed dormancy, seed maturation, cell wall growth and circadian rhythms. In particular, the relative expression analysis of sixteen genes (CYP707A, NCED, GA20ox2, GA20ox3, ABI2, PP2C, ARP3, ARP7, IAAH, IAAS, BRRK, DRM, ELF1, ELF2, SFR6, and SUS) in embryo and endosperm and at two temperatures indicated that these related genes may be candidates for clarifying the molecular basis of seed dormancy in P. polyphlla var. yunnanensis.