De novo assembly and characterization of Camelina sativa transcriptome by paired-end sequencing

Critical metabolic pathways and genes cooperate for epoxy fatty acid-enriched oil production in developing seeds of Vernonia galamensis, an industrial oleaginous plant

Background

Vernonia galamensis native to Africa is an annual oleaginous plant of Asteraceae family. As a newly established industrial oil crop, this plant produces high level (> 70%) of vernolic acid (cis-12-epoxyoctadeca-cis-9-enoic acid), which is an unusual epoxy fatty acid (EFA) with multiple industrial applications. Here, transcriptome analysis and fatty acid profiling from developing V. galamensis seeds were integrated to uncover the critical metabolic pathways responsible for high EFA accumulation, aiming to identify the target genes that could be used in the biotechnological production of high-value oils.

Results

Based on oil accumulation dynamics of V. galamensis seeds, we harvested seed samples from three stages (17, 38, and 45 days after pollination, DAP) representing the initial, fast and final EFA accumulation phases, and one mixed sample from different tissues for RNA-sequencing, with three biological replicates for each sample. Using Illumina platform, we have generated a total of 265 million raw cDNA reads. After filtering process, de novo assembly of clean reads yielded 67,114 unigenes with an N50 length of 1316 nt. Functional annotation resulted in the identification of almost all genes involved in diverse lipid-metabolic pathways, including the novel fatty acid desaturase/epoxygenase, diacylglycerol acyltransferases, and phospholipid:diacylglycerol acyltransferases. Expression profiling revealed that various genes associated with acyl editing, fatty acid β-oxidation, triacylglycerol assembly and oil-body formation had greater expression levels at middle developmental stage (38 DAP), which were consistent with the fast accumulation of EFA in V. galamensis developing seed, these genes were detected to play fundamental roles in EFA production. In addition, we isolated some transcription factors (such as WRI1, FUS3 and ABI4), which putatively regulated the production of V. galamensis seed oils. The transient expression of the selected genes resulted in a synergistic increase of EFA-enriched TAG accumulation in tobacco leaves. Transcriptome data were further confirmed by quantitative real-time PCR for twelve key genes in EFA biosynthesis. Finally, a comprehensive network for high EFA accumulation in V. galamensis seed was established.

Conclusions

Our findings provide new insights into molecular mechanisms underlying the natural epoxy oil production in V. galamensis. A set of genes identified here could be used as the targets to develop other oilseeds highly accumulating valued epoxy oils for commercial production.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02120-2.

Genetic variation and structural diversity in major seed proteins among and within Camelina species

Genetic variation in seed protein composition, seed protein gene expression and predictions of seed protein physiochemical properties were documented in C. sativa and other Camelina species. Seed protein diversity was examined in six Camelina species (C. hispida, C. laxa, C. microcarpa, C. neglecta, C. rumelica and C. sativa). Differences were observed in seed protein electrophoretic profiles, total seed protein content and amino acid composition between the species. Genes encoding major seed proteins (cruciferins, napins, oleosins and vicilins) were catalogued for C. sativa and RNA-Seq analysis established the expression patterns of these and other genes in developing seed from anthesis through to maturation. Examination of 187 C. sativa accessions revealed limited variation in seed protein electrophoretic profiles, though sufficient to group the majority into classes based on high MW protein profiles corresponding to the cruciferin region. C. sativa possessed four distinct types of cruciferins, named CsCRA, CsCRB, CsCRC and CsCRD, which corresponded to orthologues in Arabidopsis thaliana with members of each type encoded by homeologous genes on the three C. sativa sub-genomes. Total protein content and amino acid composition varied only slightly; however, RNA-Seq analysis revealed that CsCRA and CsCRB genes contributed > 95% of the cruciferin transcripts in most lines, whereas CsCRC genes were the most highly expressed cruciferin genes in others, including the type cultivar DH55. This was confirmed by proteomics analyses. Cruciferin is the most abundant seed protein and contributes the most to functionality. Modelling of the C. sativa cruciferins indicated that each type possesses different physiochemical attributes that were predicted to impart unique functional properties. As such, opportunities exist to create C. sativa cultivars with seed protein profiles tailored to specific technical applications.

Exploring Camelina sativa lipid metabolism regulation by combining gene co-expression and DNA affinity purification analyses

Camelina (Camelina sativa) is an annual oilseed plant that is gaining momentum as a biofuel cover crop. Understanding gene regulatory networks is essential to deciphering plant metabolic pathways, including lipid metabolism. Here, we take advantage of a growing collection of gene expression datasets to predict transcription factors (TFs) associated with the control of Camelina lipid metabolism. We identified approximately 350 TFs highly co-expressed with lipid-related genes (LRGs). These TFs are highly represented in the MYB, AP2/ERF, bZIP, and bHLH families, including a significant number of homologs of well-known Arabidopsis lipid and seed developmental regulators. After prioritizing the top 22 TFs for further validation, we identified DNA-binding sites and predicted target genes for 16 out of the 22 TFs tested using DNA affinity purification followed by sequencing (DAP-seq). Enrichment analyses of targets supported the co-expression prediction for most TF candidates, and the comparison to Arabidopsis revealed some common themes, but also aspects unique to Camelina. Within the top potential lipid regulators, we identified CsaMYB1, CsaABI3AVP1-2, CsaHB1, CsaNAC2, CsaMYB3, and CsaNAC1 as likely involved in the control of seed fatty acid elongation and CsaABI3AVP1-2 and CsabZIP1 as potential regulators of the synthesis and degradation of triacylglycerols (TAGs), respectively. Altogether, the integration of co-expression data and DNA-binding assays permitted us to generate a high-confidence and short list of Camelina TFs involved in the control of lipid metabolism during seed development.

Realizing the Potential of Camelina sativa as a Bioenergy Crop for a Changing Global Climate

Camelina sativa (L.) Crantz. is an annual oilseed crop within the Brassicaceae family. C. sativa has been grown since as early as 4000 BCE. In recent years, C. sativa received increased attention as a climate-resilient oilseed, seed meal, and biofuel (biodiesel and renewable or green diesel) crop. This renewed interest is reflected in the rapid rise in the number of peer-reviewed publications (>2300) containing “camelina” from 1997 to 2021. An overview of the origins of this ancient crop and its genetic diversity and its yield potential under hot and dry growing conditions is provided. The major biotic barriers that limit C. sativa production are summarized, including weed control, insect pests, and fungal, bacterial, and viral pathogens. Ecosystem services provided by C. sativa are also discussed. The profiles of seed oil and fatty acid composition and the many uses of seed meal and oil are discussed, including food, fodder, fuel, industrial, and medical benefits. Lastly, we outline strategies for improving this important and versatile crop to enhance its production globally in the face of a rapidly changing climate using molecular breeding, rhizosphere microbiota, genetic engineering, and genome editing approaches.

Overlapping Functions of the Paralogous Proteins AtPAP2 and AtPAP9 in Arabidopsis thaliana

Arabidopsis thaliana purple acid phosphatase 2 (AtPAP2), which is anchored to the outer membranes of chloroplasts and mitochondria, affects carbon metabolism by modulating the import of some preproteins into chloroplasts and mitochondria. AtPAP9 bears a 72% amino acid sequence identity with AtPAP2, and both proteins carry a hydrophobic motif at their C-termini. Here, we show that AtPAP9 is a tail-anchored protein targeted to the outer membrane of chloroplasts. Yeast two-hybrid and bimolecular fluorescence complementation experiments demonstrated that both AtPAP9 and AtPAP2 bind to a small subunit of rubisco 1B (AtSSU1B) and a number of chloroplast proteins. Chloroplast import assays using [35S]-labeled AtSSU1B showed that like AtPAP2, AtPAP9 also plays a role in AtSSU1B import into chloroplasts. Based on these data, we propose that AtPAP9 and AtPAP2 perform overlapping roles in modulating the import of specific proteins into chloroplasts. Most plant genomes contain only one PAP-like sequence encoding a protein with a hydrophobic motif at the C-terminus. The presence of both AtPAP2 and AtPAP9 in the Arabidopsis genome may have arisen from genome duplication in Brassicaceae. Unlike AtPAP2 overexpression lines, the AtPAP9 overexpression lines did not exhibit early-bolting or high-seed-yield phenotypes. Their differential growth phenotypes could be due to the inability of AtPAP9 to be targeted to mitochondria, as the overexpression of AtPAP2 on mitochondria enhances the capacity of mitochondria to consume reducing equivalents.

Characterizing the Leaf Transcriptome of Chrysanthemum rhombifolium (Ling et C. Shih), a Drought Resistant, Endemic Plant From China

Chrysanthemum rhombifolium (Ling et C. Shih), an endemic plant that is extremely well-adapted to harsh environments. However, little is known about its molecular biology of the plant's resistant traits against stress, or even its molecular biology of overall plant. To investigate the molecular biology of C. rhombifolium and mechanism of stress adaptation, we performed transcriptome sequencing of its leaves using an Illumina platform. A total of 130,891 unigenes were obtained, and 97,496 (~74.5%) unigenes were annotated in the public protein database. The similarity search indicated that 40,878 and 74,084 unigenes showed significant similarities to known proteins from NCBI non-redundant and Swissprot protein databases, respectively. Of these, 56,213 and 42,005 unigenes were assigned to the Gene Ontology (GO) database and Cluster of Orthologous Groups (COG), respectively, and 38,918 unigenes were mapped into five main categories, including 18 KEGG pathways. Metabolism was the largest category (23,128, 59.4%) among the main KEGG categories, suggesting active metabolic processes in C. rhombifolium. About 2,459 unigenes were annotated to have a role in defense mechanism or stress tolerance. Transcriptome analysis of C. rhombifolium revealed the presence of 12,925 microsatellites in 10,524 unigenes and mono, trip, and dinucleotides having higher polymorphism rates. The phylogenetic analysis based on GME gene among related species confirmed the reliability of the transcriptomic data. This work is the first genetic study of C. rhombifolium as a new plant resource of stress-tolerant genes. This large number of transcriptome sequences enabled us to comprehensively understand the basic genetics of C. rhombifolium and discover novel genes that will be helpful in the molecular improvement of chrysanthemums.

