Gene expression profiling during seed-filling process in peanut with emphasis on oil biosynthesis networks

Omics-driven advances in the understanding of regulatory landscape of peanut seed development

Peanuts (Arachis hypogaea) are an essential oilseed crop known for their unique developmental process, characterized by aerial flowering followed by subterranean fruit development. This crop is polyploid, consisting of A and B subgenomes, which complicates its genetic analysis. The advent and progression of omics technologies-encompassing genomics, transcriptomics, proteomics, epigenomics, and metabolomics-have significantly advanced our understanding of peanut biology, particularly in the context of seed development and the regulation of seed-associated traits. Following the completion of the peanut reference genome, research has utilized omics data to elucidate the quantitative trait loci (QTL) associated with seed weight, oil content, protein content, fatty acid composition, sucrose content, and seed coat color as well as the regulatory mechanisms governing seed development. This review aims to summarize the advancements in peanut seed development regulation and trait analysis based on reference genome-guided omics studies. It provides an overview of the significant progress made in understanding the molecular basis of peanut seed development, offering insights into the complex genetic and epigenetic mechanisms that influence key agronomic traits. These studies highlight the significance of omics data in profoundly elucidating the regulatory mechanisms of peanut seed development. Furthermore, they lay a foundational basis for future research on trait-related functional genes, highlighting the pivotal role of comprehensive genomic analysis in advancing our understanding of plant biology.

Genome-wide identification and characterization of lncRNAs in sunflower endosperm

BACKGROUND

Long non-coding RNAs (lncRNAs), as important regulators, play important roles in plant growth and development. The expression and epigenetic regulation of lncRNAs remain uncharacterized generally in plant seeds, especially in the transient endosperm of the dicotyledons.

RESULTS

In this study, we identified 11,840 candidate lncRNAs in 12 day-after-pollination sunflower endosperm by analyzing RNA-seq data. These lncRNAs were evenly distributed in all chromosomes and had specific features that were distinct from mRNAs including tissue-specificity expression, shorter and fewer exons. By GO analysis of protein coding genes showing strong correlation with the lncRNAs, we revealed that these lncRNAs potential function in many biological processes of seed development. Additionally, genome-wide DNA methylation analyses revealed that the level of DNA methylation at the transcription start sites was negatively correlated with gene expression levels in lncRNAs. Finally, 36 imprinted lncRNAs were identified including 32 maternally expressed lncRNAs and four paternally expressed lncRNAs. In CG and CHG context, DNA methylation levels of imprinted lncRNAs in the upstream and gene body regions were slightly lower in the endosperm than that in embryo tissues, which indicated that the maternal demethylation potentially induce the paternally bias expression of imprinted lncRNAs in sunflower endosperm.

CONCLUSION

Our findings not only identified and characterized lncRNAs on a genome-wide scale in the development of sunflower endosperm, but also provide novel insights into the parental effects and epigenetic regulation of lncRNAs in dicotyledonous seeds.

A first insight into the genetics of maturity trait in Runner × Virginia types peanut background

'Runner' and 'Virginia', the two main market types of Arachis hypogaea subspecies hypogaea, differ in several agricultural and industrial characteristics. One such trait is time to maturation (TTM), contributing to the specific environmental adaptability of each subspecies. However, little is known regarding TTM's genetic and molecular control in peanut in general, and particularly in the Runner/Virginia background. Here, a recombinant inbred line population, originating from a cross between an early-maturing Virginia and a late-maturing Runner type, was used to detect quantitative trait loci (QTL) for maturity. An Arachis SNP-array was used for genotyping, and a genetic map with 1425 SNP loci spanning 24 linkage groups was constructed. Six significant QTLs were identified for the maturity index (MI) trait on chromosomes A04, A08, B02 and B04. Two sets of stable QTLs in the same loci were identified, namely qMIA04a,b and qMIA08_2a,b with 11.5%, 8.1% and 7.3%, 8.2% of phenotypic variation explained respectively in two environments. Interestingly, one consistent QTL, qMIA04a,b, overlapped with the previously reported QTL in a Virginia × Virginia population having the same early-maturing parent ('Harari') in common. The information and materials generated here can promote informed targeting of peanut idiotypes by indirect marker-assisted selection.

