# Overexpression of MYB115, AAD2, or AAD3 in Arabidopsis thaliana seeds yields contrasting omega-7 contents

Sébastien Baud, Hasna Ettaki, Manuel Adrián Troncoso-Ponce, Alexandra A.

To, Guillaume Barthole, Loic Lepiniec, Sebastien S. Baud

► To cite this version:

Sébastien Baud, Hasna Ettaki, Manuel Adrián Troncoso-Ponce, Alexandra A. To, Guillaume Barthole, et al.. Overexpression of MYB115, AAD2, or AAD3 in Arabidopsis thaliana seeds yields contrasting omega-7 contents. PLoS ONE, 2018, 13 (1), pp.e0192156. 10.1371/journal.pone.0192156. hal-02356207

HAL Id: hal-02356207

https://hal.science/hal-02356207

Submitted on 27 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Overexpression of *MYB115*, *AAD2*, or *AAD3* in *Arabidopsis thaliana* seeds yields contrasting omega-7 contents

Hasna Ettaki¹*, Manuel Adrián Troncoso-Ponce²*, Alexandra To¹, Guillaume Barthole¹,³, Loïc Lepiniec¹, Sébastien Baud¹*

1 Institut Jean-Pierre Bourgin, Institut National de la Recherche Agronomique, AgroParisTech, Centre National de la Recherche Scientifique, Université Paris-Saclay, Versailles, France, 2 Sorbonne Universités, Génie Enzymatique et Cellulaire, Formation de Recherche en Evolution Centre National de la Recherche Scientifique 3580, Université de Technologie de Compiègne, Compiègne, France, 3 Université Paris-Sud, Université Paris-Saclay, Orsay, France

* These authors contributed equally to this work.

* sebastien.baud@inra.fr

<figure><img src="image_2.png" /></figure>

## Abstract

Omega-7 monoenoic fatty acids (ω-7 FAs) are increasingly exploited both for their positive effects on health and for their industrial potential. Some plant species produce fruits or seeds with high amounts of ω-7 FAs. However, the low yields and poor agronomic properties of these plants preclude their commercial use. As an alternative, the metabolic engineering of oilseed crops for sustainable ω-7 FA production has been proposed. Two palmitoyl-ACP desaturases (PADs) catalyzing ω-7 FA biosynthesis were recently identified and characterized in *Arabidopsis thaliana*, together with MYB115 and MYB118, two transcription factors that positively control the expression of the corresponding PAD genes. In the present research, we examine the biotechnological potential of these new actors of ω-7 metabolism for the metabolic engineering of plant-based production of ω-7 FAs. We placed the PAD and MYB115 coding sequences under the control of a promoter strongly induced in seeds and evaluated these different constructs in *A. thaliana*. Seeds were obtained that exhibit ω-7 FA contents ranging from 10 to >50% of the total FAs, and these major compositional changes have no detrimental effect on seed germination.

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Ettaki H, Troncoso-Ponce MA, To A, Barthole G, Lepiniec L, Baud S (2018) Overexpression of *MYB115*, *AAD2*, or *AAD3* in *Arabidopsis thaliana* seeds yields contrasting omega-7 contents. PLoS ONE 13(1): e0192156.
https://doi.org/10.1371/journal.pone.0192156

**Editor:** Lam-Son Phan Tran, RIKEN Center for Sustainable Resource Science, JAPAN

**Received:** November 14, 2017

**Accepted:** January 17, 2018

**Published:** January 30, 2018

**Copyright:** © 2018 Ettaki et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** The IJPB benefits from the support of the Labex Saclay Plant Sciences-SPS (ANR-10-LABX-0040-SPS). The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

## Introduction

Vegetable oils enriched in omega-7 (ω-7) monounsaturated fatty acids (FAs) such as palmitoleic acid (cis-ω-7 C16:1) and its elongation product vaccenic acid (cis-ω-7 C18:1) potentially have uses for a number of applications. First, they are attractive for biodiesel formulations since biodiesel containing FA methyl esters produced from oils rich in ω-7s has superior functional properties [1]. In addition, ω-7s from plant oils can offer the chemist promising starting material for catalytic conversions to renewable platform chemicals. Olefin metathesis constitutes a powerful tool for polymer chemistry, and ethenolytic metathesis of ω-7 FAs from plant
oils could potentially provide a competitive source of 1-octene to make linear low-density polyethylene [2]. Finally, after the fast growth of ω-3 polyunsaturated FAs in the nutrition market [3], suppliers and marketers are today calling attention to ω-7 FAs since vegetable oils enriched in ω-7s have been ascribed a number of beneficial health properties. Palmitoleic acid in particular is considered a lipokine used by adipose tissues to communicate with distant organs and regulate systemic metabolic homeostasis [4]. Increasing evidence suggests that palmitoleic acid plays a key role in the physiopathology of insulin resistance in humans, increasing muscle response to insulin [5]. Aside from these nutritional functions, palmitoleic acid is attractive for nonfood uses in the skin care industry because of its antioxidant and antimicrobial properties [6].

Natural plant oils containing high levels of ω-7 FAs are infrequent. Foods and nutraceuticals enriched in ω-7s for health purpose are often sourced from sea buckthorn berries (*Hippophae rhamnoides*) that contain up to 50% ω-7s in pulp oil [7]. Seed oils enriched in ω-7 FAs have been described in milkweed (*Asclepias syriaca*; ~25% ω-7s), macadamia (*Macadamia integrifolia*; ~35% ω-7s), cat's claw vine (*Doxantha unguis-cati*; ~75% ω-7s), *Roureopsis obliquifoliata* (~30% ω-7s), and *Entandrophragma cylindricum* (~55% ω-7s), among others (phylofadb.bch.msu.edu). These plants exhibit low yields and poor agronomic properties. Oilseed crops contain very low amounts of ω-7 monounsaturated FAs. To a certain extent, rapeseed (*Brassica napus*) is an exception to this rule. In seeds of *B. napus*, as in seeds of its close relative *Arabidopsis thaliana*, the proportion of ω-7s accounts only for a small percentage of total seed FAs [8], but ω-7 FAs are highly concentrated in the endosperm, a seed compartment of reduced size that comprises a single cell layer surrounding the embryo in mature dry seeds [9,10]. Vaccenic acid and paullinic acid (*cis*-ω-7 C20:1) thus represent more than 20 mol% of total FAs in the endosperm fraction of *A. thaliana* seeds, and approximately 35 mol% in rapeseed [11].

