Skip to main content

Intron-mediated enhancement of DIACYLGLYCEROL ACYLTRANSFERASE1 expression in energycane promotes a step change for lipid accumulation in vegetative tissues

Abstract

Background

Metabolic engineering for hyperaccumulation of lipids in vegetative tissues is a novel strategy for enhancing energy density and biofuel production from biomass crops. Energycane is a prime feedstock for this approach due to its high biomass production and resilience under marginal conditions. DIACYLGLYCEROL ACYLTRANSFERASE (DGAT) catalyzes the last and only committed step in the biosynthesis of triacylglycerol (TAG) and can be a rate-limiting enzyme for the production of TAG.

Results

In this study, we explored the effect of intron-mediated enhancement (IME) on the expression of DGAT1 and resulting accumulation of TAG and total fatty acid (TFA) in leaf and stem tissues of energycane. To maximize lipid accumulation these evaluations were carried out by co-expressing the lipogenic transcription factor WRINKLED1 (WRI1) and the TAG protect factor oleosin (OLE1). Including an intron in the codon-optimized TmDGAT1 elevated the accumulation of its transcript in leaves by seven times on average based on 5 transgenic lines for each construct. Plants with WRI1 (W), DGAT1 with intron (Di), and OLE1 (O) expression (WDiO) accumulated TAG up to a 3.85% of leaf dry weight (DW), a 192-fold increase compared to non-modified energycane (WT) and a 3.8-fold increase compared to the highest accumulation under the intron-less gene combination (WDO). This corresponded to TFA accumulation of up to 8.4% of leaf dry weight, a 2.8-fold or 6.1-fold increase compared to WDO or WT, respectively. Co-expression of WDiO resulted in stem accumulations of TAG up to 1.14% of DW or TFA up to 2.08% of DW that exceeded WT by 57-fold or 12-fold and WDO more than twofold, respectively. Constitutive expression of these lipogenic “push pull and protect” factors correlated with biomass reduction.

Conclusions

Intron-mediated enhancement (IME) of the expression of DGAT resulted in a step change in lipid accumulation of energycane and confirmed that under our experimental conditions it is rate limiting for lipid accumulation. IME should be applied to other lipogenic factors and metabolic engineering strategies. The findings from this study may be valuable in developing a high biomass feedstock for commercial production of lipids and advanced biofuels.

Graphical abstract

Background

Production of biofuel from perennial crops will increase energy security by reducing the need for fossil fuels while increasing carbon sequestration for combating climate change [1, 2]. Therefore, many countries have implemented government policies and subsidies promoting the production and consumption of bioenergy [3].

Energycane, is one of the most promising bioenergy feedstocks. Both energycane and sugarcane are highly polyploid, interspecific hybrids between Saccharum officinarum and Saccharum spontaneum, but energycane has a higher proportion of Saccharum spontaneum in its genome. This contributes to greater biomass yields, elevated resilience under marginal conditions, lower sugar, and higher fiber content in the stems of energycane relative to sugarcane [4,5,6,7]. Energycane’s superior resilience and biomass accumulation are supported by a more vigorous root system and higher light interception efficiency than sugarcane. The elevated light interception efficiency is associated with a higher tiller density and leaf inclination angles, contributing to the greatest photosynthetic efficiency among crops [8]. Photosynthesis, which transforms sunlight into stored chemical energy creates biomass that can be converted to renewable biofuel. Plant oil, consisting mainly of triacylglycerol [TAG; (98%)] and other fatty acids [9, 10] has more than twice the energy content of carbohydrates [11, 12]. Oil is typically produced from seeds of palm and oilseed crops such as soybean (Glycine max), sunflower (Helianthus annus) and canola (Brassica napus) [13]. Vegetative plant tissues including leaves, stems and roots are able to produce lipids and TAG in small quantities [14, 15]. The genes involved in the fatty acid and TAG synthesis, and hydrolysis have been identified [14, 16]. This enables metabolic engineering of high biomass crops including sugarcane and energycane for the hyperaccumulation of lipids in their vegetative biomass and has the potential to exceed oil yields of traditional oilseed crops [12, 17,18,19,20]. Ectopic expression of individual lipogenic factors led to modest lipid accumulation in vegetative tissues. A step change in accumulation of vegetative lipids has been reported following multigene engineering termed “push–pull-protect strategy” including genes involved in fatty acid biosynthesis (‘push’), TAG assembly (‘pull’), and suppression of lipid turnover (‘protect’) [21]. This approach optimized the flux of carbon into TAG at multiple metabolic levels, leading to hyperaccumulation of TAG in the vegetative tissues of model [22,23,24,25] and high biomass crops such as ryegrass [26], sugarcane [18, 20], sorghum [27], and potato [28]. WRI1 is a transcription factor of the family of APETALA2/ethylene-responsive element binding proteins involved in the regulation of fatty acid biosynthesis [2933]. DGAT1 catalyzes the only committed step in the biosynthesis of TAG by catalyzing the transfer of an acyl group from acyl-CoA to diacylglycerol [34]. DGAT1 has been classified into three types; two endoplasmic reticulum (ER)-localized forms of DGAT1 (type-1) and DGAT2 (type-2) and a cytosolic form of the DGAT family (Type-3) [35, 3638]. Overexpression of DGAT1 alone has been reported to increase the TAG content in the transgenic leaves by 20-fold compared to control [39]. Protection of TAG from turnover by lipases has been achieved by coating oil bodies with oleosin [22]. OLEOSIN1 (OLE1) is a structural protein that protects oil bodies from coalescing and reduces lipid breakdown by lipases [40, 41]. Previous research has shown increased fatty acid content of Arabidopsis vegetative tissues through the engineering of OLE1 (CysOLE1) by the substitution of six cysteine residues for their native equivalents within the N- and C-terminal hydrophilic domains [42].

We earlier reported that in the monocot Lemna japonica, co-expression of WRI1, DGAT1 and OLE1 had a much larger effect on TAG accumulation than co-expression of WRI1 and OLE1 without DGAT1, suggesting that the elevated FA synthesis due to ectopic WRI1 expression exceeds the capacity of the native DGAT to transfer FA to TAG [43]. Elevating the expression of DGAT is therefore a promising approach for hyperaccumulation of TAG in vegetative tissues. Various optimizations including stronger promoters or enhanced UTR’s, codon optimization or the presence of intron(s) have been described to elevate transgene expression in plants [44]. Introns can have a surprisingly large effect on gene expression, revealing a gap in our understanding of the molecular basis underlying this phenomenon [4547].

In this study, we found that the addition of an intron from the predicted 5-methyltetrahydropteroyltriglutamate-homocysteine methyltransferase 1 of Sorghum bicolor into a codon-optimized TmDAGT1 from nasturtium dramatically increased the gene expression under control of the strong constitutive Panicum virgatum ubiquitin promoter (pPvUbiII). The elevated TmDGAT1 expression displayed strong synergy in TAG and TFA accumulation following the co-expression of WRI1, DGAT1, and OLE1 in energycane.

Results

Cis-regulatory elements in the 110-bp intron

A putative GT1 consensus and CAAT box were identified at the position of 52 bp and 93 bp, respectively, in the 110-bp-long intron that we inserted into DGAT1 and which originated from the predicted 5-methyltetrahydropteroyltriglutamate-homocysteine methyltransferase 1 of Sorghum bicolor (NCBI accession number: NC_012877.2) (Figure S1). Sequences with similarity to motifs that were earlier reported in IME including YATCN, the reverse complement of NGATY and similar sequence to TTNGATYTG and its reverse complement CAYATCNAA were also identified.

Overexpression of WRI1, DGAT1 and OLE1 altered TAG and TFA production in leaves and stems of transgenic energycane

WRI1, DGAT1, and OLE1 expression was analyzed from leaf and stem tissues by qRT-PCR with primers shown in Table S1. To determine potential intron-mediated regulation of gene expression, we analyzed the expression in five PCR-positive lines for DGAT1 with intron (DGAT1 (In)) and compared this with the expression of in five PCR-positive lines for DGAT1 without intron (DGAT1 (W)). DGAT1 (In) elevated expression in leaves by seven times on average of 5 transgenic lines compared to DGAT1 (W) (Table 1). Expression in stem was two times higher for DGAT1 (In) compared to DGAT1 (W) (Table 2).

Table 1 Summary of TAG and TFA content and transgene expression in leaves of energycane
Table 2 Summary of TAG and TFA content and transgene expression in stems of energycane