Transcriptome analysis and identification of genes associated with oil accumulation in upland cotton

Cotton is not only the most important fiber crop but also the fifth most important oilseed crop in the world because of its oil-rich seeds as a byproduct of fiber production. By comparative transcriptome analysis between two germplasms with diverse oil accumulation, we reveal pieces of the gene expression network involved in the process of oil synthesis in cottonseeds. Approximately, 197.16 Gb of raw data from 30 RNA sequencing samples with 3 biological replicates were generated. Comparison of the high-oil and low-oil transcriptomes enabled the identification of 7682 differentially expressed genes (DEGs). Based on gene expression profiles relevant to triacylglycerol (TAG) biosynthesis, we proposed that the Kennedy pathway (diacylglycerol acyltransferase-catalyzed diacylglycerol to TAG) is the main pathway for oil production, rather than the phospholipid diacylglycerol acyltransferase-mediated pathway. Using weighted gene co-expression network analysis, 5312 DEGs were obtained and classified into 14 co-expression modules, including the MEblack module containing 10 genes involved in lipid metabolism. Among the DEGs in the MEblack module, GhCYSD1 was identified as a potential key player in oil biosynthesis. The overexpression of GhCYSD1 in yeast resulted in increased oil content and altered fatty acid composition. This study may not only shed more light on the underlying molecular mechanism of oil accumulation in cottonseed oil, but also provide a set of new gene for potential enhancement of oil content in cottonseeds.

Combined Transcriptome and Lipidomic Analyses of Lipid Biosynthesis in Macadamia ternifolia Nuts

Macadamia nuts are considered a high-quality oil crop worldwide. To date, the lipid diversity and the genetic factors that mediate storage lipid biosynthesis in Macadamia ternifolia are poorly known. Here, we performed a comprehensive transcriptomic and lipidomic data analysis to understand the mechanism of lipid biosynthesis by using young, medium-aged, and mature fruit kernels. Our lipidomic analysis showed that the M. ternifolia kernel was a rich source of unsaturated fatty acids. Moreover, different species of triacylglycerols, diacylglycerol, ceramides, phosphatidylethanolamine, and phosphatidic acid had altered accumulations during the developmental stages. The transcriptome analysis revealed a large percentage of differently expressed genes during the different stages of macadamia growth. Most of the genes with significant differential expression performed functional activity of oxidoreductase and were enriched in the secondary metabolite pathway. The integration of lipidomic and transcriptomic data allowed for the identification of glycerol-3-phosphate acyltransferase, diacylglycerol kinase, phosphatidylinositols, nonspecific phospholipase C, pyruvate kinase 2, 3-ketoacyl-acyl carrier protein reductase, and linoleate 9S-lipoxygenase as putative candidate genes involved in lipid biosynthesis, storage, and oil quality. Our study found comprehensive datasets of lipidomic and transcriptomic changes in the developing kernel of M. ternifolia. In addition, the identification of candidate genes provides essential prerequisites to understand the molecular mechanism of lipid biosynthesis in the kernel of M. ternifolia.

Deep transcriptomic study reveals the role of cell wall biosynthesis and organization networks in the developing shell of peanut pod

BACKGROUND

Peanut (Arachis hypogaea L.) belongs to an exceptional group of legume plants, wherein the flowers are produced aerially, but the pods develop under the ground. In such a unique environment, the pod's outer shell plays a vital role as a barrier against mechanical damage and soilborne pathogens. Recent studies have reported the uniqueness and importance of gene expression patterns that accompany peanut pods' biogenesis. These studies focused on biogenesis and pod development during the early stages, but the late developmental stages and disease resistance aspects still have gaps. To extend this information, we analyzed the transcriptome generated from four pod developmental stages of two genotypes, Hanoch (Virginia-type) and IGC53 (Peruvian-type), which differs significantly in their pod shell characteristics and pathogen resistance.

RESULTS