CamRegBase: a gene regulation database for the biofuel crop, Camelina sativa

Camelina is an annual oilseed plant from the Brassicaceae family that is gaining momentum as a biofuel winter cover crop. However, a significant limitation in further enhancing its utility as a producer of oils that can be used as biofuels, jet fuels or bio-based products is the absence of a repository for all the gene expression and regulatory information that is being rapidly generated by the community. Here, we provide CamRegBase (https://camregbase.org/) as a one-stop resource to access Camelina information on gene expression and co-expression, transcription factors, lipid associated genes and genome-wide orthologs in the close-relative reference plant Arabidopsis. We envision this as a resource of curated information for users, as well as a repository of new gene regulation information.

Genome-wide identification and functional characterization of the Camelina sativa WRKY gene family in response to abiotic stress

Background

WRKY transcription factors are a superfamily of regulators involved in diverse biological processes and stress responses in plants. However, there is limited knowledge about the WRKY family in camelina (Camelina sativa), an important Brassicaceae oil crop with strong tolerance for various stresses. Here, a genome-wide characterization of WRKY proteins is performed to examine their gene structures, phylogenetics, expression, conserved motif organizations, and functional annotation to identify candidate WRKYs that mediate stress resistance regulation in camelinas.

Results

A total of 242 CsWRKY proteins encoded by 224 gene loci distributed unevenly over the chromosomes were identified, and they were classified into three groups by phylogenetic analysis according to their WRKY domains and zinc finger motifs. The 15 CsWRKY gene loci generated 33 spliced variants. Orthologous WRKY gene pairs were identified, with 173 pairs in the C. sativa and Arabidopsis genomes as well as 282 pairs in the C. sativa and B. napus genomes, respectively. A total of 137 segmental duplication events were observed, but there was no tandem duplication in the camelina genome. Ten major conserved motifs were examined, with WRKYGQK being the most conserved, and several variants were present in many CsWRKYs. Expression analysis revealed that 50% more CsWRKY genes were expressed constitutively, and a set of them displayed tissue-specific expression. Notably, 11 CsWRKY genes exhibited significant expression changes in seedlings under cold, salt, and drought stresses, showing a preferentially inducible expression pattern in response to the stress.

Conclusions

The present article describes a detailed analysis of the CsWRKY gene family and its expression profiles in 12 tissues and under several stress conditions. Segmental duplication is the major force underlying the broad expansion of this gene family, and a strong purifying pressure occurred for CsWRKY proteins during their evolution. CsWRKY proteins play important roles in plant development, with differential functions in different tissues. Exceptionally, eleven CsWRKYs, particularly five alternative spliced isoforms, were found to be the possible key players in mediating plant responses to various stresses. Overall, our results provide a foundation for understanding the roles of CsWRKYs and the precise mechanism through which CsWRKYs regulate high stress resistance as well as the development of stress tolerance cultivars among Cruciferae crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07189-3.

Comparative transcriptome analysis reveals heat stress-responsive genes and their signalling pathways in lilies (Lilium longiflorum vs. Lilium distichum)

The lily, a famous bulbous flower, is seriously affected by high temperatures, which affect their growth and production. To date, the signalling pathways and the molecular mechanisms related to heat response in Lilium have not been elucidated. In this study, a comparative transcriptome analysis was performed in an important thermo-tolerant flower, L. longiflorum, and a thermo-sensitive flower, L. distichum. Lily seedlings were first exposed to heat stress at 42°C for different lengths of time, and the optimal time-points (2 h and 24 h) were selected for RNA sequencing (RNA-seq). Approximately 66.51, 66.21, and 65.36 Mb clean reads were identified from three libraries of L. longiflorum (LL_CK, LL_T2h and LL_T24h, respectively) and 66.18, 66.03, and 65.16 Mb clean reads were obtained from three libraries of L. distichum (LD_CK, LD_T2h and LD_T24h, respectively) after rRNA removing. A total of 34,301 unigenes showed similarity to known proteins in the database NCBI non-redundant protein (NR), Swiss-Prot proteins, InterPro proteins, Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG). In addition, 1,621 genes were differentially expressed in the overlapping libraries between LL_DEGs and LD_DEGs; of these genes, 352 DEGs were obviously upregulated in L. longiflorum and downregulated in L. distichum during heat stress, including 4-coumarate, CoA ligase (4CL), caffeoyl-CoA O-methyltransferase (CCoAOMT), peroxidase, pathogenesis-related protein 10 family genes (PR10s), 14-3-3 protein, leucine-rich repeat receptor-like protein kinase, and glycine-rich cell wall structural protein-like. These genes were mainly involved in metabolic pathways, phenylpropanoid biosynthesis, plant-pathogen interactions, plant hormone signal transduction, and kinase signalling pathways. Quantitative RT-PCR was performed to validate the expression profiling of these DEGs in RNA-seq data. Taken together, the results obtained in the present study provide a comprehensive sequence resource for the discovery of heat-resistance genes and reveal potential key components that are responsive to heat stress in lilies, which may help to elucidate the heat signal transcription networks and facilitate heat-resistance breeding in lily.

Transcriptomic Analysis and Specific Expression of Transcription Factor Genes in the Root and Sporophyll of Dryopteris fragrans (L.) Schott

Background: Dryopteris fragrans, which is densely covered with glandular trichomes, is considered to be one of the ferns with the most medicinal potential. The transcriptomes from selected tissues of D. fragrans were collected and analyzed for functional and comparative genomic studies. The aim of this study was to determine the transcriptomic characteristics of wild D. fragrans sporangium in tissues from the SR (root), SL (sporophyll), and TRL (sporophyll with glandular trichomes removed). Results: Cluster analysis identified genes that were highly expressed in an organ-specific manner according to read mapping, feature counting, and normalization. The functional map identified gene clusters that can uniquely describe the function of each tissue. We identified a group of three tissue-specific transcription factors targeting the SL, SR, and TRL. In addition, highly expressed transcription factors (TFs) were found in each tissue-specific gene cluster, where ERF and bHLH transcription factors were the two types showing the most distinct expression patterns between the three different tissues. The specific expression of transcription factor genes varied between the different types of tissues. The numbers of transcription factors specifically expressed in the roots and sporophylls were 60 and 30, respectively, while only seven were found for the sporophylls with glandular trichomes removed. The expression of genes known to be associated with the development of glandular trichomes in flowering plants, including MIXTA, ATML1, and MYB106, were also validated and are discussed. In particular, a unigene encoding MIXTA was identified and exhibited the highest expression level in SL in D. fragrans. Conclusions: This study is the first report of global transcriptomic analysis in different tissues of D. fragrans, and the first to discuss these findings in the context of the development of homologous glandular trichomes. These results set the stage for further research on the development, stress resistance, and secondary metabolism of D. fragrans glandular trichomes.

Divergent receptor proteins confer responses to different karrikins in two ephemeral weeds

Wildfires can encourage the establishment of invasive plants by releasing potent germination stimulants, such as karrikins. Seed germination of Brassica tournefortii, a noxious weed of Mediterranean climates, is strongly stimulated by KAR₁, the archetypal karrikin produced from burning vegetation. In contrast, the closely-related yet non-fire-associated ephemeral Arabidopsisthaliana is unusual because it responds preferentially to KAR₂. The α/β-hydrolase KARRIKIN INSENSITIVE 2 (KAI2) is the putative karrikin receptor identified in Arabidopsis. Here we show that B. tournefortii expresses three KAI2 homologues, and the most highly-expressed homologue is sufficient to confer enhanced responses to KAR₁ relative to KAR₂ when expressed in Arabidopsis. We identify two amino acid residues near the KAI2 active site that explain the ligand selectivity, and show that this combination has arisen independently multiple times within dicots. Our results suggest that duplication and diversification of KAI2 proteins could confer differential responses to chemical cues produced by environmental disturbance, including fire.

Karrikins are germination stimulants perceived by KAI2 in Arabidopsis. Here the authors show that Brassica tournefortii, a close relative to Arabidopsis, has multiple copies of KAI2 with amino acid substitutions that confer responsiveness to the specific karrikin compounds found in wildfire smoke.

High Flux Through the Oxidative Pentose Phosphate Pathway Lowers Efficiency in Developing Camelina Seeds

Many seeds are green during development, and light has been shown to play a role in the efficiency with which maternally supplied substrates are converted into storage compounds. However, the effects of light on the fluxes through central metabolism that determine this efficiency are poorly understood. Here, we used metabolic flux analysis to determine the effects of light on central metabolism in developing embryos of false flax (Camelina sativa). Metabolic efficiency in C. sativa is of interest because, despite its growing importance as a model oilseed and engineering target and its potential as a biofuel crop, its yields are lower than other major oilseed species. Culture conditions under which steady-state growth and composition of developing embryos match those in planta were used to quantify substrate uptake and respiration rates. The carbon conversion efficiency (CCE) was 21% ± 3% in the dark and 42% ± 4% under high light. Under physiological illumination, the CCE (32% ± 2%) was substantially lower than in green and nongreen oilseeds studied previously. ¹³C and ¹⁴C isotopic labeling experiments were used together with computer-aided modeling to map fluxes through central metabolism. Fluxes through the oxidative pentose phosphate pathway (OPPP) were the principal source of CO₂ production and strongly negatively correlated with CCE across light levels. OPPP fluxes were greatly in excess of demand for NAD(P)H for biosynthesis and larger than those measured in other systems. Excess reductant appears to be dissipated via cyanide-insensitive respiration. OPPP enzymes therefore represent a potential target for increasing efficiency and yield in C. sativa.

