- Open Access
PLDα1-knockdown soybean seeds display higher unsaturated glycerolipid contents and seed vigor in high temperature and humidity environments
© The Author(s) 2019
- Received: 23 August 2018
- Accepted: 13 December 2018
- Published: 4 January 2019
Soybean oil constitutes an important source of vegetable oil and biofuel. However, high temperature and humidity adversely impacts soybean seed development, yield, and quality during plant development and after harvest. Genetic improvement of soybean tolerance to stress environments is highly desirable.
Transgenic soybean lines with knockdown of phospholipase Dα1 (PLDα1KD) were generated to study PLDα1′s effects on lipid metabolism and seed vigor under high temperature and humidity conditions. Under such stress, as compared with normal growth conditions, PLDα1KD lines showed an attenuated stress-induced deterioration during soybean seed development, which was associated with elevated expression of reactive oxygen species-scavenging genes when compared with wild-type control. The developing seeds of PLDα1KD had higher levels of unsaturation in triacylglycerol (TAG) and major membrane phospholipids, but lower levels of phosphatidic acid and lysophospholipids compared with control cultivar. Lipid metabolite and gene expression profiling indicates that the increased unsaturation on phosphatidylcholine (PC) and enhanced conversion between PC and diacylglycerol (DAG) by PC:DAG acyltransferase underlie a basis for increased TAG unsaturation in PLDα1KD seeds. Meanwhile, the turnover of PC and phosphatidylethanolamine (PE) into lysoPC and lysoPE was suppressed in PLDα1KD seeds under high temperature and humidity conditions. PLDα1KD developing seeds suffered lighter oxidative stresses than did wild-type developing seeds in the stressful environments. PLDα1KD seeds contain higher oil contents and maintained higher germination rates than the wild-type seeds.
The study provides insights into the roles of PLDα1 in developing soybean seeds under high temperature and humidity stress. PLDα1KD decreases pre-harvest deterioration and enhances acyl editing in phospholipids and TAGs. The results indicate a way towards improving production of quality soybean seeds as foods and biofuels under increasing environmental stress.
- Phospholipase D
- Oxidative stress
- Acyl editing
- High temperature and humidity
- Oil content
- Seed vigor
Developing soybean seeds suffer from high temperature and humidity stress conditions that often occur in Southern China and cause oxidative stress and seed pre-harvest deterioration. Such soybean seed deterioration can be alleviated by seed-specific knockdown of phospholipase Dα1 (PLDα1KD). Thus, PLDα1KD soybeans have improved nutrition quality and seed vigor.
PLDα1KD developing seeds have higher levels of unsaturation in triacylglycerol and phospholipids due to enhanced expression of desaturase genes and enhanced acyl editing and PC–DAG conversion activities in PLDα1KD developing seeds compared with the wild-type seeds.
The finding that PLDα1KD seeds had improved seed vigor, nutrition quality, and tolerance to high temperature and humidity provides a molecular tool for genetic improvement of soybean for adaptation to wider growth conditions.
With global temperature increasing steadily in recent decades, high temperature conditions, accompanied either by drought or by humidity in different areas, caused damages and losses on crops and yield [1–4]. However, how high temperature and humidity affect crop growth, seed development, as well as yield is not well understood [2, 4]. The high temperature and humidity conditions adversely affect membrane lipids and storage triacylglycerol (TAG) in wheat [3, 5] and oilseed crops, such as soybean (Glycine max), particularly during seed development . High temperatures impact developing soybean seeds’ sensitivity and vulnerability to stresses by causing seeds with poor germination, increased incidence of pathogen infection, and decreased economic value . Previous studies indicated that phospholipase Dα (PLDα1) is involved in seed naturally and artificially aging, by affecting both phospholipids and TAG of mature Arabidopsis and soybean seeds [7–9]. However, how PLDα1 affects lipid metabolism and storage nutrition of developing soybean seeds grown under high temperature and humidity conditions was unknown [8, 10].
Assembly of phospholipids and TAG occurs primarily in the endoplasmic reticulum (ER) and shares a common biosynthetic precursor phosphatidic acid (PA) . Membrane phospholipid synthesis is active in young and green tissues . However, in developing oilseeds, phospholipid metabolism is overwhelmingly directed to TAG accumulation [13, 14], through two major pathways: the diacylglycerol acyltransferase (DGAT)-mediated Kennedy pathway that uses acyl-CoA and the diacylglycerol (DAG) to generate TAG, and phospholipid:diacylglycerol acyltransferase (PDAT)-mediated pathway that uses phosphatidylcholine (PC) and DAG to produce TAG . Overexpression of an A-type phospholipase, pPLAIIIδ, enhanced TAG production in Arabidopsis and Camelina seeds [12, 15]. A-type phospholipase, PLA catalyzing the hydrolysis of PC to generate a free fatty acid and lysoPC (LPC), can also affect TAG production [12, 15]. Acyl-CoA:LPC acyltransferase (LPCAT) can modify PC saturation by introducing an acyl-CoA into a new PC [16, 17]. De novo PC biosynthesis from choline by the actions of choline/ethanolamine kinase (CEK), choline-phosphate cytidylyltransferase (CCT), and DAG cholinephosphotransferase (DAG:CPT) is also important for phospholipid and TAG metabolism . The interconversion between PC and DAG by PC:DAG cholinephosphotransferase (PDCT, also ROD1) can significantly affect TAG biosynthesis and unsaturation through above mentioned pathways [16, 19, 20]. PDCT transfers phosphocholine from PC to DAG actively during oil seed development, and edits TAG composition using PC that is also extensively modified by fatty acid desaturases (FADs) on their acyl chains [19, 21, 22]. In the ER, the partitioning of PA, PC, and DAG precursors for TAG or phospholipid biosynthesis may be controlled by a unexplored complex network. For instance, it has been estimated that more than 70% PC-derived DAG is used to synthesize TAG in flax seeds . Thus, PC plays multiple roles in TAG and phospholipid biosynthesis by recycling or incorporation of the newly synthesized fatty acids in TAG acyl editing .
PLD may contribute to diurnal cycling of PA, and PC acyl editing and significantly affect PC pools and TAG production [24–28]. A study demonstrated that PC-derived DAG is the major source for TAG synthesis in Camelina seeds overexpressing PLDζ1 and PLDζ2 . Recently, a genome-wide association study suggested a PLDα gene as a key locus affecting oil biosynthesis in soybean [29, 30]. Studies indicated that PLDα1 is involved in the mature seed aging and deterioration stored under high temperature and humidity [7, 9]. Two PLDα-RNAi knockdown (PLDα1KD) soybean lines had had altered unsaturation fatty acids in both phospholipids and TAG when grown in Kansas in the United States . However, it is still not understood why PLDα1KD seeds have such changes in fatty acyl chains saturation.
Although the anti-deterioration effects of PLDα1 mutation on naturally or artificially aging Arabidopsis and soybean seeds have been reported [7–9], the underlying molecular mechanism by which PLDα1 mutation affects lipids metabolism and seed quality has not been explored [8, 9]. Especially, how PLDα1 mutation affects soybean developing seeds under heat and humidity stresses is unknown. High temperature and humidity stresses occur frequently in southern China amid increasing global temperature in recent years and cause soybean pre-harvest deterioration and significant losses . High temperature and humidity are major problems and limiting factors in soybean production area, and they not only caused soybean yield loss, but also affected adversely on seed storage and reduced nutrition . In the Mid-Yangzi River region of China, soybean seed development and ripening occur during the season of high precipitations (~ more than 150 mm), high temperature (~ 34–38 °C), and high humidity (~ 75–80%) usually from the late June to the early September (Fig. 3a). The soybean seeds grown in these regions usually have pre-harvest deterioration, rapidly losing seed vigor, reducing nutrition quality severely, and are more vulnerable to pathogens during storage [6, 31]. A better understanding of the mechanism by which high temperature and humidity impact developing soybean seeds will help design effective genetic strategies to improve soybean tolerance to these environmental stress conditions. Here, we generated PLDα1KD soybean lines and investigated the PLDα1KD soybean performance in such stressful environments. The data indicate a critical role of PLDα1 in lipid metabolism and stress response in developing soybean seeds under high temperature and humidity stress conditions, as compared in normal environments.
Generation of transgenic soybean plants with seed-specific knockdown of PLDα1
Effects of PLDα1KD on the expression of other PLDs
PLDα1KD developing seeds had elevated levels of ROS-scavenging genes
Increased unsaturated fatty acids of TAG and phospholipids in PLDα1KD seeds
Up-regulation of FADs in PLDα1KD developing seeds as compared with wild-type
Quantitative RT-PCR results showed that the expression of FAD2-1B was much higher in developing seeds of PLDα1KD than in wild-type under both normal and stress conditions (Fig. 5b). Microarray data showed that FAD2s were highly expressed in seeds during development . FAD2-1B and FAD2-1A were mainly expressed at the late stages of developing seeds, and their transcripts were higher in PLDα1KD than wild-type under both conditions. The transcript levels of GmFAD3A, GmFAD3B and GmFAD3C were generally high in developing seeds (Fig. 5b, Additional file 1: Figure S5). Under high temperature, FAD3B expression levels increased to the highest level at the stage 4, and then decreased. However, under normal conditions, the FAD3B transcript level was higher at all developmental stages compared with that under stressed conditions. Furthermore, the FAD3B transcript level was much higher in PLDα1KD1 than wild-type seeds under normal conditions (Fig. 5b). The expression levels of the major FADs in seeds were higher under normal conditions than high temperature conditions, except for GmFAD2-2A and GmFAD3C, whose expression levels were generally low under normal conditions, and increased in response to elevated temperatures. The expression of chloroplast localized GmFAD6A was also higher in PLDα1KDs than in wild-type and up-regulated under stress conditions (Fig. 5b).
Increased contents of PC and PE in PLDα1KD seeds
Decreased PA pools in PLDα1KD seeds
Three PAH genes, homologous to AtPAH1 and AtPAH2, are present in the soybean genome. Two of them, Glyma.13G134500 (GmPAHβ1) and Glyma.10G046400 (GmPAHβ2), were highly expressed in developing seeds, in a trend coincident with seed oil accumulation (Additional file 1: Figure S7) . The transcript levels for both GmPAHβ1 and GmPAHβ2 in PLDα1KD were one and twofold higher, respectively, than these in wild-type seeds at stages 4 under both conditions and the expression of GmPAHs was not affected at all by high temperature except GmPAHβ1 at stage 2 and GmPAHβ2 at stage 3 (Fig. 8b). The combined effects of suppressed PLDα1KD, and a markedly higher PAH expression level contributed to the decreased PA levels, which was confirmed by mature and developing seeds (Figs. 6, 7).
Altered transcript levels of DGATs and PDATs for TAG biosynthesis in PLDα1KD seeds
To explore how PLDα1KD affected TAG biosynthesis and phospholipid metabolism in soybean seeds, we examined several major genes involved in the Kennedy pathway (Fig. 8b). DGAT synthesizes TAG by transferring an acyl group to DAG from newly synthesized or recycled acyl-CoA (Fig. 8a). The DGAT family in the soybean genome has 10 members. Type 1 DGATs, Glyma.13G106100, Glyma.09G065300, and Glyma.17G053300, were highly expressed in seeds. Type 3 DGAT Glyma.17G041600 was also highly expressed in seeds. Compared with type 1 and type 3 DGATs, type 2 DGAT, such as Glyma.16G115700 and Glyma.09G195400, were expressed at a lower level in seeds . The transcript level of these genes increased steadily over seed development (Additional file 1: Figure S8) . PLDα1KD lines have lower transcript levels for several seed-specific DGATs, such as GmDGAT1A (Glyma.13G106100), GmDGAT1C (Glyma.09G065300) and GmDGAT3B (Glyma.17G041600), over all developmental stages under both conditions, indicating that PLDα1KD lines have reduced contributions through DGAT pathway towards TAG synthesis (Fig. 8b). Meanwhile, the expression of GmDGAT1A and GmDGAT3B was increased whereas the expression of GmDGAT1C was suppressed in PLDα1KD lines and wild-type under high temperature conditions compared with normal conditions.
PDAT transfers the sn-2 acyl group from phosphatidylcholine or phosphatidylethanolamine to DAG for TAG production in plants and yeast (Fig. 8a) [45, 46], and in soybean seeds, DAG from PC is primarily used for TAG biosynthesis. PDAT and DGAT were shown to have overlapping functions in TAG biosynthesis . The soybean genome contains 6 putative PDAT genes, and among them, Glyma.12G084000, Glyma.11G190400, and Glyma.13G108100, as well as Glyma.07G036400, were highly expressed in seeds. Transcripts of these PDATs increased steadily during seed development except Glyma.07G036400 (Additional file 1: Figure S9) . The transcript level of PDAT genes in PLDα1KD seeds was more than higher than those in wild-type seeds from stages 3 to 5 under both conditions, suggesting that PLDα1KD seeds have increased PDAT-mediated, DAG-PC dependent TAG biosynthesis (Fig. 8b). Meanwhile, high temperature and humidity suppressed the expression of GmPDAT1B in both PLDα1KD and wild-type seeds.
Enhanced PC–DAG conversion and acyl editing in PLDα1KD soybean seeds
To test whether the active PC–DAG–PDAT pathway contributed to more TAG biosynthesis in PLDα1KD than in wild-type, we compared the transcript levels of relevant genes. For PC synthesis, choline/ethanolamine kinase (CEK) produces phosphocholine that is used by CTP: phosphocholine cytidylyltransferase (CCT) to synthesize CDP- choline (Additional file 1: Figures S10, S11). DAG: cholinephosphotransferase (CPT) then transfers choline form CDP-choline to DAG to generate PC (Fig. 9a). The soybean genome has two CCT genes, CCT1 (Glyma.09G051200) and CCT2 (Glyma.15G157500), and their transcript levels fluctuated during seed development (Fig. 9b, Additional file 1: Figure S11) . However, transcripts of CCTs in PLDα1KD lines were down-regulated by approximately 18% at stage 3 under both conditions. The expression of GmCCTs was suppressed in both PLDα1KD lines and wild-type seeds under high temperature and humidity conditions.
DAG-CPT and PDCT form an important PC–DAG exchange/conversion cycle to enforce the acyl editing of TAGs (Fig. 9a). In PLDα1KD developing seeds, DAG:CPTs (also called AAPTs), DAG:CPT1 and 2 (Glyma.12G081900 and Glyma.02G128300, respectively), were up-regulated as compared with those in developing seeds of wild-type under both conditions (Fig. 9b, Additional file 1: Figure S12). The soybean genome contains two PDCT genes, GmPDCT1 (Glyma.07G029800) and 2 (Glyma.08G213100). The two genes were highly expressed in developing soybean seeds (Additional file 1: Figure S13) . The transcript of GmPDCT2 was up-regulated at early developmental stages and then decreased during late seed stages under both conditions. GmPDCT1 and 2 were significantly up-regulated at stages 2-3 in PLDα1KD developing seeds under both conditions (Fig. 9b). The expression of GmPDCT1 and GmDAG-CPT2 was suppressed at stages 3-5 in PLDα1KD developing seeds under high temperature and humidity conditions. Other two genes GmPDCT2 and GmDAG-CPT1 displayed complicated expression patterns in wild-type and PLDα1KD developing seeds in both environments. These data suggest that the activity of PC and DAG interconversion is increased when PLDα1 was suppressed in developing soybean seeds.
Reduced transcript levels of PLAs but increased levels of LPCAT in PLDα1KD seeds
Since pPLA affects TAG biosynthesis , we examined the expression of pPLAs that were either specifically or highly expressed in the developing soybean seeds. The soybean genome contains a large pPLA gene family, and several pPLAs were highly expressed in developing seeds, such as Glyma.08g028800, Glyma.11g036900, and Glyma.17G145900. The pPLAs, Glyma.18g251500 and Glyma.09g243100, were expressed only in seeds (Additional file 1: Figure S14) [32, 33]. The transcripts of pPLAs (Glyma01G.002400 and Glyma08G.028800) in PLDα1KD soybean seeds were, on average, 44% lower than those in wild-type seeds at stage-4 under both conditions, which was consistent with higher levels of PCs in PLDα1KD soybean seeds and lower levels of lysophospholipids (Figs. 6, 7, 10b). Meanwhile, high temperature and humidity suppressed the expression of GmPLA2s in both PLDα1KD lines and WT seeds. The down-regulation of both PLD and PLA expression may explain the higher level of PCs in PLDα1KD seeds.
As LPC and LPE content decreased significantly in PLDα1KD than wild-type seeds by more than tenfold at stage 3, we examined the expression of genes involved in the LPC and PC cycle. LPC acyltransferases (LPCAT) catalyzes the synthesis of PC from LPC using a new fatty acyl-CoA (Fig. 10a). AtLPCAT1 and 2 in Arabidopsis control the acyl editing process by acting as the main entry of unsaturated FAs into PC . The lpcat1/lpcat2 mutant showed decreased PUFA in seed TAG [16, 17]. Soybean genome contains two LPCATs, Glyma.17G131500 (GmLPCAT1) and Glyma.05G049500 (GmLPCAT2). Transcript levels of two genes increased to the highest levels at middle stages and then decreased at the later stage of seed development in both PLDα1KD lines under stress conditions (Fig. 10b, Additional file 1: Figure S15) . Both GmLPCAT transcripts were 17% higher in PLDα1KD seeds than in wild-type seeds at later stages under stress conditions (Fig. 10b). Meanwhile, Both genes were up-regulated in both PLDα1KD lines and wild-type under stress conditions. The down-regulation of GmPLAs and up-regulation of GmLPCATs in PLDα1KD soybean seeds could lead to reduced contents of LPC and LPE, as compared with wild-type (Figs. 6, 7, 10).
Higher germination rate of PLDα1KD seeds
PLDα1 knockdown promoted fatty acid unsaturation in both TAGs and phospholipids
There are increases in the total contents of both TAGs and unsaturation fatty acids in the developing seeds of PLDα1KD than wild-type cultivar under high temperature and humidity conditions. The desaturation of PCs is primarily attributable to FAD2 and FAD3 in the ER [48–50]. The increased transcripts of FAD2s and FAD3s in PLDα1KD seeds explain the increased unsaturated fatty acids in PCs and PEs, and TAG species, together with the decreased levels of PA, LPC, and LPE species in the PLDα1KD developing seeds. The difference in lipid contents between PLDα1KD and wild-type seeds under normal conditions was similar to that reported previously . The high temperature and humidity decreased the contents of total lipids, especially fatty acid contents in seeds compared with these under normal conditions. However, our data showed that the differences in lipid contents between GmPLDα1KD and wild-type seeds under the stress conditions became much bigger than those under normal conditions. Meanwhile, the high temperature and humidity conditions affected the expression of genes involved in PLs and TAG synthesis pathways, which may eventually result in decreased lipid contents. Consistently, under these stress conditions, developing seeds of GmPLDα1KD displayed higher gene expression levels and total lipid contents than those of wild-type seeds, e.g., GmFAD2-1A and GmFAD2-2A. Mature GmPLDα1KD seeds also had better germination rates than did wild-type seeds. Therefore, suppression of GmPLDα1 improved the expression of genes involved in PLs and TAG synthesis pathways under both high temperature and humidity stress and normal conditions.
PLDα1KD enhanced the PDAT pathway and DAG-PC conversion toward TAG biosynthesis
It has been proposed that the DAG:CPT-catalyzed reaction provides an acyl editing mechanism for the production of polyunsaturated TAGs containing 18:2 and 18:3 through PCs . The significant increases in unsaturation acyl chains of PCs in developing soybean seeds of PLDα1KD plants were attributable to the higher expression levels of FADs. Since both PC and DAG play essential roles in acyl editing on glycerolipids, their interconversion is important for TAG synthesis [14, 28, 52]. The knockdown of PLDα1 in the developing soybean seeds affected expression of all genes involved in the acyl editing in soybean developing seeds (Fig. 12). The up-regulation of both PDCTs and DAG:CPTs in consistent with an increased TAG contents in soybean developing seeds, since PCs are the major source of DAGs that further flux into TAGs by the action of PDATs . Since PA–DAG–PC–DAG metabolic pathway primarily takes place in soybean developing seeds, PDCT thus becomes essential in determining the unsaturation of TAG. PDCT mutation results in a 40% decrease in polyunsaturated FAs in seed TAG without disrupting overall TAG levels . PDAT’s substrate specificity on PC and DAG further enhanced the PC–DAG conversion and resulted in higher levels of unsaturated fatty acids in TAGs of developing seeds. The markedly increased PC/PE accumulation coincides with up-regulated FAD2 and 3, PDCT and PDAT genes suggest PLDα1 negatively affects acyl editing and bridges the phospholipid turnover and TAG biosynthesis.
PLA catalyzes PC or PE hydrolysis to generate free fatty acid and lysoPLs, such as LPC or LPE. PLAs and LPCATs can form a PC turnover cycle (Fig. 12). The lower transcripts of PLAs in PLDα1KD are consistent with higher levels of PCs and lower levels of LPC and LPE in PLDα1KD developing seeds, which is also consistent with previous report in Arabidopsis PLDα1KD plants . The down-regulation of PLDα1 and PLA genes may together lower PA, LPC and LPE levels . In pah1pah2 mutant, PC contents dramatically increased and TAG content decreased due to the increased PA contents through up-regulation of CCT1 . Our results suggest that PLA–LPCAT, PLD, and PDCT–DAG:CTP play important roles in acyl editing. PLDα1KD promotes the PDCT and DAG:CPT-catalyzed conversion between PCs and DAGs, through which extensive acyl desaturation on PCs by FADs is eventually reflected by the increased TAG unsaturation. Thus, the developing PLDα1KD soybean seeds have higher levels of di18:2 than these non-transgenic lines.
Effects of high temperature and humidity on soybean seed development
While high temperature and humidity can significantly affect membrane compositions, altering membrane physiological properties and functions, how the stress affects phospholipid and TAG metabolism in developing seeds is largely unknown [5, 9]. High temperature stress causes significant increases in ROS-scavenging, lipid desaturating, oxidizing, and acylating genes, and in 18:3-TAG contents in wheat, while fatty acyl chain unsaturation in polar phospholipids decrease . ROS plays various roles in regulating cell growth, development, and cell survival under stress. Generally, moderate levels of ROS may function as signals to promote plant growth and survival, whereas a large increase of ROS can induce plant cell death. Under physiologic conditions, the balance between generation and elimination of ROS maintains the proper function of redox-sensitive signaling proteins. However, under adversary environments, such as extreme heat, light, and cold conditions, plants produce more ROS or xenobiotics than normal levels, and activate ROS-scavenging enzymes and xenobiotics-detoxification enzyme GST to reduce ROS damage to plants.
Possible mechanisms for the improvement of soybean seeds by knockdown of GmPLDα1
PLDα1 and its product PA and PLA and its product LysPLs, have been reported in other plants to function in plant stress response to salt and drought stresses, oxidative stresses, and hormones [10, 27, 55]. Our study here indicates that high temperature and humidity environments significantly affect phospholipid and TAG metabolism in developing soybean seeds, and increase the contents of PAs and lysoPLs in wild-type seed, but the contents of PAs and lysoPLs were significantly reduced in PLDα1KD mutant seeds (Figs. 6, 7). The anti-deterioration effects of PLDα1KD mutation on developing soybean seeds under high temperature and humidity conditions in our study were similar to the better oil storability and seed viability of these naturally or artificially aged Arabidopsis and soybean PLDα1-antisense mutant seeds [7, 8]. Those results suggest that the lowered PAs and lysoPLs in PLDα1KD or knockout mutant seeds might be the major causes. The higher PA and lysoPL contents may increase membrane permeability during seed development and storage [7, 8, 10]. The variations in lipid contents and gene expression in two independent PLDα1KD lines may result from the variations of levels of PAs and lysoPLs because of different degrees of PLD suppression. It has been observed that different knockdown lines of AtPLDα1 often showed variations in phenotypes [7, 10, 27]. Besides its metabolic function, PLDα1 is also involved in plant responses to salinity, drought, cold, and heat stresses  (Fig. 12). The distinct responses of ROS-scavenging genes in PLDα1KD and wild-type cultivars may also be associated to GmPLDα1/PLA and PA/lysPLs that could be critically involved in soybean adaptive response to high temperature and humidity [27, 34, 35, 55]. We posit that PA and lysoPLs generated by PLDs and PLAs under high temperature and humidity conditions could be the keys for the detrimental effects of high temperature and humidity [7–9]. PLDα1KD substantially decreases stress-induced production of PA and lysoPLs and thereby alleviates the negative impacts by these adverse conditions on developing soybean seeds.
Genetically modified PLDα1KD soybean seeds have improved oil content, seed vigor, and resistance against aging and deterioration under high temperature and humility environments. The mechanism behind the phenomena was explored by examination of relevant metabolites and transcripts of lipid metabolic genes in WT and PLDα1KD soybean developing seeds grown under both stress and normal conditions. The higher TAG content with higher proportion of unsaturation degree in PLDα1KD than wild-type at most developmental stages are attributable to the lower expression levels of PLDs, PLA2s, CCTs, and DGATs and higher expression levels of FADs, LPCAT, PDAT, PAH, PDCT, and DAG:CPT in PLDα1KD soybean. Higher levels of ROS-scavenging and xenobiotics-detoxifying genes (GST, SOD, CAT, POD, STI, and APX) in PLDα1KD soybean developing seeds and lower ROS level in germinated PLDα1KD seeds than those in wild-type seeds may explain their tolerance against high temperature and humility stress in preventing pre-harvest seed deterioration. This study presents a valuable model illustrating the role of PLD in TAG synthesis and provides novel insights into the mechanistic details of lipid metabolic pathway changes upon the knockdown of PLDs during seed development under stress. Despite of the biotechnological application potential in genetic improvement of soybean production under stress, the field trials may be required to validate the advantages for quality soybean production, especially under frequently occurred heatwaves.
Creation of PLDα1KD transgenic soybean
The forward and reverse sequences of 1151 bp from the conserved GmPLDα1 cDNA in Williams 82 were used for making RNAi constructs, which were different from the forward 820 bp and reverse 1300 bp from conserved GmPLDα1 sequence derived from cloning partial sequence in Fayette used previously . The RNAi cassette was assembled into the NotI site of the soybean expression vector pBetaConSoyHyg, flanked on its 5′ end by the seed-specific promoter for the soybean α’ subunit of β-conglycinin gene  and on its 3′ end by the 3′UTR of the phaseolin gene . The vector backbone, which was derived from pBluescript SK-(Stratagene), was engineered with a hygromycin B phosphotransferase gene  under control of the potato ubiquitin-3 promoter  for selection of transformed soybean embryos. Additional file 1: Figure S1 illustrates the cloning strategy and details of the organization of gene cassettes generated in the vector pBetaConSoyHyg for soybean transformation. The construct was introduced into soybean (Glycine max cv. Jack) somatic embryos using biolistic transformation as described . Embryos were maintained after bombardment in SHaM media  with hygromycin selection . Mature somatic embryos obtained following hygromycin selection were screened by immunoblotting analysis for PLDα1 levels. Embryos, representing independent transgenic events, that showed reductions in PLDα1 levels, were desiccated and used for regeneration of transgenic plants as described . Regenerated plants were advanced to homozygosity (> T3 generation) for use in the described studies. One hundred forty-three putative tissues were recovered after selection. Ten transgenic soybean events were recovered and fully grown in the greenhouse. Using enzyme activity, PCR and western blotting, two successfully repressed transgenic soybean events (#1020 and #1048) were selected from ten transgenic events for analysis.
Plant growth conditions and treatments
Wild-type Jack and transgenic soybean PLDα1KD lines were grown in a greenhouse with a normal (26 ± 3 °C/18 h day and 23 ± 3 °C/6 h night photoperiod, 45–65% humidity) or high temperature and high humidity (36 ± 3 °C/18 h day and 30 ± 3 °C/6 h night photoperiod, 70–85% humidity) conditions which were controlled by air condition and humidifier machine. The developing seeds were harvested at approximately 50 days after fertilization (DAF). The varieties for transcriptomic analyses at different developmental stages of soybean seeds were Williams 82 and Hokkaido Black 25- and 50-day seeds were collected and put into liquid nitrogen immediately after separating from plants. To test the suppression of PLDα1 in transgenic soybean seeds, qRT-PCR, western blotting with PLDα1 antibody, and PLDα1 activity assays were performed using developing staged 3 seeds of T3 and T5 lines. The presence of the transgene in T5 line seeds was confirmed by PCR. Monitoring soybean PLD proteins in the soybean extract used antibodies against Arabidopsis thaliana PLDα1.
Analysis of seed total fatty acids
Total fatty acids were extracted three times and analyzed with GC using triheptadecanoylglycerol as an internal standard using a method described previously [44, 62]. Analyses were carried out on an Agilent Technologies (USA) 7890A Network gas chromatograph, equipped with a flame-ionization detector (FID) and a split/splitless injector. Polyethylene glycerol:Agilent 19091 N-133 column (30 m × 0.25 mm i.d. 0.25 m film thickness) was used . An HP ChemStation (Hewlett-Packard, Palo Alto, CA, USA) was used for instrument control and data analysis.
Analysis of seed TAGs
Seed TAG and its fatty acid compositions were analyzed as described previously . Briefly, total seed lipids from developing seeds were separated into TAG, DAG, and free fatty acids on a TLC plate, and TAG spots were visualized with iodine vapor and identified by their migration. The spots of TAG were scraped off TLC plates for determination of fatty acid contents and compositions by GC as described above.
Immunoblotting and PLD activity assays
Proteins were extracted from immature or mature soybean seeds with extraction buffer containing 50 mm Tris–HCl (pH 7.5), 10 mm KCl, 1 mm EDTA, 0.5 mm phenyl methylsulfonyl fluoride and 2 mm DTT at 4 °C . Proteins were incubated at boiling water bath for 10 min with a SDS-PAGE loading buffer containing 100 μL of 50 mm Tris–HCl (pH 6.8), 10 mm DTT, 2% SDS, 0.01% bromophenol blue and 10% glycerol. Twenty micrograms of denatured proteins were fractionated by 10% SDS-PAGE and transferred to a membrane. The membrane was blotted with Arabidopsis PLDα1 antibodies as described by . Aliquots of 20 μg of native protein were used for the PLD enzyme assay as described previously .
Phospholipid profiling was performed on soybean developing seeds with ESI-MS/MS as described previously . In brief, five replicates of developing soybean seed samples were smashed in liquid nitrogen and transferred immediately to 3 mL of 85 °C hot isopropanol containing 0.01% butylated hydroxytoluene (BHT) for 15 min. Then, adding 1.5 mL of chloroform and 1.5 mL of methanol for mixing well and extraction on shaker for 30 min, and then 1.8 mL of H2O were added to each sample separately. After mixing and centrifugation, lower layer solvent containing lipids was transferred to a new tube. Repeat the extraction with 3 mL of chloroform/methanol (2:1) containing 0.01% BHT for five times for 10 min each time. The combined solvent extracts were washed with 1 mL of 1 M KCl and 1 mL of sterile ddH2O, successively, and then concentrated with nitrogen gas. The total extracts were transferred to 2.0 mL vial with Teflon-lined screw cap and dried completely. The methanol solved lipids were detected followed the published protocol . Five duplicates for each genotype or treatments were analyzed.
Quantitative reverse transcriptase-PCR (qRT-PCR) for gene expression analysis
qRT-PCR was done as previously described . Briefly, total RNA was isolated from seeds with RNA isolation kit (Bioteke Corporation). The first strand of cDNA was synthesized according to the supplier’s instructions (M-MLV First Strand kit, Life Technologies, Invitrogen, USA) and real-time RT-PCR was executed with primers listed as Additional file 2: Table S1 using the BioRad iQ5 (Bio-Rad, USA) with an SYBR green MIX (Premix system, NEWBIO INDUSTRY). Soybean ACTIN1 (Glyma.15g034000) was used as the internal positive control.
Seed germination assay
After storing seeds for 3 months in high temperature and humidity conditions, seeds of PLDα1KDs and wild-type control harvested at the same time in same size were selected for assay. Seeds were sterilized with chloride gas and washed in distilled water three times before germinating in plates containing wet sterile filter papers. Seeds were germinated in an incubator at 25 °C in 16 h day/8 h night. Seed germination was scored when the radical was elongated about 0.5 cm.
Analysis of seed hormones
Seed samples collected at different time during germination were immediately frozen in liquid nitrogen and stored at − 80 °C. Hormones were measured according to the method described previously, with modifications . Briefly, the liquid nitrogen-frozen samples were lyophilized, and 0.1 g of samples was extracted with 750 μL methanol:ddH2O:acetic acid (80:19:1, v/v/v) for three times in the dark and 4 °C. The supernatants were combined and dried under stream of nitrogen. Prior to UPLC/ESI-MS/MS analysis, the extracts were suspended in 250 μL methanol and sonicated for 10 min, and spun at 12,000 g for 10 min, and the supernatant was transferred into UPLC vials. Chromatographic separation was performed using the Agilent LC-20AD system (Agilent, Santa Clara, CA, USA) equipped with AB Sciex QTRAP® 5500 detector and a Zorbax × 300 SB-C18 (4.6 mm × 150 mm × 5 μm) column (Agilent, Palo Alto, CA, USA), as described previously .
Soybean proteins involved in lipid metabolism and their expression data were retrieved from Phytozome (http://phytozome.jgi.doe.gov/pz). For analyses, these genes were compared with homologues from Arabidopsis thaliana and other plant species retrieved from NCBI (http://www.ncbi.nlm.nih.gov/BLAST/). Phylogeny trees were constructed using the neighbor-joining tree method with MEGA6. The significance level of the neighbor-joining analysis was examined by bootstrap testing with 1000 repeats.
JZ planned and designed the research. GZ, SB, GW, YRZ, BC, YLZ performed the experiments and analyzed data. SB, GW, XW did soybean transformat screening. GZ, BC analyzed the data. JZ, GZ, XW wrote and edited the manuscript. All authors read and approved the final manuscript.
We thank Dr. Edgar B. Cahoon and Jamie Shipp (Center for Plant Science Innovation, University of Nebraska-Lincoln) for helps in generation of transgenic soybean plants and Dr. Ruth Welti (Kansas lipidomics center, Kansas State University) for parts of lipid analyses and critical reading the manuscript. The authors also thank lab members in Prof. Zhao’s lab for all assistances in experiments and data analyses.
The authors declare that they have no competing interests.
Availability of data and materials
All data supporting my findings can be available and found in the supplementary data; materials can be available for distribution upon request.
Consent for publication
The authors agree to the publication of this manuscript in the journal.
Ethics approval and consent to participate
No investigations were undertaken using humans/human samples in this study. No experimental animals were used to conduct any of the experiments reported in this manuscript. Our study did not involve endangered or protected species. No specific permits were required from the studies and Professor Jian Zhao should be contacted for future permissions.
This project is support in part by the Ministry of Science and Technology of China (2016YFD0100504), by Agriculture and Food Research Initiative (AFRI) Award No. [2016-67013-24429/project accession number 1007600] from the USDA National Institute of Food and Agriculture, and the National Natural Science Foundation of China (31670294).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- IPCC. Climate Change 2013: the physical science basis. Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change, Cambridge university press, Cambridge, United Kingdom and New York, NY, USA; 2013.Google Scholar
- Ray DK, Gerber JS, MacDonald GK, West PC. Climate variation explains a third of global crop yield variability. Nat Commun. 2015;6:5989.View ArticleGoogle Scholar
- Narayanan S, Prasad PVV, Fritz AK, Boyle DL, Gill BS. Impact of high night-time and high daytime temperature stress on winter wheat. J Agron Crop Sci. 2015;201:206–18.View ArticleGoogle Scholar
- Schauberger B, Archontoulis S, Arneth A, Balkovic J, Ciais P, Deryng D, et al. Consistent negative response of US crops to high temperatures in observations and crop models. Nat Commun. 2017;8:13931.View ArticleGoogle Scholar
- Narayanan S, Tamura PJ, Roth MR, Prasad PV, Welti R. Wheat leaf lipids during heat stress: I. High day and night temperatures result in major lipid alterations. Plant Cell Environ. 2016;39:787–803.View ArticleGoogle Scholar
- Shu Y, Tao Y, Wang S, Huang L, Yu X, Wang Z, et al. GmSBH1, a homeobox transcription factor gene, relates to growth and development and involves in response to high temperature and humidity stress in soybean. Plant Cell Rep. 2015;34:1927–37.View ArticleGoogle Scholar
- Devaiah SP, Pan X, Hong Y, Roth M, Welti R, Wang X. Enhancing seed quality and viability by suppressing phospholipase D in Arabidopsis. Plant J. 2007;50:950–7.View ArticleGoogle Scholar
- Lee J, Welti R, Roth M, Schapaugh WT, Li J, Trick HN. Enhanced seed viability and lipid compositional changes during natural ageing by suppressing phospholipase Dalpha in soybean seed. Plant Biotechnol J. 2012;10:164–73.View ArticleGoogle Scholar
- Lee J, Welti R, Schapaugh WT, Trick HN. Phospholipid and triacylglycerol profiles modified by PLD suppression in soybean seed. Plant Biotechnol J. 2011;9:359–72.View ArticleGoogle Scholar
- Zhao J. Phospholipase D and phosphatidic acid in plant defence response: from protein-protein and lipid-protein interactions to hormone signalling. J Exp Bot. 2015;66:1721–36.View ArticleGoogle Scholar
- Manan S, Chen B, She G, Wan X, Zhao J. Transport and transcriptional regulation of oil production in plants. Crit Rev Biotechnol. 2017;37:641–55.View ArticleGoogle Scholar
- Li M, Wei F, Tawfall A, Tang M, Saettele A, Wang X. Overexpression of patatin-related phospholipase AIIIdelta altered plant growth and increased seed oil content in camelina. Plant Biotechnol J. 2015;13:766–78.View ArticleGoogle Scholar
- Bates PD, Ohlrogge JB, Pollard M. Incorporation of newly synthesized fatty acids into cytosolic glycerolipids in pea leaves occurs via acyl editing. J Biol Chem. 2007;282:31206–16.View ArticleGoogle Scholar
- Bates PD, Durrett TP, Ohlrogge JB, Pollard M. Analysis of acyl fluxes through multiple pathways of triacylglycerol synthesis in developing soybean embryos. Plant Physiol. 2009;150:55–72.View ArticleGoogle Scholar
- Li M, Bahn SC, Fan C, Li J, Phan T, Ortiz M, et al. Patatin-related phospholipase pPLAIIIdelta increases seed oil content with long-chain fatty acids in Arabidopsis. Plant Physiol. 2013;162:39–51.View ArticleGoogle Scholar
- Bates PD, Fatihi A, Snapp AR, Carlsson AS, Browse J, Lu C. Acyl editing and headgroup exchange are the major mechanisms that direct polyunsaturated fatty acid flux into triacylglycerols. Plant Physiol. 2012;160:1530–9.View ArticleGoogle Scholar
- Wang L, Shen W, Kazachkov M, Chen G, Chen Q, Carlsson AS, et al. Metabolic interactions between the Lands cycle and the Kennedy pathway of glycerolipid synthesis in Arabidopsis developing seeds. Plant Cell. 2012;24:4652–69.View ArticleGoogle Scholar
- Lin YC, Liu YC, Nakamura Y. The choline/ethanolamine kinase family in Arabidopsis: essential role of CEK4 in phospholipid biosynthesis and embryo development. Plant Cell. 2015;27:1497–511.View ArticleGoogle Scholar
- Lu C, Xin Z, Ren Z, Miquel M, Browse J. An enzyme regulating triacylglycerol composition is encoded by the ROD1 gene of Arabidopsis. Proc Natl Acad Sci USA. 2009;106:18837–42.View ArticleGoogle Scholar
- Wang L, Kazachkov M, Shen W, Bai M, Wu H, Zou J. Deciphering the roles of Arabidopsis LPCAT and PAH in phosphatidylcholine homeostasis and pathway coordination for chloroplast lipid synthesis. Plant J. 2014;80:965–76.View ArticleGoogle Scholar
- Hu Z, Ren Z, Lu C. The phosphatidylcholine diacylglycerol cholinephosphotransferase is required for efficient hydroxy fatty acid accumulation in transgenic Arabidopsis. Plant Physiol. 2012;158:1944–54.View ArticleGoogle Scholar
- Wickramarathna AD, Siloto RMP, Mietkiewska E, Singer SD, Pan X, Weselake RJ. Heterologous expression of flax phospholipid:diacylglycerol cholinephosphotransferase (PDCT) increases polyunsaturated fatty acid content in yeast and Arabidopsis seeds. BMC Biotechnol. 2015;15:63.View ArticleGoogle Scholar
- Pan X, Siloto RM, Wickramarathna AD, Mietkiewska E, Weselake RJ. Identification of a pair of phospholipid:diacylglycerol acyltransferases from developing flax (Linum usitatissimum L.) seed catalyzing the selective production of trilinolenin. J Biol Chem. 2013;288:24173–88.View ArticleGoogle Scholar
- Maatta S, Scheu B, Roth MR, Tamura P, Li M, Williams TD, et al. Levels of Arabidopsis thaliana leaf phosphatidic acids, phosphatidylserines, and most trienoate-containing polar lipid molecular species increase during the dark period of the diurnal cycle. Front Plant Sci. 2012;3:49.View ArticleGoogle Scholar
- Wang X. Regulatory functions of phospholipase D and phosphatidic acid in plant growth, development, and stress responses. Plant Physiol. 2005;139:566–73.View ArticleGoogle Scholar
- Wang X, Chapman KD. Lipid signaling in plants. Front Plant Sci. 2013;4:216.PubMedPubMed CentralGoogle Scholar
- Hong Y, Zhao J, Guo L, Kim SC, Deng X, Wang G, et al. Plant phospholipases D and C and their diverse functions in stress responses. Prog Lipid Res. 2016;62:55–74.View ArticleGoogle Scholar
- Yang W, Wang G, Li J, Bates PD. Phospholipase Dzeta enhances diacylglycerol flux into triacylglycerol. Plant Physiol. 2017;174:110–23.View ArticleGoogle Scholar
- Fang C, Ma Y, Wu S, Liu Z, Wang Z, Yang R, et al. Genome-wide association studies dissect the genetic networks underlying agronomical traits in soybean. Genome Biol. 2017;18:161.View ArticleGoogle Scholar
- Manan S, Ahmad MZ, Zhang G, Chen B, Haq BU, Yang J, Zhao J. 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.View ArticleGoogle Scholar
- Wang L, Ma H, Song L, Shu Y, Gu W. Comparative proteomics analysis reveals the mechanism of pre-harvest seed deterioration of soybean under high temperature and humidity stress. J Proteomics. 2012;75:2109–27.View ArticleGoogle Scholar
- Libault M, Farmer A, Joshi T, Takahashi K, Langley RJ, Franklin LD, et al. An integrated transcriptome atlas of the crop model Glycine max, and its use in comparative analyses in plants. Plant J. 2010;63:86–99.PubMedGoogle Scholar
- Severin AJ, Woody JL, Bolon YT, Joseph B, Diers BW, Farmer AD, et al. RNA-Seq Atlas of Glycine max: a guide to the soybean transcriptome. BMC Plant Biol. 2010;10:160.View ArticleGoogle Scholar
- Galant A, Koester RP, Ainsworth EA, Hicks LM, Jez JM. From climate change to molecular response: redox proteomics of ozone-induced responses in soybean. New Phytol. 2012;194:220–9.View ArticleGoogle Scholar
- Oh M, Komatsu S. Characterization of proteins in soybean roots under flooding and drought stresses. J Proteomics. 2015;114:161–81.View ArticleGoogle Scholar
- Noctor G, Reichheld JP, Foyer CH. ROS-related redox regulation and signaling in plants. Semin Cell Dev Biol. 2017;80:3–12.View ArticleGoogle Scholar
- Heppard EP, Kinney AJ, Stecca KL, Miao GH. Developmental and growth temperature regulation of two different microsomal [omega]-6 desaturase genes in soybeans. Plant Physiol. 1996;110:311–9.View ArticleGoogle Scholar
- Tang GQ, Novitzky WP, Carol Griffin H, Huber SC, Dewey RE. Oleate desaturase enzymes of soybean: evidence of regulation through differential stability and phosphorylation. Plant J. 2005;44:433–46.View ArticleGoogle Scholar
- Li L, Wang X, Gai J, Yu D. Molecular cloning and characterization of a novel microsomal oleate desaturase gene from soybean. J Plant Physiol. 2007;164:1516–26.View ArticleGoogle Scholar
- Bilyeu K, Palavalli L, Sleper D, Beuselinck P. Mutations in soybean microsomal omega-3 fatty acid desaturase genes reduce linolenic acid concentration in soybean seeds. Crop Sci. 2005;45:1830–6.View ArticleGoogle Scholar
- Flores T, Karpova O, Su X, Zeng P, Bilyeu K, Sleper DA, Nguyen HT, Zhang ZJ. Silencing of GmFAD3 gene by siRNA leads to low alpha-linolenic acids (18:3) of fad3-mutant phenotype in soybean [Glycine max (Merr.)]. Transgenic Res. 2008;17:839–50.View ArticleGoogle Scholar
- Kim HU, Huang AH. Plastid lysophosphatidyl acyltransferase is essential for embryo development in Arabidopsis. Plant Physiol. 2004;134:1206–16.View ArticleGoogle Scholar
- Kim HU, Li Y, Huang AH. Ubiquitous and endoplasmic reticulum-located lysophosphatidyl acyltransferase, LPAT2, is essential for female but not male gametophyte development in Arabidopsis. Plant Cell. 2005;17:1073–89.View ArticleGoogle Scholar
- Chen B, Wang J, Zhang G, Liu J, Manan S, Hu H, Zhao J. Two types of soybean diacylglycerol acyltransferases are differentially involved in triacylglycerol biosynthesis and response to environmental stresses and hormones. Sci Rep. 2016;6:28541.View ArticleGoogle Scholar
- Dahlqvist A, Stahl U, Lenman M, Banas A, Lee M, Sandager L, et al. Phospholipid:diacylglycerol acyltransferase: an enzyme that catalyzes the acyl-CoA-independent formation of triacylglycerol in yeast and plants. Proc Natl Acad Sci USA. 2000;97:6487–92.View ArticleGoogle Scholar
- Oelkers P, Tinkelenberg A, Erdeniz N, Cromley D, Billheimer JT, Sturley SL. A lecithin cholesterol acyltransferase-like gene mediates diacylglycerol esterification in yeast. J Biol Chem. 2000;275:15609–12.View ArticleGoogle Scholar
- Zhang M, Fan J, Taylor DC, Ohlrogge JB. DGAT1 and PDAT1 acyltransferases have overlapping functions in Arabidopsis triacylglycerol biosynthesis and are essential for normal pollen and seed development. Plant Cell. 2009;21:3885–901.View ArticleGoogle Scholar
- Schlueter JA, Vasylenko-Sanders IF, Deshpande S, Yi J, Siegfried M, Roe BA, et al. FAD2 Gene family of soybean: insights into the structural and functional divergence of a paleopolyploid genome. Crop Sci. 2007;47(S1):14–26.Google Scholar
- Pham AT, Lee JD, Shannon JG, Bilyeu KD. Mutant alleles of FAD2-1A and FAD2-1B combine to produce soybeans with the high oleic acid seed oil trait. BMC Plant Biol. 2010;10:195.View ArticleGoogle Scholar
- Haun W, Coffman A, Clasen BM, Demorest ZL, Lowy A, Ray E, et al. Improved soybean oil quality by targeted mutagenesis of the fatty acid desaturase 2 gene family. Plant Biotechnol J. 2014;12:934–40.View ArticleGoogle Scholar
- Vogel G, Browse J. Cholinephosphotransferase and diacylglycerol acyltransferase (substrate specificities at a key branch point in seed lipid metabolism). Plant Physiol. 1996;110:923–31.View ArticleGoogle Scholar
- Chapman KD, Ohlrogge JB. Compartmentation of triacylglycerol accumulation in plants. J Biol Chem. 2012;287:2288–94.View ArticleGoogle Scholar
- Wang X. Plant phospholipases. Annu Rev Plant Physiol Plant Mol Biol. 2001;52:211–31.View ArticleGoogle Scholar
- Craddock CP, Adams N, Bryant FM, Kurup S, Eastmond PJ. Phosphatidic acid phosphohydrolase regulates phosphatidylcholine biosynthesis in Arabidopsis by Phosphatidic acid-mediated activation of CTP:phosphocholine cytidylyltransferase activity. Plant Cell. 2015;27:1251–64.View ArticleGoogle Scholar
- Hong Y, Pan X, Welti R, Wang X. The effect of phospholipase Dalpha3 on Arabidopsis response to hyperosmotic stress and glucose. Plant Signal Behav. 2008;3:1099–100.View ArticleGoogle Scholar
- Beachy RN, Chen ZL, Horsch RB, Rogers SG, Hoffmann NJ, Fraley RT. Accumulation and assembly of soybean beta-conglycinin in seeds of transformed petunia plants. EMBO J. 1985;4:3047–53.View ArticleGoogle Scholar
- Doyle JJ, Schuler MA, Godette WD, Zenger V, Beachy RN, Slightom JL. The glycosylated seed storage proteins of Glycine max and Phaseolus vulgaris. Structural homologies of genes and proteins. J Bio Chem. 1986;261:9228–38.Google Scholar
- Gritz L, Davies J. Plasmid-encoded hygromycin B resistance: the sequence of hygromycin B phosphotransferase gene and its expression in Escherichia coli and Saccharomyces cerevisiae. Gene. 1983;25:179–88.View ArticleGoogle Scholar
- Garbarino JE, Belknap WR. Isolation of a ubiquitin-ribosomal protein gene (ubi3) from potato and expression of its promoter in transgenic plants. Plant Mol Biol. 1994;24:119–27.View ArticleGoogle Scholar
- Finer JJ, McMullen MD. Transformation of soybean via particle bombardment of embryogenic suspension culture tissue. In Vitro Cell Dev Biol Plant. 1991;27:175–82.View ArticleGoogle Scholar
- Schmidt MA, Tucker DM, Cahoon EB, Parrott WA. Towards normalization of soybean somatic embryo maturation. Plant Cell Rep. 2005;24:383–91.View ArticleGoogle Scholar
- Zhang G, Li P, Zhang W, Zhao J. Analysis of multiple soybean phytonutrients by near-infrared reflectance spectroscopy. Anal Bioanal Chem. 2017;409:3515–25.View ArticleGoogle Scholar
- Zhao J, Wang X. Arabidopsis phospholipase Dalpha1 interacts with the heterotrimeric G-protein alpha-subunit through a motif analogous to the DRY motif in G-protein-coupled receptors. J Biol Chem. 2004;279:1794–800.View ArticleGoogle Scholar
- Welti R, Li W, Li M, Sang Y, Biesiada H, Zhou HE, et al. Profiling membrane lipids in plant stress responses: role of phospholipase Dα in freezing-induced lipid changes in Arabidopsis. J Biol Chem. 2002;277:31994–2002.View ArticleGoogle Scholar
- Li P, Dong Q, Ge S, He X, Verdier J, Li D, Zhao J. Metabolic engineering of proanthocyanidin production by repressing the isoflavone pathways and redirecting anthocyanidin precursor flux in legume. Plant Biotechnol J. 2016;14:1604–18.View ArticleGoogle Scholar
- Wang J, Hou Q, Li P, Yang L, Sun X, Benedito VA, et al. Diverse functions of multidrug and toxin extrusion (MATE) transporters in citric acid efflux and metal homeostasis in Medicago truncatula. Plant J. 2017;90:79–95.Google Scholar