The transcriptome study revealed a significant reprogramming of the number and nature of differentially expressed (DE) genes during shell development. Generally, the numbers of DE genes were higher in IGC53 than in Hanoch, and the R5-R6 transition was the most dynamic in terms of transcriptomic changes. Genes related to cell wall biosynthesis, modification and transcription factors (TFs) dominated these changes therefore, we focused on their differential, temporal and spatial expression patterns. Analysis of the cellulose synthase superfamily identified specific Cellulose synthase (CesAs) and Cellulose synthase-like (Csl) genes and their coordinated interplay with other cell wall-related genes during the peanut shell development was demonstrated. TFs were also identified as being involved in the shell development process, and their pattern of expression differed in the two peanut genotypes. The shell component analysis showed that overall crude fiber, cellulose, lignin, hemicelluloses and dry matter increased with shell development, whereas K, N, protein, and ash content decreased. Genotype IGC53 contained a higher level of crude fiber, cellulose, NDF, ADF, K, ash, and dry matter percentage, while Hanoch had higher protein and nitrogen content.

CONCLUSIONS

The comparative transcriptome analysis identified differentially expressed genes, enriched processes, and molecular processes like cell wall biosynthesis/modifications, carbohydrate metabolic process, signaling, transcription factors, transport, stress, and lignin biosynthesis during the peanut shell development between two contrasting genotypes. TFs and other genes like chitinases were also enriched in peanut shells known for pathogen resistance against soilborne major pathogens causing pod wart disease and pod damages. This study will shed new light on the biological processes involved with underground pod development in an important legume crop.

Identification of consistent QTL for time to maturation in Virginia-type Peanut (Arachis hypogaea L.)

BACKGROUND

Time-to-maturation (TTM) is an important trait contributing to adaptability, yield and quality in peanut (Arachis hypogaea L). Virginia market-type peanut belongs to the late-maturing A. hypogaea subspecies with considerable variation in TTM within this market type. Consequently, planting and harvesting schedule of peanut cultivars, including Virginia market-type, need to be optimized to maximize yield and grade. Little is known regarding the genetic control of TTM in peanut due to the challenge of phenotyping and limited DNA polymorphism. Here, we investigated the genetic control of TTM within the Virginia market-type peanut using a SNP-based high-density genetic map. A recombinant inbred line (RIL) population, derived from a cross between two Virginia-type cultivars 'Hanoch' and 'Harari' with contrasting TTM (12-15 days on multi-years observations), was phenotyped in the field for 2 years following a randomized complete block design. TTM was estimated by maturity index (MI). Other agronomic traits like harvest index (HI), branching habit (BH) and shelling percentage (SP) were recorded as well.

RESULTS

MI was highly segregated in the population, with 13.3-70.9% and 28.4-80.2% in years 2018 and 2019. The constructed genetic map included 1833 SNP markers distributed on 24 linkage groups, covering a total map distance of 1773.5 cM corresponding to 20 chromosomes on the tetraploid peanut genome with 1.6 cM mean distance between the adjacent markers. Thirty QTL were identified for all measured traits. Among the four QTL regions for MI, two consistent QTL regions (qMIA04a,b and qMIB03a,b) were identified on chromosomes A04 (118680323-125,599,371; 6.9Mbp) and B03 (2839591-4,674,238; 1.8Mbp), with LOD values of 5.33-6.45 and 5-5.35 which explained phenotypic variation of 9.9-11.9% and 9.3-9.9%, respectively. QTL for HI were found to share the same loci as MI on chromosomes B03, B05, and B06, demonstrating the possible pleiotropic effect of HI on TTM. Significant but smaller effects on MI were detected for BH, pod yield and SP.

CONCLUSIONS

This study identified consistent QTL regions conditioning TTM for Virginia market-type peanut. The information and materials generated here can be used to further develop molecular markers to select peanut idiotypes suitable for diverse growth environments.

Molecular cloning and expression characterization of flavonol synthase genes in peanut (Arachis hypogaea)