De novo transcriptome sequencing of radish (Raphanus sativus L.) fleshy roots: analysis of major genes involved in the anthocyanin synthesis pathway

Background

The HongXin radish (Raphanus sativus L.), which contains the natural red pigment (red radish pigment), is grown in the Fuling district of Chongqing City. However, the molecular mechanisms underlying anthocyanin synthesis for the formation of natural red pigment in the fleshy roots of HongXin radish are not well studied.

Results

De novo transcriptome of HX-1 radish, as well as that of the advanced inbred lines HX-2 and HX-3 were characterized using next generation sequencing (NGS) technology. In total, approximately 66.22 million paired-end reads comprising 34, 927 unigenes (N50 = 1, 621 bp) were obtained. Based on sequence similarity search with known proteins, total of 30, 127 (about 86.26%) unigenes were identified. Additionally, functional annotation and classification of these unigenes indicated that most of the unigenes were predominantly enriched in the metabolic process-related terms, especially for the biosynthetic pathways of secondary metabolites. Moreover, majority of the anthocyanin biosynthesis-related genes (ABRGs) involved in the regulation of anthocyanin biosynthesis were identified by targeted search for their annotation. Subsequently, the expression of 15 putative ABRGs involved in the anthocyanin synthesis-related pathways were validated using quantitative real-time polymerase chain reaction (qRT-PCR). Of those, RsPAL2, RsCHS-B2, RsDFR1, RsDFR2, RsFLS, RsMT3 and RsUFGT73B2-like were identified significantly associated with anthocyanin biosynthesis. Especially for RsDFR1, RsDFR2 and RsFLS, of those, RsDFR1 and RsDFR2 were highest enriched in the HX-3 and WG-3, but RsFLS were down-regulated in HX-3 and WG-3. We proposed that the transcripts of RsDFR1, RsDFR2 and RsFLS might be act as key regulators in anthocyanin biosynthesis pathway.

Conclusions

The assembled radish transcript sequences were analysed to identify the key ABRGs involved in the regulation of anthocyanin biosynthesis. Additionally, the expression patterns of candidate ABRGs involved in the anthocyanin biosynthetic pathway were validated by qRT-PCR. We proposed that the transcripts of RsDFR1, RsDFR2 and RsFLS might be acted as key regulators in anthocyanin biosynthesis pathway. This study will enhance our understanding of the biosynthesis and metabolism of anthocyanin in radish.

Transcriptional regulatory networks controlling taste and aroma quality of apricot (Prunus armeniaca L.) fruit during ripening

BACKGROUND

Taste and aroma, which are important organoleptic qualities of apricot (Prunus armeniaca L.) fruit, undergo rapid and substantial changes during ripening. However, the associated molecular mechanisms remain unclear. The goal of this study was to identify candidate genes for flavor compound metabolism and to construct a regulatory transcriptional network.

RESULTS

We characterized the transcriptome of the 'Jianali' apricot cultivar, which exhibits substantial changes in flavor during ripening, at 50 (turning), 73 (commercial maturation) and 91 (full ripe) days post anthesis (DPA) using RNA sequencing (RNA-Seq). A weighted gene co-expression network analysis (WGCNA) revealed that four of 19 modules correlated highly with flavor compound metabolism (P < 0.001). From them, we identified 1237 differentially expressed genes, with 16 intramodular hubs. A proposed pathway model for flavor compound biosynthesis is presented based on these genes. Two SUS1 genes, as well as SPS2 and INV1 were correlated with sugar biosynthesis, while NADP-ME4, two PK-like and mitochondrial energy metabolism exerted a noticeable effect on organic acid metabolism. CCD1 and FAD2 were identified as being involved in apocarotenoid aroma volatiles and lactone biosynthesis, respectively. Five sugar transporters (Sweet10, STP13, EDR6, STP5.1, STP5.2), one aluminum-activated malate transporter (ALMT9) and one ABCG transporter (ABCG11) were associated with the transport of sugars, organic acids and volatiles, respectively. Sixteen transcription factors were also highlighted that may also play regulatory roles in flavor quality development.

CONCLUSIONS

Apricot RNA-Seq data were obtained and used to generate an annotated set of predicted expressed genes, providing a platform for functional genomic research. Using network analysis and pathway mapping, putative molecular mechanisms for changes in apricot fruit taste and aroma during ripening were elucidated.

De novo transcriptome sequencing and metabolite profiling analyses reveal the complex metabolic genes involved in the terpenoid biosynthesis in Blue Anise Sage (Salvia guaranitica L.)

Many terpenoid compounds have been extracted from different tissues of Salvia guaranitica. However, the molecular genetic basis of terpene biosynthesis pathways is virtually unknown. In this study, approximately 4 Gb of raw data were generated from the transcriptome of S. guaranitica leaves using Illumina HiSeq 2000 sequencing. After filtering and removing the adapter sequences from the raw data, the number of reads reached 32 million, comprising 186 million of high-quality nucleotide bases. A total of 61,400 unigenes were assembled de novo and annotated for establishing a valid database for studying terpenoid biosynthesis. We identified 267 unigenes that are putatively involved in terpenoid metabolism (including, 198 mevalonate and methyl-erythritol phosphate (MEP) pathways, terpenoid backbone biosynthesis genes and 69 terpene synthases genes). Moreover, three terpene synthase genes were studied for their functions in terpenoid biosynthesis by using transgenic Arabidopsis; most transgenic Arabidopsis plants expressing these terpene synthetic genes produced increased amounts of terpenoids compared with wild-type control. The combined data analyses from the transcriptome and metabolome provide new insights into our understanding of the complex metabolic genes in terpenoid-rich blue anise sage, and our study paves the way for the future metabolic engineering of the biosynthesis of useful terpene compounds in S. guaranitica.

Camelina sativa, an oilseed at the nexus between model system and commercial crop

The rapid assessment of metabolic engineering strategies in plants is aided by crops that provide simple, high throughput transformation systems, a sequenced genome, and the ability to evaluate the resulting plants in field trials. Camelina sativa provides all of these attributes in a robust oilseed platform. The ability to perform field evaluation of Camelina is a useful, and in some studies essential benefit that allows researchers to evaluate how traits perform outside the strictly controlled conditions of a greenhouse. In the field the plants are subjected to higher light intensities, seasonal diurnal variations in temperature and light, competition for nutrients, and watering regimes dictated by natural weather patterns, all which may affect trait performance. There are difficulties associated with the use of Camelina. The current genetic resources available for Camelina pale in comparison to those developed for the model plant Arabidopsis thaliana; however, the sequence similarity of the Arabidopsis and Camelina genomes often allows the use of Arabidopsis as a reference when additional information is needed. Camelina's genome, an allohexaploid, is more complex than other model crops, but the diploid inheritance of its three subgenomes is straightforward. The need to navigate three copies of each gene in genome editing or mutagenesis experiments adds some complexity but also provides advantages for gene dosage experiments. The ability to quickly engineer Camelina with novel traits, advance generations, and bulk up homozygous lines for small-scale field tests in less than a year, in our opinion, far outweighs the complexities associated with the crop.

High-level accumulation of oleyl oleate in plant seed oil by abundant supply of oleic acid substrates to efficient wax ester synthesis enzymes

BACKGROUND

Biotechnology enables the production of high-valued industrial feedstocks from plant seed oil. The plant-derived wax esters with long-chain monounsaturated acyl moieties, like oleyl oleate, have favorite properties for lubrication. For biosynthesis of wax esters using acyl-CoA substrates, expressions of a fatty acyl reductase (FAR) and a wax synthase (WS) in seeds are sufficient.

RESULTS

For optimization of the enzymatic activity and subcellular localization of wax ester synthesis enzymes, two fusion proteins were created, which showed wax ester-forming activities in Saccharomyces cerevisiae. To promote the formation of oleyl oleate in seed oil, WSs from Acinetobactor baylyi (AbWSD1) and Marinobacter aquaeolei (MaWS2), as well as the two created fusion proteins were tested in Arabidopsis to evaluate their abilities and substrate preference for wax ester production. The tested seven enzyme combinations resulted in different yields and compositions of wax esters. Expression of a FAR of Marinobacter aquaeolei (MaFAR) with AbWSD1 or MaWS2 led to a high incorporation of C₁₈ substrates in wax esters. The MaFAR/TMMmAWAT2-AbWSD1 combination resulted in the incorporation of more C_18:1 alcohol and C_18:0 acyl moieties into wax esters compared with MaFAR/AbWSD1. The fusion protein of a WS from Simmondsia chinensis (ScWS) with MaFAR exhibited higher specificity toward C_20:1 substrates in preference to C_18:1 substrates. Expression of MaFAR/AbWSD1 in the Arabidopsis fad2 fae1 double mutant resulted in the accumulation of oleyl oleate (18:1/18:1) in up to 62 mol% of total wax esters in seed oil, which was much higher than the 15 mol% reached by MaFAR/AbWSD1 in Arabidopsis Col-0 background. In order to increase the level of oleyl oleate in seed oil of Camelina, lines expressing MaFAR/ScWS were crossed with a transgenic high oleate line. The resulting plants accumulated up to >40 mg g seed^-1 of wax esters, containing 27-34 mol% oleyl oleate.

CONCLUSIONS

The overall yields and the compositions of wax esters can be strongly affected by the availability of acyl-CoA substrates and to a lesser extent, by the characteristics of wax ester synthesis enzymes. For synthesis of oleyl oleate in plant seed oil, appropriate wax ester synthesis enzymes with high catalytic efficiency and desired substrate specificity should be expressed in plant cells; meanwhile, high levels of oleic acid-derived substrates need to be supplied to these enzymes by modifying the fatty acid profile of developing seeds.

De Novo Transcriptome Assembly and Characterization of the Synthesis Genes of Bioactive Constituents in Abelmoschus esculentus (L.) Moench