As shown in Table 1, 8 out of 10 transgenic lines co-expressed all three lipogenic factors with significant differences in transgene expression between the lines. Two of the transgenic lines (L4; L10, Table 1) expressed DGAT1 and OLE1 but not WRI1. Among the three transgenes, WRI1 expression showed the lowest variation ranging from 0.01 to 0.05 relative to GAPDH and averaged 0.02 across the five DGAT1 (W) or DGAT1 (In) lines (Table 1). In WDO combination, TAG accumulated up to 1.01% of leaf DW and 0.22% of stem DW, which was 50-fold and tenfold of that of WT, respectively (Tables 1, 2). The highest TAG accumulation was observed in WDiO combination with up to 3.85% of leaf DW and 1.14% of stem DW, which was 192-fold and 56-fold of that of WT, respectively (Tables 1, 2). Similarly, TFA was the highest in WDiO combination with up to 8.39% of leaf DW and 2.08% of stem DW (Tables 1, 2). In addition, Pearson’s correlation was evaluated for TAG accumulation and lipid gene expression levels in transgenic plants. As shown in Table S2, TAG accumulation exhibited the highest positive correlation with DGAT1(W) (r = 0.97) followed by OLE1 (W) (r = 0.89), DGAT1 (In) (r = 0.88), and WRI1 expression (r = 0.51).

Alteration of FA accumulation in leaves and stem of high TAG transgenic energycane

TAG FA composition from leaf (Fig. 1) and the stem (Fig. 2) were analyzed from the three highest oil accumulating lines (L11 – L13). Palmitic acid (PA, 16:0), oleic acid (OA, 18:1), stearic acid (SA, 18:0), linoleic acid (LA, 18:2) and α-linolenic acid (ALA, 18:3) were the major fractions. In both leaf and stem samples (Figs. 1 and 2), levels of unsaturated FA, OA (18:1), and LA (18:2) were significantly increased in the transgenic lines at the expense of saturated FAs, PA (16:0) for TAG in leaves (Fig. 1a) and stems (Fig. 2a), as well as for TFA in leaves (Figs. 1b) and stems (Fig. 2b). Additionally, in the leaf, ALA (18:3) accumulation increased for TAG (Fig. 1a) and was reduced for TFA (Fig. 1b), whereas SA (18:0) accumulation was marginally lowered for TAG (Fig. 1a) and increased for TFA (Fig. 1b). In the stem, the OA (18:1) was the highest for both TAG (Fig. 2a) and TFA (Fig. 2b), in contrast to the FA composition of the leaves, and the SA (18:0) was reduced in both TAG (Fig. 2a) and TFA (Fig. 2b) of the stem. Interestingly, ALA (18:3) was present in transgenic lines but not in wild-type lines (Fig. 2).

Fig. 1
figure 1

TAG FA composition and total FA composition in leaves of different transgenic energycane lines (L11; L12; L13) compared with non-modified energycane (WT). A TAG composition in leaves of different transgenic lines (L11; L12; L13) and non-modified energycane (WT). B TFA composition in leaves of different transgenic lines (L11; L12; L13) and WT. Vertical bars are the means ± standard error (n = 3). Values with different letters are significantly different at p ≤ 0.05 according to a one-way ANOVA test and the Duncan’s multiple range test (MRT). TAG: triacylglycerol. FA: fatty acid. WT: non-modified energycane (WT)

Fig. 2
figure 2

TAG FA composition and total FA composition in stems of different transgenic lines (L11; L12; L13) compared with non-modified energycane (WT). A TAG composition in stems of different transgenic lines (L11; L12; L13) and non-modified energycane (WT). B TFA compositions in stems of different transgenic lines (L11; L12; L13) and WT. Vertical bars are the means ± standard error (n = 3). Values with different letters are significantly different at p ≤ 0.05 according to a one-way ANOVA test and the Duncan’s multiple range test (MRT). TAG triacylglycerol. FA fatty acid. WT non-modified energycane (wild-type)

Accumulation of lipid droplets in leaves of transgenic energycane with hyperaccumulation of TAG in comparison to WT controls

For the purpose of visualizing lipid droplets, leaves from both the transgenic line L13 and its corresponding WT were subjected to boron-dipyrromethene (BODIPY) staining, followed by imaging using confocal microscopy. Notably, the transgenic line L13, with the overexpression of DGAT1 (In), demonstrated a substantial increase in the abundance of BODIPY-stained lipid droplets in comparison to its WT counterpart (Fig. 3). Intriguingly, within the transgenic line L13, the lipid droplets, encompassing both regular lipid droplets (LD) and guard cell lipid droplets (GCL), exhibited a widespread distribution. The most prominent droplet size, around 8 μm in diameter, was predominantly situated within the stomatal guard cells. In contrast, the wild-type leaves exclusively exhibited lipid droplets within stomatal guard cells, typically measuring less than 5 μm in diameter. These findings correlate well with the triacylglycerol (TAG) data depicted in Table 1.

Fig. 3
figure 3

Lipid droplets within the leaf tissue of energecane. The upper panel showcases fluorescence-stained lipid droplets within guard cells (GCL) using BODIPY 493/503, along with the autofluorescence emanating from chloroplasts in the wild-type (WT) strain. The lower panel displays both the lipid droplets within guard cells and those in regular cells (LD), alongside the chloroplast autofluorescence found in leaf tissue of the transgenic line L13. The lipid droplets within guard cells are indicated by white arrows, while a selection of typical lipid droplets are indicated by blue arrows. Both the lipid droplets found within the guard cells (GCL) and those present outside of the guard cells (LD) are rendered in shades of green, while autofluorescence originating from chlorophyll is is seen in red. Scale bar = 50 μm

Impact of hyperaccumulation of TAG on plant growth

To determine the effect of hyperaccumulation of TAG on plant growth, biomass related traits, including plant height, stem diameter, and tiller number were recorded. Hyperaccumulation of TAG corresponded to a reduction in plant height, tiller number, and stem diameter Table 3. Line 13 with the highest TAG accumulation had the greatest reduction in plant height, number of tillers and stem diameter with 31%, 81% and 50%, respectively. In addition, plant height, tiller number and stem diameter correlated negatively with TAG content (r = − 0.71, − 0.77 and − 0.50, respectively) (Additional file 2: Table S3).

Table 3 Height, stem diameter and number of tillers of the transgenic energycane plants (L11; L12; L13) and non-modified energycane (WT) under greenhouse conditions

Discussion

Energycane is among the most promising perennial feedstocks to contribute to the emerging bioeconomy due to its high biomass production, yield resilience and ease of biocontainment [4,5,6, 48]. Metabolic engineering for hyperaccumulation of lipids in vegetative tissues has been proposed for enhancing energy density and biofuel production from high biomass crops like energycane [12, 17, 20, 21]. Metabolic engineering of energycane is still in its infancy due to the recalcitrance of this crop to tissue culture and genetic transformation [17]. Here we report a step change in hyperaccumulation of triacylglycerol (TAG) and total fatty acid (TFA) in energycane using intron-mediated enhancement (IME) of the expression of DGAT1. Line-to-line variation in transgene expression in different transgenic events has been reported due to epigenetic changes [49, 50], position of transgene insertion in the genome [51, 52] and copy number of the transgene [5355]. Therefore, several transgenic events need to be analyzed per recombinant DNA construct to support conclusions on the impact of specific approach/construct. In this study, five lines were evaluated for each, TmDGAT1 with and without intron. Including an intron in TmDGAT1 (TmDGAT1 (In)) increased its expression sevenfold compared to the intron-less variant. Combination of TmDGAT1 (In) expression with SbWRI1 and SiOLE1 expression elevated TAG accumulation 3.8-fold compared to the intron-less combination. Our findings also confirm that DGAT1 is a rate-limiting enzyme for TAG accumulation in vegetative tissues of energycane.

The potential of introns to elevate transgene expression by up to 100-fold, termed intron-mediated enhancement (IME) has been described more than 30 years ago [45, 5661]. IME can strongly stimulate mRNA accumulation far downstream of the transcription start site, even in the absence of a promoter, revealing gaps in our understanding of the molecular mechanisms governing gene expression [62]. Identification of introns that confer IME requires systematic evaluation in transgenic cells since not all introns lead to transcription enhancement [62]. Introns can differ in the presence of enhancers [63, 64] or transcription factor binding sites [65]. In addition, splicing can influence transcription by affecting the phosphorylation state of RNA polymerase II [66]. Intron-dependent gene looping allows physical interaction of the promoter and the terminator regions accelerating mRNA processing events including capping and polyadenylation [6769]. Disposition of the exon junction complex proteins also enhances mRNA stability, export and translation [7072].

The impact of introns on transgene expression is influenced by a variety of factors including position of introns [7375], length of exons [76], specific intronic sequence elements [63, 74, 77], splicing efficiency [74, 76] and their ability to promote the formation of a multi-looped gene architecture [47, 74]. Therefore, not all introns lead to elevated gene expression [78]. Recent studies have identified specific cis-elements important for IME expression [46, 7983]. The intron that was used in this study is constitutively spliced and contains the NGATY core of the longer TTNGATYTG motif and two additional motifs that are similar to TTNGATYTG [46, 83]. These cis-elements could play a role in IME, in addition to reduced exon length and supporting transcript processing and nuclear export with the constitutive splicing mechanism. The online transgene design tool “Intronserter”: https://bibiserv.cebitec.uni-bielefeld.de/intronserter supports systematic intron insertion into transgenes for reduction of exon length, along with codon optimization and removal of undesired sequence elements. The program is fully customizable for individual eukaryotic target organisms and should support the systematic design of transgenes for IME in support of metabolic engineering strategies [84].