Flavonol is an important functional bioactive substance in peanut seeds, and plays important roles responding to abiotic stress. The flavonol content is closely related to the activity and regulation of gene expression patterns of flavonol synthase (FLS). In this study, eight FLS genes, AhFLSs were cloned and their expression characterization in different peanut organ and seedling under different abiotic stress were conducted. The results showed that the expressions levels of AhFLSs were differed in all assayed peanut organs and seedlings under abiotic stress treatments. Expression levels of AhFLS2, AhFLS3, AhFLS4, and AhFLS6 were higher than those of other AhFLSs. The flavonol contents of peanut organs and seedlings under different abiotic stress were also determined using high performance liquid chromatography (HPLC). Dried mature peanut seeds were the organ tissue with the highest flavonol content, and flavonol content increased with seed development. Under abiotic stress treatments, the types of flavonols induced differed among stress treatments. Correlation analysis results suggested that eight AhFLS genes may have different functions in peanut. Moreover, changes in the expression of the eight genes appear to has substrate preference. These results can lay the foundation for the study of improving nutritional value of peanut seed and resistance of peanut plant.

The characteristic of Arachis duranensis-specific genes and their potential function

Arachis species produce flowers aerially, and then grow into the ground, where they develop into fruits; a feature that is unique to Arachis species. We hypothesized that Arachis species evolved genes specifically involved in the control of aerial flowers and the formation of underground fruits. Arachis duranensis is more resistant to biotic and abiotic stressors. Here, we compared different legume species and identified Arachis duranensis-specific genes. We analyzed gene expression patterns, base substitution patterns and sequence features between genes that are conserved across legume plants and A. duranensis-specific genes. Furthermore, we tested the role of A. duranensis-specific genes during seed development, response to nematode Meloidogyne arenaria infection and drought stress. We found that A. duranensis-specific genes had characteristics of young genes. The gene expression level and breadth were lower in the A. duranensis-specific genes compared to conserved genes. The A. duranensis-specific genes had higher codon usage bias than conserved genes, and the polypeptide length and GC content at the three codon sites were lower compared to conserved genes. Of the A. duranensis-specific genes, single-copy and duplicated genes had different features. The RNA-seq result showed A. duranensis-specific genes were involved in seed development, as well as response to nematode infection and drought stress. In addition, we detected asymmetric functions in A. duranensis-specific duplicated genes in response to nematode infection and drought stress.

Transcriptomic Analysis Reveals the High-Oleic Acid Feedback Regulating the Homologous Gene Expression of Stearoyl-ACP Desaturase 2 (SAD2) in Peanuts

Peanuts with high oleic acid content are usually considered to be beneficial for human health and edible oil storage. In breeding practice, peanut lines with high monounsaturated fatty acids are selected using fatty acid desaturase 2 (FAD2), which is responsible for the conversion of oleic acid (C18:1) to linoleic acid (C18:2). Here, comparative transcriptomics were used to analyze the global gene expression profile of high- and normal-oleic peanut cultivars at six time points during seed development. First, the mutant type of FAD2 was determined in the high-oleic peanut (H176). The result suggested that early translation termination occurred simultaneously in the coding sequence of FAD2-A and FAD2-B, and the cultivar H176 is capable of utilizing a potential germplasm resource for future high-oleic peanut breeding. Furthermore, transcriptomic analysis identified 74 differentially expressed genes (DEGs) involved in lipid metabolism in high-oleic peanut seed, of which five DEGs encoded the fatty acid desaturase. Aradu.XM2MR belonged to the homologous gene of stearoyl-ACP (acyl carrier protein) desaturase 2 (SAD2) that converted the C18:0 into C18:1. Further subcellular localization studies indicated that FAD2 was located at the endoplasmic reticulum (ER), and Aradu.XM2MR was targeted to the plastid in Arabidopsis protoplast cells. To examine the dynamic mechanism of this finding, we focused on the peroxidase (POD)-mediated fatty acid (FA) degradation pathway. The fad2 mutant significantly increased the POD activity and H₂O₂ concentration at the early stage of seed development, implying that redox signaling likely acted as a messenger to connect the signaling transduction between the high-oleic content and Aradu.XM2MR transcription level. Taken together, transcriptome analysis revealed the feedback mechanism of SAD2 (Aradu.XM2MR) associated with FAD2 mutation during the seed developmental stage, which could provide a potential peanut breeding strategy based on identified candidate genes to improve the content of oleic acid.