Abelmoschus esculentus (okra or lady's fingers) is a vegetable with high nutritional value, as well as having certain medicinal effects. It is widely used as food, in the food industry, and in herbal medicinal products, but also as an ornamental, in animal feed, and in other commercial sectors. Okra is rich in bioactive compounds, such as flavonoids, polysaccharides, polyphenols, caffeine, and pectin. In the present study, the concentrations of total flavonoids and polysaccharides in five organs of okra were determined and compared. Transcriptome sequencing was used to explore the biosynthesis pathways associated with the active constituents in okra. Transcriptome sequencing of five organs (roots, stem, leaves, flowers, and fruits) of okra enabled us to obtain 293,971 unigenes, of which 232,490 were annotated. Unigenes related to the enzymes involved in the flavonoid biosynthetic pathway or in fructose and mannose metabolism were identified, based on Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. All of the transcriptional datasets were uploaded to Sequence Read Archive (SRA). In summary, our comprehensive analysis provides important information at the molecular level about the flavonoid and polysaccharide biosynthesis pathways in okra.

Identification of Putative Precursor Genes for the Biosynthesis of Cannabinoid-Like Compound in Radula marginata

The liverwort Radula marginata belongs to the bryophyte division of land plants and is a prospective alternate source of cannabinoid-like compounds. However, mechanistic insights into the molecular pathways directing the synthesis of these cannabinoid-like compounds have been hindered due to the lack of genetic information. This prompted us to do deep sequencing, de novo assembly and annotation of R. marginata transcriptome, which resulted in the identification and validation of the genes for cannabinoid biosynthetic pathway. In total, we have identified 11,421 putative genes encoding 1,554 enzymes from 145 biosynthetic pathways. Interestingly, we have identified all the upstream genes of the central precursor of cannabinoid biosynthesis, cannabigerolic acid (CBGA), including its two first intermediates, stilbene acid (SA) and geranyl diphosphate (GPP). Expression of all these genes was validated using quantitative real-time PCR. We have characterized the protein structure of stilbene synthase (STS), which is considered as a homolog of olivetolic acid in R. marginata. Moreover, the metabolomics approach enabled us to identify CBGA-analogous compounds using electrospray ionization mass spectrometry (ESI-MS/MS) and gas chromatography mass spectrometry (GC-MS). Transcriptomic analysis revealed 1085 transcription factors (TF) from 39 families. Comparative analysis showed that six TF families have been uniquely predicted in R. marginata. In addition, the bioinformatics analysis predicted a large number of simple sequence repeats (SSRs) and non-coding RNAs (ncRNAs). Our results collectively provide mechanistic insights into the putative precursor genes for the biosynthesis of cannabinoid-like compounds and a novel transcriptomic resource for R. marginata. The large-scale transcriptomic resource generated in this study would further serve as a reference transcriptome to explore the Radulaceae family.

Transcriptome and metabolite analyses reveal the complex metabolic genes involved in volatile terpenoid biosynthesis in garden sage (Salvia officinalis)

A large number of terpenoid compounds have been extracted from different tissues of S. officinalis. However, the molecular genetic basis of terpene biosynthesis pathways is virtually unknown. In this study, approximately 6.6 Gb of raw data were generated from the transcriptome of S. officinalis leaves using Illumina HiSeq 2000 sequencing. After filtering and removing the adapter sequences from the raw data, the number of reads reached 21 million, comprising 98 million of high-quality nucleotide bases. 48,671 unigenes were assembled de novo and annotated for establishing a valid database for studying terpenoid biosynthesis. We identified 135 unigenes that are putatively involved in terpenoid metabolism, including 70 mevalonate and methyl-erythritol phosphate pathways, terpenoid backbone biosynthesis genes, and 65 terpene synthase genes. Moreover, five terpene synthase genes were studied for their functions in terpenoid biosynthesis by using transgenic tobacco; most transgenic tobacco plants expressing these terpene synthetic genes produced increased amounts of terpenoids compared with wild-type control. The combined data analyses from the transcriptome and metabolome provide new insights into our understanding of the complex metabolic genes in terpenoid-rich sage, and our study paves the way for the future metabolic engineering of the biosynthesis of useful terpene compounds in S. officinalis.

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.

Characterization of the transcriptome and EST-SSR development in Boea clarkeana, a desiccation-tolerant plant endemic to China

BACKGROUND

Desiccation-tolerant (DT) plants can recover full metabolic competence upon rehydration after losing most of their cellular water (>95%) for extended periods of time. Functional genomic approaches such as transcriptome sequencing can help us understand how DT plants survive and respond to dehydration, which has great significance for plant biology and improving the drought tolerance of crops. Boea clarkeana Hemsl. (Gesneriaceae) is a DT dicotyledonous herb. Its genomic sequences characteristics remain unknown. Based on transcriptomic analyses, polymorphic EST-SSR (simple sequence repeats in expressed sequence tags) molecular primers can be designed, which will greatly facilitate further investigations of the population genetics and demographic histories of DT plants.

METHODS

In the present study, we used the platform Illumina HiSeq™2000 and de novo assembly technology to obtain leaf transcriptomes of B. clarkeana and conducted a BLASTX alignment of the sequencing data and protein databases for sequence classification and annotation. Then, based on the sequence information, the EST-SSR markers were developed, and the functional annotation of ESTs containing polymorphic SSRs were obtained through BLASTX.

RESULTS

A total of 91,449 unigenes were generated from the leaf cDNA library of B. clarkeana. Based on a sequence similarity search with a known protein database, 72,087 unigenes were annotated. Among the annotated unigenes, a total of 71,170 unigenes showed significant similarity to the known proteins of 463 popular model species in the Nr database, and 59,962 unigenes and 32,336 unigenes were assigned to Gene Ontology (GO) classifications and Cluster of Orthologous Groups (COG), respectively. In addition, 44,924 unigenes were mapped in 128 KEGG pathways. Furthermore, a total of 7,610 unigenes with 8,563 microsatellites were found. Seventy-four primer pairs were selected from 436 primer pairs designed for polymorphism validation. SSRs with higher polymorphism rates were concentrated on dinucleotides, pentanucleotides and hexanucleotides. Finally, 17 pairs with stable, highly polymorphic loci were selected for polymorphism screening. There was a total of 65 alleles, with 2-6 alleles at each locus. Primarily due to the unique biological characteristics of plants, the HE (0-0.196), HO (0.082-0.14) and PIC (0-0.155) per locus were very low. The functional annotation distribution centered on ESTs containing di- and tri-nucleotide SSRs, and the ESTs containing primers BC2, BC4 and BC12 were annotated to vegetative dehydration/desiccation pathways.

DISCUSSION

This work is the first genetic study of B. clarkeana as a new plant resource of DT genes. A substantial number of transcriptome sequences were generated in this study. These sequences are valuable resources for gene annotation and discovery as well as molecular marker development. These sequences could also provide a valuable basis for future molecular studies of B. clarkeana.

Two Acyltransferases Contribute Differently to Linolenic Acid Levels in Seed Oil

Acyltransferases are key contributors to triacylglycerol (TAG) synthesis and, thus, are of great importance for seed oil quality. The effects of increased or decreased expression of ACYL-COENZYME A:DIACYLGLYCEROL ACYLTRANSFERASE1 (DGAT1) or PHOSPHOLIPID:DIACYLGLYCEROL ACYLTRANSFERASE (PDAT) on seed lipid composition were assessed in several Camelina sativa lines. Furthermore, in vitro assays of acyltransferases in microsomal fractions prepared from developing seeds of some of these lines were performed. Decreased expression of DGAT1 led to an increased percentage of 18:3n-3 without any change in total lipid content of the seed. The tri-18:3 TAG increase occurred predominantly in the cotyledon, as determined with matrix-assisted laser desorption/ionization-mass spectrometry, whereas species with two 18:3n-3 acyl groups were elevated in both cotyledon and embryonal axis. PDAT overexpression led to a relative increase of 18:2n-6 at the expense of 18:3n-3, also without affecting the total lipid content. Differential distributions of TAG species also were observed in different parts of the seed. The microsomal assays revealed that C.sativa seeds have very high activity of diacylglycerol-phosphatidylcholine interconversion. The combination of analytical and biochemical data suggests that the higher 18:2n-6 content in the seed oil of the PDAT overexpressors is due to the channeling of fatty acids from phosphatidylcholine into TAG before being desaturated to 18:3n-3, caused by the high activity of PDAT in general and by PDAT specificity for 18:2n-6. The higher levels of 18:3n-3 in DGAT1-silencing lines are likely due to the compensatory activity of a TAG-synthesizing enzyme with specificity for this acyl group and more desaturation of acyl groups occurring on phosphatidylcholine.

Genome-Wide Transcriptional Profiling to Elucidate Key Candidates Involved in Bud Burst and Rattling Growth in a Subtropical Bamboo (Dendrocalamus hamiltonii)

Bamboo, one of the fastest growing plants, can be a promising model system to understand growth. The study provides an insight into the complex interplay between environmental signaling and cellular machineries governing initiation and persistence of growth in a subtropical bamboo (Dendrocalamus hamiltonii). Phenological and spatio-temporal transcriptome analysis of rhizome and shoot during the major vegetative developmental transitions of D. hamiltonii was performed to dissect factors governing growth. Our work signifies the role of environmental cues, predominantly rainfall, decreasing day length, and high humidity for activating dormant bud to produce new shoot, possibly through complex molecular interactions among phosphatidylinositol, calcium signaling pathways, phytohormones, circadian rhythm, and humidity responses. We found the coordinated regulation of auxin, cytokinin, brassinosteroid signaling and cell cycle modulators; facilitating cell proliferation, cell expansion, and cell wall biogenesis supporting persistent growth of emerging shoot. Putative master regulators among these candidates were identified using predetermined Arabidopsis thaliana protein-protein interaction network. We got clues that the growth signaling begins far back in rhizome even before it emerges out as new shoot. Putative growth candidates identified in our study can serve in devising strategies to engineer bamboos and timber trees with enhanced growth and biomass potentials.

The developmental transcriptome atlas of the biofuel crop Camelina sativa

Camelina sativa is currently being embraced as a viable industrial bio-platform crop due to a number of desirable agronomic attributes and the unique fatty acid profile of the seed oil that has applications for food, feed and biofuel. The recent completion of the reference genome sequence of C. sativa identified a young hexaploid genome. To complement this work, we have generated a genome-wide developmental transcriptome map by RNA sequencing of 12 different tissues covering major developmental stages during the life cycle of C. sativa. We have generated a digital atlas of this comprehensive transcriptome resource that enables interactive visualization of expression data through a searchable database of electronic fluorescent pictographs (eFP browser). An analysis of this dataset supported expression of 88% of the annotated genes in C. sativa and provided a global overview of the complex architecture of temporal and spatial gene expression patterns active during development. Conventional differential gene expression analysis combined with weighted gene expression network analysis uncovered similarities as well as differences in gene expression patterns between different tissues and identified tissue-specific genes and network modules. A high-quality census of transcription factors, analysis of alternative splicing and tissue-specific genome dominance provided insight into the transcriptional dynamics and sub-genome interplay among the well-preserved triplicated repertoire of homeologous loci. The comprehensive transcriptome atlas in combination with the reference genome sequence provides a powerful resource for genomics research which can be leveraged to identify functional associations between genes and understand the regulatory networks underlying developmental processes.

Transcriptome Analysis of the Signalling Networks in Coronatine-Induced Secondary Laticifer Differentiation from Vascular Cambia in Rubber Trees

The secondary laticifer in rubber tree (Hevea brasiliensis Muell. Arg.) is a specific tissue within the secondary phloem. This tissue differentiates from the vascular cambia, and its function is natural rubber biosynthesis and storage. Given that jasmonates play a pivotal role in secondary laticifer differentiation, we established an experimental system with jasmonate (JA) mimic coronatine (COR) for studying the secondary laticifer differentiation: in this system, differentiation occurs within five days of the treatment of epicormic shoots with COR. In the present study, the experimental system was used to perform transcriptome sequencing and gene expression analysis. A total of 67,873 unigenes were assembled, and 50,548 unigenes were mapped at least in one public database. Of these being annotated unigenes, 15,780 unigenes were differentially expressed early after COR treatment, and 19,824 unigenes were differentially expressed late after COR treatment. At the early stage, 8,646 unigenes were up-regulated, while 7,134 unigenes were down-regulated. At the late stage, the numbers of up- and down-regulated unigenes were 7,711 and 12,113, respectively. The annotation data and gene expression analysis of the differentially expressed unigenes suggest that JA-mediated signalling, Ca²⁺ signal transduction and the CLAVATA-MAPK-WOX signalling pathway may be involved in regulating secondary laticifer differentiation in rubber trees.

Transcriptome Analysis of the Response to NaCl in Suaeda maritima Provides an Insight into Salt Tolerance Mechanisms in Halophytes

Although salt tolerance is a feature representative of halophytes, most studies on this topic in plants have been conducted on glycophytes. Transcriptome profiles are also available for only a limited number of halophytes. Hence, the present study was conducted to understand the molecular basis of salt tolerance through the transcriptome profiling of the halophyte Suaeda maritima, which is an emerging plant model for research on salt tolerance. Illumina sequencing revealed 72,588 clustered transcripts, including 27,434 that were annotated using BLASTX. Salt application resulted in the 2-fold or greater upregulation of 647 genes and downregulation of 735 genes. Of these, 391 proteins were homologous to proteins in the COGs (cluster of orthologous groups) database, and the majorities were grouped into the poorly characterized category. Approximately 50% of the genes assigned to MapMan pathways showed homology to S. maritima. The majority of such genes represented transcription factors. Several genes also contributed to cell wall and carbohydrate metabolism, ion relation, redox responses and G protein, phosphoinositide and hormone signaling. Real-time PCR was used to validate the results of the deep sequencing for the most of the genes. This study demonstrates the expression of protein kinase C, the target of diacylglycerol in phosphoinositide signaling, for the first time in plants. This study further reveals that the biochemical and molecular responses occurring at several levels are associated with salt tolerance in S. maritima. At the structural level, adaptations to high salinity levels include the remodeling of cell walls and the modification of membrane lipids. At the cellular level, the accumulation of glycinebetaine and the sequestration and exclusion of Na+ appear to be important. Moreover, this study also shows that the processes related to salt tolerance might be highly complex, as reflected by the salt-induced enhancement of transcription factor expression, including hormone-responsive factors, and that this process might be initially triggered by G protein and phosphoinositide signaling.

De novo transcriptome analysis in radish (Raphanus sativus L.) and identification of critical genes involved in bolting and flowering

Background

The appropriate timing of bolting and flowering is pivotal for reproductive success in Brassicaceae crops including radish (Raphanus sativus L.). Although several flowering regulatory pathways had been described in some plant species, no study on genetic networks of bolting and flowering regulation was performed in radish. In this study, to generate dataset of radish unigene sequences for large-scale gene discovery and functional pathway identification, a cDNA library from mixed radish leaves at different developmental stages was subjected to high-throughput RNA sequencing (RNA-seq).

Results

A total of 54.64 million clean reads and 111,167 contigs representing 53,642 unigenes were obtained from the radish leaf transcriptome. Among these, 50,385 unigenes were successfully annotated by BLAST searching against the public protein databases. Functional classification and annotation indicated that 42,903 and 15,382 unique sequences were assigned to 55 GO terms and 25 COG categories, respectively. KEGG pathway analysis revealed that 25,973 unigenes were classified into 128 functional pathways, among which 24 candidate genes related to plant circadian rhythm were identified. Moreover, 142 potential bolting and flowering-related genes involved in various flowering pathways were identified. In addition, seven critical bolting and flowering-related genes were isolated and profiled by T-A cloning and RT-qPCR analysis. Finally, a schematic network model of bolting and flowering regulation and pathways was put forward in radish.

Conclusions

This study is the first report on systematic identification of bolting and flowering-related genes based on transcriptome sequencing and assembly in radish. These results could provide a foundation for further investigating bolting and flowering regulatory networks in radish, and facilitate dissecting molecular genetic mechanisms underlying bolting and flowering in Brassicaceae vegetable crops.

Electronic supplementary material

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

Identification of multiple lipid genes with modifications in expression and sequence associated with the evolution of hydroxy fatty acid accumulation in Physaria fendleri

Two Brassicaceae species, Physaria fendleri and Camelina sativa, are genetically very closely related to each other and to Arabidopsis thaliana. Physaria fendleri seeds contain over 50% hydroxy fatty acids (HFAs), while Camelina sativa and Arabidopsis do not accumulate HFAs. To better understand how plants evolved new biochemical pathways with the capacity to accumulate high levels of unusual fatty acids, transcript expression and protein sequences of developing seeds of Physaria fendleri, wild-type Camelina sativa, and Camelina sativa expressing a castor bean (Ricinus communis) hydroxylase were analyzed. A number of potential evolutionary adaptations within lipid metabolism that probably enhance HFA production and accumulation in Physaria fendleri, and, in their absence, limit accumulation in transgenic tissues were revealed. These adaptations occurred in at least 20 genes within several lipid pathways from the onset of fatty acid synthesis and its regulation to the assembly of triacylglycerols. Lipid genes of Physaria fendleri appear to have co-evolved through modulation of transcriptional abundances and alterations within protein sequences. Only a handful of genes showed evidence for sequence adaptation through gene duplication. Collectively, these evolutionary changes probably occurred to minimize deleterious effects of high HFA amounts and/or to enhance accumulation for physiological advantage. These results shed light on the evolution of pathways for novel fatty acid production in seeds, help explain some of the current limitations to accumulation of HFAs in transgenic plants, and may provide improved strategies for future engineering of their production.

The role of melanin pathways in extremotolerance and virulence of Fonsecaea revealed by de novo assembly transcriptomics using illumina paired-end sequencing

Melanisation has been considered to be an important virulence factor of Fonsecaea monophora. However, the biosynthetic mechanisms of melanisation remain unknown. We therefore used next generation sequencing technology to investigate the transcriptome and digital gene expression data, which are valuable resources to better understand the molecular and biological mechanisms regulating melanisation in F. monophora. We performed de novo transcriptome assembly and digital gene expression (DGE) profiling analyses of parent (CBS 122845) and albino (CBS 125194) strains using the Illumina RNA-seq system. A total of 17 352 annotated unigenes were found by BLAST search of NR, Swiss-Prot, Gene Ontology, Clusters of Orthologous Groups and Kyoto Encyclopedia of Genes and Genomes (KEGG) (E-value <1e‒5). A total of 2 283 unigenes were judged to be the differentially expressed between the two genotypes. We identified most of the genes coding for key enzymes involved in melanin biosynthesis pathways, including polyketide synthase (pks), multicopper oxidase (mco), laccase, tyrosinase and homogentisate 1,2-dioxygenase (hmgA). DEG analysis showed extensive down-regulation of key genes in the DHN pathway, while up-regulation was noted in the DOPA pathway of the albino mutant. The transcript levels of partial genes were confirmed by real time RT-PCR, while the crucial role of key enzymes was confirmed by either inhibitor or substrate tests in vitro. Meanwhile, numbers of genes involved in light sensing, cell wall synthesis, morphology and environmental stress were identified in the transcriptome of F. monophora. In addition, 3 353 SSRs (Simple Sequence Repeats) markers were identified from 21 600 consensus sequences. Blocking of the DNH pathway is the most likely reason of melanin deficiency in the albino strain, while the production of pheomelanin and pyomelanin were probably regulated by unknown transcription factors on upstream of both pathways. Most of genes involved in environmental tolerance to oxidants, irradiation and extreme temperatures were also assembled and annotated in transcriptomes of F. monophora. In addition, thousands of identified cSSR (combined SSR) markers will favour further genetic linkage studies. In conclusion, these data will contribute to understanding the regulation of melanin biosynthesis and help to improve the studies of pathogenicity of F. monophora.

De Novo Transcriptome Sequencing of Olea europaea L. to Identify Genes Involved in the Development of the Pollen Tube

In olive (Olea europaea L.), the processes controlling self-incompatibility are still unclear and the molecular basis underlying this process are still not fully characterized. In order to determine compatibility relationships, using next-generation sequencing techniques and a de novo transcriptome assembly strategy, we show that pollen tubes from different olive plants, grown in vitro in a medium containing its own pistil and in combination pollen/pistil from self-sterile and self-fertile cultivars, have a distinct gene expression profile and many of the differentially expressed sequences between the samples fall within gene families involved in the development of the pollen tube, such as lipase, carboxylesterase, pectinesterase, pectin methylesterase, and callose synthase. Moreover, different genes involved in signal transduction, transcription, and growth are overrepresented. The analysis also allowed us to identify members in actin and actin depolymerization factor and fibrin gene family and member of the Ca(2+) binding gene family related to the development and polarization of pollen apical tip. The whole transcriptomic analysis, through the identification of the differentially expressed transcripts set and an extended functional annotation analysis, will lead to a better understanding of the mechanisms of pollen germination and pollen tube growth in the olive.

Induction of Viable but Nonculturable State in Rhodococcus and Transcriptome Analysis Using RNA-seq

Viable but nonculturable (VBNC) bacteria, which maintain the viability with loss of culturability, universally exist in contaminated and non-contaminated environments. In this study, two strains, Rhodococcus sp. TG13 and TN3, which were isolated from PCB-contaminated sediment and non-contaminated sediment respectively, were investigated under low temperature and oligotrophic conditions. The results indicated that the two strains TG13 and TN3 could enter into the VBNC state with different incubation times, and could recover culturability by reversal of unfavourable factors and addition of resuscitation-promoting factor (Rpf), respectively. Furthermore, the gene expression variations in the VBNC response were clarified by Illumina high throughput RNA-sequencing. Genome-wide transcriptional analysis demonstrated that up-regulated genes in the VBNC cells of the strain TG13 related to protein modification, ATP accumulation and RNA polymerase, while all differentially expressed genes (DEGs) in the VBNC cells of the strain TN3 were down-regulated. However, the down-regulated genes in both the two strains mainly encoded NADH dehydrogenase subunit, catalase, oxidoreductase, which further verified that cold-induced loss of ability to defend oxidative stress may play an important role in induction of the VBNC state. This study further verified that the molecular mechanisms underlying the VBNC state varied with various bacterial species. Study on the VBNC state of non-pathogenic bacteria will provide new insights into the limitation of environmental micro-bioremediation and the cultivation of unculturable species.

A De novo Transcriptomic Approach to Identify Flavonoids and Anthocyanins "Switch-Off" in Olive (Olea europaea L.) Drupes at Different Stages of Maturation

Highlights A de novo transcriptome reconstruction of olive drupes was performed in two genotypesGene expression was monitored during drupe development in two olive cultivarsTranscripts involved in flavonoid and anthocyanin pathways were analyzed in Cassanese and Leucocarpa cultivarsBoth cultivar and developmental stage impact gene expression in Olea europaea fruits. During ripening, the fruits of the olive tree (Olea europaea L.) undergo a progressive chromatic change characterized by the formation of a red-brown "spot" which gradually extends on the epidermis and in the innermost part of the mesocarp. This event finds an exception in the Leucocarpa cultivar, in which we observe a destabilized equilibrium between the metabolisms of chlorophyll and other pigments, particularly the anthocyanins whose switch-off during maturation promotes the white coloration of fruits. Despite its importance, genomic information on the olive tree is still lacking. Different RNA-seq libraries were generated from drupes of "Leucocarpa" and "Cassanese" olive genotypes, sampled at 100 and 130 days after flowering (DAF), and were used in order to identify transcripts involved in the main phenotypic changes of fruits during maturation and their corresponding expression patterns. A total of 103,359 transcripts were obtained and 3792 and 3064 were differentially expressed in "Leucocarpa" and "Cassanese" genotypes, respectively, during 100-130 DAF transition. Among them flavonoid and anthocyanin related transcripts such as phenylalanine ammonia lyase (PAL), cinnamate 4-hydroxylase (C4H), 4-coumarate-CoA ligase (4CL), chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonol 3'-hydrogenase (F3'H), flavonol 3'5 '-hydrogenase (F3'5'H), flavonol synthase (FLS), dihydroflavonol 4-reductase (DFR), anthocyanidin synthase (ANS), UDP-glucose:anthocianidin: flavonoid glucosyltransferase (UFGT) were identified. These results contribute to reducing the current gap in information regarding metabolic processes, including those linked to fruit pigmentation in the olive.

