Determination of superior Pistacia chinensis accession with high-quality seed oil and biodiesel production and revelation of LEC1/WRI1-mediated high oil accumulative mechanism for better developing woody biodiesel
BMC Plant Biology volume 23, Article number: 268 (2023)
Based on our previous studied on different provenances of Pistacia chinensis, some accessions with high quality and quantity of seed oils has emerged as novel source of biodiesel. To better develop P. chinensis seed oils as woody biodiesel, a concurrent exploration of oil content, FA profile, biodiesel yield, and fuel properties was conducted on the seeds from 5 plus germplasms to determine superior genotype for ideal biodiesel production. Another vital challenge is to unravel mechanism that govern the differences in oil content and FA profile of P. chinensis seeds across different accessions. FA biosynthesis and oil accumulation of oil plants are known to be highly controlled by the transcription factors. An integrated analysis of our recent transcriptome data, qRT-PCR detection and functional identification was performed as an attempt to highlight LEC1/WRI1-mediated transcription regulatory mechanism for high-quality oil accumulation in P. chinensis seeds.
To select ideal germplasm and unravel high oil accumulative mechanism for developing P. chinensis seed oils as biodiesel, five plus trees (accession PC-BJ/PC-AH/PC-SX/PC-HN/PC-HB) with high-yield seeds were selected to assess the variabilities in weight, oil content, FA profile, biodiesel yield and fuel property, revealing a variation in the levels of seed oil (50.76–60.88%), monounsaturated FA (42.80–70.72%) and polyunsaturated FA (18.78–43.35%), and biodiesel yield (84.98–98.15%) across different accessions. PC-HN had a maximum values of seed weight (26.23 mg), oil (60.88%) and biodiesel yield (98.15%), and ideal proportions of C18:1 (69.94%), C18:2 (17.65%) and C18:3 (1.13%), implying that seed oils of accession PC-HN was the most suitable for ideal biodiesel production. To highlight molecular mechanism that govern such differences in oil content and FA profile of different accessions, a combination of our recent transcriptome data, qRT-PCR detection and protein interaction analysis was performed to identify a pivotal role of LEC1/WRI1-mediated transcription regulatory network in high oil accumulation of P. chinensis seeds from different accessions. Notably, overexpression of PcWRI1 or PcLEC1 from P. chinensis seeds in Arabidopsis could facilitate seed development and upregulate several genes relevant for carbon flux allocation (plastidic glycolysis and acetyl-CoA generation), FA synthesis, TAG assembly and oil storage, causing an increase in seed oil content and monounsaturated FA level, destined for biodiesel fuel property improvement. Our findings may present strategies for better developing P. chinensis seed oils as biodiesel feedstock and bioengineering its high oil accumulation.
This is the first report on the cross-accessions assessments of P. chinensis seed oils to determine ideal accession for high-quality biodiesel production, and an effective combination of PcWRI1 or PcLEC1 overexpression, morphological assay, oil accumulation and qRT-PCR detection was applied to unravel a role of LEC1/WRI1-mediated regulatory network for oil accumulation in P. chinensis seeds, and to highlight the potential application of PcWRI1 or PcLEC1 for increasing oil production. Our finding may provide new strategies for developing biodiesel resource and molecular breeding.
Energy crisis, climate warming, and environment deterioration are still the most severe issues for the world, and these problems should be bound to intensify with the heavy use of fossil energy, population expansion and industrialization acceleration [1,2,3]. Hence, effective exploitation and full development of novel renewable/clean energy has become an urgent issue for sustainable economic development and global energy security. Biodiesel, known as FA methyl ester (FAME), has gradually attracted global great attention as ecofriendly fuel for its nontoxicity and biodegradability , which has been widely applied in several countries such as USA, Brazil, Netherlands, Germany, France and Malaysia . Yet, more than 95% of world biodiesel is produced from edible oils (such as sunflower, peanut and soybean) [1, 6,7,8], which should inevitably lead to reduction of available arable land and shortage of food or edible oils. Thus, exploiting an alternative biodiesel feedstock of non-food plant resources is utmost vital, especially in China with huge population and low per capita arable land. In recent years, the oils from some woody plants (such as Prunus sibirica, Xanthoceras sorbifolia, Pistacia chinensis and Jatropha curcas) have been applied as raw material for biodiesel production in China [9,10,11,12,13,14,15,16]. Notably, the Chinese government has announced to achieve the peak of carbon dioxide emissions by 2030 and carbon neutrality before 2060 at the 75th session of the UN General Assembly (UNGA) in 2020. All these will speed up to develop economical non-food plant resources for bioenergy in China.
Pistacia chinensis Bunge, one perennial deciduous tree of the Anacardiaceae family, is widely distributed in wild mountain areas (geographical coordinates approximately E96°52′-123°14′, N18°09′-40°09′) in China. It was estimated that the total area of artificial plantation and natural forest of P. chinensis was about 86.6 and 133,300 ha, respectively, and the annual yield of P. chinensis seeds was about 330,000 tons in China [11, 17,18,19]. There are many germplasms of P. chinensis with different oil content and FA compositions, but some superior germplasms have not been exploited. Based on our previous studies on different germplasms of P. chinensis, some accessions have been identified with rich oil content and high proportion of oleic and linoleic acid [17,18,19], and notably, the oil content (42.5–50.0%) of mature seeds, which was higher than that of Chinese traditional oil woody plants, such as J. curcas (38.1%), L. glauca (31.6%) and Comus wilsoniana (34.5%) [13, 16, 20]. All these revealed that P. chinensis seed oil could be as potential biodiesel feedstock in China. Considering the differences in seed oil content and FA profiles of different germplasms, it is vitally important to select superior accession with high-quality seed oil for better development of woody biodiesel. However, an interesting open challenge is the mechanism that governed such differences in oil content and FA compositions of P. chinensis seeds across different accessions. Therefore, exploration of molecular regulation mechanism for high oil production in P. chinensis seeds has become another one crucial issue for developing woody biodiesel.
It is known that oil biosynthetic process is composed of carbon source supply, FA synthesis, triacylglycerol (TAG) assembly and oil storage involved in several regulatory enzymes , and their expressions have been identified to be highly regulated by a series of transcription factors (TFs) [22,23,24,25,26,27,28,29,30]. Hence, identification of some critical TFs essential for oil accumulation in P. chinensis seeds could help to unravel TF-mediated molecular regulatory mechanism of increased oil production and may present new strategy for bioengineering oil accumulation in oil plants. Recently, we performed transcriptome sequencing of P. chinensis seeds (SRP030781) to identify some TFs in relation to oil accumulation , which could make us possible to gain better insight into transcriptional mechanism for governing differences in seed oil content and FA compositions across different germplasms of P. chinensis.
This work was undertaken to determine ideal germplasm and to unravel transcriptional regulatory mechanism of high oil production for better development of P. chinensis seed oils as potential feedstock for biodiesel, To this achieve, 5 plus trees (germplasm accessions PC-BJ, PC-AH, PC-SX, PC-HN and PC-HB) with high-yield seeds were selected for the cross-accessions evaluations of seed phenotypic characteristics (size and weight), oil accumulation (content, FA profile, and oil body amount), biodiesel yield and fuel property of the seeds from different accessions. Such exploration could help to select ideal accession for developing biodiesel. Another vital focus of this work was to highlight transcriptional mechanism that govern such difference in seed oil content and FA composition across different accessions, and to gain useful information for bioengineering oil accumulation or molecular-assisted selection. To this set, based on our recent transcriptome result of P. chinensis seeds, the comparative assay of cross-accessions correlation of transcript level of genes for the TFs by qRT-PCR with oil accumulative amount, in combination with the protein interaction assay between the TFs and oil-synthesis enzymes, was described as an initial attempt to identify LEC1/WRI1-mediated network crucial for regulating transcriptions of enzymes relevant for FA synthesis and oil accumulation in P. chinensis seeds of different accessions. Finally, a concurrent exploration was made on dynamic patterns of seedling growth, seed development, oil accumulation (content and FA profile) and some vital gene transcript relevant to FA biosynthesis, TAG assembly and oil storage in transgenic Arabidopsis lines with ectopic overexpression of WRI1 or LEC1 from P. chinensis seeds, with the aim of gaining a better insight into LEC1/WRI1-mediated transcriptional regulation mechanism for high oil production in P. chinensis seeds. Our findings could facilitate the development of P. chinensis seed oils as biodiesel feedstock, and provide new molecular target (WRI1 or LEC1) for bioengineering oil accumulation in oilseed plants.
Variability in seed size and weight of 5 selected accessions of P. chinensis
To better develop the P. chinensis seed oils as novel potential biofuel feedstock, it was crucial for us to select superior germplasm resources with high quantity and quality of seed oils for gaining ideal economic benefit. To achieve this, five plus trees with high seed yield (germplasm accessions PC-BJ, PC-AH, PC-SX, PC-HN, and PC-HB)  were selected as the materials to analyze the variability in the seed size and weight. Here, the dry weights of mature seeds varied from 15.01 mg (PC-AH) to 26.23 mg (PC-HN), followed by 18.84 mg (PC-HB), 19.12 mg (PC-BJ) and 22.87 mg (PC-SX) (Fig. 1a and Additional file 1: Table S1), indicating a notable difference in seed weight across different accessions. This allowed us to explore the differences in longitudinal and transverse diameter of the seeds from different accessions, and a notable variation on longitudinal diameter (3.68–4.87 mm) and transverse diameter (2.56–3.58 mm) of the seeds was detected among different accessions (Fig. 1b, c and Additional file 1: Table S1). Of note, the maximum seed weight (26.23 mg) with the greatest longitudinal diameter (4.87 mm) and transverse diameter (3.58 mm) was noted for accession PC-HN, reflecting that PC-HN as ideal accession could be utilized potentially for the plantation.
Variability in seed oil content and oil body amount of 5 selected accessions of P. chinensis
High-quality oil content is known as one vital determinant of whether oil plant can be suitable for biodiesel production. Here, the highest content of seed oils was detected for accession PC-HN (60.88%), followed by PC-SX (55.37%), PC-HB (54.42%), PC-BJ (52.09%), and PC-AH (50.76%) (Fig. 2a and Additional file 1: Table S1), pointing to a difference in seed oil content across 5 selected accessions. In consideration of oil stored as a form of oil body, it was important to explore the relationship between oil content and oil body amount in the seeds of different accessions. In this work, the morphological difference of seed oil bodies across different accessions was observed by histochemical staining method. Of these, the larger and denser red-stained areas distributed on both longitudinal and transverse sections of the seeds were observed for accession PC-HN (Fig. 2b) compared with the seeds of other accessions (Fig. 2c-f), and the density of stained red area of seed oil bodies of PC-HN was also higher than that of other accessions (Fig. 2g, h and Additional file 1: Table S1), and thus concluded that the seeds of accession PC-HN contained more numbers of oil body cells required for oil storage, which was in line with the results of seed oil contents of different accessions (Fig. 2a).
Variability in FA profiles of the seeds among 5 selected accessions of P. chinensis
Ideal FA compositions (FAs) of plant seed oils is vital for producing high-quality biodiesel, and thus FA compounds and its content of seed oils were analyzed across 5 selected accessions. Here, 7 kinds of FAs were identified in seed oils of all accessions (Table 1), including palmitic acid (C16:0), palmitoleic acid (C16:1), stearic acid (C18:0), oleic acid (C18:1), linoleic acid (C18:2), linolenic acid (C18:3) and arachidic acid (C20:0), of which the richest composition was oleic acid (42.05–69.94%) with an average of 51.97%, followed by linoleic acid (17.65–38.94%), palmitic acid (9.10–12.44%), and linolenic acid (1.13–4.41%), but the others exhibited small amount (0.24–2.11%). Of note, the maximum value of C18:1 (69.94%) and minimum level of C18:2 (17.65%) and C18:3 (1.13%) were all detected for accession PC-HN. Additionally, the total content of C18:1 and C18:2 in seed oils varied from 80.99% to 87.59%, and two accessions (PC-HN and PC-SX) had C18:1 content of more than 55% (Table 1), emphasizing that they were ideal raw material for biodiesel.
It is worth noting that raw plant oils required for ideal biodiesel production should contain a small level (< 30%) of saturated FA (SFA), high content (> 30%) of monounsaturated FA (MUFA), and low amount (< 60%) of polyunsaturated FA (PUFA) [12, 31]. In this work, the contents of MUFA and PUFA were detected for all accessions with the values of 42.80–70.72% and 18.78–43.35%, respectively, but the SFA content of all tested accessions was less than 15% (Table 1). Of these, PC-HN had the maximum content of MUFA (70.72%) and the minimum amounts of PUFA (18.78%) and SFA (10.05%), revealing that the seed oils from accession PC-HN with ideal FA profiles could satisfy the demand of superior biodiesel production.
Another concern about low ratio of PUFA/MUFA (especially C18:2/C18:1) or C20-24/C16-18 as a critical index for assessing tribological properties of plant oils for utilization . The ratios of PUFA/MUFA (0.265–1.013), C18:2/C18:1 (0.252–0.926) and C20-24/C16-18 (0.002–0.004) varied among different accessions (Fig. 3a-c and Additional file 2: Table S2), of which the minimum ratios of PUFA/MUFA (0.265), C18:2/C18:1 (0.252) and C20-24/C16-18 (0.002) were all recorded for accession PC-HN, indicating a good tribological property for the seed oils from accession PC-HN.
Evaluation of biodiesel yield and fuel properties from seed oils among different accessions of P. chinensis
Considering a notable variation on oil content and total proportion of C18:1 and C18:2 in P. chinensis seeds of different accessions (Table 1 and Fig. 2a), it was essential to assess yield and fuel property of biodiesel derived from seed oils of different accessions for determining ideal accession as biodiesel feedstock. Here, the biodiesel yield from seed oils varied across different accessions, ranging from 84.78% (PC-AH) to 98.15% (PC-HN) with an average value of 92.17% (Table 2), of which biodiesel yield of PC-SX (96.76%) and PC-HN (98.15%) was in the standard of EN 14,214 (96.5%).
Biodiesel properties, including iodine value (IV), cetane number (CN), oxidation stability (OS), cold filter plugging point (CFPP), cloud point (CP), density (D), kinematic viscosity (KV), chain length saturated factor (LCSF), and degree of unsaturation (DU), were evaluated with a difference across all tested accessions and described as the follows.
IV, one paramount structural parameter for estimating the unsaturated degree (DU) of FA and OS of biodiesel , varied from 98.15 to 114.53 in the FAMEs from seed oils of different accessions (Table 2), all of which were less than the specified maximum limit (120) of ASTM D6751, EN 14,214 and GB/T20,828 standards. Also, CN is one vital index for determining ignition delay time and combustion quality of biofuels, and its values of all accessions ranged from 48.91 to 51.87, all of which could satisfy the standard of USA (ASTM D6751: CN > 47), and most accessions (except PC-AH) was also in concordance with China standard (GB/T 20,828: CN > 49), while only PC-HN (51.87) could meet Europe standard (EN 14,214: 51 < CN < 65) (Table 2), indicating that the biodiesels from seed oils of all accessions had ideal ignition quality.
KV is one key criteria for determining flow ability of biodiesel, and spray penetration and atomization of fuel . The range of KV value (4.21–4.43 mm2/s) of all accessions could meet the standards of EN 14,214 (3.5 < KV < 5.0), ASTM D6751 (1.9 < KV < 6.0), and GB/T20828 (1.9 < KV < 6.0) (Table 2), referring a good flow or spray ability of biodiesel fuel from seed oils of all accessions. Also, density (D), one key index for assessing combustion efficiency of biodiesel , ranged from 875.89 to 878.84 kg/m3 of all accessions (Table 2), could satisfy the standards of GB/T20,828 (820 < D < 900) and EN 14,214 (860 < D < 900), implying ideal combustion efficiency for the biodiesels from P. chinensis seed oils of all accessions.
OS is one vital technical parameter relevant for stability of biodiesel reaction with air . Our detected OS values of biodiesel fuel from the seed oils of all accessions varied from 2.63 to 3.57 h (Table 2), all of which did not reach the minimum limit (6 h) specified in the standards of EN 14,214 and GB/T 20,828, but only accession PC-HN (3.57 h) could meet the ASTM D6751 standard (OS > 3.0 h), which may be attributable to low amount of PUFA (18.78%) in seed oils of PC-HN compared with other accessions (34.59–44.25%) (Table 1).
Both CFPP and CP are known as two key low-temperature parameters to determine the maximum of filterability, but not limited by the US and European standards. The CFPP valve ranged from -8.90 °C to -6.56 °C for biodiesels across all accessions (Table 2), less than the maximum limit (0 °C) of Germany standard (DIN V51,606) in summer, but higher than the minimum limit (-10.0 °C) for spring and autumn. Also, all accessions had an ideal CP value (-6.28 °C to -3.98 °C) (Table 2) within ASTM D6751 standard (-12 °C < CP < -3 °C), implying a good cold flow property of biodiesel from seed oils of all accessions, especially PC-HN with the minimum value of CFPP (-8.90 °C) and CP (-6.28 °C), which was attributed to high level of unsaturated FA (Fig. 3a) and small amount of length chain saturated factor (LCSF) (Table 2).
Also noteworthy was the C18:3 content and the FAs with four double bonds in FAMEs. Low amount (1.13–4.41%) of C18:3 in seed oils of all accessions and no four double-bond FAs (C18:4 and C20:4) in the FAMEs (Table 1) all met the EN14,214-2008 specification (< 12% and 1%, respectively).
Another concern about DU and LCSF as two additional parameters based on the type/structure of FAs. The values of CN, IV and OS are affected greatly by the DU , but both CFPP and CP mostly depends on LCSF [31, 32]. The DU values were detected in the range of 109.40 to 133.91 across all accessions, of which PC-HN had a minimum of 109.40, with the highest values of CN (51.87) and OS (3.57 h) and the lowest value of IV (98.15) (Table 2). However, a negative correlation of LCSF with CFPP or CP was noted across all accessions (Table 2). It seems therefore that both DU and LCSF may be as vital indicators for evaluating biodiesel fuel properties of seed oils of P. chinensis from different accessions.
Global identification of potential TFs specific for high oil accumulation in P. chinensis seeds
Another one interesting question was the mechanism that govern such difference in seed oil content and FA profile of P. chinensis across all accessions (Table 1 and Fig. 2a). FA synthesis and oil accumulation of plants are highly controlled by the TFs, but the question of which TF and how to regulate seed oil accumulation of P. chinensis is still elusive. Recently, differential transcripts were annotated for 12 TFs (LEC1, ABI3, L1L, FUS3, WRI1, VAL1, GL2, AP2, VAL2, VAL3, EMF2 and DOF1) involving potentially in seed oil synthesis of P. chinensis by our transcriptome assay (SRP030781) , which promoted to unravel the potential association of their transcript levels with accumulative amount of seed oils across all accessions with the aim to identify TFs essential for high oil production of P. chinensis seeds. Here, transcript levels of WRI1 and LEC1 were much higher than those of the other 9 TFs (Fig. 4a and Additional file 3: Table S3), and importantly, only both LEC1 and WRI1 transcripts exhibited a pattern that highly correlated with seed oil content across all accessions (Figs. 2a and 4a). Hence, we further performed protein interaction assay for our all identified TFs as an important attempt to explore TF-mediated regulation mechanism of high oil production of the seeds, in which some enzymes crucial for oil accumulation, including carbon source provision (PK, GPI, ENO, E1-α, E1-β, and E3), de novo FA biosynthesis and elongation (ACC2, BCCP1/2, KASI/III, EAR, MAT and FATA/B), FA desaturation (FAD2/3), and TAG assembly and storage (TAG1, PDAT1 and GPAT9), were identified to be highly correlated with WRI1 (interaction score > 0.7), indicating a central regulator of WRI1 for oil accumulation in P. chinensis seeds. Notably, LEC1 displayed the highest confidence level (interaction score = 0.958) to interact with WRI1 among the selected TFs (Fig. 4b), and thus considered that LEC1 and its targeted WRI1, located in center position of interaction network, may contribute mostly to regulate transcriptions of enzymes in relation to oil accumulation of P. chinensis seeds.
Together, accession PC-HN with high seed oil content and biodiesel yield, and ideal fuel property, could be suitable for high-quality biodiesel production, and LEC1/WRI1-mediated transcription regulatory network may play a major role in FA synthesis and oil accumulation in P. chinensis seeds of different accessions. Hence, our following work focused on identifying the functions of WRI1 and LEC1 and exploring their molecular regulatory mechanism in P. chinensis seeds.
Isolation and structural analysis of PcWRI1 and PcLEC1 genes from P. chinensis seeds
To gain a better insight into LEC1/WRI1-mediated regulatory mechanism for high oil accumulation of P. chinensis seeds, based on our recent transcriptome sequencing results from P. chinensis seeds , the full-length cDNA sequences of WRI1 (1,377 bp) and LEC1 (988 bp) (designated as PcWRI1 and PcLEC1, respectively) were obtained from the seeds of ideal accession PC-HN, with the open reading frames (ORFs) of 1,278 bp and 741 bp, encoding proteins containing 425 and 247 amino acids, respectively (Fig. 5a, b and Additional file 4: Fig. S1).
To predict the potential biological functions and evolutionary relationships of PcWRI1 and PcLEC1, the assays were conducted for phylogenetic tree, multiple sequence alignment, conserved domain and 3D structure. Here, PcWRI1 and PcLEC1 had typical ‘VYL’ and ‘CBFD_NFYB_HMF’ motifs respectively, and showed high-similar amino acid sequence and 3D structure with other oilseed plants (such as Glycine max and Zea mays) (Additional file 5: Fig. S2), pointing to a semblable regulatory function of PcWRI1 and PcLEC1 in FA biosynthesis.
To determine the function of PcWRI1 or PcLEC1 from P. chinensis seeds, the CaMV35S promoter-driven expressing vectors of PcWRI1 (CaMV35S::PcWRI1) and PcLEC1 (CaMV35S::PcLEC1) were constructed and transformed into model plant Arabidopsis respectively. The T3 homozygous seeds were obtained by resistance selection, PCR assay (Fig. 5c-i) and chi-square (χ2) test (Additional file 6: Table S4), from which T4 transgenic plants were used for further analyses.
Morphological analyses of transgenic PcWRI1 and PcLEC1 lines
To examine whether ectopic overexpression of PcWRI1 or PcLEC1 affected morphological changes in transgenic lines, a concurrent exploration was conducted on root length, plant height, rosette and leaf size, and seed dry weight. After germination of T3 seeds on half-strength MS medium for 10 d, the root lengths of PcLEC1 and PcWRI1 transgenic lines were detected with 4.98 and 4.75 cm, respectively, which were all longer than those of wild type (WT, 3.05 cm) and null transgenic lines (NT, 3.28 cm) (Fig. 6a, h). PcLEC1 transgenic line had the earliest bolting time (Fig. 6b), about 3 and 5 d earlier than that of the control and PcWRI1 transgenic lines, respectively, and its mature plant height (42.33 cm) was higher than PcWRI transgenic plants (34.67 cm), both of which were higher than that of WT (31.33 cm) and NT (30.04 cm) (Fig. 6c, j). Compared with the controls (WT and NT), PcWRI1 and PcLEC1 transgenic lines could produce the larger and longer rosette radius and leaf, but the number of rosette leaves showed no significant change among all tested lines (Fig. 6d-g). These results implied that ectopic overexpression of PcWRI1 (especially PcLEC1) could promote growth and development of root, stem and leaf of transgenic plants. Also, the dry weights of mature 1,000 seeds from transgenic PcWRI1 and PcLEC1 lines were 20.60 mg and 16.21 mg, respectively, about 0.65- and 0.30-fold higher than those of the controls (Fig. 6i and Additional file 7: Table S5), emphasizing that ectopic overexpression of PcLEC1 (especially PcWRI1) could facilitate reproductive development and seed growth of transgenic lines.
Effect of PcWRI1 or PcLEC1 overexpression on seed oil accumulation and biodiesel fuel property
To explore potential effect of PcWRI1 or PcLEC1 overexpression on oil accumulation of transgenic seeds, we detected dynamic patterns of oil content and FA compositions during seed development of WT, NT, and transgenic PcWRI1 and PcLEC1 lines by GC–MS. Here, the oil contents of developing seeds of all tested lines gradually increased at 4–16 d after anthesis (DAA) with a rapid increase at 12–16 DAA, followed by a slight decrease at 20 DAA (full rape) (Fig. 7b and Additional file 8: Table S6). Of note, the oil contents of ripen dry seeds of 35S::PcWRI1 and 35S::PcLEC1 lines were 59.88% and 57.04%, respectively, both of which was higher than that of the WT (48.45%) and NT (47.19%), implying that overexpression of PcWRI1 or PcLEC1 increased seed oil content, but PcWRI1 overexpression showed the relatively great effect on seed oil accumulation.
A total of 7 kinds of FAs were identified in developing seed oils of all tested lines (35S::PcWRI1, 35S::PcLEC1, NT, and WT), including C16:0, C18:0, C18:1, C18:2, C18:3, C20:0 and C20:1 (Fig. 7a), all of which displayed no significant difference in oil content of developing seeds at 4–8 DAA across all lines. The contents of C16:0 and C18:0 decreased in the seeds of all tested lines during whole development (4–20 DAA), but an obvious decrease was noted for transgenic seeds. The contents of C18:1, C18:2 and C20:1 in seed oils of both NT and WT gradually increased at 4–20 DAA, but an increase of these was recorded for transgenic seeds before 16, 12 and 16 DAA, respectively. Also, a notable increase of C18:3 content in the controls (NT and WT) and transgenic seeds was observed before 16 and 8 DAA, respectively. Given that a higher amount of C18:1 and C20:1 and a lower content of other FAs (C16:0, C18:0, C18:2 and C18:3) were detected for ripen seeds of transgenic lines compared with the control, it seems likely that overexpression of PcWRI1 or PcLEC1 could result in an alteration of FA profile in seed oils of transgenic lines.
Another concern was that whether the overexpression of PcWRI1 or PcLEC1 was implicated in affecting biodiesel fuel property of seed oils from transgenic lines. To this end, the percentages of SFA, MUFA and PUFA in ripen seed oils of all tested lines (Fig. 7c) were calculated as three angular points to establish a triangular graph. It was found that all tested lines were presented in light-blue area of our constructed prediction model, but PcWRI1 and PcLEC1 transgenic lines was located at the far end of PUFA angular point (lower left vertex) and SFA angular point (lower right vertex) (Fig. 7d), implying that seed oils from transgenic lines had good fuel properties (CN, IN, CFPP, OS and CP) [12, 16]. Thus, PcWRI1 or PcLEC1 overexpression could improve biodiesel property of seed oils, which coincided with the fact that a higher MUFA mount (55.83% and 52.76%, respectively) and a lower PUFA level (30.52% and 33.15%, respectively) were detected for ripen seed oils of transgenic PcWRI1 and PcLEC1 lines compared with the WT (37.04% and 46.46%, respectively) (Fig. 7c).
Analyses of spatiotemporal expressions of PcWRI1 and PcLEC1 in transgenic lines
To gain insight into molecular regulation mechanism of PcWRI1 and PcLEC1 on FA synthesis and oil accumulation, we explored spatial–temporal expression pattern of PcWRI1 or PcLEC1 in three T3 independent transgenic lines (Fig. 8 and Additional file 9: Table S7). Here, low transcript level of PcWRI1 was identified in vegetative organs (flowers, leaves, stems and roots) of PcWRI1 transgenic lines, but a gradual and significant increase was marked for the seeds during development (Fig. 8a), emphasizing that PcWRI1 was specifically expressed in seed with development-dependently increase pattern. Yet, PcLEC1 transcript was detected with high level in the leaves, stems and roots of PcLEC1 transgenic lines, and a sharp peak in transcript was observed in developing seeds from 8 to 12 DAA, followed by a decrease (Fig. 8b). Also, transcript level of PcLEC1 in the leaves, stems and roots was higher than that of PcWRI1 in those vegetative organs of PcWRI1 transgenic lines. These results revealed a difference in spatiotemporal expression pattern between PcWRI1 and PcLEC1.
Effect of ectopic PcWRI1 or PcLEC1 on transcript of oil synthesis-related TFs in developing seeds of transgenic lines
To further examine the mechanism required for an alteration in seed oil content in PcWRI1 or PcLEC1 over-expressing transgenic lines (Fig. 7b), it was vital to explore whether ectopic expression of PcWRI1 or PcLEC1 affected transcript of endogenous oil-related TFs (such as L1L, LEC2, FUS3 and ABI3) in the transgenic seeds during development (Fig. 9 and Additional file 10: Table S8). Compared with the WT, ectopic overexpression of PcWRI1 or PcLEC1 could result in a decrease in transcript level of endogenous AtLEC1 and AtLEC2 during seed development (Fig. 9a, b), while the notably upregulated transcript was noted for AtABI3, AtWRI1, AtFUS3 and AtL1L (Fig. 9c-f), and thereby reflected that ectopic overexpression of PcLEC1 could prompt transcript of its targets (AtABI3, AtWRI1, AtFUS3 and AtL1L).
Temporal transcripts of oil accumulation-related genes in developing seeds of PcWRI1 or PcLEC1 transgenic lines
To highlight how PcWRI1 and PcLEC1 regulate FA biosynthesis and oil accumulation at molecular level, the comparative analysis of temporal transcript changes was conducted on 28 key genes encoding for regulatory enzymes in developing seeds of the WT and transgenic lines, involved specifically in carbon source supply (glycolysis and acetyl-CoA formation) and oil synthetic process (FA synthesis and modification, and TAG assembly and oil storage).
Acetyl-CoA, one key precursor for FA biosynthesis, is mostly derived from pyruvate (PYR), involved in glycolysis and PYR dehydrogenase complex (PDC). By qRT-PCR, the differentially up-regulated transcripts were marked for PYR kinase (PK) and phosphofructokinase (PFK) of plastidic and cytosolic glycolysis, and mitochondrial and plastidic PDC subunits (E1-α, E1-β, E2 and E3) relevant for acetyl-CoA formation during seed development of both the control and transgenic lines, of which plastidic PFK, PK and PDC all exhibited a higher transcript (Fig. 10a, b), emphasizing an importance of plastidic glycolysis and PDC in providing carbon sources (PYR and acetyl-CoA) for FA synthesis in developing seeds. Yet, abundant transcript of PFK, PK and PDC was all identified in developing seeds of transgenic PcWRI1 lines as compared with transgenic seeds of PcLEC1, both of which were higher than that of the WT (Fig. 10a, b), indicating that PcWRI1 or PcLEC1 overexpression could strongly increase plastidic glycolysis capacity and PDC activity via effective transcriptional regulation to provide large acetyl-CoA destined for high FA production in developing seeds of transgenic lines.
High oil accumulation in oil plants is highly dependent on high FA synthesis and TAG assembly as well as oil storage, and thus we performed qRT-PCR to explore quantitative association of oil accumulative amount with oil-synthetic gene transcript level in developing seeds of all tested lines. Compared with the WT, overexpression of PcLEC1 or PcWRI1 significantly increased transcripts of FA synthetic enzymes in developing seeds, including acetyl-CoA carboxylase (ACC), 3-ketoacyl ACP synthase I/II/III (KAS I/II/III), fatty acyl-ACP thioesterase A/B (FATA/B), malonyl-CoA-ACP transferase (MAT), 18:0-ACP desaturase 6 (SAD6) and enoyl-ACP reductase (EAR), but a higher transcript level of them was noted for PcWRI1 transgenic seeds (Fig. 10c), as was the case for key enzymes of TAG assembly [acyl-CoA:G3P acyltransferase 9 (GPAT9), acyl-CoA:LPA acyltransferase 4 (LPAAT4), PA phosphatase 2 (PAP2), DAG acyltransferase 1 (PDAT1) and acyl-CoA:DAG acyltransferase 1 (DGAT1)] (Fig. 10d) and TAG accumulation-related oil body protein [oleosin-1/-2/-3/-4/-5 (OLEO1/2/3/4/5)] (Fig. 10e). A temporal and strong correlation of their transcript levels with oil content of developing seeds of transgenic lines (Figs. 7b and 10c, d) reflected that PcLEC1 overexpression (especially PcWRI1) could activate transcript of multiple genes required for FA synthesis, TAG assembly and oil storage in developing seeds of transgenic lines. Of note, transcript level of DGAT1 was much higher than that of PDAT1 across all tested lines (Fig. 10d), implying that DGAT1 may be a crucial regulator for TAG assembly, and the most abundant transcript of OLEO1 implies (Fig. 10e) its importance in seed oil body formation. Also noteworthy was the role of long chain acyl-CoA synthetase (LACS) in the export of plastidic free FA to produce acyl-CoA pool for TAG formation . High transcript of LACS8/9 (Fig. 10d) in developing seeds of transgenic lines was paralleled to that of FA-synthesized enzymes (Fig. 10c), and showed a pattern that correlated with oil amount during seed development (Fig. 7b), and thus referred that overexpression of PcLEC1 or PcWRI1 could promote LACS8/9 transcript to export FA from plastid into ER for acyl-CoA pool generation destined to TAG assembly.
Another concerned was that the PUFA (C18:2 and C18:3) level, as a key factor affecting the quality of plant oils , was produced from desaturation of C18:1 by FAD2 and FAD3, respectively. The differentially upregulated transcripts of FAD2 (high) and FAD3 (low) were detected in developing seeds of all lines, but the relative high transcript of them was identified for the WT (Fig. 10d), reflecting that PcLEC1 or PcWRI1 overexpression could suppress transcript of FAD2/3, causing high content of C18:1 and low amount of C18:2 and C18:3 in ripen seed oils of both transgenics lines (Fig. 7a).
P. chinensis seed oil as novel and potential woody feedstock for biodiesel production in China
To determine ideal accession for developing P. chinensis seed oil as novel woody feedstock for biodiesel production, the cross-accessions comparisons were conducted on a subsequent evaluation of seed phenotypic characteristics (size and weight), oil body amount, oil content, FA profile, biodiesel yield and fuel properties in P. chinensis seeds of 5 selected accessions (Tables 1, 2; Figs. 1, 2 and 3). The seed oil content (50.76–60.88%) of P. chinensis from all accessions (Fig. 2a) was not only greater than that of other traditional woody oil plants (11.68–40.28%) [12,13,14, 37,38,39,40,41,42,43,44], but also higher than that reported for P. chinensis seeds from other germplasms (42.5–50.0%) [17,18,19], implying that P. chinensis accessions, especially for high oil-bearing PC-HN (60.88%) and PC-SX (55.37%), had great potential as biodiesel feedstock. This, integrated with high yield of biodiesel from seed oils of PC-HN (98.15%) and PC-SX (96.76%) (Table 2) could satisfy the standard of EN 14,214 (96.5%), indicated that the P. chinensis seed oils of PC-HN and PC-SX accessions may be advantage for high-yield biodiesel production. It is worth noticing that the high-quality oils for ideal biodiesel production should contain suitable amount of SFA, high percentage of MUFA and low level of PUFA [12, 16, 31,32,33, 45]. In this work, all accession seeds were detected with high content of MUFA (42.80–70.72%) and low levels of PUFA (18.78–43.35%) and SFA (10.50–14.54%) (Table 1), and importantly, the total amount of C18:1 and C18:2 in all accessions accounted for more than 81% of seed oils (Table 1), emphasizing that the seed oils from all accession may be as raw material suitable for biodiesel, corresponded to previous results of several traditional woody oil plants [12, 46,47,48,49,50]. Considering that accession PC-HN had the maximum content of MUFA (70.72%), the minimum amounts of PUFA (18.78%) and SFA (10.50%) and a small ratio (0.252) of C18:2/C18:1 (Fig. 3a, b), together with superior value of IV (98.15), CN (51.87), OS (3.57 h), CP (-6.28 °C) and CFPP (-8.90 °C) in biodiesel (Table 2), it seems certain that the seed oils from PC-NH with superior FA compositions could satisfy the requirement for high-quantity biodiesel production.
LEC1/WRI1-mediated transcription regulation for high oil production in seeds of P. chinensis and transgenic lines
In oil plants, FA synthesis and oil accumulation as a part of seed development, involving a series of gene expression, is highly regulated by several TFs [22,23,24,25,26,27,28,29, 51,52,53,54,55,56,57,58,59]. In this work, of our recent annotated 12 TFs (LEC1, ABI3, FUS3, WRI1, VAL1, L1L, GL2, AP2, VAL2, VAL3, EMF2 and DOF1) in P. chinensis seeds by transcriptome analysis , only both LEC1 and WRI1 transcripts increased and coincided with seed oil content across all accessions (Fig. 4a). Given that LEC1 and its targeted WRI1 at the center of interaction network were highly correlated with oil-synthetized enzymes (Fig. 4b), it seems clear that both WRI1 and LEC1 may act as critical positive regulators crucial for gene transcriptional expression relevant to seed oil accumulation of P. chinensis across all accessions, as was noted for several oil plants [29, 60,61,62,63]. Of note, WRI1, has been shown to control the expression of targeted genes involved in glycolysis, FA synthesis and TAG accumulation during oilseed development [29, 53, 57, 64,65,66]. This was clearly demonstrated by our findings that the abundantly coordinated transcripts of the enzymes relevant for carbon allocation (plastidic glycolysis and acetyl-CoA formation), FA synthesis, TAG assembly and oil storage in developing seeds of transgenic PcWRI1 lines (Fig. 10) were temporally associated with the increased amount of seed oils (Fig. 7b), emphasizing that these genes, as the targets of WRI1, were sequentially and transcriptionally activated by PcWRI1 overexpression, giving rise to a strong collaborative funneling of carbon and acyl flux into oil biosynthetic pathway for high oil production in developing seeds of transgenic PcWRI1 lines, destined for an 23.59% increase in oil content of ripen seed as compared with the WT (Fig. 7b). This was consistent with former results of other oilseed plants [52, 53, 65, 67,68,69,70,71], and also evidenced by the fact of an increase in FA synthesis and oil accumulation by ectopic overexpression of WRI1 from several oilseed plants (such as rapeseed, oil palm, sunflower, coconut, maize, soybean, Camelina sativa and J. curcas) [22,23,24, 27, 28, 65,66,67, 70, 72].
As a crucial TF governing plant oil accumulation, WRI1 is known to be up-regulated by LEC1/2, ABI3 or FUS3 . Here, a similar pattern of high transcript for LEC1 and WRI1 in P. chinensis seeds across all accessions (Fig. 4a) revealed WRI1 as one target gene of LEC1, which was also identified in A. thaliana and P. sibirica seeds [14, 51], but contrasted with the idea that WRI1 function may be independent of upstream TFs in controlling oil biosynthesis in oil palm . LEC1 has been shown as a central regulator in TF-mediated regulation network for oil biosynthesis by transcriptional control distinct gene sets during seed development [62, 64, 73, 74], and LEC1 overexpression increased oil content by up-regulating several genes relevant for FA and oil synthesis in oil plants [25, 54, 64, 74, 75], but lec1 loss-of-function mutation represented suppression of oil accumulation . In this study, the transcripts of endogenous oil-related TFs (AtABI3, AtWRI1, AtFUS3 and AtL1L) were significantly up-regulated in developing seeds of transgenic PcLEC1 lines (Fig. 9c-f), showing the same temporal co-transcription of oil biosynthetic genes (Fig. 10c-e), reflecting indeed that ectopic expression of PcLEC1 could promote transcript of its targets (AtABI3, AtWRI1, AtFUS3 and AtL1L). Given that a higher root length, plant height, and rosette and leaf size was detected in transgenic PcLEC1 line than in transgenic PcWRI1 line (Fig. 6), integrated with the relatively low amount of seed oil (Fig. 7b) and low transcript level of oil accumulation-related enzymes (carbon flux allocation, FA synthesis, TAG assembly and oil storage) in developing seeds of transgenic PcLEC1 line compared with the transgenic PcWRI1 lines (Fig. 10), it seems certain that overexpression of PcLEC1 may contribute mostly to plant development of transgenic line, but an increase of seed oil content may be likely attributed to PcWRI1 overexpression.
Overall, the strong correlation of seed oil accumulation (Fig. 7b) with rich transcript of the enzymes essential for oil biosynthetic process (carbon allocation, FA synthesis, TAG assembly and oil body formation) in developing seeds of transgenic lines with PcLEC1 or PcWRI1 (Fig. 10) revealed that LEC1/WRI1-mediated transcription regulatory network may play a major role in plant growth, seed development and oil accumulation.
In this work, the cross-accessions comparisons on seed weight, oil content, FA composition, biodiesel yield, and fuel properties were conducted to determine ideal accession PC-HN as novel promising feedstock for superior biodiesel production. The application of an integrated assay on the cross-accessions correlation of seed oil accumulative amount with TF transcript level, and protein interaction between the TFs and oil-synthesis enzymes, has led to identification of two key TFs (LEC1 and WRI1) involved in LEC1/WRI1-mediated regulatory network essential for high-quality seed oil accumulation of P. chinensis of different accessions. Notably, the transgenic results indicated that overexpression of PcWRI1 or PcLEC1 from P. chinensis seeds in Arabidopsis could facilitate seedling growth and seed development, and increase transcripts of several genes relevant for carbon allocation (plastidic glycolysis and PDC), FA synthesis (ACC, BCCP1/2, MAT, FATA/B, KAS I/II/III, EAR and SAD6), TAG assembly (LACS8/9, GPAT9, LPAAT4, PAP2 and DGAT1) and oil body formation (OLEO1/2/3/5), leading to an increase in ripen seed oil content (17.73%-23.59%) and MUFA amount (42.44%-50.73%), destined for biodiesel property improvement. Our findings may present strategies for developing P. chinensis seed oils as biodiesel feedstock, and highlight translational potential of PcWRI1 or PcLEC1 for further bioengineering of high oil yield and quality destined for ideal biodiesel in P. chinensis and other oilseed plants.
Materials and methods
Collection of fruit samples of P. chinensis and detections of seed size and weight
Based on our previous studies on 18 provenances of P. chinensis located in south of Qinling Mountains (geographical coordinates approximately E113°00′-119°26′, N28°26′-33°49′) [17, 19], a total of 5 plus trees (defined as superior trees) of P. chinensis (germplasm accessions PC-BJ, PC-AH, PC-SX, PC-HN and PC-HB) with high seed yield were selected and identified, and deposited in the Forest and Flower Germplasm Resource Genebank of Beijing Forestry University in China (Voucher No. 1111C0003102003374, 1111C0003102003402, 1111C0003102003408, 1111C0003102003113 and 1111C0003102003108, respectively). The ripen fresh fruits were collected from 10-year-old tree, and 4000 fruits from 10 trees (400 fruits each tree) of each accession were randomly selected each time. After removing the sarcocarp, the undamaged and plump fresh seeds were selected and immediately frozen in liquid nitrogen and stored at -80 °C until use for FA detection, qRT-PCR assay and gene cloning. The size and weight of seeds were respectively detected with vernier caliper and electronic balance. All analyses were conducted in triplicate.
Plant material and growth conditions
The wild type (WT) Arabidopsis was Col-0 ecotypes, and the seeds of WT and transgenic line with the empty vector pCAMBIA2300/CaMV35S::GUS (null transgenic, NT) were provided by our laboratory. To prepare Arabidopsis plants for genetic transformation of PcWRI1 and PcLEC1, the WT seeds were surface-sterilized and germinated on half-strength MS medium, and imbibed in the dark for 3–4 d at 4 °C, and then were placed into an incubator under a continuous artificial light period of 16 h (22 °C) at a photon flux density of 200 μmol m−2 s−1 and a dark period of 8 h (18 °C). After 20 d, the seedlings were transplanted into the pots containing a soil mixture (humus-soil/vermiculite/perlite, 3:3:1, v/v/v) and grown under a condition of 16 h photoperiod (24 °C) and 8 h dark period (22 °C) for gene transformation and seed harvest.
Section and observation of seed lipid droplet of P. chinensis
The mature seeds of each accession (PC-BJ, PC-AH, PC-SX, PC-HN and PC-HB) used for lipid droplet observation were horizontally and vertically transected, and colored by safranin with fast green. The resulting sections were scanned by PANNORAMIC MIDI (3DHISTECH, Hungary), and the density of lipid droplet was calculated by Image J software 1.8.0. Ten individual seeds were measured in each accession.
Oil extraction and trans-esterification
The analyses of seed oil content and FA compositions were conducted on the P. chinensis seeds of all accessions, and the developing seeds of WT, NT, and T3 generation of PcWRI1 and PcLEC1 transgenic Arabidopsis lines.
About 10 g of the undamaged and plump dry seeds of P. chinensis (3 samples per accession) were pulverized into powders with domestic grinder and then lyophilized by freeze dryer (LGJ-12, China) to extract oil with petroleum ether using Soxhlet apparatus at 45–50 °C for 6–8 h , and then the oils were separated from organic mixture by rotary evaporator and dried to constant in ventilated oven at 105 °C. The content of extracted oil from each accession seed was calculated as the difference between the weights of seed sample before and after extraction, and expressed as the percentage of the extracted oil weight to dry seed weight (%, g/g). To analyze FA compositions, the seed oils from each accession were trans-esterified as previously described .
A total of 20 Arabidopsis seeds (per developing period) were weighted and trans-methylated at 95 °C for 1.5 h in 1 mL of methanol (containing 50 μL H2SO4, 50 μg BHT and 25 0 μg C17:0) and 300 μL toluene, and then 1.5 mL NaCl solution (0.9%, m/v) was added. The obtained oils were recovered by three sequential extractions with 1 mL of hexane, and evaporated under a stream of nitrogen and re-dissolved in 50 μL hexane. The oil contents were calculated by the following equations : Oil content (%, g/g) = 100 × [ΣAi + 4 × (ΣAi /MWi)/3]/W. Where Ai is the weight of the ith FAME, MWi is the molecular mass of each component, W is the seed weight of A. thaliana, 4 is the Mr difference between three moles of FAME and TAG.
All the determination was conducted for three biological replicates of seed sample with three technical replicates.
Analyses of FAMEs and biodiesel yield of P. chinensis seed oils
The FA methyl esters (FAMEs) obtained from each accession seed was used to detect FA compositions by Agilent 6890 (California, USA) gas chromatograph equipped with flame ionization detector (GC-FID) . The HP-INNOWax capillary column (inner diameter 0.32 mm, filmthickness 0.5 μm, split 1:20) was used, and the temperature was programmed at 60 °C, with a rise of 4 °C min−1 to 220 °C and heated to 240 °C for 10 min. The carrier gas was helium with a flow rate of 1.0 mL min−1. The peaks of FAMEs were identified by comparing their retention time with that of the known standards, and peak integration was performed by applying HP3398A software. Each FAME assay was performed in triplicate. The biodiesel yield from P. chinensis seed oils was calculated by the previous method [12, 16], where the yield was defined as the percentage (%, g/g) of the obtained total amount of FAMEs (g) to the used amount of raw oils (g).
Evaluation of biodiesel fuel property of P. chinensis seed oils
The biodiesel fuel properties (CN, IV, CFPP, CP, OS, KV and D) of P. chinensis seed oils from 5 accessions were evaluated by the FAME compositions according to our previously method . The biodiesel fuel characteristics of raw oils were compared with the relevant biodiesel standards of EN 14,214–2008 (European, 2008), ASTM D6751-2010 (USA, 2010), DIN V51,606-1997 (Germany, 1997) and GB/T 20,828–2007 (China, 2007).
Protein interaction network analysis
Based on previous transcriptome data of P. chinensis seeds , the Retrieval of Interacting Genes/Proteins (STRING version 9.0, http://string90.embl.de/) was employed to analyze the potential interactions between all our identified TFs and functional proteins involved in FA synthesis and oil accumulation. STRING analysis was conducted using high confidence (score 0.7), and cluster analysis was performed by using k-means with a value of k = 3.
Isolation and structural analysis of PcWRI1 and PcLEC1
Total RNA isolated from P. chinensis seeds by RNeasy Plant Mini Kit (Qiagen, Inc., USA) was used to synthesize cDNA by PrimeScript™ RT reagent Kit (TaKaRa Dalian, China). The full-length sequences of PcWRI1 and PcLEC1 were amplified using the primers of PcWRI1-U/PcWRI1-D and PcLEC1-U/PcLEC1-D, respectively (Additional file 11: Table S9), and then subcloned into pGEM-T vector (Promega, USA) for sequencing. Multiple sequence alignment and phylogenetic tree were constructed by DNAMAN 7.0 with ClustalX program and RAxML-NG (https://raxmlng.vital-it.ch/) with maximum likelihood method, respectively . The 3D model for PcWRI1 and PcLEC1 as well as homologous proteins from A. thaliana, Z. may and G. max was constructed by online tool Phyre2 (http://www.sbg.bio.ic.ac.uk/phyre2/html/).
Heterologous expression of PcWRI1 and PcLEC1 in Arabidopsis thaliana
The ORFs of PcWRI1 and PcLEC1 were amplified with the primer pairs PcWRI1-F-EcoRI/PcWRI1-R-BamHI and PcLEC1-F-KpnI/PcLEC1-R-XbaI, and subcloned into pCAMBIA2300 driven by the CaMV35S promoter to construct the vectors of pCAMBIA2300/CaMV35S::PcWRI1 and pCAMBIA2300/CaMV35S::PcLEC1, and then transformed into Agrobacterium tumefaciens strain EHA105 by freeze–thaw method and delivered into wild-type Arabidopsis plants (ecotype Col-0) by floral dip method . The transgenic lines were screened on 1/2MS medium (containing 50 μg/mL kanamycin, 20 g/ L sucrose and 20 g/L agar) and detected by PCR to obtain T3 seeds. Finally, 3 independent-homogenous transgenic lines of PcWRI1 and PcLEC1 were obtained and examined by χ2 test for further analyses.
Analyses of plant phenotypic and seed traits
To explore the effect of PcWRI1 or PcLEC1 overexpression on phenotypes of transgenic lines, a comparative analysis of morphological features (root length, plant height, leaf size and seed weight) was conducted on the seedlings of WT, NT and T3 transgenic Arabidopsis lines (35S::PcWRI1 and 35S::PcLEC1). The root length of 10-day-old seedling germinated on half-strength MS medium was detected. The 20-day-old seedlings were transfer to pots and grown in separate flats in same incubator used for detections of plant height and flowering time. One thousand dry seeds from the controls (WT and NT) and T3 transgenic lines were weighed, and the results was expressed as thousand seed weight (mg/1000 seeds). All the determination was conducted in triplicate.
To harvest Arabidopsis seeds at different developing stages, the flowers on main axis were tagged on the day of anthesis, and then the seeds (10 plants of WT and transgenic lines per set) were collected at 5 developing stages from 4 to 20 DAA with an interval of 4 days. The flowers, leaves, stems and roots were collected at 0 DAA. All samples were frozen in liquid nitrogen and stored at − 80 °C for analyses of oil content, FA compositions, RNA extraction and qRT-PCR.
Construction of prediction model for biodiesel property of raw oils from ripen seeds of transgenic Arabidopsis
Triangular prediction model for biodiesel property of raw oils from ripen seeds of transgenic lines was constructed based on the influence of FA compositions [12, 16]. To predict biodiesel properties of seed oils from transgenic lines, the percentages of SFA, MUFA and PUFA in ripen seed oils from each transgenic line were calculated to outline triangular prediction graph, in which three angular points of the triangle meant the 100% of SFA, MUFA and PUFA, respectively. In triangular graph, the region existed at the far end of the polyunsaturated angular point (lower left vertex) and the saturated angular point (lower right vertex) was delineated to predict the biodiesel fuel properties, taking into account the CN, IV, CFPP and OS [12, 16].
Gene expression analysis by qRT-PCR
Total RNA was extracted from P. chinensis seeds by RNeasy Plant Mini Kits (Qiagen) and quantified by Nanodrop ND-1000 spectrophotometer (N Wilmington, USA) to be reverse-transcribed by Reverse Transcription System (Promega). qRT-PCR was conducted on 7500 Real-Time PCR System using SYBR Premix Ex Taq Kit (TaKaRa). All amplified primers (Additional file 12: Table S10) were designed by PrimerQuest (http://www.idtdna.com/PrimerQuest/Home/Index). The genes encoding large subunit ribosomal protein L32e (RPL32e) and ubiquitin-conjugating enzyme (UBC) were used as inner references [16, 79] for the detection of TF gene transcription in P. chinensis seeds, and its expression level in the seeds of different accessions was set to 1.00 for standardization. As for detections of gene expression in developing seeds of the controls (WT and NT) and transgenic lines, Arabidopsis ACTIN gene was used as endogenous reference, and expression level from seed sample at 4 DAA was set to 1.00 for standardization. Three biological replicates with three technical repetitions each were performed for all qRT-PCR assays.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its Additional files.
ABSCISIC ACID INSENSITIVE3
Cold filter plugging point
Acyl-CoA:DAG acyltransferase 1
Fatty acid (FA) desaturase
Fatty acyl-ACP thioesterase A/B
FA methyl esters
- KAS I/II/III:
3-Ketoacyl ACP synthase I/II/III
Long-chain acyl CoA synthetase
Acyl-CoA:LPA acyltransferase 2
Diacylglycerol (DAG) acyltransferase
Pyruvate (PYR) dehydrogenase complex
Quantitative real time-PCR
18:0-ACP desaturase 6
Bora AP, Gupta DP, Durbha KS. Sewage sludge to bio-fuel: A review on the sustainable approach of transforming sewage waste to alternative fuel. Fuel. 2020;259: 116262.
Musa SD, Zhonghua T, Ibrahim AO, Habib M. China’s energy status: A critical look at fossils and renewable options. Renew Sustain Energy Rev. 2018;81:2281–90.
Benti NE, Gurmesa GS, Argaw T, Aneseyee AB, Gunta S, Kassahun GB, et al. The current status, challenges and prospects of using biomass energy in Ethiopia. Biotechnol Biofuels. 2021;14:209.
Knothe G, Razon LF. Biodiesel fuels. Prog Energy Combust Sci. 2017;58:36–59.
Atabani AE, Silitonga AS, Badruddin IA, Mahlia TMI, Masjuki HH, Mekhilef S. A comprehensive review on biodiesel as an alternative energy resource and its characteristics. Renew Sustain Energy Rev. 2012;16(4):2070–93.
Monirul IM, Masjuki HH, Kalam MA, Zulkifli NWM, Rashedul HK, Rashed MM, et al. A comprehensive review on biodiesel cold flow properties and oxidation stability along with their improvement processes. RSC Adv. 2015;5(105):86631–55.
Graham-Rowe D. Agriculture: Beyond food versus fuel. Nature. 2011;474(7352):S6–8.
Park M-E, Lee K-R, Chen GQ, Kim HU. Enhanced production of hydroxy fatty acids in Arabidopsis seed through modification of multiple gene expression. Biotechnol Biofuels Bioproducts. 2022;15:66.
Yu HY, Zhou S. Preparation of biodiesel from Xanthoceras sorblfolia Bunge seed oil. China Oils Fats. 2009;3:43–5.
Li ZY, Xia F, Zhang LY, Fan CY, Wei GT. Research advances in biodiesel production on from Jatropha oil by chemical methods. China Science Paper. 2013;8(3):221–4.
Tang M, Zhang P, Zhang L, Li M, Wu L. A Potential Bioenergy Tree: Pistacia chinensis Bunge. Energy Procedia. 2012;16:737–46.
Wang J, Lin WJ, Yin ZD, Wang LB, Dong SB, An JY, et al. Comprehensive evaluation of fuel properties and complex regulation of intracellular transporters for high oil production in developing seeds of Prunus sibirica for woody biodiesel. Biotechnol Biofuels. 2019;12:6.
Guo JY, Li HY, Fan SQ, Liang TY, Yu HY, Li JR, et al. Genetic variability of biodiesel properties in some Prunus L.(Rosaceae) species collected from Inner Mongolia, China. Ind Crop Prod. 2015;76:244–8.
Niu J, An JY, Wang LB, Fang CL, Ha DL, Fu CY, et al. Transcriptomic analysis revealed the mechanism of oil dynamic accumulation during developing Siberian apricot (Prunus sibirica L.) seed kernels for the development of woody biodiesel. Biotechnol Biofuels. 2015;8:29.
Wang L, Chu J. Optimization of biodiesel production from Siberian Apricot (Prunus sibirica L.) oil using response surface methodology. Asian J Chem. 2013;25(5):2577.
Lin Z, An J, Wang J, Niu J, Ma C, Wang L, et al. 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. Biotechnol Biofuels. 2017;10:134.
Dong S, Liu Y, Xiong B, Jiang X, Zhang Z. Transcriptomic analysis of a potential bioenergy tree, Pistacia chinensis Bunge, and identification of candidate genes involved in the biosynthesis of oil. BioEnergy Research. 2016;9(3):740–9.
Li H, Zhang Z, Lin S, Li X. Research advances in the study of Pistacia chinensis Bunge, a superior tree species for biomass energy. For Stud in China. 2007;9(2):164–8.
Li XX, Qin HL, Lin SZ, Hou BN, Zhang ZX. Comparative analysis on 100-seed weight and oil content of Pistacia chinensis seed from different provenances. Modern Agricult Sci. 2008;15(6):11–3.
Wang R, Hanna MA, Zhou WW, Bhadury PS, Chen Q, Song BA, et al. Production and selected fuel properties of biodiesel from promising non-edible oils: Euphorbia lathyris L., Sapium sebiferum L. and Jatropha curcas L. Bioresour Technol. 2011;102(2):1194–9.
Bates PD, Stymne S, Ohlrogge J. Biochemical pathways in seed oil synthesis. Curr Opin Plant Biol. 2013;16(3):358–64.
Yeap Wc, Lee Fc, Shan S, Kumar D, Musa H, Appleton DR, et al. WRI1–1, ABI5, NF-YA3 and NF-YC2 increase oil biosynthesis in coordination with hormonal signaling during fruit development in oil palm. Plant J. 2017;91(1):97–113.
Adhikari ND, Bates PD, Browse J. WRINKLED1 rescues feedback inhibition of fatty acid synthesis in hydroxylase-expressing seeds. Plant Physiol. 2016;171(1):179–91.
Bhattacharya S, Das N, Maiti MK. Cumulative effect of heterologous AtWRI1 gene expression and endogenous BjAGPase gene silencing increases seed lipid content in Indian mustard Brassica juncea. Plant Physiol Biochem. 2016;107:204–13.
Elahi N, Duncan RW, Stasolla C. Modification of oil and glucosinolate content in canola seeds with altered expression of Brassica napus LEAFY COTYLEDON1. Plant Physiol Biochem. 2016;100:52–63.
Manan S, Ahmad MZ, Zhang G, Chen B, Haq BU, Yang J, et al. Soybean LEC2 regulates subsets of genes involved in controlling the biosynthesis and catabolism of seed storage substances and seed development. Front Plant Sci. 2017;8:1604.
Ye J, Wang C, Sun Y, Qu J, Mao H, Chua N-H. Overexpression of a transcription factor increases lipid content in a woody perennial Jatropha curcas. Front Plant Sci. 2018;9:1479.
Vogel PA, Bayon de Noyer S, Park H, Nguyen H, Hou L, Changa T, et al. Expression of the Arabidopsis WRINKLED 1 transcription factor leads to higher accumulation of palmitate in soybean seed. Plant Biotechnol J. 2019;17:1–11.
Jo L, Pelletier JM, Harada JJ. Central role of the LEAFY COTYLEDON1 transcription factor in seed development. J Integr Plant Biol. 2019;61(5):564–80.
Han X, Zhang Y-W, Liu J-Y, Zuo J-F, Zhang Z-C, Guo L, et al. 4D genetic networks reveal the genetic basis of metabolites and seed oil-related traits in 398 soybean RILs. Biotechnol Biofuels Bioproducts. 2022;15:92.
Wang LB, Yu HY, He XH, Liu RY. Influence of fatty acid composition of woody biodiesel plants on the fuel properties. J Fuel Chem Technol. 2012;40(4):397–404.
Yaşar F. Comparision of fuel properties of biodiesel fuels produced from different oils to determine the most suitable feedstock type. Fuel. 2020;264: 116817.
Singh D, Sharma D, Soni SL, Sharma S, Kumari D. Chemical compositions, properties, and standards for different generation biodiesels: A review. Fuel. 2019;253:60–71.
Sakthivel R, Ramesh K, Purnachandran R, Mohamed SP. A review on the properties, performance and emission aspects of the third generation biodiesels. Renew Sustain Energy Rev. 2018;82:2970–92.
Zhu K, Li N, Zheng X, Sarwar R, Li Y, Cao J, et al. Overexpression the BnLACS9 could increase the chlorophyll and oil content in Brassica napus. Biotechnol Biofuels Bioproducts. 2023;16:3.
Sun Y, Liu B, Xue J, Wang X, Cui H, Li R, et al. Critical metabolic pathways and genes cooperate for epoxy fatty acid-enriched oil production in developing seeds of Vernonia galamensis, an industrial oleaginous plant. Biotechnol Biofuels Bioproducts. 2022;15:21.
Wang LB. Evaluation of Siberian Apricot (Prunus sibirica L.) germplasm variability for biodiesel properties. J Am Oil Chem Soc. 2012;89(9):1743–7.
Chen YH, Chen JH, Chang CY, Chang CC. Biodiesel production from tung (Vernicia montana) oil and its blending properties in different fatty acid compositions. Bioresour Technol. 2010;101(24):9521–6.
Wang L, Yu H. Biodiesel from Siberian apricot (Prunus sibirica L.) seed kernel oil. Bioresour Technol. 2012;112:355–8.
Seiler GJ, Gulya T, Kong G, Thompson S, Mitchell J. Oil concentration and fatty-acid profile of naturalized Helianthus annuus populations from Australia. Genet Resour Crop Evol. 2018;65(8):2215–29.
Kaushik N, Bhardwaj D. Screening of Jatropha curcas germplasm for oil content and fatty acid composition. Biomass Bioenerg. 2013;58:210–8.
Hoseini SS, Najafi G, Ghobadian B, Mamat R, Ebadi MT, Yusaf T. Ailanthus altissima (tree of heaven) seed oil: Characterisation and optimisation of ultrasonication-assisted biodiesel production. Fuel. 2018;220:621–30.
Amalfitano C, Golubkina NA, Del Vacchio L, Russo G, Cannoniero M, Somma S, et al. Yield, antioxidant components, oil content, and composition of onion seeds are influenced by planting time and density. Plants. 2019;8(8):293.
Liu P, Zhang LN, Wang XS, Gao JY, Yi JP, Deng RX. Characterization of Paeonia ostii seed and oil sourced from different cultivation areas in China. Ind Crop Prod. 2019;133:63–71.
Ramos MJ, Fernández CM, Casas A, Rodríguez L, Pérez Á. Influence of fatty acid composition of raw materials on biodiesel properties. Bioresour Technol. 2009;100(1):261–8.
Ma Y, Wang S, Liu X, Yu H, Yu D, Li G, et al. Oil content, fatty acid composition and biodiesel properties among natural provenances of Siberian apricot (Prunus sibirica L.) from China. GCB Bioenergy. 2021;13(1):112–32.
Mohan MR, Jala RCR, Kaki SS, Prasad R, Rao B. Swietenia mahagoni seed oil: A new source for biodiesel production. Ind Crop Prod. 2016;90:28–31.
Lovato L, Pelegrini BL, Rodrigues J, de Oliveira AJB, Ferreira ICP. Seed oil of Sapindus saponaria L.(Sapindaceae) as potential C16 to C22 fatty acids resource. Biomass Bioenerg. 2014;60:247–51.
Yu HY, Fan SQ, Bi QX, Wang SX, Hu XY, Chen MY, et al. Seed morphology, oil content and fatty acid composition variability assessment in yellow horn (Xanthoceras sorbifolium Bunge) germplasm for optimum biodiesel production. Ind Crop Prod. 2017;97:425–30.
Rodríguez-Rodríguez MF, Sánchez-García A, Salas JJ, Garcés R, Martínez-Force E. Characterization of the morphological changes and fatty acid profile of developing Camelina sativa seeds. Ind Crop Prod. 2013;50:673–9.
Baud S, Mendoza MS, To A, Harscoët E, Lepiniec L, Dubreucq B. WRINKLED1 specifies the regulatory action of LEAFY COTYLEDON2 towards fatty acid metabolism during seed maturation in Arabidopsis. Plant J. 2007;50(5):825–38.
Baud S, Wuillème S, To A, Rochat C, Lepiniec L. Role of WRINKLED1 in the transcriptional regulation of glycolytic and fatty acid biosynthetic genes in Arabidopsis. Plant J. 2009;60(6):933–47.
Maeo K, Tokuda T, Ayame A, Mitsui N, Kawai T, Tsukagoshi H, et al. An AP2-type transcription factor, WRINKLED1, of Arabidopsis thaliana binds to the AW-box sequence conserved among proximal upstream regions of genes involved in fatty acid synthesis. Plant J. 2009;60(3):476–87.
Shen B, Allen WB, Zheng P, Li C, Glassman K, Ranch J, et al. Expression of ZmLEC1 and ZmWRI1 increases seed oil production in maize. Plant Physiol. 2010;153(3):980–7.
Pouvreau B, Baud S, Vernoud V, Morin V, Py C, Gendrot G, et al. Duplicate maize Wrinkled1 transcription factors activate target genes involved in seed oil biosynthesis. Plant Physiol. 2011;156(2):674.
Tan HL, Yang XH, Zhang FX, Zheng X, Qu CM, Mu JY, et al. Enhanced seed oil production in canola by conditional expression of Brassica napus LEAFY COTYLEDON1 and LEC1-LIKE in developing seeds. Plant Physiol. 2011;156(3):1577–88.
To A, Joubès J, Barthole G, Lécureuil A, Scagnelli A, Jasinski S, et al. WRINKLED transcription factors orchestrate tissue-specific regulation of fatty acid biosynthesis in Arabidopsis. Plant Cell. 2012;24(12):5007.
Ivarson E, Leiva-Eriksson N, Ahlman A, Kanagarajan S, Bülow L, Zhu L-H. Effects of overexpression of WRI1 and hemoglobin genes on the seed oil content of Lepidium campestre. Front Plant Sci. 2017;7:2032.
Cernac A, Benning C. WRINKLED1 encodes an AP2/EREB domain protein involved in the control of storage compound biosynthesis in Arabidopsis. Plant J. 2004;40(4):575–85.
Niu J, Hou XY, Fang CL, An JY, Ha DL, Qiu L, et al. Transcriptome analysis of distinct Lindera glauca tissues revealed the differences in the unigenes related to terpenoid biosynthesis. Gene. 2015;599:22–30.
Lin ZX, An JY, Wang J, Niu J, Ma C, Wang LB, et al. 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. Biotechnol Biofuels. 2017;10(1):134.
Jo L, Pelletier JM, Hsu S-W, Baden R, Goldberg RB, Harada JJ. Combinatorial interactions of the LEC1 transcription factor specify diverse developmental programs during soybean seed development. Proc Natl Acad Sci USA. 2020;117(2):1223.
Wan X, Liu Q, Dong B, Vibhakaran Pillai S, Huang F-H, Singh SP, et al. Molecular and biochemical analysis of the castor caruncle reveals a set of unique genes involved in oil accumulation in non-seed tissues. Biotechnol Biofuels. 2019;12:158.
Pelletier JM, Kwong RW, Park S, Le BH, Baden R, Cagliari A, et al. LEC1 sequentially regulates the transcription of genes involved in diverse developmental processes during seed development. Proc Natl Acad Sci USA. 2017;114(32):E6710–9.
Chen B, Zhang G, Li P, Yang J, Guo L, Benning C, et al. Multiple GmWRI1s are redundantly involved in seed filling and nodulation by regulating plastidic glycolysis, lipid biosynthesis and hormone signalling in soybean (Glycine max). Plant Biotechnol J. 2020;18(1):155–71.
Sun J, Chen T, Liu M, Zhao D, Tao J. Analysis and functional verification of PoWRI1 gene associated with oil accumulation process in Paeonia ostii. Int J Mol Sci. 2021;22(13):6996.
Baud S, Lepiniec L. Physiological and developmental regulation of seed oil production. Prog Lipid Res. 2010;49(3):235–49.
Bourgis F, Kilaru A, Cao X, Ngando-Ebongue G-F, Drira N, Ohlrogge JB, et al. Comparative transcriptome and metabolite analysis of oil palm and date palm mesocarp that differ dramatically in carbon partitioning. Proc Natl Acad Sci USA. 2011;108(30):12527–32.
Troncoso-Ponce MA, Kilaru A, Cao X, Durrett TP, Fan J, Jensen JK, et al. Comparative deep transcriptional profiling of four developing oilseeds. Plant J. 2011;68(6):1014–27.
Hofvander P, Ischebeck T, Turesson H, Kushwaha SK, Feussner I, Carlsson AS, et al. Potato tuber expression of Arabidopsis WRINKLED1 increase triacylglycerol and membrane lipids while affecting central carbohydrate metabolism. Plant Biotechnol J. 2016;14(9):1883–98.
Dussert S, Guerin C, Andersson M, Joët T, Tranbarger TJ, Pizot M, et al. Comparative transcriptome analysis of three oil palm fruit and seed tissues that differ in oil content and fatty acid composition. Plant Physiol. 2013;162(3):1337–58.
Lim AR, Kong Q, Singh SK, Guo L, Yuan L, Ma W. Sunflower WRINKLED1 plays a key role in transcriptional regulation of oil biosynthesis. Int J Mol Sci. 2022;23(6):3054.
Lotan T, Ohto M-a, Yee KM, West MAL, Lo R, Kwong RW, et al. Arabidopsis LEAFY COTYLEDON1 is sufficient to induce embryo development in vegetative cells. Cell. 1998;93(7):1195–205.
Mu J, Tan H, Zheng Q, Fu F, Liang Y, Zhang J, et al. LEAFY COTYLEDON1 is a key regulator of fatty acid biosynthesis in Arabidopsis. Plant Physiol. 2008;148(2):1042.
Tang G, Xu P, Ma W, Wang F, Liu Z, Wan S, et al. Seed-specific expression of AtLEC1 increased oil content and altered fatty acid composition in seeds of peanut (Arachis hypogaea L.). Frontiers in Plant Science. 2018;9:260.
Li Y, Beisson F, Pollard M, Ohlrogge J. Oil content of Arabidopsis seeds: The influence of seed anatomy, light and plant-to-plant variation. Phytochemistry. 2006;67(9):904–15.
Kozlov AM, Darriba D, Flouri T, Morel B, Stamatakis A. RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics. 2019;35(21):4453–5.
Clough SJ, Bent AF. Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 1998;16(6):735–43.
Niu J, Chen Y, An J, Hou X, Cai J, Wang J, et al. Integrated transcriptome sequencing and dynamic analysis reveal carbon source partitioning between terpenoid and oil accumulation in developing Lindera glauca fruits. Sci Rep. 2015;5(1):1–12.
We thank the National Natural Sciences Foundation of China (Grant No. 31972952 and 31570653) for financial support. We also like to thank the reviewers for positive criticism to improve the quality of the manuscript.
Research involving plants
We confirm that the experimental material (Pistacia chinensis) is owned by the authors, and deposited in the Forest and Flower Germplasm Resources Genebark of Beijing Forestry University, China.
This work was financially supported by The National Natural Science Foundation of China (Grant No. 31972952 and 31570653).
Ethics approval and consent to participate
The authors declare that all the experimental methods were carried out in accordance with relevant guidelines and regulations.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
“Yu Xiu” made a great contribution (organized and planned this research) to this article.
About this article
Cite this article
Chen, F., Lin, W., Li, W. et al. Determination of superior Pistacia chinensis accession with high-quality seed oil and biodiesel production and revelation of LEC1/WRI1-mediated high oil accumulative mechanism for better developing woody biodiesel. BMC Plant Biol 23, 268 (2023). https://doi.org/10.1186/s12870-023-04267-y