Gene network of oil accumulation reveals expression profiles in developing embryos and fatty acid composition in Upland cotton

Cottonseed oil accumulated dramatically from 20 days post-anthesis (DPA) to 30 DPA in Upland cotton (Gossypium hirsutum L.). To reveal the gene network of oil accumulation and fatty acid composition in developing embryos, embryos at 10, 20 and 30 DPA in cottonseed were sampled and used for transcriptome sequencing (RNA-Seq). In total, 8629, 7891, and 12,555 differentially expressed genes (DEGs) were identified in the comparison sets of '20 DPA vs 10 DPA', '30 DPA vs 20 DPA', and '30 DPA vs 10 DPA', respectively. The gene network highlighted the dynamic expression profiles of oil accumulation in fatty acid (FA) synthesis, FA desaturation, and triacylglycerol (TAG) biosynthesis. WRI1 and NF-YB6 were suggested elite transcription factors in regulating lipid metabolism. Compared with the gene expression levels in developing seeds, GhPDAT was highly expressed and might play a more important role than GhDGAT in transforming diacylglycerol to TAG in cotton. Expression patterns of 12 FA-biosynthesis-related genes were validated by quantitative real-time PCR (qRT-PCR) method. To reveal the reason for the high content of linoleic acid (C18:2) in cottonseed oil, we carried out a comparative analysis of gene expression levels in Upland cotton, rapeseed (Brassica napus), and oleaster (Olea europaea). Compared with in rapeseed and oleaster, GhFAD2 genes were up-regulated and GhFAD3 genes down-regulated in cottonseed, taking into account the relative high amount of C18:2 but low content of linolenic acid (C18:3) in Upland cotton. The present study offers new information to interpret the mechanism of the FA biosynthesis network and to alter FA composition in cotton breeding projects.

Homoeolog expression bias in allopolyploid oleaginous marine diatom Fistulifera solaris

Background

Allopolyploidy is a genomic structure wherein two or more sets of chromosomes derived from divergent parental species coexist within an organism. It is a prevalent genomic configuration in plants, as an important source of genetic variation, and also frequently confers environmental adaptability and increased crop productivity. We previously reported the oleaginous marine diatom Fistulifera solaris JPCC DA0580 to be a promising host for biofuel production and that its genome is allopolyploid, which had never previously been reported in eukaryotic microalgae. However, the study of allopolyploidy in F. solaris was hindered by the difficulty in classifying the homoeologous genes based on their progenitor origins, owing to the shortage of diatom genomic references.

Results

In this study, the allopolyploid genome of F. solaris was tentatively classified into two pseudo-parental subgenomes using sequence analysis based on GC content and codon frequency in each homoeologous gene pair. This approach clearly separated the genome into two distinct fractions, subgenome Fso_h and Fso_l, which also showed the potency of codon usage analysis to differentiate the allopolyploid subgenome. Subsequent homoeolog expression bias analysis revealed that, although both subgenomes appear to contribute to global transcription, there were subgenomic preferences in approximately 61% of homoeologous gene pairs, and the majority of these genes showed continuous bias towards a specific subgenome during lipid accumulation. Additional promoter analysis indicated the possibility of promoter motifs involved in biased transcription of homoeologous genes. Among these subgenomic preferences, genes involved in lipid metabolic pathways showed interesting patterns in that biosynthetic and degradative pathways showed opposite subgenomic preferences, suggesting the possibility that the oleaginous characteristics of F. solaris derived from one of its progenitors.

Conclusions

We report the detailed genomic structure and expression patterns in the allopolyploid eukaryotic microalga F. solaris. The allele-specific patterns reported may contribute to the oleaginous characteristics of F. solaris and also suggest the robust oleaginous characteristics of one of its progenitors. Our data reveal novel aspects of allopolyploidy in a diatom that is not only important for evolutionary studies but may also be advantageous for biofuel production in microalgae.

Electronic supplementary material

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

Identification of the Candidate Proteins Related to Oleic Acid Accumulation during Peanut (Arachis hypogaea L.) Seed Development through Comparative Proteome Analysis