Transcriptome profiling of Camelina sativa to identify genes involved in triacylglycerol biosynthesis and accumulation in the developing seeds

BACKGROUND

Camelina sativa is an emerging dedicated oilseed crop designed for biofuel and biodiesel applications as well as a source for edible and general-purpose oils. Such valuable oilseed crop is subjected to plant breeding programs and is suggested for large-scale production of better seed and oil quality. To accomplish this objective and to further enhance its oil content, a better understanding of lipid metabolism at the molecular level in this plant is critical. Here, we applied tissue transcriptomics and lipid composition analysis to identify and profile the genes and gene networks associated with triacylglycerol (TAG) biosynthesis, and to investigate how those genes are interacting to determine the quantity and quality of Camelina oil during seed development.

RESULTS

Our Camelina transcriptome data analysis revealed an approximate of 57,854 and 57,973 genes actively expressing in developing seeds (RPKM ≥ 0.1) at 10-15 (Cs-14) and 16-21 (Cs-21) days after flowering (DAF), respectively. Of these, 7932 genes showed temporal and differential gene expression during the seed development (log2 fold change ≥1.5 or ≤-1.5; P ≤ 0.05). The differentially expressed genes (DEGs) were annotated and were found to be involved in distinct functional categories and metabolic pathways. Furthermore, performing quantitative real-time PCR for selected candidate genes associated with TAG biosynthesis validated RNA-seq data. Our results showed strong positive correlations between the expression abundance measured using both qPCR and RNA-Seq technologies. Furthermore, the analysis of fatty-acid content and composition revealed major changes throughout seed development, with the amount of oil accumulate rapidly at early mid seed development stages (from 16-28 DAF onwards), while no important changes were observed in the fatty-acid profile between seeds at 28 DAF and mature seeds.

