Triacylglycerol synthesis by PDAT1 in the absence of DGAT1 activity is dependent on re-acylation of LPC by LPCAT2
© Xu et al; licensee BioMed Central Ltd. 2012
Received: 21 July 2011
Accepted: 10 January 2012
Published: 10 January 2012
The Arabidopsis thaliana dgat1 mutant, AS11, has an oil content which is decreased by 30%, and a strongly increased ratio of 18:3/20:1, compared to wild type. Despite lacking a functional DGAT1, AS11 still manages to make 70% of WT seed oil levels. Recently, it was demonstrated that in the absence of DGAT1, PDAT1 was essential for normal seed development, and is a dominant determinant in Arabidopsis TAG biosynthesis.
Biochemical, metabolic and gene expression studies combined with genetic crossing of selected Arabidopsis mutants have been carried out to demonstrate the contribution of Arabidopsis PDAT1 and LPCAT2 in the absence of DGAT1 activity.
Through microarray and RT-PCR gene expression analyses of AS11 vs. WT mid-developing siliques, we observed consistent trends between the two methods. FAD2 and FAD3 were up-regulated and FAE1 down-regulated, consistent with the AS11 acyl phenotype. PDAT1 expression was up-regulated by ca 65% while PDAT2 expression was up-regulated only 15%, reinforcing the dominant role of PDAT1 in AS11 TAG biosynthesis. The expression of LPCAT2 was up-regulated by 50-75%, while LPCAT1 expression was not significantly affected. In vitro LPCAT activity was enhanced by 75-125% in microsomal protein preparations from mid-developing AS11 seed vs WT. Co-incident homozygous knockout lines of dgat1/lpcat2 exhibited severe penalties on TAG biosynthesis, delayed plant development and seed set, even with a functional PDAT1; the double mutant dgat1/lpcat1 showed only marginally lower oil content than AS11.
Collectively, the data strongly support that in AS11 it is LPCAT2 up-regulation which is primarily responsible for assisting in PDAT1-catalyzed TAG biosynthesis, maintaining a supply of PC as co-substrate to transfer sn-2 moieties to the sn-3 position of the enlarged AS11 DAG pool.
Triacylglycerols (TAGs) are the major storage lipids which accumulate in developing seeds, flower petals, anthers, pollen grains, and fruit mesocarp of a number of plant species [1, 2]. TAGs are thought to be not only the major energy source for seed germination but also essential for pollen development and sexual reproduction in many plants [3, 4]. In oilseeds, TAG bioassembly is catalyzed by the membrane-bound enzymes of the Kennedy pathway that operate in the endoplasmic reticulum . The biosynthesis of TAGs is catalyzed by the sequential acylation of the glycerol backbone involving three acyltransferases: glycerol-3-phosphate acyltransferase (GPAT), lyso-phosphatidic acid acyltransferase (LPAAT) and diacylglycerol acyltransferase (DGAT). DGAT catalyses the final acylation of DAG to give TAG, which has been suggested to be the rate-limiting step in plant lipid accumulation.
In the traditional Kennedy pathway DGAT was thought to be the only enzyme that is exclusively committed to TAG biosynthesis using acyl-CoA as its acyl donor. The first DGAT gene was cloned from mouse and is a member of the DGAT1 family, which had high sequence similarity with sterol: acyl-CoA acyltransferase .
We had previously characterized an EMS-induced mutant of Arabidopsis, designated AS11, which displayed a decrease in stored seed TAG, delayed seed development, and an altered fatty acid composition . We analyzed WT vs. AS11 lipid pools and Kennedy pathway enzyme activities in fractions isolated from green mid-developing seed, and performed parallel labeling of intact seeds at this developmental stage, with [14 C] acetate. We found that compared to WT, there was an increase in all fatty acids in the DAG pool of AS11 seeds at mid-development, and, to a lesser extent, an associated backup of fatty acids in the PC pool. DAG was elevated from 1% in WT to 10-12% in AS11 and PC pools were elevated from about 2% in WT, to 8-12% in AS11. Cell-free fractions from WT and AS11 green seeds at mid-development were compared for their ability to incorporate [14 C]-18:1-CoA into glycerolipids in the presence of G-3-P. Proportions of labeled LPA and PA formed during the incubation period were similar in WT vs AS11, indicating that the activities of the Kennedy pathway enzymes GPAT and LPAAT (EC 22.214.171.124) were relatively unaffected in the AS11 mutant. However, the proportion of labeled TAG was much lower and that of DAG was much higher in AS11. The TAG/DAG ratio was therefore consistently 3-5-fold lower in AS11 compared to WT at all developmental stages (early-, mid- and late development). Furthermore, the ratio of 18:3/20:1 dramatically increased about 7-10 fold .
Cumulatively, this data suggested a lesion in DGAT1 which was subsequently demonstrated upon cloning the mutated gene from AS11. There is an 81 bp in-frame insertion consisting entirely of exon 2 in the transcript from AS11. The exon 2 in the repeat is properly spliced, thus the alteration of the transcript does not disturb the reading frame. However, this additional exon 2 sequence in the AS11 transcript would result in an altered DGAT protein with a 27 amino acid insertion (131SHAGLFNLCVVVLIAVNSRLIIENLMK157) . It is important to note that two other labs independently and simultaneously cloned the A. thaliana DGAT1 [8, 9].
Studies manipulating the expression of DGAT1 followed: We demonstrated that expression of the Arabidopsis DGAT1 cDNA in a seed specific manner in the AS11 mutant restored wild type levels of TAG and VLCFA content. The acyl distribution, specifically, the sn-3 composition of the TAGs, was also restored to WT proportions. Furthermore, overexpression of the Arabidopsis DGAT1 in wild type plants led to an increase in seed oil content and seed weight . Over the past 10 years, DGAT1 expression has been genetically manipulated to produce Brassica napus prototypes containing increased seed oil [11, 12].
A second family of DGAT genes (DGAT2), first identified in the oleaginous fungus Morteriella ramanniana, has no sequence similarity with DGAT1 . A human DGAT2 and plant DGAT2s from tung and castor were subsequently identified by Cases et al. , Shockey et al.  and Kroon et al. , respectively. The putative DGAT2 from Arabidopsis has been studied by several labs including ours; functional expression in yeast has not been successful, and therefore whether it is a true functioning DGAT is still in question. Notably, an Arabidopsis dgat2 knockout mutant has a wild-type seed oil content and fatty acid composition .
TAG can also be formed by an acyl-CoA-independent enzyme, phosphatidylcholine: diacylglycerol acyltransferase (PDAT), in which the transfer of an acyl group from the sn-2 position of PC to the sn-3 position of DAG yields TAG and sn-1 lyso-PC [18, 19]. During the exponential growth phase in yeast, PDAT1 is a major determinant in TAG synthesis. In Arabidopsis, two close homologs to the yeast PDAT gene have been identified: PDAT1 At5g13640 and PDAT2 At3g44830 . Mhaske et al.  isolated and characterized a knockout mutant of Arabidopsis which has a T-DNA insertion in the PDAT1 locus At5g13640 (PDAT1, EC 126.96.36.199). Lipid analyses were conducted on this mutant to assess the contribution of PDAT1 to seed lipid biosynthesis; surprisingly, and in contrast to the situation in yeast, the oil content and composition in seeds did not show significant changes in the mutant. At the time, these results were interpreted to indicate that PDAT1 activity as encoded by At5g13640 is not a major determining factor for TAG synthesis in Arabidopsis seeds.
Nonetheless, because the Arabidopsis DGAT1 mutant AS11 shows only a 30-35% decrease in oil content [6, 9], it was apparent that other enzymes must contribute to oil synthesis in the developing seed . Thus, an examination of the contribution of DGAT2, PDAT2 or PDAT1 to oil deposition in an AS11 background was studied by performing double mutant crosses with AS11 . While the dgat2-ko line has no oil phenotype, homozygous double mutants from cross of AS11 with dgat2-ko mutant showed an oil fatty acid profile similar to AS11. The same pattern was observed with the pdat2-ko mutant alone and in crosses of the pdat2-ko mutant with AS11. In contrast, while the pdat1-ko has no oil or fatty acid composition phenotype, crosses of the pdat1-ko with AS11 were embryo-lethal in the double homozygous condition; only heterozygous lines produced by having expression of the pdat1 or dgat1 gene only partially inhibited using RNAi, allowed an examination of the double mutants. These detailed studies resulted in the finding that DGAT1 and PDAT1 have overlapping functions in both embryo development and TAG biosynthesis in the developing seed and in pollen. When DGAT1 is compromised in AS11, it is PDAT1, and not DGAT2 or PDAT2, that is responsible for the remaining 65-70% of TAG which is synthesized. This finding suggested a major, perhaps dominant role of PDAT1 in this process .
Recently, a castor bean-specific PDAT, PDAT1-2, was cloned and found to be highly expressed in developing seeds and localized in the ER, similar to the castor FAH12 hydroxylase. Transgenic Arabidopsis co-expressing the castor PDAT1-2 and FAH12 showed enhanced ricinoleate accumulation to up to 25% in TAGs (compared to 17% in FAH12 -only transgenics) [23, 24]. This study suggests that specialized PDATs may play a significant role in channeling PC-synthesized unusual fatty acids such as ricinoleic (from castor), or epoxy fatty acids (e.g. from Vernonia galamensis), into TAGs.
The discovery of the critical role PDAT1 has in TAG synthesis  suggested the importance of re-acylation of lysophospholipids and especially LPC, since it is produced as a result of PDAT1 activity. We hypothesized that in a situation where the responsibility for TAG synthesis is shifted to PDAT1 such as in the AS11 dgat1 mutant, profound changes in activity of re-acylation enzymes would be evident. As candidates of genes coding such enzymes, we proposed the two Arabidopsis genes LPLAT1 (At1g12640) and LPLAT2 (At1g63050), characterized by Ståhl et al.  as possessing broad specificity lysophospholipid acyltransferase activity. Although both of the genes to some extent could acylate a range of different lysophospholipids, they were shown to have a strong preference towards LPC. For this reason we chose to name the two genes LPCAT1 (At1g12640) and LPCAT2 (At1g63050), in the current study.
Here we report the further genetic and biochemical characterization of the AS11 dgat1 mutant. During the course of microarray and qRT-PCR studies of AS11 vs WT gene expression in mid-developing siliques, we found that LPCAT2, encoding acyl-CoA:lysophosphatidylcholine acyltransferase 2 (EC 188.8.131.52), was up-regulated while LPCAT1 was not affected. By a series of biochemical studies and key crosses of AS11 with either lpcat1 or lpcat2 knockout mutants, we determined that LPCAT2 (and much less so LPCAT1) is most critical for TAG synthesis in the AS11 mutant, primarily to maintain the PC pool for TAG assembly primarily catalyzed by PDAT1.
Results and discussion
Summary of the AS11mutant developmental and oil phenotypes
The AS11 mutant line was about one week behind WT in bolting and entering the generative phase and thus, under our growing conditions, AS11 seed set was also delayed to four weeks post-anthesis instead of three, as typically observed in WT. For comparison, we studied another DGAT1 mutant, which we designated A7, which is a homozygous SALK line (Salk 039456) with a T-DNA insertion in the last exon of the same DGAT1 gene (At2g19450). A7 shows a developmental delay similar to that exhibited by AS11 (Additional file 1: Figure S1).
Heterologous expression of mutated Arabidopsis DGAT1s from both A7 and AS11 in Yeast Strain H1246MATα
Thus our results strongly support that DGAT1 is not functioning in the developing seed of the AS11 mutant line. Clearly, the (radiolabeled) TAG formation observed in protein fractions from AS11 was coming from another path and not via reduced DGAT1 catalysis.
Given the importance of PDAT1 in oil biosynthesis in Arabidopsis as we earlier defined , and combining this new information with the 14C TAG biosynthesis results in assays of WT vs AS11 protein fractions reported above, it raised the question that, if not from DGAT1, DGAT2 nor PDAT2, what biochemical steps besides PDAT1 may be critical for TAG biosynthesis?
AS11 microarray and qRT-PCR analyses
Some general observations from these combined gene expression studies follow: In AS11, FAE1 is down-regulated, FADs 2 & 3 are up-regulated; this is consistent with the AS11 acyl composition profile (reduced 20:1 and elevated 18:3). DGAT1 expression was not significantly affected. Interestingly, LPCAT2 was up-regulated by an average of 65%, while LPCAT1 expression was indifferent. PDAT2 was only marginally affected, but PDAT1 expression was up-regulated by an average of 62% compared to WT.
Given the critical role of PDAT1 in embryo development and TAG deposition in both pollen and seeds , we were interested in the relative differences in LPCAT1 and LPCAT2 expression between the AS11 mutant and WT. We hypothesized that an acyl-CoA-dependent LPCAT may be critical to maintain the PC pool as one of the co-substrates for PDAT1-catalyzed TAG synthesis, particularly in the absence of a functional DGAT1, and performed a series of genetic and metabolic studies to examine these relationships.
RT-PCR Characterization of the lpcat1, lpcat2 and dgat1 (AS11) mutants
Based on these trends we were confident that in addition to DGAT1 and PDAT1, LPCAT plays a significant role in TAG biosynthesis, and postulated that it is important not only for membrane development and acyl turnover therein, but also to replenish the supply of PC for PDAT1-catalyzed TAG biosynthesis in the developing seed. To resolve this question we needed to study what occurs when both DGAT1 and LPCAT expression are co-disrupted.
It had previously been shown, following expression in a yeast lca1Δ mutant background, that both Arabidopsis LPCAT1 (At1g12640) and LPCAT2 (At1g63050) are functional LPCATs. In the presence of sn-1 [14 C] 16:0-LPC and 18:1-CoA (the same combination of substrates used in the current study), expressed LPCAT1 and LPCAT2 proteins had activities of 12 nmol/min/mg/mg protein and 18 nmol/min/mg protein, respectively .
As confirmed by RT-PCR (discussed above), lpact1 (At1g12640) and lpcat2 (At1g63050) are T-DNA insertion knockouts for their respective genes, while AS11 encodes a mutated and non-functional (Figure 2) dgat1. To examine the relationship between these two LPCATs and DGAT1 in TAG biosynthesis, we performed crosses of the AS11 dgat1 with the lpcat1 or the lpcat2 T-DNA insertion knock-out mutants and characterized the hemizygous/homozygous and double knockout seed oil profiles.
dgat1 × lpcat1
dgat1 × lpcat2
The acyl composition of the dgat1/lpcat2 He/H lines was again essentially like the LPCAT2 NS (Figure 12). In contrast, the dgat1/lpcat2 H/He lines have an AS11-like reduction in oil content which is about 70% of that found in the WT and the LPCAT2 null segregant (Figure 11). Interestingly, there is no developmental penalty in this heterozygous combination and seed set is normal as in the AS11 background, provided expression of LPCAT2 is in the heterozygous state (Additional file 1: Figure S3). Not unexpectedly, the acyl composition of the oil in these lines is consistent with the AS11 profile (Figure 12).
Summarizing the mutant crossing experiments, in lpcat2 × dgat1, oil synthesis, vegetative growth and seed set are severely penalized. Compared to the relative oil content in dgat1 alone (75% of WT), there is only a 5% further relative decrease in dgat1/lpcat1 oil (70% of WT) which suggests that LPCAT2 is the major determinant and can almost fully compensate for the loss of LPCAT1. In contrast, in the co-incident absence of LPCAT2 and DGAT1 (lpcat2 × dgat1), relative oil content drops to 30% of WT, which suggests that LPCAT1 can only partially compensate for the loss of LPCAT2 with respect to oil. Metabolic and microarray studies in AS11 alone support this hypothesis as LPCAT2 expression and LPCAT activity are up-regulated in the absence of DGAT1. In addition, the acyl profiles indicate that any altered oil composition phenotype observed in the lpcat2 x dgat1 is predominantly influenced by the dgat1 lesion. Collectively, the data from the crosses of AS11 with lpcat1 or lpcat2 mutants strongly suggest that in an AS11 background, while both LPCAT2 and LPCAT1 can contribute to PDAT1-catalyzed TAG synthesis, it is LPCAT2 up-regulation (and not LPCAT1), which is primarily responsible for assisting in PDAT1-catalyzed TAG biosynthesis by supplying an enhanced pool of PC as co-substrate to transfer sn-2 moieties to the sn-3 position of the enlarged AS11 DAG pool. Nonetheless, while LPCAT2 is the most important for maintaining PDAT1-catalyzed TAG synthesis in the AS11 background, the partial over-lapping function of LPCATs 1 and 2 identified in this study require further investigation to determine their relative roles in oil synthesis when DGAT1 is functional in a WT background.
The recent discovery of a phosphatidylcholine: diacylglycerol cholinephosphotransferase (PDCT), encoded by the Arabidopsis ROD1 gene, indicates that the issue of DAG-PC interconversion during TAG biosynthesis is considerably more complex than previously thought . The PDCT enzyme catalyzes transfer of the phosphocholine headgroup from PC to DAG and is thought to be a major reaction for the transfer of 18:1 into PC for desaturation by FAD2 and FAD3, as well as enabling the reverse transfer of 18:2 and 18:3 into the TAG synthesis pathway. In the rod1 mutant, the accumulation of PUFAs in TAGs is reduced by 40%. However, in AS11 where 18:3 is strongly increased, and notably, shows a strong stereospecific shift to the sn-3 position of TAGs , it is unclear whether there is a strong involvement of PDCT; were it so, one would expect a more even distribution of 18:3 at the sn-2 and sn-3 positions. However, the strong re-distribution of 18:3 proportions to the sn-3 and away from the sn-2 position in AS11 suggests that the major shift in positional composition has to do with PDAT1 catalysis. Thus there are still many questions to address with respect to the roles ROD1 and other PC-related metabolic steps such as CPTase might play in oil synthesis in Arabidopsis and in AS11 in particular.
This adjustment in the TAG assembly route in AS11 is aided by enhanced LPCAT2 activity to supply the additional PC required by PDAT1. Given these findings and the implicit plasticity of Arabidopsis TAG assembly mechanisms to overcome critical bottlenecks, it will be interesting to determine the relative roles played by DGAT1 and [PDAT1 + LPCAT2] during TAG assembly in higher oilseeds (e.g. canola, soybean, sunflower, safflower, flax).
Plant materials and growth conditions
Arabidopsis lines, WT (ecotype Columbia-0), and AS11 (CS3861) EMS-induced mutant, and A7 (SALK_039456), lpcat1 (SALK_123480) and lpcat2 (SAIL_357_H01) T-DNA insertion mutant lines (all in a Columbia background) were obtained from the Salk Institute via the Nottingham Arabidopsis Stock Centre (University of Nottingham, UK). Seeds of these lines and progeny from genetic crosses were grown in a growth chamber at 22°C with photoperiod of 16 h light (120 μE·m-2·s-1) and 8 h dark.
Arabidopsis siliques containing mid-green developing embryos (pooled siliques at stages 6- 8, as designated on the University of Toronto Arabidopsis eFP Expression Browser at: http://bar.utoronto.ca/affydb/cgi-bin/affy_db_exprss_browser_in.cgi; see e.g. seed staging in LPCAT2 developmental map at http://bbc.botany.utoronto.ca/efp/cgi-bin/efpWeb.cgi?primaryGene=AT1G63050&modeInput=Absolute), were harvested and frozen at -80°C for lipid analyses, enzyme assays and DNA and RNA extraction.
Preparation of total lipid extracts (TLEs) and study of lipid classes, determination of oil content and acyl composition in seeds of WT and the AS11, lpcat1 and lpcat 2 mutant lines and progeny from crosses were performed as described previously [26, 36]. In all cases, the data represent the averages of three-five determinations.
Preparation of Arabidopsis protein fractions
In general, enzyme preparations were made from 200 Arabidopsis siliques of AS11 and WT containing mid-green developing seeds (stages 7-8 as defined above), and immediately powdered with liquid nitrogen in a mortar and pestle. Grinding medium (100 mM HEPES-KOH, pH 7.4 containing 0.32 M Sucrose, 1 mM EDTA, and 1 mM dithiothreitol; 8 mL/50 siliques) along with 65 mg polyvinylpolypyrrolidone was immediately added, and grinding continued on ice for 5 min. The slurried cell-free homogenate was filtered through two layers of Miracloth (Calbiochem, La Jolla, CA), centrifuged at 3,000 × g for 5 min, the pellet discarded and the supernatant re-centrifuged at 15,000 × g for 30 min. The supernatant was re-centrifuged at 100,000 × g for 1 h and the resultant pellet was resuspended in 2 mL of grinding medium, probe-sonicated on ice for 30 s and protein concentrations were determined using BioRad™ reagent based on the method of Bradford . Protein concentrations were normalized to the same value for WT and AS11 in each experiment.
Assay of TAG synthesis activity in Arabidopsis mid-developing seeds
TAG assembly assays were conducted at pH 7.4, with shaking at 100 rev/min in a water bath at 30°C for 60 min. Assay mixtures (500 μL final volume) contained 100-300 μg protein normalized as described above, 90 mM HEPES-NaOH, 0.5 mM ATP, 0.5 mM CoASH, 1 mM MgCl2 in the presence of 100 μM [1-14 C] sn-1,2 diolein in 0.02% Tween-20 (specific activity 10 nCi/nmol; pre-purified by TLC on 5% borate silica G plates) and 18 μM unlabeled 18:1-CoA. Reactions were stopped using Isopropanol:CH2Cl2 (2:1:) v/v, and the TLE prepared as described previously [26, 36]. The 14 C-labeled products were resolved by TLC on silica gel G plates developed in hexane: diethyl ether: acetic acid (70:30:1 v/v/v), the 14 C-triolein band visualized on a Bioscan AR-2000 radio-TLC scanner using Win-Scan 2D© software (Bioscan Inc., Washington DC) and the band scraped and quantified on a scintillation counter.
Expression of mutated AtDGAT1 in yeast
The mutated AtDGAT1 (from AS11 as well as A7) in pYES2.1/NT B plasmid were transformed into a quadruple yeast mutant H1246MATα  using the S.c. EasyComp™Transformation Kit (Invitrogen, San Francisco, CA). Yeast cells transformed with pYES2.1/NT B plasmid containing the WT DGAT1 or with empty plasmid were used as positive or negative controls, respectively. Transformants were selected by growth on synthetic complete medium lacking uracil (SC-ura), supplemented with 2% (w/v) glucose. The colonies were transferred into liquid SC-ura with 2% (w/v) glucose and grown at 30°C overnight. The overnight culture was diluted to an OD 0.4 in induction medium (SC-ura + 2% Galactose + 1% Raffinose), and were induced for 24 h at 30°C. The yeast cells were collected and broken with glass beads using a Beadbeater™. The protein concentrations of the yeast cell lysates were normalized using the Biorad™ assay and assayed for DGAT activity. DGAT assays were conducted at pH 7.4, with shaking at 100 rev/min in a water bath at 30°C for 60 min. Assay mixtures (500 μL final volume) contained 100 μg of lysate protein, 90 mM HEPES-NaOH, 100 μM sn-1,2 diolein or sn-1,2 dierucin (pre-purified by TLC on 10% borate silica H plates and emulsified in 0.02% Tween-20), and 18 μM 1-14 C 18:1-CoA (specific activity 10 nCi/nmol) as the acyl donor. Reactions were stopped using isopropanol:CH2Cl2 (2:1:) v/v, and extraction of the lipid fraction performed as described previously [26, 36]. The 14 C-labeled TAGs were resolved by TLC on silica gel G plates developed in hexane: diethyl ether: acetic acid (70:30:1 v/v/v), the 14 C-triolein band visualized on a Bioscan AR-2000 radio-TLC scanner using Win-Scan 2D© software (Bioscan Inc., Washington DC) and the band scraped and quantified on a scintillation counter.
The yeast cell lysates were run on a 10% Tris-HCl SDS-PAGE gel, the proteins were then transferred to a nitrocellulose membrane (Nitrobind, Thermo Fisher Scientific, Waltham, MA). The membrane was blocked in PBST (phosphate buffered saline containing 0.5% Tween 20) containing 4% skim milk for 60 min, and then incubated with the primary antibody, Anti-Xpress (epitope-tagged) antibody (Invitrogen, San Francisco, CA) diluted to 1:5000 with PBST containing 2% skim milk, for 60 min. The membrane was submitted to three washes with PBST followed by three washes with PBS to remove any unbound antibody. Next, the membrane was incubated with a goat anti-mouse IgG peroxidase antibody (Sigma-Aldrich, A2554, St. Louis, MO), diluted to 1:5000 with PBST containing 2% skim milk, for 60 min. The membrane was washed three times with PBST followed by three times with PBS, then the proteins were detected using the Amersham ECL Plus Western Blotting Detection Kit (GE Healthcare Life Sciences) with a 10 min exposure time.
LPCAT assays were conducted at pH 7.4, with shaking at 100 rev/min in a water bath at 30°C in a 20-minute time course. Assay mixtures (500 μL final volume) contained 300 μg protein normalized as described above, 90 mM HEPES-NaOH, 0.5 mM ATP, 0.5 mM CoASH, 1 mM MgCl2 in the presence of 6 μM L-α-palmitoyl- [1- 3H methyl] lyso-3-phosphatidylcholine (specific activity 60 μCi/nmol; 1.48-2.22 TBq/mmol) and 18 μM [1-14 C] 18:1-CoA (specific activity 10 nCi/nmol; 0.37 GBq/mmol). Reactions were stopped at each time point by adding isopropanol: CH2Cl2 (2:1:) v/v, and the TLE prepared as described previously [26, 36]. The 14 C and 3H-labeled products were resolved by 3D TLC performed as follows: The first 2 dimensions were run as described by Yokoyama et al. ; 1st D: CHCl3: MeOH: Formic Acid (88%):H2O 60: 30: 9 : 2; 2nd D: CHCl3: MeOH: Ammonia solution (28%): H2O 50: 40: 7: 3; 3D: The 3 rd D: 100% ethyl ether. Radiolabeled spots corresponding to standards of LPC and PC were identified using a radio-TLC scanner and the bands scraped and counted on a scintillation counter with a dual 3H/14 C isotope measurement program.
In vivo feeding experiments
Twenty bolted stems with mid-developing siliques of both WT and AS11 plants were harvested under water and then immediately placed in a solution containing 1 μCi of 14 C sodium acetate in 100 μL of water and incubated at room temperature in a fume hood. Once most of this solution was taken up, the plants were supplied with equal aliquots of distilled water during the chase period. At 0, 7, 24 and 30 h, 5 bolted stems were harvested, siliques counted and weighed and then ground with a polytron and a TLE extraction performed as described by Taylor et al. . The lipid extract was resuspended in 1 mL CHCl3:MeOH 2:1. 10 μL of the solution was counted and the remainder spotted on TLC plates, 3D TLC performed, radiolabeled spots resolved on a radio-TLC scanner, silica bands scraped and radioactivity measured with a scintillation counter.
qRT-PCR comparison of AS11 and WT expression of key lipid genes
Triplicate biological samples of total RNA were extracted using the method of Wang and Vodkin  from Arabidopsis siliques of AS11 and WT containing mid-green developing seeds (stages 7-8 as defined above) and used for qRT-PCR . One microgram of total RNA was reverse-transcribed using SuperScript® II Reverse Transcriptase (Invitrogen) according to the manufacturer's protocol. The resulting cDNA was then amplified by PCR using Taq DNA polymerase (Invitrogen). PCR conditions comprised an initial cycle of 94°C for 3 min; followed by 30 cycles of 94°C for 30 s, 56°C for 30 s, and 72°C for 1 min; then 72°C for 10 min to complete the reaction. Concentrations for the individual samples within each group (AS11 and WT) were normalized using 18S rRNA levels. The gene-specific primers (listed in Additional file 2: Table S1) were designed according to the target gene sequences as annotated in Genbank™ (http://www.ncbi.nlm.nih.gov/genbank/). Amplified RT-PCR products were resolved by electrophoresis on 1% agarose gels and gel photos taken. The densities of the PCR bands were analyzed using ImageJ image analysis software (National Institutes of Health, USA http://rsbweb.nih.gov/nih-image/) and quantified relative to the "housekeeping gene" 18S rRNA transcript signal. Then the ratio of AS11:WT gene expression level was calculated. Specifically, for each gene of interest, first we calculated the ratio of AS11 density/18S rRNA density (= "AS11 ratio") and WT density/18S rRNA density (= "WT ratio"). The difference in specific gene expression in AS11 vs. WT, was then calculated as ("AS11 ratio")/("WT ratio") = V. If V >1.0, then the gene is up-regulated in AS11 relative to WT; if V <1.0, then the gene is down-regulated in AS11 relative to WT. The % relative increase or decrease is then = ([V × 100]-100).
Affymetrix microarray analysis
The same triplicate biological samples prepared for the qRT-PCR study (above) were sub-sampled for microarray analysis. Affymetrix microarray hybridizations using the Ath1 whole genome array, containing probe sets representing ~22,800 genes, were performed using three biological replicate samples for each genotype. Labeling, hybridization, and scanning were performed by the Affymetrix Gene Chip Facility at the University of Toronto (http://www.csb.utoronto.ca/resources/facilities/affymetrix-genechip). Data analysis was performed using GeneSpring™ software version 7.2. To identify key lipid genes that were differentially expressed between the two genotypes, a per-gene normalization was applied to the values and a parametric test was performed. Genes that exhibited a false discovery rate of p < 0.05 and passed the minimum signal and fold-change threshold were determined to be differentially expressed. Reported candidate genes were selected by comparing those listed by Beisson et al.  with our microarray and qRT-PCR data. The "fold change" was converted to % change. For the group of genes detailed in Figure 3, the final changes ranged from -60% to + 150%.
Confirmation of lpcat1 and lpcat2 T-DNA insertion mutation lines
Arabidopsis insertion mutation lines  SALK_123480 putative for LPCAT1 (At1g12640) and SAIL_357_H01 putative for LPCAT2 (At1g63050) were identified in the Salk Institute T-DNA insertion library database (http://signal.salk.edu/cgi-bin/tdnaexpress), and seeds were obtained from the Nottingham Arabidopsis Stock Centre (University of Nottingham, UK). According to annotation in the database SALK_123480 contains a T-DNA insertion in the middle of the 7th exon of the LPCAT1 gene and SAIL_357_H01 a T-DNA insertion in the 8th exon of the LPCAT2 gene. Individual plants homozygous for a T-DNA insertion in each the LPCAT1 or LPCAT2 genes were identified by PCR screening using primers SALK_123480LP, SALK_123480RP and SALK_LBb1 (Additional file 2: Table S2) for LPCAT1 and SAIL_357_H01LP, SAIL_357_H01RP and SAIL_LB2 (Additional file 2: Table S2) for LPCAT2. Individual plants from each mutant line lacking a T-DNA insertion in the LPCAT1 or LPCAT2 genes (null segregants) were also identified in the same PCR primer set. Annotation of lines from each set of crosses is as designated in Additional file 2: Table S3.
Transcript analyses of lpcat1, lpcat2, dgat1 mutant lines and crosses
lpcat1 (See Figure 4)
Total RNA was isolated from developing siliques from control and SALK_12340 T-DNA insertion mutation plants using CONCERT Plant RNA Reagent (Invitrogen, San Francisco, CA) according to the technical specifications given by the manufacturer. The RNA was then DNAse-treated with TURBO™ DNase (Ambion, Huntingdon, UK)) and used (1 mg) for synthesize first strand cDNA with Maxima First Strand cDNA Synthesis Kit (Fermentas.Thermo Fisher Scientific, St. Leon-Rot, Germany). The cDNA was amplified by PCR using Taq DNA Polymerase (Sigma-Aldrich, St. Louis, MO) and the primers SALK_12640-1 s and SALK_12640-6r as listed in Additional file 2: Table S2. PCR conditions comprised an initial cycle of 94°C for 1 min; followed by 40 cycles of 94°C for 20 s, 64°C for 20 s, and 72°C for 35 s; then 72°C for 8 min to complete the reaction. PCR samples were separated on 1% agarose gel using λ-DNA digested with PstI as size marker.
lpcat2 and AS11 dgat1 (See Figure 5)
Total RNA was isolated from developing siliques from WT control, AS11(dgat1), SAIL_357_H01 insertion mutation (lpcat2) and AS11/lpcat2 plants by the method of Wang and Vodkin . One microgram of total RNA was reverse-transcribed using SuperScript® II Reverse Transcriptase (Invitrogen) according to the manufacturer's protocol. The resulting cDNA was then amplified by PCR using Taq DNA polymerase (Invitrogen). PCR conditions comprised an initial cycle of 94°C for 3 min; followed by 30 cycles of 94°C for 30 s, 56°C for 30 s, and 72°C for 1 min; then 72°C for 10 min to complete the reaction. Concentrations for the individual samples within each group (AS11 and WT) were normalized using 18S rRNA levels. Gene specific primers spanning the T-DNA insertional region were used to check for gene transcription. For DGAT1 transcript analysis in lines of WT, the AS11 (dgat1) mutant, lpcat2 mutant and crosses of lpcat2 with dgat1, the primer pairs were TAG1-mut-primerA and TAG1-mut-primerB. For LPCAT2 transcript analysis in these same lines, the primer pairs were LPCAT2-F2 and LPCAT2-R2 as listed in Additional file 2: Table S2.
Creating double mutants of dgat1 × lpcat1 and dgat1 × lpcat2
Crosses between the AS11 mutant and lpcat1 and lpcat2 T-DNA insertion mutation lines, respectively, were made and F1 plants heterozygous for AS11 and the insertion mutations were identified by PCR using primers listed in Additional file 2: Table S2. F2 seeds segregating for the mutations were planted and screened by PCR. All PCR screening was done as above for LPCAT1 and LPCAT2 and for the AS11 mutation vs WT using primers dgat1-mut-primerA and dgat1-mut-primerB (Additional file 2: Table S2), as designed by Zou et al. . Identification of individuals homozygous for both the AS11 mutation and lpcat2 insertion mutation was only possible after growing a segregating F2 seed population (homozygous for AS11 and heterozygous for lpcat2 mutant) on agar media containing a 1/3 strength MS and 1% sucrose (Additional file 2: Table 3).
This work was partially supported by a grant to D.C.T. from Genome Canada: Genome Prairie Phase II and by funding to A.S.C. from FORMAS (The Swedish Research Council for Environment, Agricultural Sciences and Spatial Planning). The authors also acknowledge, with thanks, the kind assistance of Drs. Mark Smith and Alison Ferrie for their careful reviews of the manuscript and Dr. John Ohlrogge for helpful discussions during the course of this work. This is National Research Council of Canada Publication No. 50191. We thank the Salk Institute Genomic Analysis Laboratory for providing the sequence-indexed Arabidopsis T-DNA insertion mutant and the Nottingham Arabidopsis Stock Centre for providing seeds. Funding for the SIGNAL indexed insertion mutant collection was provided by the National Science Foundation.
- Stymne S, Stobart AK: Triacylglycerol biosynthesis. The biochemistry of plants. Lipids: structure and funchtion. Volume 9. Edited by: Stumpf PK. New York: Academic Press; 1987:175-214.Google Scholar
- Murphy DJ, Vance J: Mechanisms of lipid-body formation. Trends Biochem Sci. 1999, 24: 109-115. 10.1016/S0968-0004(98)01349-8.PubMedView ArticleGoogle Scholar
- Wolters-Arts M, Lush WM, Mariani C: Lipids are required for directional pollen-tube growth. Nature. 1998, 392: 818-821. 10.1038/33929.PubMedView ArticleGoogle Scholar
- Zheng Z, Xia Q, Dauk M, Shen W, Selvaraj G, Zou J: Arabidopsis AtGPAT1, a member of the membrane-bound glycerol-3-phosphate acyltransferase gene family, is essential for tapetum differentiation and male fertility. Plant Cell. 2003, 15: 1872-1887. 10.1105/tpc.012427.PubMedPubMed CentralView ArticleGoogle Scholar
- Cases S, Smith SJ, Zheng YW, Myers HM, Lear SR, Sande E, Novak S, Collins C, Welch CB, Lusis AJ, Erickson S, Farese RV: Identification of a gene encoding an acyl CoA: diacylglycerol acyltransferase, a key enzyme in triacylglycerol synthesis. Proc Natl Acad Sci USA. 1998, 95: 13018-13023. 10.1073/pnas.95.22.13018.PubMedPubMed CentralView ArticleGoogle Scholar
- Katavic V, Reed DW, Taylor DC, Giblin EM, Barton DL, Zou J, MacKenzie SL, Covello PS, Kunst L: Alteration of seed fatty acid composition by an ethyl methanesulfonateinduced mutation in arabidopsis thaliana affecting diacylglycerol acyltransferase activity. Plant Physiol. 1995, 108: 399-409. 10.1104/pp.108.1.399.PubMedPubMed CentralView ArticleGoogle Scholar
- 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-653. 10.1046/j.1365-313x.1999.00555.x.PubMedView ArticleGoogle Scholar
- Hobbs DH, Lu C, Hills MJ: Cloning of a cDNA encoding diacylglycerol acyltransferase from Arabidopsis thaliana and its functional expression. FEBS. 1999, 452: 145-149. 10.1016/S0014-5793(99)00646-8.View ArticleGoogle Scholar
- Routaboul JM, Benning C, Bechtold N, Caboche M, Lepiniec L: The TAG1 locus of Arabidopsis encodes for a diacylglycerol acyltransferase. Plant Physiol and Biochem. 1999, 37: 831-840. 10.1016/S0981-9428(99)00115-1.View ArticleGoogle Scholar
- Jako C, Kumar A, Wei Y, Zou J, Barton DL, Giblin EM, Covello PS, Taylor DC: Seedspecific over-expression of an Arabidopsis cDNA encoding a diacylglycerol acyltransferase enhances seed oil content and seed weight. Plant Physiol. 2001, 126: 861-874. 10.1104/pp.126.2.861.PubMedPubMed CentralView ArticleGoogle Scholar
- Weselake RJ, Shah S, Tang M, Quant PA, Snyder CL, Furukawa-Stoffer TL, Zhu W, Taylor DC, Zou J, Kumar A, Hall L, Laroche A, Rakow G, Raney P, Moloney MM, Harwood JL: Metabolic control analysis is helpful for informed genetic manipulation of oilseed rape (Brassica napus) to increase seed oil content. J Exp Bot. 2008, 59: 3543-3549. 10.1093/jxb/ern206.PubMedPubMed CentralView ArticleGoogle Scholar
- Taylor DC, Zhang Y, Kumar A, Francis T, Giblin EM, Barton DL, Ferrie JR, Laroche A, Shah S, Zhu W, Snyder CL, Hall L, Rakow G, Harwood JL, Weselake RJ: Molecular modification of triacylglycerol accumulation under field conditions to produce canola with increased seed oil content. Botany. 2009, 87: 533-543. 10.1139/B08-101.View ArticleGoogle Scholar
- Lardizabal KD, Mai JT, Wagner NW, Wyrick A, Voelker T, Hawkins DJ: DGAT2 is a new diacylglycerol acyltransferase gene family: purification, cloning, and expression in insect cells of two polypeptides from Mortierella ramanniana with diacylglycerol acyltransferase activity. J Biol Chem. 2001, 276: 38862-38869. 10.1074/jbc.M106168200.PubMedView ArticleGoogle Scholar
- Cases S, Stone SJ, Zhou P, Yen E, Tom B, Lardizabal KD, Farese RV: Cloning of DGAT2, a second mammalian diacylglycerol acyltransferase, and related family members. J Biol Chem. 2001, 276: 38870-38876. 10.1074/jbc.M106219200.PubMedView ArticleGoogle Scholar
- Shockey JM, Gidda SK, Chapital DC, Kuan JC, Dhanoa PK, Bland JM, Rothstein SJ, Mullen RT, Dyer JM: Tung Tree DGAT1 and DGAT2 have nonredundant functions in triacylglycerol biosynthesis and are localized to different subdomains of the endoplasmic reticulum. Plant Cell. 2006, 15: 1872-1887.Google Scholar
- Kroon JT, Wei W, Simon WJ, Slabas AR: Identification and functional expression of a type 2 acyl-CoA: diacylglycerol acyltransferase (DGAT2) in developing castor bean seeds which has high homology to the major triglyceride biosynthetic enzyme of fungi and animals. Phytochemistry. 2006, 67: 2541-2549. 10.1016/j.phytochem.2006.09.020.PubMedView 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-3901. 10.1105/tpc.109.071795.PubMedPubMed CentralView ArticleGoogle Scholar
- Dahlqvist A, Stahl U, Lenman M, Banas A, Lee M, Sandager L, Ronne H, Stymne S: 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-6492. 10.1073/pnas.120067297.PubMedPubMed CentralView ArticleGoogle Scholar
- Banas A, Dahlqvist A, Stahl U, Lenman M, Stymne S: The involvement of phospholipid: diacylglycerol acyltransferases in triacylglycerol production. Biochem Soc Trans. 2000, 28: 703-705. 10.1042/BST0280703.PubMedView ArticleGoogle Scholar
- Stahl U, Carlsson AS, Lenman M, Dahlqvist A, Huang B, Banas W, Banas A, Stymne S: Cloning and functional characterization of a phospholipid: diacylglycerol acyltransferase from Arabidopsis. Plant Physiol. 2004, 135: 1324-1335. 10.1104/pp.104.044354.PubMedPubMed CentralView ArticleGoogle Scholar
- Mhaske V, Beldjilali K, Ohlrogge J, Pollard M: Isolation and characterization of an Arabidopsis thaliana knockout line for phospholipid: diacylglycerol transacylase gene (At5g13640). Plant Physiol Biochem. 2005, 43: 413-417. 10.1016/j.plaphy.2005.01.013.PubMedView ArticleGoogle Scholar
- Lu C, de Noyer SB, Hobbs DH, Kang J, Wen Y, Krachtus D, Hills MJ: Expression pattern of diacylglycerol acyltransferase-1, an enzyme involved in triacylglycerol biosynthesis, in Arabidopsis thaliana. Plant Mol Biol. 2003, 52: 31-41. 10.1023/A:1023935605864.PubMedView ArticleGoogle Scholar
- Kim HU, Lee K-R, Go YS, Jung JH, Suh M-C, Kim JB: Endoplasmic reticulum-located PDAT1-2 from castor bean enhances hydroxy fatty acid accumulation in transgenic plants. Plant Cell Physiol. 2011, 52: 983-993. 10.1093/pcp/pcr051.PubMedView ArticleGoogle Scholar
- Van Erp H: Castor phospholipid: diacylglycerol acyltransferase facilitates efficient metabolism of hydroxy fatty acids in transgenic Arabidopsis. Plant Physiol. 2011, 155: 683-693. 10.1104/pp.110.167239.PubMedPubMed CentralView ArticleGoogle Scholar
- Stahl U, Stalberg K, Stymne S, Ronne H: A family of eukaryotic lysophospholipid acyltransferases with broad specificity. FEBS Lett. 2008, 582: 305-309. 10.1016/j.febslet.2007.12.020.PubMedView ArticleGoogle Scholar
- Taylor DC, Weber N, Hogge LR, Underhill EW, Pomeroy MK: Formation of Trierucoylglycerol (Trierucin) from 1,2-Dierucoylglycerol by a homogenate of microspore-derived embryos of Brassica napus L. J Am Oil Chem Soc. 1992, 69: 355-358. 10.1007/BF02636067.View ArticleGoogle Scholar
- Sandager L, Gustavsson MH, Stahl U, Dahlqvist A, Wiberg E, Banas A, Lenman M, Ronne H, Stymne S: Storage lipid synthesis is non-essential in yeast. J Biol Chem. 2002, 277: 6478-6482. 10.1074/jbc.M109109200.PubMedView ArticleGoogle Scholar
- Xu J, Francis T, Mietkiewska E, Giblin EM, Barton DL, Zhang Y, Zhang M, Taylor DC: Cloning and characterization of an acyl-CoA-dependent diacylglycerol acyltransferase 1 (DGAT1) gene from Tropaeolum majus, and a study of the functional motifs of the DGAT protein using site-directed mutagenesis to modify enzyme activity and oil content. Plant Biotechnol J. 2008, 6: 799-818. 10.1111/j.1467-7652.2008.00358.x.PubMedView 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. The Plant Cell. 2009, 21: 3885-3901. 10.1105/tpc.109.071795.PubMedPubMed CentralView ArticleGoogle Scholar
- Mañas-Fernández A, Vilches-Ferrón M, Garrido-Cárdenas J, Belarbi E-H, Alonso D, García-Maroto F: Cloning and molecular characterization of the acyl-CoA: diacylglycerol acyltransferase 1 (DGAT1) Gene from Echium. Lipids. 2009, 44: 555-568. 10.1007/s11745-009-3303-9.PubMedView ArticleGoogle Scholar
- Guihéneuf F, Leu S, Zarka A, Khozin-Goldberg I, Khalilov I, Boussibab S: Cloning and molecular characterization of a novel acyl-CoA: diacylglycerol acyltransferase 1-like gene (PtDGAT1) from the diatom Phaeodactylum tricornutum. FEBS J. 2011, 278: 3651-3666. 10.1111/j.1742-4658.2011.08284.x.PubMedView ArticleGoogle Scholar
- Turchetto-Zolet AC, Maraschin FS, de Morais GL, Cagliari A, Andrade CMB, Margis-Pinheiro M, Margis R: Evolutionary view of acyl-CoA diacylglycerol acyltransferase (DGAT), a key enzyme in neutral lipid biosynthesis. BMC Evolutionary Biology. 2011, 11: 263-10.1186/1471-2148-11-263. doi:10.1186/1471-2148-11-263PubMedPubMed CentralView ArticleGoogle Scholar
- Beisson F, Koo AJK, Ruuska S, Schwender J, Pollard M, Thelen JJ, Paddock T, Salas JJ, Savage L, Milcamps A, Mhaske VB, Cho Y, Ohlrogge JB: Arabidopsis genes involved in acyl lipid metabolism a 2003 census of the candidates, a study of the distribution of expressed sequence tags in organs, and a web-based database. Plant Physiol. 2003, 132: 681-697. 10.1104/pp.103.022988.PubMedPubMed CentralView ArticleGoogle Scholar
- Shen W, Chen Q, Wang L, Zou J: Functional Study of Lysophosphatidylcholine Acyltransferase in Arabidopsis. Proceedings of the 18th International Symposium on Plant Lipids. Bordeaux; 2008:102.Google Scholar
- Lu C, Xin Z, Rin Z, Miquel M, Browse J: An enzyme regulating triacylglycerol composition is encoded by the ROD1 gene of Arabidopsis. Proc Nat'l Acad Sci USA. 2009, 106: 18837-18842. 10.1073/pnas.0908848106.View ArticleGoogle Scholar
- Taylor DC, Barton DL, Rioux KP, MacKenzie SL, Reed DW, Underhill EW, Pomeroy MK, Weber N: Biosynthesis of acyl lipids containing very long chain fatty acids in microspore-derived and zygotic embryos of Brassica napus L cv Reston. Plant Physiol. 1992, 99: 1609-1618. 10.1104/pp.99.4.1609.PubMedPubMed CentralView ArticleGoogle Scholar
- Bradford MM: A rapid and sensitive method for the quantification of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976, 72: 248-254. 10.1016/0003-2697(76)90527-3.PubMedView ArticleGoogle Scholar
- Yokoyama K, Shimiza F, Setaka M: Simultaneous separation of lysophospholipids from the total lipid fraction of crude biological samples using two-dimensional thinlayer chromatography. J Lip Res. 2000, 41: 142-147.Google Scholar
- Wang CS, Vodkin LO: Extraction of RNA from tissues containing high levels of procyanidins that bind RNA. Plant Mol Bio Rep. 1994, 12: 132-145. 10.1007/BF02668374.View ArticleGoogle Scholar
- Meadus WJ: A semi-quantitative RT-PCR method to measure the in vivo effect of dietary conjugated linoleic acid on porcine muscle PPAR gene expression. Biol Prod Online. 2003, 5: 20-28. 10.1251/bpo43.View ArticleGoogle Scholar
- Alonso JM, Stepanova AN, Leisse TJ, Kim CJ, Chen H, Shinn P, Stevenson DK, Zimmerman J, Barajas P, Cheuk R, Gadrinab C, Heller C, Jeske A, Koesema E, Meyers CC, Parker H, Prednis L, Ansari Y, Choy N, Deen H, Geralt M, Hazari N, Hom E, Karnes M, Mulholland C, Ndubaku R, Schmidt I, Guzman P, Aguilar-Henonin L, Schmid M, Weigel D, Carter DE, Marchand T, Risseeuw E, Brogden D, Zeko A, Crosby WL, Berry CC, Ecker JR: Genome-wide insertional mutagenesis of Arabidopsis thaliana. Science. 2003, 301: 653-657. 10.1126/science.1086391.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution ArticleCopyright: License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.