Peanuts (Arachis hypogaea L.) are an important oilseed crop, containing high contents of protein and fatty acids (FA). The major components of FA found in peanut oil are unsaturated FAs, including oleic acid (OA, C18:1) and linoleic acid (LOA, C18:2). Moreover, the high content of OA in peanut oil is beneficial for human health and long-term storage due to its antioxidant activity. However, the dynamic changes in proteomics related to OA accumulation during seed development still remain largely unexplored. In the present study, a comparative proteome analysis based on iTRAQ (isobaric Tags for Relative and Absolute Quantification) was performed to identify the critical candidate factors involved in OA formation. A total of 389 differentially expressed proteins (DEPs) were identified between high-oleate cultivar Kainong176 and low-oleate cultivar Kainong70. Among these DEPs, 201 and 188 proteins were upregulated and downregulated, respectively. In addition, these DEPs were categorized into biosynthesis pathways of unsaturated FAs at the early stage during the high-oleic peanut seed development, and several DEPs involved in lipid oxidation pathway were found at the stage of seed maturation. Meanwhile, 28 DEPs were sporadically distributed in distinct stages of seed formation, and their molecular functions were directly correlated to FA biosynthesis and degradation. Fortunately, the expression of FAB2 (stearoyl-acyl carrier protein desaturase), the rate-limiting enzyme in the upstream biosynthesis process of OA, was significantly increased in the early stage and then decreased in the late stage of seed development in the high-oleate cultivar Kainong176. Furthermore, real-time PCR verified the expression pattern of FAB2 at the mRNA level, which was consistent with its protein abundance. However, opposite results were found for the low-oleate cultivar Kainong70. Overall, the comparative proteome analysis provided valuable insight into the molecular dynamics of OA accumulation during peanut seed development.

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.

Integrated analysis of 454 and Illumina transcriptomic sequencing characterizes carbon flux and energy source for fatty acid synthesis in developing Lindera glauca fruits for woody biodiesel

BACKGROUND

Lindera glauca fruit with high quality and quantity of oil has emerged as a novel potential source of biodiesel in China, but the molecular regulatory mechanism of carbon flux and energy source for oil biosynthesis in developing fruits is still unknown. To better develop fruit oils of L. glauca as woody biodiesel, a combination of two different sequencing platforms (454 and Illumina) and qRT-PCR analysis was used to define a minimal reference transcriptome of developing L. glauca fruits, and to construct carbon and energy metabolic model for regulation of carbon partitioning and energy supply for FA biosynthesis and oil accumulation.

RESULTS

We first analyzed the dynamic patterns of growth tendency, oil content, FA compositions, biodiesel properties, and the contents of ATP and pyridine nucleotide of L. glauca fruits from seven different developing stages. Comprehensive characterization of transcriptome of the developing L. glauca fruit was performed using a combination of two different next-generation sequencing platforms, of which three representative fruit samples (50, 125, and 150 DAF) and one mixed sample from seven developing stages were selected for Illumina and 454 sequencing, respectively. The unigenes separately obtained from long and short reads (201, and 259, respectively, in total) were reconciled using TGICL software, resulting in a total of 60,031 unigenes (mean length = 1061.95 bp) to describe a transcriptome for developing L. glauca fruits. Notably, 198 genes were annotated for photosynthesis, sucrose cleavage, carbon allocation, metabolite transport, acetyl-CoA formation, oil synthesis, and energy metabolism, among which some specific transporters, transcription factors, and enzymes were identified to be implicated in carbon partitioning and energy source for oil synthesis by an integrated analysis of transcriptomic sequencing and qRT-PCR. Importantly, the carbon and energy metabolic model was well established for oil biosynthesis of developing L. glauca fruits, which could help to reveal the molecular regulatory mechanism of the increased oil production in developing fruits.

CONCLUSIONS

This study presents for the first time the application of an integrated two different sequencing analyses (Illumina and 454) and qRT-PCR detection to define a minimal reference transcriptome for developing L. glauca fruits, and to elucidate the molecular regulatory mechanism of carbon flux control and energy provision for oil synthesis. Our results will provide a valuable resource for future fundamental and applied research on the woody biodiesel plants.