CONCLUSIONS

This study is highly useful for understanding the regulation of TAG biosynthesis and identifying the rate-limiting steps in TAG pathways at seed development stages, providing a precise selection of candidate genes for developing Camelina varieties with improved seed and oil yields.

Synthesis of oleyl oleate wax esters in Arabidopsis thaliana and Camelina sativa seed oil

Seed oil composed of wax esters with long-chain monoenoic acyl moieties represents a high-value commodity for industry. Such plant-derived sperm oil-like liquid wax esters are biodegradable and can have excellent properties for lubrication. In addition, wax ester oil may represent a superior substrate for biodiesel production. In this study, we demonstrate that the low-input oil seed crop Camelina sativa can serve as a biotechnological platform for environmentally benign wax ester production. Two biosynthetic steps catalysed by a fatty alcohol-forming acyl-CoA reductase (FAR) and a wax ester synthase (WS) are sufficient to achieve wax ester accumulation from acyl-CoA substrates. To produce plant-derived sperm oil-like liquid wax esters, the WS from Mus musculus (MmWS) or Simmondsia chinensis (ScWS) were expressed in combination with the FAR from Mus musculus (MmFAR1) or Marinobacter aquaeolei (MaFAR) in seeds of Arabidopsis thaliana and Camelina sativa. The three analysed enzyme combinations Oleo3:mCherry:MmFAR1∆c/Oleo3:EYFP:MmWS, Oleo3:mCherry:MmFAR1∆c/ScWS and MaFAR/ScWS showed differences in the wax ester molecular species profiles and overall biosynthetic performance. By expressing MaFAR/ScWS in Arabidopsis or Camelina up to 59% or 21% of the seed oil TAGs were replaced by wax esters, respectively. This combination also yielded wax ester molecular species with highest content of monounsaturated acyl moieties. Expression of the enzyme combinations in the Arabidopsis fae1 fad2 mutant background high in oleic acid resulted in wax ester accumulation enriched in oleyl oleate (18:1/18:1 > 60%), suggesting that similar values may be obtained with a Camelina high oleic acid line.

De-novo assembly and characterization of Chlorella minutissima UTEX2341 transcriptome by paired-end sequencing and the identification of genes related to the biosynthesis of lipids for biodiesel

Chlorella minutissima is considered to be one of the promising feedstocks for biofuels in the future. In this study, the transcriptome from the oil-rich strain UTEX2341 of C. minutissima was generated based on Illumina paired-end sequencing. Through de-novo assembly, a total of 14,905 isogenes were obtained and compacted into 6216 unigenes. A total of 80% of the unigenes were assigned with GO terms and were further subdivided into 55 sub-categories. KEGG analysis demonstrated that 37.2% of the unigenes could be accessed and mapped into 278 pathways. Interestingly, the genes that encoded key enzymes that are involved in the biosynthesis, elongation, and metabolism of fatty acids were identified, including malonyl-CoA-ACP transacylase, 3-ketoacyl-ACP synthase, 3-ketoacyl-ACP reductase, and others. Moreover, the genes that are involved in triacylglycerol (TAG) biosynthesis and metabolism were also observed. Therefore, the transcriptome analysis of C. minutissima UTEX2341 not only supplies comprehensive insight into the molecular pathway that is involved in the biosynthesis of biofuel precursors but also provides substantial valuable genomic resources to accelerate the further development and utilization of biofuels.