Production of saturated FAs is performed in a stepwise manner, in plastids, by the FA synthase of type II. Saturated acyl chains, bound to acyl-carrier proteins (ACPs), can then be desaturated by stromal soluble acyl-ACP desaturases (AADs) to form *cis*-monoenes. Depending on both their substrate specificity and their regiospecificity, activities of AAD isoforms can give rise to a range of *cis*-monoenes [12]. $\Delta^9$ Stearoyl-ACP desaturases (SADs) efficiently desaturate C18:0 (stearic acid) to form *cis*-ω-9 C18:1 ($\Delta^9$ 18:1; oleic acid) and represent the archetypal AAD [13]. Some SADs also exhibit low levels of $\Delta^9$ palmitoyl-ACP desaturase (PAD) activity, thus forming *cis*-ω-7 C16:1 ($\Delta^9$ 16:1; palmitoleic acid) [14,15]. However, production of ω-7 monoenes in tissues storing oils enriched in ω-7 FAs relies on variant specialized PAD isoforms that display a distinct substrate specificity, with a preference for C16:0-ACP. Whereas SAD isoforms have been cloned and characterized in a relatively wide range of plant species, a limited number of PADs have been identified so far. A PAD-encoding cDNA was first cloned in cat's claw vine [16]. More recently, two PAD isoforms named AAD2 and AAD3 were identified in *A. thaliana* [11,17]. Three-dimensional structures of these related, but functionally divergent, soluble desaturases have provided clues on how they recognize the chain length of substrates [18,19]. A hydrophobic channel within the enzymes accommodates FA substrates bound to ACP, and the side chains of amino acids lining the lower portion of this channel set constraints on the chain length of FA substrates. The structure of the channel of the archetypelike $\Delta^9$ SAD is highly conserved among plant species. This channel is deep enough to accommodate C18:0-ACP substrates. On the contrary, the bottom part of the substrate channel of $\Delta^9$ PADs is lined by hydrophobic residues with bulky lateral chains that reduce the depth of their substrate pockets and favor the binding of shorter C16:0-ACP substrates [17,19]. Despite similar enzymatic activities, $\Delta^9$ PADs from cat's claw vine and *A. thaliana* do not have identical
substrate channels, suggesting that these activities may have appeared independently during
the evolution of AADs.