In the current study, by co-expression of WRI1, DGAT1 (In), and OLE1, we showed an average TAG hyperaccumulation of 3.85% of leaf DW. This TAG content was 192-fold higher than in WT energycane and 2.5-fold higher than reported recently for energycane [17]. The hyperaccumulation of total fatty acid was also elevated to 8.4% of leaf DW which is almost 56 times that of WT and 1.7-fold than in the first report of metabolic engineering of energycane [17]. High level of TAG and TFA accumulation in vegetative tissues of plants typically requires co-expression of several lipogenic factors involved in upregulation of fatty acid biosynthesis, TAG biosynthesis and its protection from degradation [18, 20, 23, 27]. Constitutive co-expression of WRI1, TmDGAT1 (In), and OLE1 in this study led to 2.5-fold higher TAG and 1.9-fold higher TFA accumulation in energycane than using ZmDGAT1-2 without intron and OLE1 in combination with RNAi co-suppression of the TAG lipase sdp1 and tgd1, a lipid transporter from the endoplasmic reticulum (ER) to chloroplasts [17]. We expect the optimization of expression cassettes using IME for all lipogenic factors, codon optimization, stacking of additional lipid genes in combination with gene editing approaches [8589] will further elevate oil accumulation in energycane.

Reported TFA contents following ectopic co-expression of WRI1, DGAT1, and OLE1, in grain sorghum [27] (up to 9.7% of DW) or tobacco [24] (up to 17.7% of DW) exceed the levels of TFA reported here for energycane (8.4% of leaf DW) by 1.16- to 2.11-fold, respectively. However, biomass yield for energycane are typically in the range of 41 to 49 t ha−1 in the Southeastern USA [90, 91] which is 3- to 12-fold higher than typical biomass yields from grain sorghum (4 to 12 t ha−1) [92, 93] and 3- to 9-fold higher than typical tobacco yields (5.5 to 13 t ha−1) [94, 95]. This should lead to higher or similar oil yields from the energycane described here compared to the alternative crops. Biodiesel derived from lipid accumulating energycane is an attractive approach to minimize the use of food and feed crops in the production of renewable transportation fuels. Additional feedstocks and engineering strategies will be needed to establish a robust pipeline for production of sustainable aviation fuel [96].

Coordinating the expression of WRI1 with stem specific or developmentally regulated promoters active late in plant development [97, 98] or inducible promoters [43, 99] may mitigate the negative impact on biomass accumulation due to constitutive expression of WRI1 and other lipogenic factors described in this and several other studies [18, 25, 27].

The overall TAG and TFA compositions were significantly altered by the co-expression of the lipogenic factors in energycane. Most notable was the increase in oleic acid (18:1) and linoleic acid (18:2) in both, TAG and TFA of transgenic leaves and stems at expense of saturated fatty acids. Similar findings were earlier reported for sugarcane by Parajuli et al. [18]. A significant reduction of alpha linoleic acid (18:3) was only observed in the TFA of leaves. Interestingly in stem tissues oleic acid (18:1) linoleic acid (18:2) and alfa-linolenic acids (18:3) were significantly elevated in both TAG and TFA which differs from the observation in sugarcane using different lipogenic factors [18]. The elevated amount of unsaturated FA in vegetative tissues of energycane should lead to improved fuel stability and cold flow [100].

The presence of leaf GCL droplets as seen in guard cells of WT have been reported for more than three decades [101], where they serve as an energy source to power stomatal opening [102]. Subsequent work on Arabidopsis and the lycophyte Selaginella provide evidence that TAG breakdown is an evolutionarily conserved mechanism in light-induced stomatal opening [103]. In the oil accumulating L13, constitutive expression of WRI1, DGAT1 and OLE1 resulted in a substantial boost in GCL abundance, and LDs were generally visible in surrounding tissue, that were absent in WT controls. The pattern of widespread LD accumulation in transgenic energycane leaves is consistent with reports from transgenic sugarcane [18], tobacco [21], and the duckweed Lemna japonica [43].

Conclusions

Energycane is a favored feedstock for the production of bioproducts and advanced renewable fuels because of its high biomass production and resilience under marginal conditions. In this study, we demonstrated that DGAT1 is the principal rate-limiting step for TAG and TFA accumulation in energycane leaves. IME of DGAT1 caused a step change in its expression and its co-expression with WRI1 and OLE1 elevated TAG or TFA content by 192- or 56-fold relative to levels found in the leaves of unmodified energycane, respectively. These results lay the foundation for the commercial production of biodiesel and other FA derivatives and reveal opportunities for further optimizations of high biomass energycane.

Materials and methods

Design and synthesis of the optimized DGAT1 gene

For codon optimization, the native DGAT1 gene from Tropaeolum majus (Gene bank AY084052) was optimized by replacing the preferred codons with synonymous codons of Sorghum bicolor, a close relative of sugarcane, without changing the encoded amino acids (Additional file 1: Figure S1). Two optimized DGAT1, including DGAT1(W) without intron and DGAT1 (In) with the 110-bp-long intron from the predicted 5-methyltetrahydropteroyltriglutamate-homocysteine methyltransferase 1 of Sorghum bicolor (NCBI Accession number: NC_012877.2), were synthesized by GenScript (Piscataway NJ, USA). The specific intron was selected based on its compact size and the presence of cis-regulatory elements that may impact IME (Additional file 1: Figure S1).

Construction of the multigene expression vectors

Two different optimized DGAT1 genes obtained from Tropaeolum majus, including DGAT1(W) without intron and DGAT1 (In) with a 110-bp-long intron were compared in this study. Using the Golden Gate cloning strategy [104, 105], three multigene expression vectors were created (Fig. 4). The first vector, shown in Fig. 4A, contains WRI1 from Sorghum bicolor controlled by the Brachypodium distachyon elongation factor 1-α promoter (pBdEF1α) and Panicum virgatum ubiquitin (tPvUbiII) terminator, and the selectable marker gene neomycin phosphotransferase (nptII) driven by the maize ubiquitin promoter (pZmUbi) and Sorghum bicolor heat shock protein terminator (tSbHSP; Fig. 4A). Vectors shown in Fig. 4B and C contained expression cassettes for two genes: OLE1, obtained from Sesamum indicum, was under the control of the Brachypodium distachyon ubiquitin promoter (pBdUbi10) and Nicotiana benthamiana actin 3’ UTR terminator (tNbACT3), and DGAT1 from Tropaeolum majus, with or without intron was cloned under the control of the Panicum virgatum ubiquitin promoter (pPvUbiII) and Nicotiana tabacum extensin terminator (tNtEU; Fig. 4B and C).

Fig. 4
figure 4

Recombinant DNA constructs used in energy cane transformation. A Expression cassettes of WRI1 from Sorghum bicolor controlled by the Brachypodium distachyon elongation factor 1-α promoter (pBdEF1α) and Panicum virgatum ubiquitin (tPvUbiII) terminator, and the selectable marker gene neomycin phosphotransferase (nptII) driven by the maize ubiquitin promoter (pZmUbi) and Sorghum bicolor heat shock protein terminator (tSbHSP). B Expression cassettes for OLE1, obtained from Sesamum indicum, was under the control of the Brachypodium distachyon ubiquitin promoter (pBdUbi10) and Nicotiana benthamiana actin 3ʹ UTR terminator (tNbACT3), and DGAT1 from Tropaeolum majus, without intron was cloned under the control of the Panicum virgatum ubiquitin promoter (pPvUbiII) and Nicotiana tabacum extensin terminator (tNtEU). C Expression cassettes for OLE1, obtained from Sesamum indicum, was under the control of the Brachypodium distachyon ubiquitin promoter (pBdUbi10) and Nicotiana benthamiana actin 3ʹ UTR terminator (tNbACT3), and DGAT1 from Tropaeolum majus, including a 110-bp-long intron from the predicted Sorghum bicolor 5-methyltetrahydropteroyltriglutamate-homocysteine methyltransferase 1 was cloned under the control of the Panicum virgatum ubiquitin promoter (pPvUbiII) and Nicotiana tabacum extensin terminator (tNtEU) W: no intron. In: intron

Tissue culture and genetic transformation of energycane

For this study, the energy cane genotype UFCP 84-1047 (WT) was chosen for its superior biomass and resilience [106]. Immature leaf whorl cross-sections (Fig. 5A) were cultured to induce embryogenic callus (Fig. 5B) as described by Fouad et al., [107]. The media for callus induction, transformation, selection, and shoot and root regeneration were prepared according to Fouad et al. [107].