Lipid labeling from acetate or glycerol in cultured embryos of Camelina sativa seeds: A tale of two substrates

Studies on the metabolism of lipids in seeds frequently use radiolabeled acetate and glycerol supplied to excised developing seeds to track the biosynthesis of acyl and lipid head groups, respectively. Such experiments are generally restricted to shorter time periods and the results may not quantitatively reflect in planta rates. These limitations can be removed by using cultured embryos, provided they mimic growth and lipid deposition observed for embryos in planta. Mid-maturation embryos from Camelina sativa were cultured in vitro to assess the use of sufficient acetate or glycerol concentrations and labeling periods for stable isotope labeling and mass spectrometric detection. Maximum incorporation of exogenous acetate into fatty acids occurred at 1mM and above. This provides about 5% of the total carbon flux entering fatty acids, enough for (13)C isotopomer analysis while maintaining normal biosynthetic rates for over 24h. Labeling analysis indicates that acetate reports lipid metabolism uniformly across the embryo. At higher acetate concentrations with longer incubations, the rate of fatty acid synthesis is reduced and the composition of newly synthesized fatty acids changes. While the mole fractions of oleate that undergo Δ12-desaturation or elongation are independent of biosynthetic flux, Δ15-desaturation shows a bimodal dependence. These observations are consistent with changes occurring in planta over seed development. Incorporation rates of the glyceryl moiety into lipids saturates at about 0.5mM exogenous glycerol. At saturation, the exogenous glycerol almost completely replaces the endogenous supply of glycerol-3-phosphate without affecting net lipid accumulation or fatty acid composition. It is concluded that acetate and glycerol labeling of cultured C. sativa embryos can provide an accurate representation of lipid metabolism in embryos in vivo, and that in Camelina embryos glycerol-3-phosphate levels do not co-limit triacylglycerol synthesis.

Lipid analysis of developing Camelina sativa seeds and cultured embryos

Camelina sativa is a cultivated oilseed rich in triacylglycerols containing oleic, linoleic, α-linolenic and eicosenoic acids. As it holds promise as a model species, its lipid synthesis was characterized in vivo and in culture. Lipid accumulates at a maximum rate of about 26 μg/day/seed (11.5 mg lipid/day/g fresh seed weight), a rate comparable with other oilseeds. Noteworthy is a late stage surge in α-linolenic acid accumulation. Small amounts of unusual epoxy and hydroxy fatty acids are also present in the triacylglycerols. These include 15,16-epoxy- and 15-hydroxy-octadecadienoic acids and homologous series of ω7-hydroxy-alk-ω9-enoic and ω9/10-hydroxy-alkanoic acids. Mid-maturation embryos cultured in vitro have growth and lipid deposition rates and fatty acid compositions that closely match that of seeds, but extended culture periods allow these rates to rise and surpass those observed in planta. Optimized thin layer chromatography systems for characterization of labeled products from acetate or glycerol labeling are described. Glycerol label is only found in acylglycerols, largely as the intact glyceryl backbone, but acetate can label acyl groups and sterols, the latter to a much higher relative specific activity. This presumably occurs because mevalonic acid precursor is derived from the non-plastid pool of acetyl-CoA that is also the source for malonyl-CoA to drive FAE1-dependent chain elongation. Particular attention has been paid to the separation of sterols and diacylglycerols, and to hydrogenation of triacylglycerols to simplify their analysis. These improved methods will allow more accurate analyses of the fluxes of lipid metabolism in cultured plant embryos.

De novo assembly of the Japanese lawngrass (Zoysia japonica Steud.) root transcriptome and identification of candidate unigenes related to early responses under salt stress

Japanese lawngrass (Zoysia japonica Steud.) is an important warm-season turfgrass that is able to survive in a range of soils, from infertile sands to clays, and to grow well under saline conditions. However, little is known about the molecular mechanisms involved in its resistance to salt stress. Here, we used high-throughput RNA sequencing (RNA-seq) to investigate the changes in gene expression of Zoysia grass at high NaCl concentrations. We first constructed two sequencing libraries, including control and NaCl-treated samples, and sequenced them using the Illumina HiSeq™ 2000 platform. Approximately 157.20 million paired-end reads with a total length of 68.68 Mb were obtained. Subsequently, 100,800 unigenes with an N50 length of 1104 bp were assembled using Trinity, among which 70,127 unigenes were functionally annotated (E ≤ 10(-5)) in the non-redundant protein (NR) database. Furthermore, three public databases, the Kyoto Encyclopedia of Genes and Genomes (KEGG), Swiss-prot, and Clusters of Orthologous Groups (COGs), were used for gene function analysis and enrichment. The annotated genes included 46 Gene Ontology (GO) terms, 120 KEGG pathways, and 25 COGs. Compared with the control, 6035 genes were significantly different (false discovery rate ≤0.01, |log2Ratio|≥1) in the NaCl-treated samples. These genes were enriched in 10 KEGG pathways and 58 GO terms, and subjected to 25 COG categories. Using high-throughput next-generation sequencing, we built a database as a global transcript resource for Z. japonica Steud. roots. The results of this study will advance our understanding of the early salt response in Japanese lawngrass roots.

De novo assembly and characterisation of the field pea transcriptome using RNA-Seq

BACKGROUND

Field pea (Pisum sativum L.) is a cool-season grain legume that is cultivated world-wide for both human consumption and stock-feed purposes. Enhancement of genetic and genomic resources for field pea will permit improved understanding of the control of traits relevant to crop productivity and quality. Advances in second-generation sequencing and associated bioinformatics analysis now provide unprecedented opportunities for the development of such resources. The objective of this study was to perform transcriptome sequencing and characterisation from two genotypes of field pea that differ in terms of seed and plant morphological characteristics.

RESULTS

Transcriptome sequencing was performed with RNA templates from multiple tissues of the field pea genotypes Kaspa and Parafield. Tissue samples were collected at various growth stages, and a total of 23 cDNA libraries were sequenced using Illumina high-throughput sequencing platforms. A total of 407 and 352 million paired-end reads from the Kaspa and Parafield transcriptomes, respectively were assembled into 129,282 and 149,272 contigs, which were filtered on the basis of known gene annotations, presence of open reading frames (ORFs), reciprocal matches and degree of coverage. Totals of 126,335 contigs from Kaspa and 145,730 from Parafield were subsequently selected as the reference set. Reciprocal sequence analysis revealed that c. 87% of contigs were expressed in both cultivars, while a small proportion were unique to each genotype. Reads from different libraries were aligned to the genotype-specific assemblies in order to identify and characterise expression of contigs on a tissue-specific basis, of which 87% were expressed in more than one tissue, while others showed distinct expression patterns in specific tissues, providing unique transcriptome signatures.

CONCLUSION

This study provided a comprehensive assembled and annotated transcriptome set for field pea that can be used for development of genetic markers, in order to assess genetic diversity, construct linkage maps, perform trait-dissection and implement whole-genome selection strategies in varietal improvement programs, as well to identify target genes for genetic modification approaches on the basis of annotation and expression analysis. In addition, the reference field pea transcriptome will prove highly valuable for comparative genomics studies and construction of a finalised genome sequence.

Mining and identification of polyunsaturated fatty acid synthesis genes active during camelina seed development using 454 pyrosequencing

BACKGROUND

Camelina (Camelina sativa L.) is well known for its high unsaturated fatty acid content and great resistance to environmental stress. However, little is known about the molecular mechanisms of unsaturated fatty acid biosynthesis in this annual oilseed crop. To gain greater insight into this mechanism, the transcriptome profiles of seeds at different developmental stages were analyzed by 454 pyrosequencing.

RESULTS

Sequencing of two normalized 454 libraries produced 831,632 clean reads. A total of 32,759 unigenes with an average length of 642 bp were obtained by de novo assembly, and 12,476 up-regulated and 12,390 down-regulated unigenes were identified in the 20 DAF (days after flowering) library compared with the 10 DAF library. Functional annotations showed that 220 genes annotated as fatty acid biosynthesis genes were up-regulated in 20 DAF sample. Among them, 47 candidate unigenes were characterized as responsible for polyunsaturated fatty acid synthesis. To verify unigene expression levels calculated from the transcriptome analysis results, quantitative real-time PCR was performed on 11 randomly selected genes from the 220 up-regulated genes; 10 showed consistency between qRT-PCR and 454 pyrosequencing results.

CONCLUSIONS

Investigation of gene expression levels revealed 32,759 genes involved in seed development, many of which showed significant changes in the 20 DAF sample compared with the 10 DAF sample. Our 454 pyrosequencing data for the camelina transcriptome provide an insight into the molecular mechanisms and regulatory pathways of polyunsaturated fatty acid biosynthesis in camelina. The genes characterized in our research will provide candidate genes for the genetic modification of crops.

De novo assembly and transcriptomic profiling of the grazing response in Stipa grandis

Background

Stipa grandis (Poaceae) is one of the dominant species in a typical steppe of the Inner Mongolian Plateau. However, primarily due to heavy grazing, the grasslands have become seriously degraded, and S. grandis has developed a special growth-inhibition phenotype against the stressful habitat. Because of the lack of transcriptomic and genomic information, the understanding of the molecular mechanisms underlying the grazing response of S. grandis has been prohibited.

Results