Due to the high economic value of plant ω-7s, there is an increasing interest in either find-
ing new natural plant sources with potentially better agronomical performances or improving
existing oilseed crops for better ω-7 production. Bryant et al. [11] determined the ω-7 content
of endosperm oil in 10 inbred rapeseed varieties and reported significant variation, with the
relative proportion of ω-7s reaching 47% in the endosperm of the Swede (Rutabaga) variety
‘Huguenot’. This enrichment in ω-7s concerns only the aleurone, where approximately 5% of
the FA in rapeseed is stored. The authors suggested that extraction of ω-7s from such a variety
could be economically viable provided that rapeseed is dehulled before oil extraction. Alterna-
tively, approaches of metabolic engineering can be implemented for enhancing production of
ω-7s in the main oil-storing compartment of oilseeds. This has been successfully achieved in
*A. thaliana* and in crops such as camelina (Camelina sativa), soybean (Glycine max), and rape-
seed [1,20,21]. These approaches all relied on a strong and seed-specific expression of PAD
sequences either cloned from species naturally accumulating high ω-7 levels (cat's claw vine)
or obtained through the reengineering of archetypal SADs to achieve desired substrate speci-
ficity while retaining the stability and turnover characteristics of a paralog [22]. Additional
metabolic modifications were then introduced in some of these systems to efficiently redirect
fluxes for high levels of ω-7 accumulation in the host seed (reviewed in [23]). Examples include
reducing the elongation and the export of C16:0 from the plastid to enhance substrate avail-
ability for the PADs. Ultimately, iterative optimization of their strategy led Nguyen et al. [1] to
use an assembly of six transgenes yielding 60–65% ω-7 FAs in the oil of transgenic camelina
seeds.

The evolutionary history of the AADs has led to the emergence of isoforms catalyzing the
biosynthesis of ‘exotic’ monoenes such as the ω-7s, therefore allowing the specialization of oil
metabolism in certain oleaginous species. Interestingly, recent studies have underlined that
specialization of ω-7 metabolism also relied on the complex setup of dedicated transcriptional
machinery able to precisely control the spatiotemporal expression of PAD genes. In *A. thali-
ana*, two closely related members of the MYB family of transcription factors (TFs), *MYB115*
and *MYB118*, are transcriptionally induced at the onset of the maturation phase in the endo-
sperm. In this compartment, they activate a set of common transcriptional targets such as the
two PAD-coding genes *AAD2* and *AAD3*. Even though the molecular mechanism underlying
this transcriptional activation has not been fully elucidated, the two MYBs were shown to be
necessary for endosperm-specific activation of *AAD2* and *AAD3* and the subsequent accumu-
lation of ω-7s in this seed compartment. Accordingly, the relative proportion of ω-7s is drasti-
cally reduced in the endosperm oil of *myb115 myb118* double mutants, just like in that of the
*aad2 aad3* double mutants [17].

In the present report, we systematically examine the biotechnological potential of the actors
of ω-7 metabolism newly isolated in *A. thaliana* for increasing ω-7 FA content in seeds. Our
objective is to identify new gene sequences useful for metabolic engineering for efficient plant-
based production of ω-7 FAs. We therefore placed coding sequences of *A. thaliana* PADs or
their transcriptional activator MYB115 under the control of a promoter strongly induced in
seeds and evaluated these different constructs in *A. thaliana*. Oils were obtained that exhibited
ω-7 FA contents ranging from 10 to >50% of the total FAs. Even though all the constructs
tested efficiently yielded expected compositional changes, direct overexpression of the *PADs*
appeared more efficient than that of *MYB115* for stimulating ω-7 FA biosynthesis. Remark-
ably, these major compositional changes had no detrimental defect on seed germination.
# Materials and methods

## Plant material and growth conditions

Plants were cultured as described in Baud et al. [24]. To sample embryos and endosperm/seed coats, seeds were imbibed and dissected under an optical glass binocular magnifier with a scalpel and dissecting tweezers. Samples used for RNA extraction were frozen in liquid nitrogen immediately after dissection and then stored at -80°C. Weight determinations of seed samples were realized on a M2P balance (Sartorius, Aubagne, France). For germination assays, seeds were sown in triplicate in Petri dishes containing 0.5% (w/v) solidified agarose. After stratification, plates were kept in a growth cabinet (with continuous light at 25°C), and germination was scored based on radicle emergence.

## Constructs and plant transformation

The sequences of primers used for DNA amplification are indicated in S1 Table.

For construction of the *ProAT2S2:uidA*, *ProAT2S2:MYB115*, *ProAT2S2:AAD2*, and *ProAT2S2:AAD3* transgenes, cDNA was amplified with the proofreading Pfu Ultra DNA polymerase (Stratagene, Les Ulis, France) from the pGWB3 vector ([25]; *uidA*) or from a mixture of seed cDNA (Col-0 accession; *MYB115*, *AAD2*, and *AAD3*). The PCR products thus obtained were introduced by BP recombination into the pDONR207 entry vector (Invitrogen, Paris, France) and transferred into the ProAT2S2-R1R2-HYGRO destination vector [26] by LR recombination. The resulting binary vectors were electroporated into the *Agrobacterium tumefaciens* C58C¹ strain and used for agroinfiltration of flower buds of *A. thaliana* [27]. Primary transformants were selected on Murashige and Skoog medium containing hygromycin (50 mg.l⁻¹). They were then transferred to soil for further characterization. The progeny of these primary transformants (T2 seeds) was subjected to segregation analyses, and lines segregating 3:1 for hygromycin resistance were selected (heterozygous lines, one insertion locus). T2 lines were then grown in a greenhouse, and their progeny (T3 seeds) was subjected to segregation analyses. Lines producing 100%-resistant plantlets were selected (homozygous lines, single insertion locus) and used for further analyses.

## RNA analyses

RNA extraction, reverse transcription, and real-time RT quantitative PCR were performed as previously described [28]. The sequences of primers used for real-time RT-qPCR are indicated in S2 Table. Purity of the different seed fractions sampled was assessed as described in Barthole et al. [29]. Briefly, marker genes for each of the seed fractions sampled, namely *ZHOUPI* (endosperm-specific) and *At2g23230* (embryo-specific), were quantified on cDNA prepared from these fractions, thus demonstrating that no significant contamination occurred between fractions.

## Fatty acid analyses

Gas chromatography analyses of total FAs were performed as previously described [10] on pools of *A. thaliana* seeds or seed fractions. The endosperm was analyzed with the seed coat attached as described in Penfield et al. [9]. The procedure did not bias our evaluation of the endosperm oil content since the seed coat does not accumulate storage lipids and undergo programmed cell death during seed maturation [10,30].

## Microscopy

Bright-field microscope observations and histochemical detection of GUS activity were carried out as previously described in Baud et al. [24].
## Accession numbers

Sequence data can be found in the GenBank/EMBL data libraries under the following accession numbers: AAD2, At4g02610; AAD3, At5g16230; AT2S2, At4g27140; MYB115, At5g40360; MYB118, At3g27785.

## Results

### Validation of a promoter sequence suitable for overexpression of genes of interest in zygotic tissues of the seed

Because of the importance of the transcriptional activation of PAD-coding genes by MYB115 and MYB118 for ω-7 production in *A. thaliana* seeds [17], these MYBs appeared as potential new targets for the engineering of seed oils enriched in ω-7 monoenes. Before evaluating the biotechnological potential of these TFs, we looked for a promoter sequence allowing their overexpression in the oil-storing compartments of seeds. The use of a 35S promoter sequence for high and ubiquitous expression was not possible since the *Pro35S dual*:MYB115 and *Pro35S dual*:MYB118 cassettes dramatically affect plant growth and reproduction, yielding dwarfism and partial sterility [17,29]. To overcome the negative impact of a strong expression of *MYB115* or *MYB118* in non-seed tissues, we instead chose the promoter sequence of the *AT2S2* gene (At4g27150) coding for a seed storage protein. RNA gel blot analyses, *in situ* hybridization approaches, and laser-capture microdissection of maturing seeds followed by mRNA quantification using stringent analyses of Affymetrix ATH1 GeneChip hybridization data previously showed that *AT2S2* is highly induced in the endosperm and in the embryo of maturing seeds [31,32]. We verified this expression pattern by RT-qPCR on a set of cDNA prepared from a range of plant organs of the wild-type accession Columbia-0 (Col-0). To begin with, a range of plant organs were considered, and *AT2S2* appeared to be highly expressed in siliques (Fig 1A). To further characterize the expression pattern of *AT2S2*, a time-course analysis of *AT2S2* mRNA abundance was carried out in developing seeds excised from the siliques, which revealed a peak of transcript accumulation during maturation (Fig 1B). To gain further insight into the tissue specificity of *AT2S2* expression, maturing seeds were dissected. The two fractions thus obtained, embryo and endosperm/seed coat, were independently analyzed. A strong induction of *AT2S2* during the course of maturation was observed in the two seed fractions under study (Fig 1C). These data confirmed that the *AT2S2* promoter was well adapted for overexpressing genes of interest in maturing seeds.

A 984-bp *AT2S2* gene promoter fragment exhibiting a close association of RY and G-box elements essential for the activation of maturation genes by the master regulators of seed maturation [33] was then cloned (S3 Table) and transcriptionally fused to the *uidA* reporter gene for validating its activity. The corresponding construct was assayed for the resulting *uidA* expression pattern in transgenic *A. thaliana* lines. Intense GUS activity was observed in the embryos (Fig 1D–1G) and in the endosperm of maturing seeds (Fig 1H–1J). This promoter fragment was therefore used for driving the expression of *MYB115* and *MYB118* in seeds.

### Overexpression of *MYB115* in seeds

Both *ProAT2S2:MYB115* and *ProAT2S2:MYB118* cassettes were prepared (see Methods). Unfortunately, we failed to recover *ProAT2S2:MYB118* transformants. The corresponding antibiotic-resistant plantlets selected *in vitro* did not survive when transferred to soil. On the contrary, transgenic lines carrying the *ProAT2S2:MYB115* construct could be generated. Five independent transformants with a single insertion locus were characterized. Homozygous lines were grown under controlled conditions, and neither vegetative growth nor fertility was
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

**Fig 1. Validation of the AT2S2 promoter sequence.** (A-C) Analysis of relative mRNA accumulation of AT2S2 was performed in different plant organs (A), in developing seeds (B), and in developmental series of endosperm (Endo.) and embryo fractions (C). The results obtained were standardized to the *EF1αA4 (EF)* gene expression level. Values are the means and SE of three to six replicates carried out on cDNA dilutions obtained from three independent mRNA extractions. DAA, days after anthesis; Fl, flowers; RI, rosette leaves; Ro, roots; St, stems. (D-J) Pattern of activity of the
*ProAT2S2:uidA* cassette in maturing embryos harvested 10 (D), 12 (E), 14 (F), or 16 days after anthesis (DAA) (G) and in endosperm fractions harvested 14 (H) or 16 DAA (I). A close-up of a peeled endosperm layer aged 14 DAA is presented in (J). For histochemical detection of GUS activity, tissues were incubated 4 hours in a buffer containing 2 mM each of potassium ferrocyanide and potassium ferricyanide. Microscopy observations were performed using Nomarski optics. Bars = 100 µm in (D-I), 10 µm in (J).

https://doi.org/10.1371/journal.pone.0192156.g001

affected in these lines, despite moderate activity of the *ProAT2S2* promoter used in vegetative and in reproductive organs (S1 Fig). Likewise, whole-mount clearing of developing seeds showed that the structure and early development of the three tissues composing the seeds were unmodified in a *ProAT2S2:MYB115* background (S2 Fig).

Using a RT-qPCR strategy, *MYB115* mRNA level was quantified in maturing seeds of the transgenic lines (Fig 2A). Expression levels ranging from 14% to 25% of the *EF1αA4* (At5g60390) mRNA level (constitutive expression; [34]) were detected 14 days after anthesis (DAA), indicating that *MYB115* was efficiently overexpressed in these lines (from 660- to 1,170-fold compared to expression in the wild type).

To evaluate the effect of *MYB115* overexpression on seed filling, whole mature dry seeds were first analyzed. Seed dry weight and total FA content were significantly increased in the transgenic lines (from 10% to 20%) as a consequence of increased seed dimensions (Fig 3A-3E). However, FA concentration in these seeds remained unchanged. Mature seeds were then dissected. The two fractions obtained, embryo and endosperm/seed coat, were collected separately before total FA quantification by gas chromatography. A significant increase in total FAs was measured in transgenic embryos (Fig 3F). On the contrary, the total amount of FAs stored in the endosperm fraction was unmodified (Fig 3G). The FA composition of the oil stored in transgenic seeds was then examined, with a particular focus placed on monoenes of the ω-7 series. At the whole seed level, overexpression of *MYB115* yielded a strong enrichment in ω-7 FAs (*cis*-ω-7 C16:1, *cis*-ω-7 C18:1, and *cis*-ω-7 C20:1), the relative proportion of which was more than doubled with respect to the wild type (Fig 2B; S4 Table). This phenotype hid contrasting results in the two zygotic compartments of the seed, with the concentration of ω-7 being unchanged in endosperm oil, whereas that of embryo oil was more than tripled (Fig 2C and 2D; S5 and S6 Tables).

To test whether the increase in ω-7 content measured in transgenic seeds correlated with higher *AAD2* and *AAD3* transcript abundance, we first quantified the relative accumulation of these transcripts on cDNA prepared from whole seeds harvested 14 DAA. There was no significant variation of the *AAD* transcript levels in transgenic seeds with respect to the wild type (S3 Fig). However, a tendency for increased *AAD3* mRNA accumulation in transgenic seeds was observed. FA analyses carried out on dissected seeds (see above) led us to hypothesize that statistically significant variation of the *AAD* transcript levels occurring in transgenic embryos might be hidden by the naturally elevated but unmodified abundance of these transcripts in corresponding albumens. For a subset of lines, we therefore prepared cDNA from dissected embryos and verified that increased embryo ω-7 levels correlated with enhanced transcriptional activation of both *AAD2* and *AAD3* in this seed compartment (Fig 4). Unfortunately, the very small size of the endosperm fraction prevented us from collecting sufficient material on the different lines considered to perform similar analyses on this fraction. As a consequence, we ignore whether the *AADs* were efficiently overexpressed in the albumen of transgenic seeds.

## Overexpression of *AAD2* or *AAD3* in seeds

To compare the effect of *MYB115* overexpression on ω-7 accumulation with that of one or the other of the *PADs* induced by *MYB115*, transgenic plants carrying either the *ProAT2S2:AAD2*
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

**Fig 2.** Increased omega-7 fatty acid content of *ProAT2S2:MYB115* seeds. The five independent transformants considered, TCR4, TAR3, TXR2, TGR5, and TJR2 (in order from left to right), are displayed in the same order in the four graphs in the figure. (A) RT-qPCR analysis of transcript abundance in cDNA prepared from maturing seeds aged 14 DAA to assess efficient overexpression of the transgene. Values are the means and SE of three to 12 replicates performed on cDNA dilutions obtained from three independent mRNA extractions. (B-D) Relative proportion of ω-7
fatty acids (*cis*-ω-7 C16:1, *cis*-ω-7 C18:1, and *cis*-ω-7 C20:1) in mature dry seeds (B), in embryos (C), and in endosperm fractions (D) dissected from mature seeds. Values are the means and SE of five replicates performed on batches of 20 seeds from five plants. Asterisks indicate significant differences from the wild type according to *t*-test at *** P<0.001 and ** P<0.01, respectively.

https://doi.org/10.1371/journal.pone.0192156.g002

construct or the *ProAT2S2:AAD3* construct were then obtained. For each construct, five independent transformants were considered. With a RT-qPCR strategy, *AAD* mRNA levels were quantified in maturing seeds 14 DAA (**Fig 5A**). Significantly increased *AAD2* mRNA levels were detected in maturing *ProAT2S2:AAD2* seeds, indicating that the transgene was efficiently overexpressed (from 20- to 60-fold compared to the wild type). Similarly, a significant overexpression of *AAD3* (from 10- to 14-fold compared to the wild type) was measured in the *ProAT2S2:AAD3* lines. To evaluate the effect of *AAD2* and *AAD3* overexpression on seed filling, whole mature dry seeds were then subjected to fine biochemical analysis. Dry weight of the transgenic mature seeds was unmodified (**S4A Fig**). Likewise, the overall oil content of the two zygotic compartments of the seeds was unchanged. The FA composition of the oil stored in transgenic seeds was then examined. At the seed level, overexpression of *AAD2* and *AAD3* yielded a strong enrichment in ω-7s, the relative proportion of which was multiplied by eight (*AAD3*) to 10 (*AAD2*) with respect to the wild type (**Fig 5B; S4 Table**). The ω-7 FA contents measured at the seed level faithfully reflected the relative proportions of ω-7 FA species observed both in the embryo and in the endosperm fractions of transgenic seeds (**Fig 5C and 5D; S5 and S6 Tables**).

**Increased seed omega-7 content does not affect germination**

Since drastic alterations of FA composition of seed oils sometimes affect triacylglycerol lipolysis during germination and seed vigor [35], germination of the engineered seeds was then examined (**Fig 6**). Under standard light conditions, *ProAT2S2:MYB115* seeds exhibited a transient but reproducible increase of the germination rate 24 hours after the first exposure to the light. This phenotype was apparently not linked with the increased ω-7 content of these seeds, given that the germination rate of seeds overexpressing *AAD2* or *AAD3* was similar to that of wild-type seeds. We then indirectly evaluated the ability of transgenic seedlings to remobilize triacylglycerols enriched in ω-7 FAs. For this purpose, cold-stratified seeds were submitted to a brief light impulse aimed at triggering the germination process and then placed in the dark, where carbon fixation by photosynthesis is prevented. After one week, the germination rate of engineered seeds was similar to that of wild-type seeds. The hypocotyl length of the transgenic etiolated seedlings thus obtained was unmodified, suggesting that triacylglycerols enriched in ω-7 FAs were efficiently remobilized.

**Discussion**

Iterative optimization of unusual FA production in seeds not only requires the stacking of multiple traits [36], but also implies the evaluation and comparison of different gene versions for a given trait and that of regulatory elements for driving the expression of selected transgenes. The results presented in this study demonstrate that overexpression of either of the two *A. thaliana* PADs in seeds is sufficient to partially shift monoene production from ω-9 to ω-7 FAs. This result can be attained both by overexpression of the desaturase-coding sequences placed under the control of the *AT2S2* promoter and by overexpression of one of their transcriptional activators, such as *MYB115*. With regard to the efficiency of the different strategies evaluated, the latter appears to have lower efficiency. The relative ω-7 FA content of *ProAT2S2:MYB115* seeds is five time less than that of *ProAT2S2:AAD2* seeds on average. Further,
<figure><img src="image_13.png" /><figcaption>**Fig 3. Characterization of transgenic seeds overexpressing MYB115.** The five independent transformants considered, TCR4, TAR3, TXR2, TGR5, and TJR2 (in order from left to right), are displayed in the same order in every graph in the figure. (A) Mature seed dry weight. Values are the means and SE of five replicates carried out on batches of 20 seeds from five plants. (B) Mature seed length. Values are the means and SE of 100 measurements carried out on seeds from five plants. (C) Mature seed width. Values are the means and SE of 100 measurements carried out on seeds from five plants. (D, E) Total fatty acid content of mature dry seeds, expressed in µg·mg⁻¹ DW (D) or in µg per seed (E). Values are the means and SE of five replicates carried out on batches of 20 seeds from five plants. (F) Total seed fatty acid content of embryos dissected from mature dry seeds. Values are the means and SE of five replicates carried out on batches of 20 embryos from five plants. (G) Total seed fatty acid content of endosperm fractions dissected from mature dry seeds. Values are the means and SE of five replicates carried out on batches of 20 endosperm fractions from five plants. Asterisks indicate significant differences from the wild type according to *t*-test at *** P<0.001, ** P<0.01, and * P<0.05, respectively.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0192156.g003
<figure><img src="image_14.png" /><figcaption>Fig 4. Analysis of *MYB115*, *AAD2*, and *AAD3* transcript levels by quantitative RT-PCR in *ProAT2S2:MYB115* embryos. The four independent transformants considered, TAR3, TXR2, TGR5, and TJR2 (in order from left to right), are displayed in the same order in the three graphs in the figure. RT-qPCR analysis of transcript abundance in cDNA prepared from excised embryos aged 14 DAA was carried out to assess efficient overexpression of the transgene (*MYB115*) and of its targets (*AAD2* and *AAD3*). Values are the means and SE of three replicates performed on cDNA</figcaption></figure>
dilutions obtained from three independent mRNA extractions. Asterisks indicate significant differences from the wild type according to *t*-test at \*\*\* P<0.001, \*\*P<0.01, and \*P<0.05, respectively.

https://doi.org/10.1371/journal.pone.0192156.g004

increased ω-7 levels are observed only in embryos of *MYB115* overexpressing lines when the production of ω-7 monoenes is increased both in the embryo and in the endosperm of *ProAT2S2:AAD2* and *ProAT2S2:AAD3* seeds. These results suggest that transcriptional activation of PAD genes constitutes the main limiting factor in the control of ω-7 accumulation in *A. thaliana* seeds. This is in agreement with the QTL approach described by Bryant et al. [11], which led to the identification of a single major QTL determining the ω-7 FA content in seeds, the peak of which colocalized with *AAD3*. Other actors involved in the biosynthesis (acyl-carrier proteins and ferredoxins), the export, or the acylation of ω-7 monoenes seem to be expressed in seeds at a level sufficient to support an important shift from ω-9 to ω-7 FA production. It is tempting to speculate that the substrate range of the enzymes of lipid metabolism using ω-9 monoenes in a wild-type context is broad enough to efficiently metabolize ω-7 monoenes in a transgenic background. This is probably due to the fact that the structures of ω-9 and ω-7 FAs are similar enough. The structure of *cis*-ω-12 C18:1 (petroselenic acid) or *cis*-ω-10 C16:1 (sapienic acid) is more divergent. These monoenes are accumulated at high levels in seeds of *Coriandrum sativum* and *Thunbergia alata*, respectively. Specialized AADs catalyzing the synthesis of these unusual monoenes have been identified and characterized [37,38]. Expression of the corresponding cDNAs in seeds of *A. thaliana* yielded limited accumulation (<15 mol%) of the corresponding FAs [39], suggesting that specialized actors other than the AADs [40–42] may be required for efficient channeling of these unusual monoene species.

The contrasted seed ω-7 contents measured in lines overexpressing *AAD2* and *AAD3* (Fig 5; Table 1) demonstrate that *AAD2* constitutes a better candidate for ω-7 biosynthesis in an *A. thaliana* context. A first reason for this may reside in the improved stability of *AAD2* mRNAs, as denoted by the higher transcript levels measured in the seeds of several *ProAT2S2:AAD2* lines (Fig 5). Even in comparisons of *ProAT2S2:AAD2* and *ProAT2S2:AAD3* lines exhibiting roughly equivalent levels of transgene expression, the lines overexpressing *AAD2* always yield the highest seed ω-7 contents, suggesting that post-translational mechanisms also contribute to favor the activity of *AAD2* over that of *AAD3*. Distinct characteristics of the two desaturases, such as protein stability, kinetic properties, or ability to interact with protein partners may account for these differences. In a wild-type context, however, the contribution of *AAD3* to seed ω-7 FA biosynthesis is more important than that of *AAD2*, as indicated by the phenotypes of the corresponding single mutants [17] and the results of the QTL approach carried out by Bryant et al. [11]. The contrasted expression levels of the two desaturases in a wild-type context help explain these apparent discrepancies: in maturing seeds of the wild type, *AAD3* mRNAs accumulate at significantly higher levels, suggesting that the activity of the *AAD3* promoter is higher than that of the *AAD2* promoter [15,17]. A functional analysis of the two promoters will be required to elucidate the molecular bases of their contrasted activities.

Comparing strategies of metabolic engineering implemented by different groups to accumulate ω-7 monoenes in *A. thaliana* seeds by means of the overexpression of PAD-coding sequences originating from different organisms is rendered difficult by the different promoters and accessions used (Table 1). Despite the fact that caution should be exercised with this comparison, the *ProAT2S2:AAD2* construct assessed in this study appears unambiguously to be a particularly efficient tool to stimulate ω-7 biosynthesis in seeds of *A. thaliana*. It would now be interesting to test the efficiency of this construct in other Brassicaceae species, such as camellina, before making use of it in a final assembly of several transgenes for further improvement of the seed ω-7 content. Nguyen and colleagues [1,21] indeed established that additional
<figure><img src="image_15.png" /><figcaption>**Fig 5.** Increased ω-7 fatty acid content in seeds of transgenic *A. thaliana* lines overexpressing *AAD2* or *AAD3*. Five independent lines are presented for each construct tested. The five *ProAT2S2:AAD2* transformants considered, T2R3, T19R1, T7R3, T17R1, and T13R8 (in order from left to right); and the five *ProAT2S2:AAD3* transformants considered, T7R1, T19R5, T16R1, T20R5, and T13R2 (in order from left to right), are always displayed in the same order in the graphs in the figure. (A) RT-qPCR analysis of transcript abundance in cDNA prepared from maturing seeds aged 14 DAA to assess efficient overexpression of the transgenes (noted in italics). Values are the means and SE of three to six replicates performed on cDNA dilutions obtained from three independent mRNA extractions. (B-D) Relative proportion of ω-7 fatty acids (*cis*-ω-7 C16:1, *cis*-ω-7 C18:1, and *cis*-ω-7 C20:1) in mature dry seeds (B), in embryos (C), and in endosperm fractions (D) dissected from mature seeds. Values are the means and SE of five replicates performed on batches of 20 seeds from five plants. Asterisks indicate significant differences from the wild type according to *t*-test at *** P<0.001 and **P<0.01, respectively.*</figcaption></figure>
<figure><img src="image_16.png" /><figcaption>**Fig 6.** Effect of increased omega-7 content on seed germination. Five independent transgenic lines are presented for each construct tested. (A) Germination time-courses of unstratified seeds. White circles represent wild-type seeds, and filled circles represent transgenic seeds. (B) Germination time-courses of seeds after 5 days of cold stratification at 4°C in the dark. White circles represent wild-type seeds, and filled circles represent transgenic seeds. (C) Germination rates of seeds cold-stratified 5 days at 4°C, then exposed to a 1-hour light impulse, and finally germinated in the dark at 25°C for 7 days. The five independent *ProAT2S2:MYB115* lines considered, TAR3, TCR4, TGR5, TJR2, and TXR2 (in order from left to right); the five independent *ProAT2S2:AAD2* transformants considered, T2R3, T19R1, T7R3, T17R1, and T13R8 (in order from left to right); and the five independent *ProAT2S2:AAD3* transformants considered, T7R1, T19R5, T16R1, T20R5, and T13R2 (in order from left to right), are always displayed in the same order in the graphs in the figure.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0192156.g006
<table><caption>Table 1. Comparison of the omega-7 fatty acid content of seeds from engineered lines of *A. thaliana* obtained by different groups.</caption><thead><tr><th rowspan="2">Genotype</th><th colspan="4">Fatty acid composition (mol% of total fatty acids)</th><th rowspan="2">Reference</th></tr><tr><th>C16:1<br>omega-7</th><th>C18:1<br>omega-7</th><th>C20:1<br>omega-7</th><th>Total<br>omega-7</th></tr></thead><tbody><tr><td>ProAT2S2:AAD2 (Col-0)</td><td>4.26</td><td>30.21</td><td>17.03</td><td>51.50</td><td>This study</td></tr><tr><td>ProAT2S2:AAD3 (Col-0)</td><td>3.18</td><td>20.76</td><td>12.31</td><td>36.45</td><td>This study</td></tr><tr><td>ProNAP:DuPAD (WS)</td><td>2.38</td><td>14.02</td><td>11.50</td><td>27.90</td><td>[20]</td></tr><tr><td>ProPHAS:COM25</td><td>1.6</td><td>12.8</td><td>n.d.</td><td>n.d.</td><td>[21]</td></tr></tbody></table>

Bondaruk et al. [20] used a PAD cDNA from cat's claw (*Doxantha unguis-cati* L.) cloned downstream of a NAPIN promoter (from *Brassica napus*). Nguyen et al. [21] used a variant of the castor (*Ricinus communis*) Δ⁹-stearoyl-ACP desaturase named Com25 and obtained from enzyme evolution experiments to enhance the Δ⁹-palmitoyl-ACP desaturase activity of the enzyme. The corresponding cDNA was cloned under the control of a PHASEOLIN promoter (from *Phaseolus vulgaris*). The fatty acid composition of the best lines characterized (when available) was reported in this table. The accession of the engineered lines is indicated, when known, between brackets. n.d., not determined.

https://doi.org/10.1371/journal.pone.0192156.t001

increases in ω-7 monoene production could be achieved by appropriate redirection of metabolic flux in maturing seeds.

The *ProAT2S2:MYB115* construct, though less efficient for promoting ω-7 FA biosynthesis in *A. thaliana* seeds, contributes to improve other seed traits of agronomic importance. The size of the *ProAT2S2:MYB115* seeds is first increased by 10–20% on average without any reduction of the total FA concentration so that the total lipid content of mature dry seeds is also raised by 10–20% (Fig 3). When examining the germination of the *ProAT2S2:MYB115* seeds, a transient but significantly improved rate of germination is then observed 24 hours after the first exposure to light (Fig 6). These phenotypes are independent of the seed ω-7 content since they are not observed in seeds of the *ProAT2S2:AAD2* and *ProAT2S2:AAD3* lines. Target genes of MYB115 other than *AAD2* and *AAD3* may contribute to control these important seed traits. *In vivo* chromatin immunoprecipitation experiments would be useful to characterize the spectrum of putative targets of MYB115 and identify candidate genes that may take part in the control of seed development and physiology. Since MYB115 also participates in the transcriptional activation of PAD-coding genes in *A. thaliana* seeds, this transcription factor constituted a promising candidate for stimulating ω-7 production by biotechnological means. Unfortunately, transgenic *ProAT2S2:MYB115* plantlets selected *in vitro* did not survive when transferred to soil. This suggests that the *in planta* functions of the closely related MYB115 and MYB115 TFs are not completely overlapping, even though they share several transcriptional targets.

In summary, we assessed the efficiency of new molecular tools for metabolically engineering the production of ω-7 monoenes in plants, and we showed that high ω-7 levels (>50%) could be attained in seeds of *A. thaliana* by overexpressing an endogenous PAD isoform, AAD2. The corresponding construct may now be tested, alone or in combination with others, in oleaginous species of agronomic importance. More importantly, this research invites us to further characterize and exploit, when available, the genetic resources (gene versions and promoter sequences) originating from the species to be transformed. If numerous reports have put forward the advantages of the introgression of exogenous genetic material (e.g. limitation of gene silencing and of negative post-translational regulations of enzymes), it seems that in the particular case of ω-7 FA biosynthesis, further research might benefit by taking advantage of endogenous genetic resources. Therefore, this study may set the foundation for implementing new strategies to develop plant platforms dedicated to the production of ω-7 monoenes.
**Supporting information**

**S1 Fig. Complementary results for the characterization of AT2S2 promoter activity.** Pattern of activity of the *ProAT2S2:uidA* cassette in rosette (A) and cauline leaves (B), in inflorescences (C), and in flowers (D). For histochemical detection of GUS activity, tissues were incubated overnight in a buffer containing 2 mM each of potassium ferrocyanide and potassium ferricyanide. Microscopy observations were performed using Nomarski optics. Bars = 5 mm in (A-C), 500 µm in (D).
(TIF)

**S2 Fig. Characterization of early seed development in *ProAT2S2:MYB115* lines.** Five independent *ProAT2S2:MYB115* lines are presented: TAR3, TCR4, TGR5, TJR2, and TXR2. (A) Observation of seed development. Whole mounts of early developing seeds (from 4 to 8 DAA) and of maturing embryos (10 and 12 DAA) were observed with Nomarski optics. Bars = 50 µm. (B) Observation of maturing endosperm. Whole mounts of peeled endosperms (14 DAA) were observed with Nomarski optics. Cross-sections (upper panel) and lateral views (lower panel) are presented. Bars = 20 µm. DAA, days after anthesis; WT, wild type.
(TIF)

**S3 Fig. Analysis of *AAD2* and *AAD3* transcript levels by quantitative RT-PCR in *ProAT2S2:MYB115* seeds.** The five independent transformants considered are TCR4, TAR3, TXR2, TGR5, and TJR2 (in order from left to right). RT-qPCR analysis of transcript abundance in cDNA prepared from excised embryos aged 14 DAA was carried out to assess efficient overexpression of *AAD2* and *AAD3*. Values are the means and SE of six replicates performed on cDNA dilutions obtained from three independent mRNA extractions.
(TIF)

**S4 Fig. Complementary results for the characterization of *A. thaliana* lines overexpressing *AAD2* or *AAD3*.** The five independent *ProAT2S2:AAD2* transformants considered, T2R3, T19R1, T7R3, T17R1, and T13R8 (in order from left to right) and the five independent *ProAT2S2:AAD3* transformants considered, T7R1, T19R5, T16R1, T20R5, and T13R2 (in order from left to right), are always displayed in the same order in the graphs in the figure. (A) Mature seed dry weight. (B) Total fatty acid content of mature dry seeds, expressed in µg.mg⁻¹ DW. (C) Total seed fatty acid content of embryos dissected from mature dry seeds. (D) Total seed fatty acid content of endosperm fractions dissected from mature dry seeds. Values are the means and SE of five replicates carried out on batches of 20 individuals from five plants. Asterisks indicate significant differences from the wild type according to *t*-test at **P**<0.01 and \*P\<0.05, respectively.
(TIF)

**S1 Table. Primers used for construct preparation.**
(PDF)

**S2 Table. Primers used for quantitative RT-PCR.**
(PDF)

**S3 Table. At4g27140/AT2S2 promoter sequence (5'→3').**
(PDF)

**S4 Table. Total fatty acid composition (in mol%) of seeds from engineered lines of *A. thaliana*.**
(PDF)
**S5 Table.** Total fatty acid composition (in mol%) of endosperm fractions dissected from seeds of engineered lines of *A. thaliana*.
(PDF)

**S6 Table.** Total fatty acid composition (in Mol%) of embryos dissected from seeds of engineered lines of Arabidopsis.
(PDF)

**Acknowledgments**

We thank G. Tremblais and C. Boulard for their technical assistance, J. Hultstrand and M. Miquel for their careful and critical reading of the manuscript.

**Author Contributions**

**Conceptualization:** Sébastien Baud.

**Investigation:** Hasna Ettaki, Manuel Adrián Troncoso-Ponce, Alexandra To, Guillaume Barthole.

**Methodology:** Guillaume Barthole.

**Supervision:** Loïc Lepiniec, Sébastien Baud.

**Writing – original draft:** Sébastien Baud.

**Writing – review & editing:** Loïc Lepiniec, Sébastien Baud.

**References**

1. Nguyen HT, Park H, Koster KL, Cahoon RE, Nguyen HTM, Shanklin J, et al. Redirection of metabolic flux for higher levels of omega-7 monounsaturated fatty acid accumulation in camelina seeds. Plant Bio- tech J. 2015; 13: 38–50.

2. Rybak A, Fokou PA, Meier MAR. Metathesis as a versatile tool in oleochemistry. Eur J Lipid Sci Tech- nol. 2008; 110: 797–804.

3. Rasti B, Erfanian A, Selamat J. Novel nanoliposomal encapsulated omega-3 fatty acids and their appli- cations in food. Food Chem. 2017; 230: 690–696. https://doi.org/10.1016/j.foodchem.2017.03.089 PMID: 28407968

4. Cao HM, Gerhold K, Mayers JR, Wiest MM, Watkins SM, Hotamisligil GS. Identification of a lipokine, a lipid hormone linking adipose tissue to systemic metabolism. Cell. 2008; 134: 933–944. https://doi.org/ 10.1016/j.cell.2008.07.048 PMID: 18805087

5. Stefan N, Kantartzis K, Celebi N, Staiger H, Machann J, Schick F, et al. Circulating palmitoleate strongly and independently predicts insulin sensitivity in humans. Diabetes Care. 2009; 33: 405–407. https://doi. org/10.2337/dc09-0544 PMID: 19889804

6. Wu Y, Li R, Hildebrand DF. Biosynthesis and metabolic engineering of palmitoleate production, an important contributor to human health and sustainable industry. Prog Lipid Res. 2012; 51: 240–249.

7. Fatima T, Snyder CL, Schroeder WR, Cram D, Datla R, Wishart D, et al. Fatty acid composition of developing sea buckthorn (*Hippophae rhamnoides* L.) berry and the transcriptome of the mature seed. PLoS ONE. 2012; 7: e34099. https://doi.org/10.1371/journal.pone.0034099 PMID: 22558083

8. Barthet VJ. (n-7) and (n-9) cis-monounsaturated fatty acid contents of 12 Brassica species. Phytochem- istry. 2008; 69: 411–417. https://doi.org/10.1016/j.phytochem.2007.08.016 PMID: 17889912

9. Penfield S, Rylott EL, Gilday AD, Graham S, Larson TR, Graham IA. Reserve mobilization in the Arabi- dopsis endosperm fuels hypocotyl elongation in the dark, is dependent of abscisic acid, and requires PHOSPHOENOLPYRUVATE CARBOXYKINASE1. Plant Cell. 2004; 16: 2705–2718. https://doi.org/ 10.1105/tpc.104.024711 PMID: 15367715

10. Li Y, Beisson F, Ohlrogge J. Oil content of Arabidopsis seeds: the influence of seed anatomy, light and plant-to-plant variation. Phytochemistry. 2006; 67: 904–915. https://doi.org/10.1016/j.phytochem.2006. 02.015 PMID: 16600316
11. Bryant FM, Munoz-Azcarate O, Kelly AA, Beaudoin F, Kurup S, Eastmond PJ. ACYL-ACYLCARRIER PROTEIN DESATURASE2 and 3 are responsible for making omega-7 fatty acids in the Arabidopsis Aleurone. Plant Physiol. 2016; 172: 154-162. https://doi.org/10.1104/pp.16.00836 PMID: 27462083

12. Shanklin J, Cahoon EB. Desaturation and related modifications of fatty acids. Annu Rev Plant Physiol Plant Mol Biol. 1998; 49: 611-641. https://doi.org/10.1146/annurev.arplant.49.1.611 PMID: 15012248

13. Shanklin J, Somerville C. Stearoyl-acyl-carrier-protein desaturase from higher plants is structurally unrelated to the animal and fungal homologs. Proc Natl Acad Sci USA. 1991; 88: 2510-2514. PMID: 2006187

14. McKeon TA, Stumpf PK. Purification and characterization of the stearoyl-acyl carrier protein desaturase and the acyl-acyl carrier protein thioesterase from maturing seeds of safflower. J Biol Chem. 1982; 257: 12141-12147. PMID: 7118934

15. Kachroo A, Shanklin J, Whittle E, Lapchyk L, Hildebrand D, Kachroo P. The Arabidopsis stearoyl-acyl carrier protein-desaturase family and the accumulation of leaf isoforms to oleic acid synthesis. Plant Mol Biol. 2007; 63: 257-271. https://doi.org/10.1007/s11103-006-9086-y PMID: 17072561

16. Cahoon EB, Shah S, Shanklin J, Browse J. A determinant of substrate specificity predicted from the acyl-acyl carrier protein desaturase of developing cat's claw seed. Plant Physiol. 1998; 117: 593-598. PMID: 9625712

17. Troncoso-Ponce MA, Barthole G, Tremblais G, To A, Miquel M, Lepiniec L, et al. Transcriptional activation of two delta-9 palmitoyl-ACP desaturase genes by MYB115 and MYB118 is critical for biosynthesis of omega-7 monounsaturated fatty acids in the endosperm of Arabidopsis seeds. Plant Cell. 2016; 28: 2666-2682. https://doi.org/10.1105/tpc.16.00612 PMID: 27681170

18. Lindqvist Y, Huang W, Schneider G, Shanklin J. Crystal structure of Δ⁹ stearoyl-acyl carrier protein desaturase from castor seed and its relationship to other di-iron proteins. EMBO J. 1996; 15: 4081-4092. PMID: 8861937

19. Cahoon EB, Lindqvist Y, Schneider G, Shanklin J. Redesign of soluble fatty acid desaturases from plants for altered substrate specificity and double bond position. Proc Natl Acad Sci USA. 1997; 94: 4872-4877. PMID: 9144157

20. Bondaruk M, Johnson S, Degafu A, Boora P, Bilodeau P, Morris J, et al. Expression of a cDNA encoding palmitoyl-acyl carrier protein desaturase from cat's claw (*Doxantha unguis-cati*) L. in *Arabidopsis thaliana* and *Brassica napus* leads to accumulation of unusual unsaturated fatty acids and increased stearic acid content in the seed oil. Plant Breed. 2007; 126: 186-194.

21. Nguyen FT, Mishra G, Whittle E, Pidkowich MS, Bevan SA, Owens Merlo A, et al. Metabolic engineering of seeds can achieve levels of ω-7 fatty acids comparable with the highest levels found in natural plant sources. Plant Physiol. 2010; 154: 1897-1904. https://doi.org/10.1104/pp.110.165340 PMID: 20943853

22. Whittle E, Shanklin J. Engineering Δ⁹-16:0-acyl carrier protein (ACP) desaturase specificity based on combinatorial saturation mutagenesis and logical redesign of the castor Δ⁹-18:0-ACP desaturase. J Biol Chem. 2001; 276: 21500-21505.

23. Haslam RP, Sayanova O, Jin Kim H, Cahoon EB, Napier JA. Synthetic redesign of plant lipid metabolism. Plant J. 2016; 87: 76-86. https://doi.org/10.1111/tpj.13172 PMID: 27483205

24. Baud S, Wuillème S, Dubreucq B, de Almeida A, Vuagnat C, Lepiniec L, et al. Function of plastidial pyruvate kinase in seeds of *Arabidopsis thaliana*. Plant J. 2007; 52: 405-419. https://doi.org/10.1111/j.1365-313X.2007.03232.x PMID: 17892448

25. Nakagawa T, Kurose T, Hino T, Tanaka K, Kawamukai M, Niwa Y, et al. Development of series of gate-way binary vectors, pGWBs, for realizing efficient construction of fusion genes for plant transformation. J Biosci Bioeng. 2007; 104: 34-41. https://doi.org/10.1263/jbb.104.34 PMID: 17697981

26. Pouvreau B, Baud S, Vernoud V, Morin V, Py C, Gendrot G, et al. Duplicate maize Wrinkled1 transcription factors activate target genes involved in seed oil biosynthesis. Plant Physiol. 2011; 156: 674-686. https://doi.org/10.1104/pp.111.173641 PMID: 21474435

27. Bechtold N, Ellis J, Pelletier G. In planta infiltration of adult Arabidopsis plants. C. R. Acad. Sci. Paris Life Sci. 1993; 316: 1194-1199.

28. Baud S, Vaultier M-N, Rochat C. Structure and expression profile of the sucrose synthase multigene family in Arabidopsis. J Exp Bot. 2004; 55: 397-409. https://doi.org/10.1093/jxb/erh047 PMID: 14739263

29. Barthole G, To A, Marchive C, Brunaud V, Soubigou-Taconnat L, Berger N, et al. MYB118 represses endosperm maturation in seeds of Arabidopsis. Plant Cell. 2014; 26: 3519-3537. https://doi.org/10.1105/tpc.114.130021 PMID: 25194028

30. Beeckman T, De Rycke R, Viane R, Inze D. Histological study of seed coat development in *Arabidopsis thaliana*. J Plant Res. 2000; 113: 139-148.
31. Guerche P, Tire C, Grossi de Sa F, De Clercq A, Van Montagu M, Krebbers E. Differential expression of the Arabidopsis 2S albumin genes and the effect of increasing gene family size. Plant Cell. 1990; 2: 469-478. https://doi.org/10.1105/tpc.2.5.469 PMID: 12354963

32. Le BH, Cheng C, Bui AQ, Wagmaister JA, Henry K, Pelletier J, et al. Global analysis of gene activity during Arabidopsis seed development and identification of seed-specific transcription factors. Proc Natl Acad Sci USA. 2010; 107: 8063-8070. https://doi.org/10.1073/pnas.1003530107 PMID: 20385809

33. Baud S, Kelemen Z, Thévenin J, Boulard C, Blanchet S, To A, et al. Deciphering the molecular mechanisms underpinning the transcriptional control of gene expression by master transcriptional regulators in Arabidopsis seed. Plant Physiol. 2016; 171: 1099-1112. https://doi.org/10.1104/pp.16.00034 PMID: 27208266

34. Nesi N, Debeaujon I, Jond C, Pelletier G, Caboche M, Lepiniec L. The *TT8* gene encodes a basic helix-loop-helix domain protein required for expression of *DFR* and *BAN* genes in Arabidopsis siliques. Plant Cell. 2000; 12: 1863-1878. PMID: 11041882

35. Shrestha P, Callahan DL, Singh SP, Petrie JR, Zhou X-RZ. Reduced triacylglycerol mobilization during seed germination and early seedling growth in Arabidopsis containing nutritionally important polyunsaturated fatty acids. Front Plant Sci. 2016; 7: 1402. https://doi.org/10.3389/fpls.2016.01402 PMID: 27725822

36. Napier JA. 2007. The production of unusual fatty acids in transgenic plants. Annu Rev Plant Biol. 2007; 58: 295-319. https://doi.org/10.1146/annurev.arplant.58.032806.103811 PMID: 17472567

37. Cahoon EB, Shanklin J, Ohlrogge JB. Expression of a coriander desaturase results in petroselenic acid production in transgenic tobacco. Proc Natl Acad Sci USA. 1992; 89: 11184-11188. PMID: 1454797

38. Cahoon EB, Becker CK, Shanklin J, Ohlrogge JB. cDNAs for isoforms of the Λ⁹-stearoyl-acyl carrier protein desaturase from *Thunbergia alata* endosperm. Plant Physiol. 1994; 106: 807-808. PMID: 7991701

39. Suh MC, Schultz DJ, Ohlrogge JB. What limits production of unusual monoenoic fatty acids in transgenic plants? Planta. 2002; 215: 584-595. https://doi.org/10.1007/s00425-002-0751-3 PMID: 12172841

40. Cahoon EB, Dörmann P, Ohlrogge JB. Petroselenic acid biosynthesis and production in transgenic plants. Prog Lipid Res. 1994; 33: 155-163. PMID: 8190735

41. Dörmann P, Frentzen M, Ohlrogge JB. Specificities of the acyl-acyl carrier protein (ACP) thioesterase and glycerol-3-phosphate acyltransferases for octadecenoyl-ACP isomers. Plant Physiol. 1994; 104: 839-844.

42. Suh MC, Schultz DJ, Ohlrogge JB. Isoforms of acyl carrier protein involved in seed-specific fatty acid synthesis. Plant J. 1999; 17: 679-688. PMID: 10366274