Fig. 5
figure 5

Generation of transgenic energy cane. A Leaf whorl cross-sections were used as explants. B Callus induced from leaf whorl cross-section was used as the target for biolistic gene transfer. C Calli regenerating on medium with geneticin for selection of nptII expression. D Regenerated plants before transferring to soil. E Transgenic energy cane lines growing in greenhouse. Bar = 5 mm

Prior gene transfer, plasmids were digested using the restriction enzyme AscI to remove the backbone. Linearized fragments were gel extracted and purified according to Fouad et al. [107]. Two different combinations including WDO (WRI1 + nptII: DGAT1(W) + OLE1), and WDiO (WRI1 + nptII: DGAT1(In) + OLE1) were delivered to embryogenic callus by biolistic gene transfer. For each combination, a total of 10 shots of DNA were delivered to callus in a 1:2 molar ratio (WRI1 + nptII: DGAT1 + OLE1) using a Biolistic PDS-1000/He particle delivery system (Biorad, Hercules, California, USA) as described earlier [108]. Putative transgenic plants (Fig. 5C, D) were regenerated following selection on geneticin-containing culture medium as described by Fouad et al., 2015 [107]. Rooted plantlets were transferred to soil and grown under greenhouse conditions under natural photoperiod and 25-29ºC day and 20 to 24 °C night temperature (Fig. 5E).

Cis-regulatory elements in the 110-bp intron

PLACE and PlantCARE software were used to evaluate the presence of potential cis-regulatory elements in the 110-bp-long intron from the predicted 5-methyltetrahydropteroyltriglutamate-homocysteine methyltransferase 1 of Sorghum bicolor (NCBI accession number: NC_012877.2). We also searched for sequence motifs that were earlier reported for their contribution to IME including sequence TTNGATYTG and its reverse complement CAYATCNAA as well as NGATY (or its reverse complement YATCN) [76].

PCR analysis of genomic DNA extracts

Genomic DNA was extracted from leaves using the cetyl trimethyl ammonium bromide (CTAB) method [109]. 100 ng of genomic DNA was used as a template in a 20-μL PCR reaction. PCR amplification was conducted using Hot Start Taq DNA Polymerase (NEB, Ipswich, Massachusetts, USA) for WRI1, nptII, OLE1, and DGAT1 individually following the manufacturer’s protocol. Primers for each target gene are listed in Additional file 2: Table S1. Five transgenic lines in each combination were confirmed by PCR for WRI1, DGAT1, OLE1, and nptII (Additional file 1: Figure S2).

Greenhouse propagation of transgenic, lipid accumulating energycane

In the greenhouse, transgenic energycane (Fig. 5E) was propagated by nodal stem cuttings to obtain biological replicates and three plants per transgenic line and non-transgenic plants were each planted in a pot with a 15 cm diameter containing potting mix (Jolly Gardener C/G, Oldcastle Lawn and Garden, Atlanta GA). Plants were irrigated and fertilized daily with a drip fertigation system. In the greenhouse, the temperature was controlled by evaporation cooling to 25–29 °C during the day and 20–24 °C during the night using natural photoperiod with a maximum daily light intensity of 1000 to 1500 μmol m−2 s−1.

Sampling of leaves and stem for quantitative real-time PCR (qRT-PCR), TAG and FA analysis

Leaves from 15 transgenic lines were sliced in half at the midrib; one half was used for qRT-PCR analysis, while the other half was used for TAG and FA analyses. Stem internodes from the three highest TAG lines were also sampled and divided into two halves before grinding in a Retsch cryo bead mill (Verder Scientific, Haan, Germany). For qRT-PCR analysis, samples were flash-frozen in liquid nitrogen and maintained at − 80 °C for later RNA extraction. For TAG and FA analyses, each sample consisted of at least 100 mg of leaf or stem fresh weight and was freeze-dried for 72 h in a lyophilizer (Labconco, MO, USA). For TAG and FA analyses samples were shipped to Brookhaven National Laboratory on dry ice.

qRT-PCR analysis for transgene expression

TRIzol reagent was used to extract total RNA from leaf/stem samples (Life Technologies, Thermo Fisher Scientific, Waltham, MA, USA). A High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, USA) was used to synthesize cDNA from approximately 1.0 µg of total RNA from each sample. DGAT1, OLE1, and WRI1 expressions were measured using primers shown in Additional file 2: Table S1. The glyceraldehyde 3-phosphate dehydrogenase (GAPDH) gene was used as a housekeeping gene for transcript normalization [110]. qRT-PCR was performed in a Real-Time PCR Detection System (Biorad, Hercules, CA, USA) using SsoAdvanced Universal SYBR green supermix (Biorad, Hercules, CA, USA). The formula 2{Ct (GAPDH) – Ct (transgene)}, was used to determine the relative transcription levels of WRI1, DGAT1, and OLE1.

Analysis of TAG and total FA composition

TAG and FA analyses were conducted according to Zale et al., 2016 [20]. In brief, 10.0 mg of lyophilized leaf/stem tissue was extracted with 700 µL of extraction solution containing 2:1:0.1 by volume methanol:chloroform:formic acid. After 1 h of mixing on a vortex mixer and then standing at RT overnight, extract total lipids were as the lower phase by adding 350 μL of 1 M KCl and 0.2 M H3PO4 followed by briefly vortexing and centrifuging at 2000 g for 10 min Total lipids were separated by thin layer chromatography (TLC) using a mixture of hexane, diethyl ether and acetic acid (70:30:1, by volume) as mobile phase and TLC Silica gel 60 (EMD Millipore, Cat. No. 1056260001). By spraying 0.05% primuline (in 80% by volume of acetone), lipids were visualized under UV light. TAG fractions were scraped from the TLC plate after identifying by comparing mobility with that of authentic TAG standards. After being incubated in 1.0 mL of boron trichloride–methanol for 40 min at 80–85 °C, TAG fractions were transmethylated into FA methyl esters (FAMEs). FAMES were extracted into hexane and dried under a nitrogen stream before being dissolved in 100 μl hexane and analyzed by GC–MS with an Agilent Technologies (Santa Clara, CA, USA), 7890A GC System equipped with an Agilent 60 m DB23 capillary column (ID 250-μm), and a 5975C mass selective detector in full mass scan mode. The oven temperature was ramped from 100 °C to 240 °C at a rate of 15 °C min−1 and held at 240 °C for 3 min with a flow rate of 1.2 ml min−1. Following incubation in 1.0 mL of boron trichloride–methanol, total lipid extracts were immediately transmethylated into FAMES for total FA analysis. Using gas chromatography–mass spectrometry, FAMES was dissolved in 100 µL hexane and measured. An internal standard of 5.0 g of C17:0 was employed.

Visualization of lipid droplets with confocal microscopy

For the visualization of lipid droplets, the third leaf from the non-necrotic tip region of both the transgenic line L13 and its corresponding WT was carefully excised. Subsequently, these leaf samples underwent fixation using a FAA buffer (4% formalin, 5% glacial acetic acid, 50% ethanol, 41% water, V/V) under vacuum conditions for a duration of 1 h. Following fixation, the samples were stained with a solution containing 12 μg/mL of BODIPY 493/503 (Invitrogen, Eugene, OR, USA) and 0.1% Triton X-100, employing vacuum assistance for 20 min. The imaging process was carried out using a Leica TCS SP5 laser scanning confocal microscope, employing an excitation wavelength of 496 nm for BODIPY and capturing emissions within the range of 505–583 nm. Additionally, chlorophyll autofluorescence was captured within the wavelength range of 661–800 nm.

Analysis of sequence elements