Using the Illumina HiSeq 2000 platform, two libraries prepared from non-grazing (FS) and overgrazing samples (OS) were sequenced. De novo assembly produced 94,674 unigenes, of which 65,047 unigenes had BLAST hits in the National Center for Biotechnology Information (NCBI) non-redundant (nr) database (E-value < 10^-5). In total, 47,747, 26,156 and 40,842 unigenes were assigned to the Gene Ontology (GO), Clusters of Orthologous Group (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. A total of 13,221 unigenes showed significant differences in expression under the overgrazing condition, with a threshold false discovery rate ≤ 0.001 and an absolute value of log₂Ratio ≥ 1. These differentially expressed genes (DEGs) were assigned to 43,257 GO terms and were significantly enriched in 32 KEGG pathways (q-value ≤ 0.05). The alterations in the wound-, drought- and defense-related genes indicate that stressors have an additive effect on the growth inhibition of this species.

Conclusions

This first large-scale transcriptome study will provide important information for further gene expression and functional genomics studies, and it facilitated our investigation of the molecular mechanisms of the S. grandis grazing response and the associated morphological and physiological characteristics.

Identification of hydroxy fatty acid and triacylglycerol metabolism-related genes in lesquerella through seed transcriptome analysis

Background

Castor oil is the only commercial source of hydroxy fatty acid that has industrial value. The production of castor oil is hampered by the presence of the toxin ricin in its seed. Lesquerella seed also accumulates hydroxy fatty acid and is free of ricin, and thus it is being developed as a new crop for hydroxy fatty acid production. A high-throughput, large-scale sequencing of transcripts from developing lesquerella seeds was carried out by 454 pyrosequencing to generate a database for quality improvement of seed oil and other agronomic traits. Deep mining and characterization of acyl-lipid genes were conducted to uncover candidate genes for further studies of mechanisms underlying hydroxy fatty acid and seed oil synthesis.

Results

A total of 651 megabases of raw sequences from an mRNA sample of developing seeds was acquired. Bioinformatic analysis of these sequences revealed 59,914 transcripts representing 26,995 unique genes that include nearly all known seed expressed genes. Based on sequence similarity with known plant proteins, about 74% (19,861) genes matched with annotated coding genes. Among them, 95% (18,868) showed highest sequence homology with Arabidopsis genes, which will allow translation of genomics and genetics findings from Arabidopsis to lesquerella. Using Arabidopsis acyl-lipid genes as queries, we searched the transcriptome assembly and identified 615 lesquerella genes involved in all known pathways of acyl-lipid metabolism. Further deep mining the transcriptome assembly led to identification of almost all lesquerella genes involved in fatty acid and triacylglycerol synthesis. Moreover, we characterized the spatial and temporal expression profiles of 15 key genes using the quantitative PCR assay.

Conclusions

We have built a lesquerella seed transcriptome that provides a valuable reference in addition to the castor database for discovering genes involved in the synthesis of triacylglycerols enriched with hydroxy fatty acids. The information obtained from data mining and gene expression profiling will provide a resource not only for the study of hydroxy fatty acid metabolism, but also for the biotechnological production of hydroxy fatty acids in existing oilseed crops.

Electronic supplementary material

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

De novo transcriptome assembly of pummelo and molecular marker development

Pummelo (Citrus grandis) is an important fruit crop worldwide because of its nutritional value. To accelerate the pummelo breeding program, it is essential to obtain extensive genetic information and develop relative molecular markers. Here, we obtained a 12-Gb transcriptome dataset of pummelo through a mixture of RNA from seven tissues using Illumina pair-end sequencing, assembled into 57,212 unigenes with an average length of 1010 bp. The annotation and classification results showed that a total of 39,584 unigenes had similar hits to the known proteins of four public databases, and 31,501 were classified into 55 Gene Ontology (GO) functional sub-categories. The search for putative molecular markers among 57,212 unigenes identified 10,276 simple sequence repeats (SSRs) and 64,720 single nucleotide polymorphisms (SNPs). High-quality primers of 1174 SSR loci were designed, of which 88.16% were localized to nine chromosomes of sweet orange. Of 100 SSR primers that were randomly selected for testing, 87 successfully amplified clear banding patterns. Of these primers, 29 with a mean PIC (polymorphic information content) value of 0.52 were effectively applied for phylogenetic analysis. Of the 20 SNP primers, 14 primers, including 54 potential SNPs, yielded target amplifications, and 46 loci were verified via Sanger sequencing. This new dataset will be a valuable resource for molecular biology studies of pummelo and provides reliable information regarding SNP and SSR marker development, thus expediting the breeding program of pummelo.

Optimal assembly strategies of transcriptome related to ploidies of eukaryotic organisms

BACKGROUND

Several de novo transcriptome assemblers have been developed recently to assemble the short reads generated from the next-generation sequencing platforms and different strategies were employed for assembling transcriptomes of various eukaryotes without genome sequences. Though there are some comparisons among these de novo assembly tools for assembling transcriptomes of different eukaryotic organisms, there is no report about the relationship between assembly strategies and ploidies of the organisms.

RESULTS

When we de novo assembled transcriptomes of sweet potato (hexaploid), Trametes gallica (a diploid fungus), Oryza meyeriana (a diploid wild rice), five assemblers, including Edena, Oases, Soaptrans, IDBA-tran and Trinity, were used in different strategies (Single-Assembler Single-Parameter, SASP; Single-Assembler Multiple-Parameters, SAMP; Combined De novo Transcriptome Assembly, CDTA, that is multiple assembler multiple parameter). It was found that CDTA strategy has the best performance compared with other two strategies for assembling transcriptome of the hexaploid sweet potato, whereas SAMP strategy with assembler Oases is better than other strategies for assembling transcriptomes of diploid fungus and the wild rice transcriptomes.

CONCLUSION

Based on the results from ours and others, it is suggested that CDTA strategy is better used for transcriptome assembly of polyploidy organisms and SAMP strategy of Oases is outperformed for those diploid organisms without genome sequences.

A photorespiratory bypass increases plant growth and seed yield in biofuel crop Camelina sativa

BACKGROUND

Camelina sativa is an oilseed crop with great potential for biofuel production on marginal land. The seed oil from camelina has been converted to jet fuel and improved fuel efficiency in commercial and military test flights. Hydrogenation-derived renewable diesel from camelina is environmentally superior to that from canola due to lower agricultural inputs, and the seed meal is FDA approved for animal consumption. However, relatively low yield makes its farming less profitable. Our study is aimed at increasing camelina seed yield by reducing carbon loss from photorespiration via a photorespiratory bypass. Genes encoding three enzymes of the Escherichia coli glycolate catabolic pathway were introduced: glycolate dehydrogenase (GDH), glyoxylate carboxyligase (GCL) and tartronic semialdehyde reductase (TSR). These enzymes compete for the photorespiratory substrate, glycolate, convert it to glycerate within the chloroplasts, and reduce photorespiration. As a by-product of the reaction, CO2 is released in the chloroplast, which increases photosynthesis. Camelina plants were transformed with either partial bypass (GDH), or full bypass (GDH, GCL and TSR) genes. Transgenic plants were evaluated for physiological and metabolic traits.

RESULTS

Expressing the photorespiratory bypass genes in camelina reduced photorespiration and increased photosynthesis in both partial and full bypass expressing lines. Expression of partial bypass increased seed yield by 50-57 %, while expression of full bypass increased seed yield by 57-73 %, with no loss in seed quality. The transgenic plants also showed increased vegetative biomass and faster development; they flowered, set seed and reached seed maturity about 1 week earlier than WT. At the transcriptional level, transgenic plants showed differential expression in categories such as respiration, amino acid biosynthesis and fatty acid metabolism. The increased growth of the bypass transgenics compared to WT was only observed in ambient or low CO2 conditions, but not in elevated CO2 conditions.

CONCLUSIONS

The photorespiratory bypass is an effective approach to increase photosynthetic productivity in camelina. By reducing photorespiratory losses and increasing photosynthetic CO2 fixation rates, transgenic plants show dramatic increases in seed yield. Because photorespiration causes losses in productivity of most C3 plants, the bypass approach may have significant impact on increasing agricultural productivity for C3 crops.

A photorespiratory bypass increases plant growth and seed yield in biofuel crop Camelina sativa

BACKGROUND

Camelina sativa is an oilseed crop with great potential for biofuel production on marginal land. The seed oil from camelina has been converted to jet fuel and improved fuel efficiency in commercial and military test flights. Hydrogenation-derived renewable diesel from camelina is environmentally superior to that from canola due to lower agricultural inputs, and the seed meal is FDA approved for animal consumption. However, relatively low yield makes its farming less profitable. Our study is aimed at increasing camelina seed yield by reducing carbon loss from photorespiration via a photorespiratory bypass. Genes encoding three enzymes of the Escherichia coli glycolate catabolic pathway were introduced: glycolate dehydrogenase (GDH), glyoxylate carboxyligase (GCL) and tartronic semialdehyde reductase (TSR). These enzymes compete for the photorespiratory substrate, glycolate, convert it to glycerate within the chloroplasts, and reduce photorespiration. As a by-product of the reaction, CO2 is released in the chloroplast, which increases photosynthesis. Camelina plants were transformed with either partial bypass (GDH), or full bypass (GDH, GCL and TSR) genes. Transgenic plants were evaluated for physiological and metabolic traits.

RESULTS

Expressing the photorespiratory bypass genes in camelina reduced photorespiration and increased photosynthesis in both partial and full bypass expressing lines. Expression of partial bypass increased seed yield by 50-57 %, while expression of full bypass increased seed yield by 57-73 %, with no loss in seed quality. The transgenic plants also showed increased vegetative biomass and faster development; they flowered, set seed and reached seed maturity about 1 week earlier than WT. At the transcriptional level, transgenic plants showed differential expression in categories such as respiration, amino acid biosynthesis and fatty acid metabolism. The increased growth of the bypass transgenics compared to WT was only observed in ambient or low CO2 conditions, but not in elevated CO2 conditions.

CONCLUSIONS

The photorespiratory bypass is an effective approach to increase photosynthetic productivity in camelina. By reducing photorespiratory losses and increasing photosynthetic CO2 fixation rates, transgenic plants show dramatic increases in seed yield. Because photorespiration causes losses in productivity of most C3 plants, the bypass approach may have significant impact on increasing agricultural productivity for C3 crops.