The intron sequence was analyzed in PLACE database (htpp://www.dna.affrc.go.jp/PLACE/signalup.html/, accessed on 14th August) [111] and PlantCARE database (https://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 14th August 2022) [112].

Statistical analysis

The TAG content, total FA content, biomass, and qRT-PCR analysis data were expressed as means ± standard error (SE). To conduct statistical analysis of the means, one-way ANOVA was completed using SPSS version 20.0 (Armonk, NY, USA). A value of p ≤ 0.05 was considered statistically significant. Pearson’s correlation coefficient was also calculated using SPSS. For each statistical analysis, three independent biological replications were used.

Availability of data and materials

All data generated or analyzed during this study are included in this published article (and its Additional information files).

References

  1. Khounani Z, Nazemi F, Shafiei M, Aghbashlo M, Tabatabaei M. Techno-economic aspects of a safflower-based biorefinery plant co-producing bioethanol and biodiesel. Energy Convers Manag. 2019;201:112184. https://doi.org/10.1016/j.enconman.2019.112184.

    Article  CAS  Google Scholar 

  2. Field JL, Richard TL, Smithwick EAH, Cai H, Laser MS, LeBauer DS, et al. Robust paths to net greenhouse gas mitigation and negative emissions via advanced biofuels. PNAS. 2020;117:21968–77. https://doi.org/10.1073/pnas.1920877117.

    Article  CAS  PubMed Central  Google Scholar 

  3. Kilian L, Zhou X. The impact of rising oil prices on US inflation and inflation expectations in 2020–23. Energy Econ. 2022;113:106228. https://doi.org/10.2139/ssrn.3980337.

    Article  Google Scholar 

  4. Alexander AG. The energy cane alternative. 1985. https://doi.org/10.1016/0144-4565(86)90046-6.

  5. Carvalho-Netto OV, Bressiani JA, Soriano HL, Fiori CS, Santos JM, Barbosa GVS, et al. The potential of the energy cane as the main biomass crop for the cellulosic industry. Chem Biol Technol Agric. 2014;1:1–8. https://doi.org/10.1186/s40538-014-0020-2.

    Article  CAS  Google Scholar 

  6. Diniz AL, Ferreira SS, Ten-Caten F, Margarido GRA, Dos Santos JM, de Barbosa GVS, et al. Genomic resources for energy cane breeding in the post genomics era. Comput Struct Biotechnol J. 2019;17:1404–14. https://doi.org/10.1016/j.csbj.2019.10.006.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Knoll JE, Anderson WF, Missaoui A, Hale A, Hanna WW. Biomass production and stability of five energycane cultivars at two latitudes in Georgia. Agrosyst Geosci Environ. 2021;4:e20146. https://doi.org/10.1002/agg2.20146.

    Article  CAS  Google Scholar 

  8. Cruz LP, Pacheco VS, Silva LM, Almeida RL, Miranda MT, Pissolato MD, et al. Morpho-physiological bases of biomass production by energy cane and sugarcane: a comparative study. Ind Crops Prod. 2021;171:113884. https://doi.org/10.1016/j.indcrop.2021.113884.

    Article  CAS  Google Scholar 

  9. Gurr MI. The biosynthesis of triacylglycerols. Lipids: structure and function. Elsevier. 1980; 205–48. https://doi.org/10.1016/b978-0-12-675404-9.50014-x.

  10. Williams KA. Oils, fats and fatty foods. Their practical examination.(Bolton & Revis). 1949. https://doi.org/10.1002/jctb.5000471829.

  11. Miller R, Durrett TP, Kosma DK, Lydic TA, Muthan B, Koo AJK, et al. Altered lipid composition and enhanced nutritional value of Arabidopsis leaves following introduction of an algal diacylglycerol acyltransferase 2. Plant Cell. 2013;25:677–93. https://doi.org/10.1105/tpc.112.104752.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Ohlrogge J, Chapman K. The seeds of green energy: expanding the contribution of plant oils as biofuels. Biochem. 2011;33:34–8. https://doi.org/10.1042/bio03302034.

    Article  Google Scholar 

  13. Harwood JL, Ramli US, Tang M, Quant PA, Weselake RJ, Fawcett T, et al. Regulation and enhancement of lipid accumulation in oil crops: the use of metabolic control analysis for informed genetic manipulation. Eur J Lipid Sci Technol. 2013;115:1239–46. https://doi.org/10.1002/ejlt.201300257.

    Article  CAS  Google Scholar 

  14. Durrett TP, Benning C, Ohlrogge J. Plant triacylglycerols as feedstocks for the production of biofuels. Plant J. 2008;54:593–607. https://doi.org/10.1111/j.1365-313x.2008.03442.x.

    Article  CAS  PubMed  Google Scholar 

  15. Yang Z, Ohlrogge JB. Turnover of fatty acids during natural senescence of Arabidopsis, Brachypodium, and switchgrass and in Arabidopsis β-oxidation mutants. Plant Physiol. 2009;150:1981–9. https://doi.org/10.1104/pp.109.140491.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Murphy DJ. The biogenesis and functions of lipid bodies in animals, plants and microorganisms. Prog Lipid Res. 2001;40:325–438. https://doi.org/10.1016/s0163-7827(01)00013-3.

    Article  CAS  PubMed  Google Scholar 

  17. Luo G, Cao VD, Kannan B, Liu H, Shanklin J, Altpeter F. Metabolic engineering of energycane to hyperaccumulate lipids in vegetative biomass. BMC Biotechnol. 2022;22:24. https://doi.org/10.1186/s12896-022-00753-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Parajuli S, Kannan B, Karan R, Sanahuja G, Liu H, Garcia-Ruiz E, et al. Towards oilcane: engineering hyperaccumulation of triacylglycerol into sugarcane stems. GCB Bioenergy. 2020;12:476–90. https://doi.org/10.1111/gcbb.12684.

    Article  CAS  Google Scholar 

  19. Vanhercke T, Dyer JM, Mullen RT, Kilaru A, Rahman MM, Petrie JR, et al. Metabolic engineering for enhanced oil in biomass. Prog Lipid Res. 2019;74:103–29. https://doi.org/10.1016/j.plipres.2019.02.002.

    Article  CAS  PubMed  Google Scholar 

  20. Zale J, Jung JH, Kim JY, Pathak B, Karan R, Liu H, et al. Metabolic engineering of sugarcane to accumulate energy-dense triacylglycerols in vegetative biomass. Plant Biotechnol J. 2016;14:661–9. https://doi.org/10.1111/pbi.12411.

    Article  CAS  PubMed  Google Scholar 

  21. Vanhercke T, El Tahchy A, Liu Q, Zhou X, Shrestha P, Divi UK, et al. Metabolic engineering of biomass for high energy density: oilseed-like triacylglycerol yields from plant leaves. Plant Biotechnol J. 2014;12:231–9. https://doi.org/10.1111/pbi.12131.

    Article  CAS  PubMed  Google Scholar 

  22. Fan J, Yan C, Zhang X, Xu C. Dual role for phospholipid: diacylglycerol acyltransferase: enhancing fatty acid synthesis and diverting fatty acids from membrane lipids to triacylglycerol in Arabidopsis leaves. Plant Cell. 2013;25:3506–18. https://doi.org/10.1105/tpc.113.117358.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Vanhercke T, El Tahchy A, Shrestha P, Zhou X-R, Singh SP, Petrie JR. Synergistic effect of WRI1 and DGAT1 coexpression on triacylglycerol biosynthesis in plants. FEBS Lett. 2013;587:364–9. https://doi.org/10.1016/j.febslet.2012.12.018.

    Article  CAS  PubMed  Google Scholar 

  24. Vanhercke T, Petrie JR, Singh SP. Energy densification in vegetative biomass through metabolic engineering. Biocatal Agric Biotechnol. 2014;3:75–80. https://doi.org/10.1016/j.bcab.2013.11.010.

    Article  Google Scholar 

  25. Yang Y, Munz J, Cass C, Zienkiewicz A, Kong Q, Ma W, et al. Ectopic expression of WRINKLED1 affects fatty acid homeostasis in Brachypodium distachyon vegetative tissues. Plant Physiol. 2015;169:1836–47. https://doi.org/10.1104/pp.15.01236.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Beechey-Gradwell Z, Cooney L, Winichayakul S, Andrews M, Hea SY, Crowther T, et al. Storing carbon in leaf lipid sinks enhances perennial ryegrass carbon capture especially under high N and elevated CO2. J Exp Bot. 2020;71:2351–61. https://doi.org/10.1093/jxb/erz494.

    Article  CAS  PubMed  Google Scholar 

  27. Vanhercke T, Belide S, Taylor MC, El Tahchy A, Okada S, Rolland V, et al. Up-regulation of lipid biosynthesis increases the oil content in leaves of Sorghum bicolor. Plant Biotechnol J. 2019;17:220–32. https://doi.org/10.1111/pbi.12959.

    Article  CAS  PubMed  Google Scholar 

  28. Liu Q, Guo Q, Akbar S, Zhi Y, El Tahchy A, Mitchell M, et al. Genetic enhancement of oil content in potato tuber (Solanum tuberosum L.) through an integrated metabolic engineering strategy. Plant Biotechnol J. 2017;15:56–67. https://doi.org/10.1111/pbi.12590.

    Article  CAS  PubMed  Google Scholar 

  29. Focks N, Benning C. wrinkled1: a novel, low-seed-oil mutant of Arabidopsis with a deficiency in the seed-specific regulation of carbohydrate metabolism. Plant Physiol. 1998;118:91–101. https://doi.org/10.1104/pp.118.1.91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. 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:575–85. https://doi.org/10.1111/j.1365-313x.2004.02235.x.

    Article  CAS  PubMed  Google Scholar 

  31. Baud S, Wuilleme 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:933–47. https://doi.org/10.1111/j.1365-313x.2009.04011.x.

    Article  CAS  PubMed  Google Scholar 

  32. Baud S, Lepiniec L. Physiological and developmental regulation of seed oil production. Prog Lipid Res. 2010;49:235–49. https://doi.org/10.1016/j.plipres.2010.01.001.

    Article  CAS  PubMed  Google Scholar 

  33. 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:476–87. https://doi.org/10.1111/j.1365-313x.2009.03967.x.

    Article  CAS  Google Scholar 

  34. Lung S, Weselake RJ. Diacylglycerol acyltransferase: a key mediator of plant triacylglycerol synthesis. Lipids. 2006;41:1073–88. https://doi.org/10.1007/s11745-006-5057-y.

    Article  CAS  PubMed  Google Scholar 

  35. Saha S, Enugutti B, Rajakumari S, Rajasekharan R. Cytosolic triacylglycerol biosynthetic pathway in oilseeds. Molecular cloning and expression of peanut cytosolic diacylglycerol acyltransferase. Plant Physiol. 2006;141:1533–43. https://doi.org/10.1104/pp.106.082198.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Milcamps A, Tumaney AW, Paddock T, Pan DA, Ohlrogge J, Pollard M. Isolation of a gene encoding a 1, 2-diacylglycerol-sn-acetyl-CoA acetyltransferase from developing seeds of Euonymus alatus. J Bio Chem. 2005;280:5370–7. https://doi.org/10.1074/jbc.m410276200.

    Article  CAS  Google Scholar 

  37. Zou J, Wei Y, Jako C, Kumar A, Selvaraj G, Taylor DC. The Arabidopsis thaliana TAG1 mutant has a mutation in a diacylglycerol acyltransferase gene. Plant J. 1999;19:645–53. https://doi.org/10.1046/j.1365-313x.1999.00555.x.

    Article  CAS  PubMed  Google Scholar 

  38. Routaboul J-M, Benning C, Bechtold N, Caboche M, Lepiniec L. The TAG1 locus of Arabidopsis encodes for a diacylglycerol acyltransferase. Plant Physiol Biochem. 1999;37:831–40. https://doi.org/10.1016/s0981-9428(99)00115-1.

    Article  CAS  Google Scholar 

  39. Andrianov V, Borisjuk N, Pogrebnyak N, Brinker A, Dixon J, Spitsin S, et al. Tobacco as a production platform for biofuel: overexpression of Arabidopsis DGAT and LEC2 genes increases accumulation and shifts the composition of lipids in green biomass. Plant Biotechnol J. 2010;8:277–87. https://doi.org/10.1111/j.1467-7652.2009.00458.x.

    Article  CAS  PubMed  Google Scholar 

  40. Capuano F, Beaudoin F, Napier JA, Shewry PR. Properties and exploitation of oleosins. Biotechnol Adv. 2007;25:203–6. https://doi.org/10.1016/j.biotechadv.2006.11.006.

    Article  CAS  Google Scholar 

  41. Parthibane V, Rajakumari S, Venkateshwari V, Iyappan R, Rajasekharan R. Oleosin is bifunctional enzyme that has both monoacylglycerol acyltransferase and phospholipase activities. J Biol Chem. 2012;287:1946–54. https://doi.org/10.1074/jbc.m111.309955.

    Article  CAS  PubMed  Google Scholar 

  42. Winichayakul S, Scott RW, Roldan M, Hatier J-HB, Livingston S, Cookson R, et al. In vivo packaging of triacylglycerols enhances Arabidopsis leaf biomass and energy density. Plant Physiol. 2013;162:626–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Liang Y, Yu X, Anaokar S, Shi H, Dahl WB, Cai Y, et al. Engineering triacylglycerol accumulation in duckweed (Lemna japonica). Plant Biotechnol J. 2023;21:317–30. https://doi.org/10.1111/pbi.13943.

    Article  CAS  PubMed  Google Scholar 

  44. Rozov SM, Deineko EV. Strategies for optimizing recombinant protein synthesis in plant cells: classical approaches and new directions. Mol Biol. 2019;53:157–75. https://doi.org/10.1134/S0026893319020146.

    Article  CAS  Google Scholar 

  45. Callis J, Fromm M, Walbot V. Introns increase gene expression in cultured maize cells. Genes Dev. 1987;1:1183–200. https://doi.org/10.1101/gad.1.10.1183.

    Article  CAS  PubMed  Google Scholar 

  46. Gallegos JE, Rose AB. An intron-derived motif strongly increases gene expression from transcribed sequences through a splicing independent mechanism in Arabidopsis thaliana. Sci Rep. 2019;9:1–9. https://doi.org/10.1038/s41598-019-50389-5.

    Article  CAS  Google Scholar 

  47. Dwyer K, Agarwal N, Pile L, Ansari A. Gene architecture facilitates intron-mediated enhancement of transcription. Front Mol Biosci. 2021;8:669004. https://doi.org/10.3389/fmolb.2021.669004.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Richard TL. Challenges in scaling up biofuels infrastructure. Science. 1979;2010(329):793–6. https://doi.org/10.1126/science.1189139.

    Article  CAS  Google Scholar 

  49. Qin H, Dong Y, von Arnim AG. Epigenetic interactions between Arabidopsis transgenes: characterization in light of transgene integration sites. Plant Mol Biol. 2003;52:217–31. https://doi.org/10.1023/A:1023941123149.

    Article  CAS  PubMed  Google Scholar 

  50. Gibney ER, Nolan CM. Epigenetics and gene expression. Heredity. 2010;105(1):4–13. https://doi.org/10.1038/hdy.2010.54.

    Article  CAS  PubMed  Google Scholar 

  51. Hobbs SL, Kpodar P, DeLong CM. The effect of T-DNA copy number, position and methylation on reporter gene expression in tobacco transformants. Plant Mol Biol. 1990;15:851–64. https://doi.org/10.1007/BF00039425.

    Article  CAS  PubMed  Google Scholar 

  52. Matzke AJ, Matzke MA. Position effects and epigenetic silencing of plant transgenes. Curr Opin Plant Biol. 1998;1(2):142–8. https://doi.org/10.1016/S1369-5266(98)80016-2.

    Article  CAS  PubMed  Google Scholar 

  53. Assaad FF, Tucker KL, Signer ER. Epigenetic repeat-induced gene silencing (RIGS) in Arabidopsis. Plant Mol Biol. 1993;22:1067–85. https://doi.org/10.1007/BF00028978.

    Article  CAS  PubMed  Google Scholar 

  54. Hobbs SL, Warkentin TD, DeLong CM. Transgene copy number can be positively or negatively associated with transgene expression. Plant Mol Biol. 1993;21:17–26. https://doi.org/10.1007/BF00039614.

    Article  CAS  PubMed  Google Scholar 

  55. Xu-Gang LI, Song-Biao CH, Zi-Xian LU, Tuan-Jie CH, Qian-Chun ZE, Zhen ZH. Impact of copy number on transgene expression in tobacco. J Integr Plant Biol. 2002;44(1):120.

    Google Scholar 

  56. Dean C, Favreau M, Bond-Nutter D, Bedbrook J, Dunsmuir P. Sequences downstream of translation start regulate quantitative expression of two petunia rbcS genes. Plant Cell. 1989;1:201–8. https://doi.org/10.1105/tpc.1.2.201.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Mascarenhas D, Mettler IJ, Pierce DA, Lowe HW. Intron-mediated enhancement of heterologous gene expression in maize. Plant Mol Biol. 1990;15:913–20. https://doi.org/10.1007/bf00039430.

    Article  CAS  PubMed  Google Scholar 

  58. McElroy D, Zhang W, Cao J, Wu R. Isolation of an efficient actin promoter for use in rice transformation. Plant Cell. 1990;2:163–71. https://doi.org/10.1105/tpc.2.2.163.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Oard JH, Paige D, Dvorak J. Chimeric gene expression using maize intron in cultured cells of breadwheat. Plant Cell Rep. 1989;8:156–60. https://doi.org/10.1007/bf00716830.

    Article  CAS  PubMed  Google Scholar 

  60. Silva EM, Mettler IJ, Dietrich PS, Sinibaldi RM. Enhanced transient expression in maize protoplasts. Genome. 1988;30:72. https://doi.org/10.1007/978-1-4612-2694-9_109.

    Article  Google Scholar 

  61. Vasil V, Clancy M, Ferl RJ, Vasil IK, Hannah LC. Increased gene expression by the first intron of maize shrunken-1 locus in grass species. Plant Physiol. 1989;91:1575–9. https://doi.org/10.1104/pp.91.4.1575.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Rose AB. Introns as gene regulators: a brick on the accelerator. Front Genet. 2019;9:672. https://doi.org/10.3389/fgene.2018.00672.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Kim MJ, Kim H, Shin JS, Chung C-H, Ohlrogge JB, Suh MC. Seed-specific expression of sesame microsomal oleic acid desaturase is controlled by combinatorial properties between negative cis-regulatory elements in the SeFAD2 promoter and enhancers in the 5′-UTR intron. Mol Genet Genom. 2006;276:351–68. https://doi.org/10.1007/s00438-006-0148-2.

    Article  CAS  Google Scholar 

  64. Ott CJ, Suszko M, Blackledge NP, Wright JE, Crawford GE, Harris A. A complex intronic enhancer regulates expression of the CFTR gene by direct interaction with the promoter. J Cell Mol Med. 2009;13:680–92. https://doi.org/10.1111/j.1582-4934.2008.00621.x.

    Article  CAS  PubMed  Google Scholar 

  65. Morello L, Bardini M, Sala F, Breviario D. A long leader intron of the Ostub16 rice β-tubulin gene is required for high-level gene expression and can autonomously promote transcription both in vivo and in vitro. Plant J. 2002;29:33–44. https://doi.org/10.1046/j.0960-7412.2001.01192.x.

    Article  CAS  PubMed  Google Scholar 

  66. Koga M, Hayashi M, Kaida D. Splicing inhibition decreases phosphorylation level of Ser2 in Pol II CTD. Nucleic Acids Res. 2015;43:8258–67. https://doi.org/10.1093/nar/gkv740.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Dahan O, Gingold H, Pilpel Y. Regulatory mechanisms and networks couple the different phases of gene expression. Trends Genet. 2011;27:316–22. https://doi.org/10.1016/j.tig.2011.05.008.

    Article  CAS  PubMed  Google Scholar 

  68. Moabbi AM, Agarwal N, El Kaderi B, Ansari A. Role for gene looping in intron-mediated enhancement of transcription. PNAS. 2012;109:8505–10. https://doi.org/10.1073/pnas.1112400109.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Al-Husini N, Medler S, Ansari A. Crosstalk of promoter and terminator during RNA polymerase II transcription cycle. BBA-Gene Regul Mech. 2020;1863:194657. https://doi.org/10.1016/j.bbagrm.2020.194657.

    Article  CAS  Google Scholar 

  70. Le Hir H, Gatfield D, Izaurralde E, Moore MJ. The exon–exon junction complex provides a binding platform for factors involved in mRNA export and nonsense-mediated mRNA decay. EMBO J. 2001;20:4987–97. https://doi.org/10.1093/emboj/20.17.4987.

    Article  PubMed  PubMed Central  Google Scholar 

  71. Akua T, Shaul O. The Arabidopsis thaliana MHX gene includes an intronic element that boosts translation when localized in a 5′ UTR intron. J Exp Bot. 2013;64:4255–70. https://doi.org/10.1093/jxb/ert235.

    Article  CAS  PubMed Central  Google Scholar 

  72. Nott A, Le Hir H, Moore MJ. Splicing enhances translation in mammalian cells: an additional function of the exon junction complex. Genes Dev. 2004;18:210–22. https://doi.org/10.1101/gad.1163204.

    Article  CAS  PubMed Central  Google Scholar 

  73. Bourdon V, Harvey A, Lonsdale DM. Introns and their positions affect the translational activity of mRNA in plant cells. EMBO Rep. 2001;2:394–8. https://doi.org/10.1093/embo-reports/kve090.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Laxa M. Intron-mediated enhancement: a tool for heterologous gene expression in plants? Front Plant Sci. 2017;7:1977. https://doi.org/10.3389/fpls.2016.01977.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Snowden KC, Buchhholz WG, Hall TC. Intron position affects expression from the tpi promoter in rice. Plant Mol Biol. 1996;31:689–92. https://doi.org/10.1007/bf00042241.

    Article  CAS  PubMed  Google Scholar 

  76. Ding F, Elowitz MB. Constitutive splicing and economies of scale in gene expression. Nat Struct Mol Biol. 2019;26:424–32. https://doi.org/10.1038/s41594-019-0226-x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Xiao G, Zhang ZQ, Yin CF, Liu RY, Wu XM, Tan TL, et al. Characterization of the promoter and 5′-UTR intron of oleic acid desaturase (FAD2) gene in Brassica napus. Gene. 2014;545:45–55. https://doi.org/10.1016/j.gene.2014.05.008.

    Article  CAS  PubMed  Google Scholar 

  78. Ikeda C, Taku K, Miyazaki T, Shirai R, Nelson RS, Nyunoya H, et al. Cooperative roles of introns 1 and 2 of tobacco resistance gene N in enhanced N transcript expression and antiviral defense responses. Sci Rep. 2021;11:1–14. https://doi.org/10.1038/s41598-021-94713-4.

    Article  CAS  Google Scholar 

  79. Andersson R, Sandelin A. Determinants of enhancer and promoter activities of regulatory elements. Nat Rev Genet. 2020;21:71–87. https://doi.org/10.1038/s41576-019-0173-8.

    Article  CAS  PubMed  Google Scholar 

  80. Della Rosa M, Spivakov M. Silencers in the spotlight. Nat Genet. 2020;52:244–5. https://doi.org/10.1038/s41588-020-0583-8.

    Article  CAS  PubMed  Google Scholar 

  81. Kopp F, Mendell JT. Functional classification and experimental dissection of long noncoding RNAs. Cell. 2018;172:393–407. https://doi.org/10.1016/j.cell.2018.01.011.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Shlyueva D, Stampfel G, Stark A. Transcriptional enhancers: from properties to genome-wide predictions. Nat Rev Genet. 2014;15:272–86. https://doi.org/10.1038/nrg3682.

    Article  CAS  PubMed  Google Scholar 

  83. Rose AB, Carter A, Korf I, Kojima N. Intron sequences that stimulate gene expression in Arabidopsis. Plant Mol Biol. 2016;92:337–46. https://doi.org/10.1007/s11103-016-0516-1.

    Article  CAS  PubMed  Google Scholar 

  84. Jaeger D, Baier T, Lauersen KJ. Intronserter, an advanced online tool for design of intron containing transgenes. Algal Res. 2019;42:101588. https://doi.org/10.1016/j.algal.2019.101588.

    Article  Google Scholar 

  85. Jung JH, Altpeter F. TALEN mediated targeted mutagenesis of the caffeic acid O-methyltransferase in highly polyploid sugarcane improves cell wall composition for production of bioethanol. Plant Mol Biol. 2016;92:131–42. https://doi.org/10.1007/s11103-016-0499-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Kannan B, Jung JH, Moxley GW, Lee S, Altpeter F. TALEN-mediated targeted mutagenesis of more than 100 COMT copies/alleles in highly polyploid sugarcane improves saccharification efficiency without compromising biomass yield. Plant Biotechnol J. 2018;16:856–66. https://doi.org/10.1111/pbi.12833.

    Article  CAS  PubMed  Google Scholar 

  87. Eid A, Mohan C, Sanchez S, Wang D, Altpeter F. Multiallelic, targeted mutagenesis of magnesium chelatase with CRISPR/Cas9 provides a rapidly scorable phenotype in highly polyploid sugarcane. FGEED. 2021;3:654996. https://doi.org/10.3389/fgeed.2021.654996.

    Article  Google Scholar 

  88. Oz MT, Altpeter A, Karan R, Merotto A, Altpeter F. CRISPR/Cas9-mediated multi-allelic gene targeting in sugarcane confers herbicide tolerance. FGEED. 2021;3:673566. https://doi.org/10.3389/fgeed.2021.673566.

    Article  Google Scholar 

  89. Zhao Y, Karan R, Altpeter F. Error-free recombination in sugarcane mediated by only 30 nucleotides of homology and CRISPR/Cas9 induced DNA breaks or cre-recombinase. Biotechnol J. 2021;16:2000650. https://doi.org/10.1002/biot.202000650.

    Article  CAS  Google Scholar 

  90. Woodard KR, Prine GM. Dry matter accumulation of elephantgrass, energycane, and elephantmillet in a subtropical climate. Crop Sci. 1993;33(4):818–24. https://doi.org/10.2135/cropsci1993.0011183x003300040038x.

    Article  Google Scholar 

  91. Erickson JE, Soikaew A, Sollenberger LE, Bennett JM. Water use and water-use efficiency of three perennial bioenergy grass crops in Florida. Agriculture. 2012;2(4):325–38. https://doi.org/10.3390/agriculture2040325.

    Article  Google Scholar 

  92. Sigua GC, Stone KC, Bauer PJ, Szogi AA. Biomass and nitrogen use efficiency of grain sorghum with nitrogen and supplemental irrigation. Agron J. 2018;110(3):1119–27. https://doi.org/10.2134/agronj2017.09.0533

    Article  CAS  Google Scholar 

  93. Mahama GY, Prasad PV, Mengel DB, Tesso TT. Influence of nitrogen fertilizer on growth and yield of grain sorghum hybrids and inbred lines. Agron J. 2014;106(5):1623–30. https://doi.org/10.2134/agronj14.0092

    Article  Google Scholar 

  94. Perry J. Field evaluation of tobacco engineered for high leaf-oil accumulation. Theses and Dissertations—Plant and Soil Sciences. 2019. https://doi.org/10.13023/etd.2019.283.

  95. Berbeć AK, Matyka M. Biomass characteristics and energy yields of tobacco (Nicotiana tabacum L.) cultivated in eastern Poland. Agriculture. 2020;10(11):551.

    Article  Google Scholar 

  96. Van Gelder K, Oliveira-Filho ER, Messina CD, Venado RE, Wilker J, Rajasekar S, Ané JM, Amthor JS, Hanson AD. Running the numbers on plant synthetic biology solutions to global problems. Plant Sci. 2023 335:111815. https://doi.org/10.1016/j.plantsci.2023.111815

    Article  CAS  Google Scholar 

  97. Mudge SR, Basnayake SWV, Moyle RL, Osabe K, Graham MW, Morgan TE, et al. Mature-stem expression of a silencing-resistant sucrose isomerase gene drives isomaltulose accumulation to high levels in sugarcane. Plant Biotechnol J. 2013;11:502–9. https://doi.org/10.1111/pbi.12038.

    Article  CAS  PubMed  Google Scholar 

  98. Wang J, Li Y, Wai CM, Beuchat G, Chen L. Identification and analysis of stem-specific promoters from sugarcane and energy cane for oil accumulation in their stems. GCB Bioenergy. 2021;13:1515–27. https://doi.org/10.1111/gcbb.12872.

    Article  CAS  Google Scholar 

  99. Vanhercke T, Divi UK, El Tahchy A, Liu Q, Mitchell M, Taylor MC, et al. Step changes in leaf oil accumulation via iterative metabolic engineering. Metab Eng. 2017;39:237–46. https://doi.org/10.1016/j.ymben.2016.12.007.

    Article  CAS  PubMed  Google Scholar 

  100. Graef G, LaVallee BJ, Tenopir P, Tat M, Schweiger B, Kinney AJ, et al. A high-oleic-acid and low-palmitic-acid soybean: agronomic performance and evaluation as a feedstock for biodiesel. Plant Biotechnol J. 2009;7:411–21. https://doi.org/10.1111/j.1467-7652.2009.00408.x.

    Article  CAS  PubMed  Google Scholar 

  101. Sato N. Lipid biosynthesis in epidermal, guard and mesophyll cell protoplasts from leaves of Vicia faba L. Plant Cell Physiol. 1985;26(5):805–11. https://doi.org/10.1093/oxfordjournals.pcp.a076974.

    Article  CAS  Google Scholar 

  102. Daloso DM, Antunes WC, Pinheiro DP, Waquim JP, Araújo WL, Loureiro ME, Fernie AR, Williams TC. Tobacco guard cells fix CO 2 by both Rubisco and PEP case while sucrose acts as a substrate during light-induced stomatal opening. Plant, Cell Environ. 2015;38(11):2353–71. https://doi.org/10.1111/pce.12555.

    Article  CAS  PubMed  Google Scholar 

  103. McLachlan DH, Lan J, Geilfus CM, Dodd AN, Larson T, Baker A, Horak H, Kollist H, He Z, Graham I, Mickelbart MV. The breakdown of stored triacylglycerols is required during light-induced stomatal opening. Curr Biol. 2016;26(5):707–12. https://doi.org/10.1016/j.cub.2016.01.019.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  104. Lebedenko EN, Birikh KR, Plutalov OV, Berlin YA. Method of artificial DNA splicing by directed ligation (SDL). Nucleic Acids Res. 1991;19:6757–61. https://doi.org/10.1093/nar/19.24.6757.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  105. Szybalski W, Kim SC, Hasan N, Podhajska AJ. Class-IIS restriction enzymes—a review. Gene. 1991;100:13–26. https://doi.org/10.1016/0378-1119(91)90345-c.

    Article  CAS  PubMed  Google Scholar 

  106. Gordon VS, Comstock JC, Sandhu HS, Gilbert RA, Sood S, Korndorfer P, et al. Registration of ‘UFCP 84–1047’sugarcane for use as a biofuel feedstock. J Plant Regist. 2016;10:251–7. https://doi.org/10.3198/jpr2015.03.0021crc.

    Article  Google Scholar 

  107. Fouad MW, Hao W, Xiong Y, Steeves C, Sandhu KS, Altpeter F. Generation of transgenic energy cane plants with integration of minimal transgene expression cassette. Curr Pharm Biotechnol. 2015;16:407–13. https://doi.org/10.2174/1389201016666150303151559.

    Article  CAS  PubMed  Google Scholar 

  108. Sandhu S, Altpeter F. Co-integration, co-expression and inheritance of unlinked minimal transgene expression cassettes in an apomictic turf and forage grass (Paspalum notatum Flugge). Plant Cell Rep. 2008;27:1755–65. https://doi.org/10.1007/s00299-008-0599-5.

    Article  CAS  PubMed  Google Scholar 

  109. Murray MG, Thompson W. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 1980;8:4321–6. https://doi.org/10.1093/nar/8.19.4321.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  110. Iskandar HM, Simpson RS, Casu RE, Bonnett GD, Maclean DJ, Manners JM. Comparison of reference genes for quantitative real-time polymerase chain reaction analysis of gene expression in sugarcane. Plant Mol Biol Report. 2004;22:325–37. https://doi.org/10.1007/bf02772676.

    Article  CAS  Google Scholar 

  111. Higo K, Ugawa Y, Iwamoto M, Korenaga T. Plant cis-acting regulatory DNA elements (PLACE) database: 1999. Nucleic Acids Res. 1999;27:297–300. https://doi.org/10.1093/nar/27.1.297.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  112. Lescot M, Déhais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, et al. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002;30:325–7. https://doi.org/10.1093/nar/30.1.325.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Dr. Hardev Sandhu, Everglades Research and Education Center, UF-IFAS for providing energycane tops.

Funding

This research was co-funded by the Biological and Environmental Research (BER) program, U.S. Department of Energy, under Award Number DE-SC0018254 and by the DOE Center for Advanced Bioenergy and Bioproducts Innovation (U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research) under award number DE-SC0018420. Award number DE-SC0018254 supported vector construction, generation of transgenic plants and their molecular analysis. Award Number DE-SC0018420 supported the lipid analysis of transgenic plants. The confocal microscopic study at the Center for Functional Nanomaterials at Brookhaven National Laboratory was supported under Office of Basic Energy Sciences of DOE with contract no. DE-SC0012704 (to CFN). The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. This work was also supported by the USDA National Institute of Food and Agriculture, Hatch project 1020425.

Author information

Authors and Affiliations

Authors

Contributions

FA and JS conceived and designed the experiments. FA designed the recombinant DNA constructs and managed the research project. GL and VDC carried out the tissue culture and transformation of energycane, grew the transgenic plants in the greenhouse and sampled for TAG, FA and qRT-PCR analysis. VDC and SK carried out qRT-PCR analysis. HL conducted analyses of TAG and total fatty acid content and composition. YL carried out the experiments leading to visualization of lipid droplets in leaf tissues. VDC analyzed data. VDC and FA wrote the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to John Shanklin or Fredy Altpeter.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1. Comparison of a native

DGAT1 gene (I) from Tropaeolum majus and an optimized DGAT1 gene (II). The underlined letters indicate nucleotides that were modified for codon optimization. The lowercase letters indicate the inserted intron. The TTNGATYTG-like motif is highlighted in purple with nucleotides that deviated from the motif are marked in green. Extra nucleotide compared to the TTNGATYTG is highlighted in gray. The GT1-consensus, the CAAT box and the NGATY core of the longer TTNGATYTG motif are highlighted in blue, yellow, and dark yellow, respectively. Figure S2. PCR analysis of transgenic plants. A to F. PCR amplification of DGAT1(W), DGAT1(In), OLE1, WRI1, or nptII, from genomic DNA of transgenic plants, respectively. PC. Positive control (plasmid used for transformation), Genomic DNA extracts of non-transgenic energycane plant (WT) were used as the negative control. W: no intron. In: intron. Arrows indicate target amplicon.

Additional file 2: Table S1.

List of primers used for gene expression analysis. Table S2. Correlation of total FA content with TAG content in transgenic energycane during plant development. Table S3. Correlation of TAG content with total FA, biomass yield DW, height, circumference, stem diameter and tiller number in transgenic energycane.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cao, V.D., Luo, G., Korynta, S. et al. Intron-mediated enhancement of DIACYLGLYCEROL ACYLTRANSFERASE1 expression in energycane promotes a step change for lipid accumulation in vegetative tissues. Biotechnol Biofuels 16, 153 (2023). https://doi.org/10.1186/s13068-023-02393-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-023-02393-1

Keywords