Insect pest management with sex pheromone precursors from engineered oilseed plants

Pheromones have become an environmentally friendly alternative to conventional insecticides for pest control. Most current pheromone-based pest control products target lepidopteran pests of high-value crops, as today’s manufacturing processes cannot yet produce pheromones at low enough costs to enable their use for lower-value crops, especially commodity crops. Camelina sativa seeds genetically modified to express (Z)-11-hexadecenoic acid, a sex pheromone precursor of several moth species, provided the oil from which the precursor was isolated, purified and transformed into the final pheromone. Trap lures containing this pheromone were then assessed for their capacity to manage moth pests in the field. Plant-derived pheromone lures proved equally effective as synthetic pheromone lures in monitoring the diamondback moth, Plutella xylostella, in cabbage and disrupting mating of cotton bollworm, Helicoverpa armigera, in common bean fields. Our study demonstrates the biological efficacy and economic feasibility of pheromone production in plant factories by metabolic engineering of an oilseed crop. There is an urgent need to control food pests more sustainably. This study assesses and confirms the potential and feasibility of growing pest insect sex pheromones from high-productivity oil-crop seeds.

H erbivorous insects cause losses of more than one-fifth of the world's total crop production annually 1 , losses projected to increase due to global warming 2 . Conventional chemical insecticides are extensively used for crop protection, with 400,000 tonnes of active ingredients (AIs) applied per year globally 3,4 , causing severe adverse impacts on ecosystems and public health 5 . Rapid evolution of resistance to insecticides makes pest control ever more challenging, as it requires the use of increasing doses of pesticides per treated area to achieve pest suppression 5,6 . For example, the Arthropod Pesticide Resistance Database (http:// www.pesticideresistance.org) reports 879 cases of resistance to 52 different insecticide AIs by cotton bollworm, Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae), a worldwide key pest of row crops, and 980 cases of resistance to 101 AIs by diamondback moth, Plutella xylostella (L) (Lepidoptera: Plutellidae), a global key pest of brassicas. Public concerns about insecticide residues in food and adverse impacts of pesticides on fragile ecosystems have increased market demand for safer, more environmentally friendly bio-based pesticides 1 .
To address these concerns, many countries have imposed regulations aimed at reducing reliance on neuroactive insecticides for pest control 7 , creating an opportunity for alternative pest control technologies to gain a commercial foothold. One particularly promising strategy is the use of sex pheromones to prevent or mitigate damage by insect pests, through monitoring, mass trapping or mating disruption [8][9][10] . Sex pheromones are relatively species-specific chemical signals that elicit behavioural or physiological reactions in individuals of the same species. Their species specificity, often achieved by a combination of structurally similar molecules in a precise ratio, is related to their roles in mate communication and reproductive isolation 11 . To date, mating disruption using sex pheromones has found its readiest market in high-value horticultural fruit and nut crops but has also been put to use in forestry and stored products 10 .
The global market for insect pheromones in agriculture is projected to reach US$5.7 billion by 2025 and US$7.2 billion by 2027 (https://www.marketresearch.com; https://www.reportlinker.com). Adoption of mating disruption solutions for pest control in high-value specialty crops, such as fruits, vegetables and nuts, has steadily increased over the past few decades 7,12,13 . However, pheromone-based pest control is still rare in many economically important crops, especially row crops that are cultivated in large expanses, such as soybean, maize and cotton. This is largely due to the high cost of conventional chemical synthesis of the pheromone AIs, which puts mating disruption solutions beyond the financial reach of growers of lower-value crops 12,14,15 .
To lower the cost of pheromone synthesis and promote their wider use in agricultural pest management, several alternative bio-based methods of pheromone production have been evaluated, starting with proof-of-concept metabolic engineering studies conducted in bakerʼs yeast, Saccharomyces cerevisiae 14,16 , to modify it to produce certain moth pheromones. Use of transgenic plants to, in essence, 'grow' pheromone precursors or components represents an additional green chemistry alternative to synthetic pheromone production. This not only reduces the cost of these production processes but also eliminates the need for the petroleum-based and other chemical feedstocks used in conventional pheromone synthesis. In tobacco, Nicotiana benthamiana, Ding et al. transiently introduced up to four genes coding for consecutive biosynthetic steps, successfully producing multi-component sex pheromones used by two small ermine moths, Yponomeuta evonymella and Y. padella 17 . Mixtures of acetylated fatty alcohols derived from the tobacco leaves trapped male moths of the target species as effectively as synthetically produced Insect pest management with sex pheromone precursors from engineered oilseed plants pheromones 17 . This study clearly demonstrated the possibility of transiently producing a variety of moth pheromone components in genetically modified plants. Stable transformation of tobacco species (Nicotiana spp.) to obtain the pheromone precursor (Z)-11-hexadecenoic acid over generations in tobacco leaves was subsequently shown by Xia et al. 18 .
Most lepidopteran sex pheromones are mono-or di-unsaturated long straight chain (C 10 -C 18 ) fatty acid derivatives with primary alcohols, acetates and aldehydes as functional groups. Many oilseed crops have prodigious biosynthetic machinery to produce high levels of fatty acids, which are stored in the seed tissues as triacylglycerols. Over the past ten years, the authors have explored the feasibility of producing moth pheromone compounds in plants at commercially relevant levels by tailoring the fatty acid composition of the oilseed crop, Camelina sativa (hereafter, Camelina), a short life cycle Brassicaceae crop that has spring and winter forms and can be cultivated in diverse climates and soil types 19 . Camelina produces fatty acids in concentrations of ≤40% of its seed weight and can be easily genetically transformed 20 , making it an attractive metabolic engineering platform for specialty seed oils for applications such as omega-3 fish oil functionality 21 and renewable jet fuels 22 .
Here we report the metabolic engineering of Camelina for seed oil production of (Z)-11-hexadecenoic acid, the immediate fatty acid precursor of Z11-16:OH and its aldehyde and acetate derivatives. Co-expression of genes for a plant fatty acyl-acyl carrier protein thioesterase and an insect desaturase yields high production of a pheromone precursor acid in seed triacylglycerols with desired carbon-chain length and specific double bond position and configuration. We also show the feasibility of downstream seed oil processing, including a combination of simple techniques for extraction, purification and transformation of target compounds at low cost, paving the way for viable industrial-scale production of insect pheromones from plant-produced precursors. Finally, field bioassays targeting the worldwide pests, P. xylostella and H. armigera, in cabbage and common bean crops proved that Camelina seed oil-derived pheromones are as effective as conventionally synthesized pheromones in monitoring and mating disruption.

Results
Our overall workflow, shown in Fig. 1, involved generation of Camelina transgenic lines by Agrobacterium tumefaciens-mediated floral dip transformation, identification and advance of the top-producing line, field propagation of the selected line, seed oil processing and chemical conversion of the precursor to pheromone and demonstration of the bioactivity of the Camelina-derived pheromone formulations in field monitoring and trap shutdown experiments (Fig. 1). Trap shutdown experiments measure the capacity of formulations containing sex pheromone to disrupt the ability of males to locate pheromone-lured traps in the same area. The assumption is that if males cannot locate traps containing optimized pheromone lures, which are more attractive than a calling female, then males cannot locate calling females in that same field either, indicating successful mating disruption.

Generation of pheromone precursor-producing Camelina lines.
A previously reported Camelina line, engineered to produce ~43% palmitic acid (16:0) by seed-specific expression of a Cuphea FatB cDNA of a 16:0-acyl-acyl carrier protein thioesterase 22,23 , was used as the host for seed-specific expression of the pheromone gland-specific fatty acyl-CoA desaturase Atr∆11 cDNA from the navel orangeworm, Amyelois transitella. For these experiments, the Atr∆11 cDNA was placed under control of the soybean seed-specific glycinin-1 promoter and assembled into a binary vector, with selection conferred by constitutive expression of a resistance gene for the herbicide phosphinothricin (or Basta) ( Supplementary Fig. 1). Following Agrobacterium tumefaciens-mediated transformation of the high palmitic acid Camelina line with the T-DNA from the Atr∆11 cDNA-containing binary vector, phosphinothricin-resistant independent lines (namely T 1 ) were identified, which produced T 2 seeds. The T 2 seeds were screened for (Z)-11-hexadecenoic acid production (Figs. 1 and 2a,b and Table 1). The average production of the target acid, Z11-16:Acid in T 2 seeds from 16 positively transformed plants was 10.5 ± 0.5% (mean ± standard error of the mean (SE), weight percent throughout the text) (Fig. 3a). Large individual variation of the target acid production among T 2 seeds was observed ( Supplementary Fig. 2). For example, from one of the promising transformants, CpuFatB+AtrD11_#37, which produced a relatively high content of Z11-16:Acid, the percentage of Z11-16:Acid ranged from 3.5% to 24.4% among 25 individually analysed seeds. Seeds from this plant together with two other transformants, namely CpuFatB+AtrD11_#33 and CpuFatB+AtrD11_#40, were grown and selfed to obtain T 3 seeds (Fig. 3b). Through this process, we obtained homozygous T 4 seeds with an average of 22.1 ± 0.7% of the target acid and greatly reduced individual variation ( Fig. 3c and Table 1).
Seed propagation in open fields. Depending on the region and mode of cultivation, two to three generations of Camelina may be sown and harvested in a single year. An initial propagation trial of C. sativa in open fields was performed in eastern Nebraska (USA) in 2016, sowing 15.5 g of T 3 seeds from a positive transformant line. Approximately 0.3 kg of T 4 seeds were produced from this field trial, containing ~9% of the target Z11-16:Acid of the total fatty acids. Separately, T 4 plants originated from CpuFatB+AtrD11_#37 were propagated in a greenhouse in Lund (Sweden) to produce T 5 seeds (containing 12.7 ± 0.6% of the target acid, Z11-16:Acid) and T 6 seeds (15.8 ± 1.2% Z11-16:Acid) for larger open field propagation. In 2018, 62.5 g of T 5 seeds were sown in Nebraska and 121 g of T 6 seeds in Borgeby, Skåne (Sweden), respectively. Propagation in Nebraska produced ~1.8 kg of T 6 seeds (12.7 ± 1.7% Z11-16:Acid), whereas propagation in Borgeby produced 5.5 kg of T 7 seeds (19.2 ± 0.4% Z11-16:Acid).
Seed oil processing and transformation of acid to pheromones. The field-harvested seeds containing the target acid were segregated into batches corresponding to each open field propagation trial, and the oil of each batch was extracted and purified. The fatty acid components of the crude seed oil were first converted into their corresponding methyl esters and separated based on carbon-chain length via short-path distillation (Fig. 4a). This procedure removed the majority of C 18 and longer compounds, especially those with full saturation and monounsaturation. The short-path distillates were then subjected to urea complexation 24 to remove most of the remaining saturated components (Fig. 4b).
Batch 1. From the 0.3 kg of T 4 seed material from the Nebraska 2016 open field propagation, ~60 g fatty acid methyl esters (FAMEs) were obtained, with 9.0% Z11-16:Acid. After short-path distillation and urea complexation, AgNO 3 -silica gel column chromatography was applied to produce a fraction containing 38% Z11- 16:Acid. An aliquot of this fraction was further fractionated on argentation thin layer chromatography (TLC) to produce a second fraction containing 83% Z11-16:Acid. The samples were subsequently reduced to corresponding alcohols and then converted to aldehydes and acetates ( Supplementary Fig. 3) for trapping/monitoring activity field assay. Batch 2. Another batch of seed oil was obtained from a combination of 1.8 kg T 6 seeds from Nebraska and 2.5 kg T 7 seeds from Skåne. Approximately 1 kg of FAMEs was obtained after base methanolysis, and subsequently, 240 g of the short-path distillates containing most of the C 16 components were subjected to urea complexation. After urea complexation, a total of 182 g FAMEs containing ~82 g of Z11-16:Acid was obtained in two purity levels, 40% and 65%, with some non-target common fatty acids as the major impurities. From this batch, ~144 g FAMEs were used for reduction, resulting in 126 g of corresponding alcohols, and out of that, 73 g alcohols were converted into aldehydes, which resulted in a total of 30 g of the AI Z11-16:Ald at two purity levels, 45% and 75%, for the trap shutdown activity field assay ( Supplementary Fig. 4).
Field activity evaluation. Pheromone from Batch 1 was assessed for its field monitoring trap activity against P. xylostella, and pheromone from Batch 2 was tested in a field trap shutdown trial with H. armigera. The activity of plant-derived pheromone for trapping of P. xylostella, a global pest of the genus Brassica, was tested in early 2017 in a choy sum field (Brassica rapa var. parachinensis) in Guangzhou, China. The sex pheromone of this species is a mixture of (Z)-11-hexadecenal (Z11-16:Al), (Z)-11-hexadecenyl acetate (Z11-16: OAc) and (Z)-11-hexadecenyl alcohol (Z11-16:OH), with a variable blend ratio among different geographical populations of the insect [25][26][27] . The plant-derived pheromones from Batch 1 were maintained at two purity levels, high (83%) and low (38%). Each purity level was adjusted to replicate the sex pheromone  27 , and formulated into rubber septa (100 µg pheromone per septum). In the same experiment, we also used conventionally synthesized pheromone at the same ratio and dose formulated into rubber septa as a positive control, and the solvent n-heptane alone formulated into rubber septa as a negative control. This study assessed the capacity of lures containing plant-derived pheromone at two levels of purity to capture male P. xylostella compared with lures using high-purity synthetic pheromone (Supplementary Fig. 3 and Fig. 5a). No significant differences were observed in average number of male moths captured per trap per week for weeks 1, 2, 4, 6 and 7 for all pheromone-lured traps, independent of the origin of the pheromone. However, catches by synthetic pheromone-lure traps were numerically slightly larger. Except for week 8, all pheromone-baited traps had significantly higher male trap catches than those baited with a solvent control. Interestingly, traps baited with the lower purity of plant-derived pheromone (38%) performed as well as those baited with synthetic pheromone lures throughout the first seven weeks of the monitoring experiment. The data from this B. rapa field monitoring study indicated that P. xylostella males do not discriminate between plant-derived and synthetic-derived pheromone lures. It also suggests that the impurities found in plant-derived pheromone lures are not detrimental to the attraction and capture of P. xylostella males in monitoring traps.
To verify the mating disruption field activity of formulations containing the plant-derived pheromone, Z11-16:Ald, a trap shutdown experiment was performed in a field of common bean, Phaseolus vulgaris, in São Paulo, Brazil, targeting the noctuid species, H. armigera, a global agricultural pest of row crops that uses Z11-16:Ald as its major sex pheromone component 28 . Two purity levels of the plant-derived pheromone were assessed, 45% and 75% (from Batch 2), compared with that of high-purity synthetic pheromone (≥98%). The results showed that when compared to those of the negative control treatment, the presence of the tested mating disruption formulation had a clear suppressive effect on male H. armigera captures in every one of the three pheromone treatments. Efficacy did not vary significantly between mating disruption treatments using pheromone derived from either source (synthetically produced or plant-derived), nor between those with higher or lower purity of plant-derived pheromones (45% or 75%) (Fig. 5b). The latter finding suggests that impurities found in plant-derived pheromone do not decrease efficacy of mating disruption of H. armigera.

Discussion
With growing concerns over the sustainability of chemical pesticide use, safe and eco-friendly alternatives are needed to preserve crop yield, quality and food security. A promising alternative to toxic pesticides is the use of pheromones to suppress pest populations in the field. Despite decades of research demonstrating the efficacy of pheromone-based pest control, their adoption has been limited to high-value niche markets because of the high cost of the pheromone AIs 29 but also due to availability and cost of labour involved in deploying the traditional pheromone dispensers, which are devices that must be manually deployed in the field. As described here, we established a pipeline from Camelina engineering for pheromone fatty acid production, field cultivation of engineered lines, pheromone precursor fatty acid extraction and enrichment from engineered Camelina seed oil, pheromone synthesis from Camelina-derived fatty acids, and finally, effective insect pest monitoring and mating disruption using sustainably produced pheromone formulations that are flowable and thus amenable to field application using conventional farm equipment. Notably, we found that plant-derived insect pheromones are as effective as conventional synthetic pheromones for monitoring and mating disruption, regardless of the product purities. Hence, we demonstrated that engineered Camelina can be used as a plant factory for scalable production of moth sex pheromone precursors of high market value.
Insects have evolved conserved biosynthetic pathways employing a limited number of key enzymes controlling fatty acid synthesis, desaturation, limited chain shortening or elongation and the downstream reactions to modify the oxygen-containing functional groups 30 . In an engineered heterologous platform, it is possible to reconstruct the insect biosynthetic pathway or to design novel pathways that lead to the production of the same compound. Camelina, which has high yield potential and oil content and simple genetic transformation methodology 19 , provides an ideal platform for production of moth pheromones. Monounsaturated C 18 , C 16 and C 14 pheromone precursor acids with double bonds in ∆9 or ∆11 position can be produced by only slight modification of the available pool of saturated fatty acids. Our studies show the feasibility of introducing different thioesterases for export of optimal chain length fatty acids from the plastid and co-expression of desaturases with the desired substrate and product specificity to generate pheromone precursor fatty acids with limited impacts on plant fitness. Introducing desaturases with different fatty acid chain length and regiospecificities into a Camelina production platform could expand the portfolio of monounsaturated pheromone precursor acids. Production of pheromone compounds with shorter chain length or with multiple double bonds involves additional challenges introducing the necessary genes and controlling their action. Notably, the individual variation among T 2 seeds was observed, but the selection during the following generation was aimed to take care of this, generating stable production of the target acids in the seed material. After selection over three generations of Camelina, the level of the target precursor Z11-16:Acid in the seed oil was stabilized around 20%. On the basis of variation in the target acid production between parental seeds multiplied in greenhouses and those harvested in the field, it is likely that careful selection of field production sites based on environmental conditions, including growth temperatures, will be useful for optimizing pheromone precursor fatty acid production levels.
Field data from the monitoring trap experiment demonstrated that the plant-derived pheromone lures are attractive to P. xylostella  male moths, independent of their purity. We also found that traps with plant-derived pheromone lures were as efficient as traps with high-purity synthetic pheromone lures in attracting and trapping these moths. Traps with plant-derived pheromone lures, independent of their purity, captured moths that oriented upwind, located the pheromone point source and landed on it. Plant-derived pheromone, therefore, can be used in lures for monitoring traps (Fig. 5a) and in baits for attract-and-kill formulations containing minute amounts of pesticide to kill attracted insects 15 . The mating disruption trial, using trap shutdown as an indication of suppression of orientation of males to mates, also demonstrated that the plant-derived pheromone mating disruption formulations were as effective as the high-purity synthetic pheromone mating disruption formulation. Because plant-derived pheromone formulations disrupted the males' ability to locate sex pheromone sources in the field, a key behavioural mechanism underlying mating  disruption 12 , it is conceivable that these plant-derived pheromones will be useful in the establishment of commercial pest population suppression through mating disruption. It is much more difficult to induce impairment of source finding in small test plots than larger treated areas; so our data also suggest that mating disruption in areawide programmes might be achieved with lower doses of pheromone per treated area. Together, the results of these field trials suggest that Camelina-derived pheromone could effectively replace petroleum-based synthetic pheromone when formulating pest management products. Conventional pheromone synthesis requires complex dedicated manufacturing infrastructures. Conventional pheromone synthesis is frequently a long and complex process, involving the import and transport of raw materials, chemical modification using large amounts of organic solvents and expensive catalysts and usually, production of toxic by-products, which are expensive and difficult to dispose of. Furthermore, this long process requires large financial investments over many months to obtain the target AIs for formulation and commercial use. All of these challenges combine to put pheromone-based control methods beyond the financial reach of growers of anything but high-value crops, such as orchard fruit and vegetables, with most pheromone AIs costing between US$1,000-3,500 kg −1 (ref. 15 ). Usually, pheromone formulations developed for mating disruption require 40-120 g ha −1 of AI 15 . In pheromone AI alone, formulators will spend from US$40 ha −1 for the lowest dose of the most inexpensive pheromones, up to US$400 ha −1 for the highest dose of the most expensive pheromones. Pheromone AI usually accounts for 70% of the cost of the formulation and ~30% of the cost of the final product to the grower. As a result, the area covered by crop protectants using sex pheromones as AIs has been relegated almost exclusively to high-value crops over the past five decades, reaching ~1 million ha globally 13 , which is a very small fraction of the ~1.5 billion ha of total cultivated area (Food and Agriculture Organization of the United Nations: https://www.fao. org/sustainability/news/detail/en/c/1274219).
New bio-based methods of pheromone production have the potential to reduce the cost of these AIs. In pheromone manufacturing, production of Z11-16:Acid using a conventional synthetic chemistry path is estimated to cost US$150-US$400 kg −1 , depending on the pathway, whereas bioproduction of Z11-16:Acid through the Camelina pathway costs between US$10 kg −1 and US$25 kg −1 , depending on the agronomical production of the crop and industrial extraction and isolation protocols utilized (A.M.-N., unpublished observations). The transformation of the acid precursor into the final sex pheromone adds another US$60 kg −1 to US$100 kg −1 to the cost, depending on the final AI product. Bioproduction of pheromones also reduces the dependence of pheromone manufacturers on supply chains of petroleum-based raw materials (A.M.-N., unpublished observations). The remarkable high-frequency occurrence of low-probability events have substantial negative impact on global production and commerce over the past two decades and increased the uncertainty and complexity of global supply chains to unprecedented levels 31 (McKinsey.com: https://www.mckinsey.com/business-functions/ risk/our-insights/covid-19-implications-for-business and https:// www.mckinsey.com/business-functions/operations/our-insights/ why-now-is-the-time-to-stress-test-your-industrial-supply-chain).
Here we developed plants that accumulate sex pheromone precursors for downstream processing. A different approach would be to design genetically modified plants to release the volatiles in the field for the purpose of attraction or mating disruption in a push-and-pull approach 32 , but this strategy has not yet been proven to work in the field. For example, a wheat variety was stably transformed to release (E)-β-farnesene, the alarm pheromone of many aphid pests, but the technology has not yet demonstrated to suppress or control aphid populations in the field 33 . Xia et al 34 . transiently modified the tobacco, Nicotiana benthamiana, to release a mixture of noctuid moth pheromone components, more specifically Z11-16:OH, Z11-16:OAc and Z11-16:Ald, in the laboratory, but the engineering of plants to release specific moth pheromone blends at efficacious ratios and biologically relevant doses in the field for successful mating disruption remains an outstanding challenge 34 . In contrast, pheromones derived from downstream processing of precursors from oil crops could easily replace conventionally produced synthetic pheromone used by formulators for pest control tools used in mating disruption and trapping programmes. We anticipate initially growing up to 4,000 acres (1,618 ha) of engineered Camelina under a Plant Made Pharmaceutical and Industrial permit, obtained in consultation with United States Department of Agriculture Animal and Plant Health Inspection Service (USDA APHIS). This will provide the opportunity to optimize field and bioprocessing production and logistic methods at commercial scale. Given that the engineered Camelina may produce 150 kg to 350 kg of precursor per hectare, this production level would yield ≥566 metric tons of precursor or enough to generate mating disruption formulation for ≥14.5 million acres of row crops. Plots under USDA APHIS permit are considered experimental production and require a high degree of grower stewardship to be compliant with permit conditions. Because of the cost and risks associated with permit production, we may consider deregulation of these lines, a process that may cost > US$3 million and take several years. Deregulation of these lines, which should streamline production, will ultimately be an economic decision that balances market demands, cost and risks. Genetically modified Camelina, engineered to produce insect pheromone precursors, can pave the way for adoption of pheromone-based control methods in lower-value but higher-volume row crops, such as maize, cotton and soybean, where pheromone-based solutions currently have very little penetration. Efficient downstream processing and scale up of the production of sex pheromone products, based on pheromone precursors produced by genetically modified Camelina seeds, is essential to achieve the economic threshold needed for the adoption of pheromone tools in row crop pest protection. Costs of biological production relative to conventionally produced pheromones will depend on the purities needed for successful application. Whereas attraction of male moths to traps may depend on high purity and specific ratios of pheromone components, it is likely that successful mating disruption may be achieved with lower purity AIs. We estimate that US$30 ha −1 per application is the upper threshold cost to the grower for the adoption of pheromone-based pest control solutions in row crops, which may be achieved if the cost of AI production can be reduced to US$100 kg −1 or less. This would allow pheromone-based pest control products to finally compete with conventional pesticides, reducing global reliance on these environmentally hazardous chemicals. Additional positive impacts of bio-based pheromones include increased capacity to control pests that have developed resistance to traditional pesticides, promotion of integrated pest management and organic farming, production of agricultural products with no insecticide residues and improved crop protection due to high specificity of pheromones, preserving non-target species and preventing secondary pest problems.

Generation of transgene constructs and transformation.
Camelina (Camelina sativa) var. Suneson, previously engineered for seed-specific expression of a Cuphea FatB thioesterase, was used as the transformation background 22,23 . Plants were grown on a 14 h day/10 h night, with day lengths maintained as needed with supplemental lighting (400-500 μmoles m −2 s −1 ) in a greenhouse. Day temperatures were maintained at ~25 °C, and night temperatures were held at ~19 °C. The pheromone gland-specific fatty acyl-CoA desaturase Atr∆11 cDNA from the navel orangeworm, Amyelois transitella, was introduced into the binary vector pBinGlyBar1 35 using Gateway cloning. The Atr∆11 cDNA was flanked on its 5′ end by the soybean glycinin-1 promoter and on its 3′ end by the soybean glycinin-1 3′UTR. The pBinGlyBar1 binary vector contains a constitutively expressed resistance gene for the herbicide phosphinothricin (or Basta). The resulting binary vector was introduced into Agrobacterium tumefaciens C58.
A. tumefaciens cells harbouring the binary vector were transformed by floral infiltration into the high palmitic acid Camelina background as previously described 36 . Engineered lines from collected T 1 seeds were screened for Basta resistance as described 35 . This selection identified >18 independent T 1 lines. T 2 seeds from these lines were screened for (Z)-11-hexadecenoic acid production by gas chromatography-flame ionization detection as described below. Top (Z)-11-hexadecenoic acid-producing lines were advanced under greenhouse conditions to homozygosity in the T 3 generation.

Characterization of positively transformed plants.
The seeds harvested from the florally dipped plants were sown in the soil for Basta selection. The surviving T 1 plants were considered potential positive transformants and were grown until fully mature. Seeds produced from the T 1 plants were sampled for chemical analysis to confirm the transformation.
Fatty acid composition analysis. Fatty acid methyl esters (FAMEs) were generated either by grinding pooled 25 individual seeds (for production analysis of one transformant) or by grinding each of 15 individual seeds per plant separately (for variation analysis within one transformant) in 1 ml 2% H 2 SO 4 in methanol for 1 h at 90 °C in a 4 ml glass vial. After cooling, 1 ml water and 1 ml heptane were added and vortexed. Then, the heptane phase containing the FAMEs was transferred to 1.5 ml autosampler vials for chemical analysis.
An Agilent 5975 mass-selective detector coupled to an Agilent 6890 series gas chromatograph (GC/MS) equipped with a polar column (HP-INNOWax, 30 m length × 0.25 mm ID, 0.25 μm film thickness) was used for seed fatty acid composition analysis. Helium was used as carrier gas, and samples were injected under a splitless mode, with a constant flow of 1 ml min −1 corresponding to a linear velocity of 36 cm s −1 . The inlet temperature was set at 260 °C, and the oven temperature was set at 80 °C for 1 min, then increased to 230 °C at a rate of 10 °C min −1 and held for 10 min. The temperature of transfer line and MS source were set at 280 °C and 230 °C, respectively. Reference synthetic fatty acid methyl esters (Supelco 37 FAME mix, Sigma-Aldrich, or individual compound from our lab stock) were used to confirm the identity of each fatty acid component in the seed lipid extracts in the form of their methyl esters, by comparing their retention times and mass spectra obtained from both polar and non-polar columns. Dimethyl disulfide derivatization was performed with the samples of seed FAMEs to confirm the double bond position in the target compound 37 . The samples of Dimethyl disulfide-adducts were analysed by an Agilent 5975 mass-selective detector coupled to an Agilent 7890 series gas chromatograph equipped with a non-polar column (HP-5MS, 30 m × 0.25 mm, 0.25 μm), for which the oven temperature was programmed at 80 °C for 2 min, then increased at a rate of 15 °C min −1 to 140 °C, and then increased at a rate of 5 °C min −1 to 260 °C and held for 30 min. Seed oil process, target compound extraction and purification. Approximately 0.3 kg and 4.3 kg (including 1.8 kg from Nebraska and 2.5 kg from Skåne) seed material, harvested from field trials in 2016 and 2018, respectively, were extracted in batches at room temperature. The seeds were ground and immersed in 2 volumes of n-heptane in a 2.5 l flask and stirred overnight. After extraction, the organic phase was filtered through a Buchner funnel, and the solvent was removed by a rotary evaporator. The fatty acid content in the resulting seed oil was transformed into corresponding methyl esters by base methanolysis in a 2 l flask. For each batch of the crude oil, 2 volumes of 0.5 M KOH in methanol was added and the mixture was stirred overnight at room temperature. The reaction was monitored by TLC until the band of triacylglycerol disappeared, and then an equal volume of 0.5 M HCl was added, and a volume of n-heptane equal to the neutralized mixture was added to extract the FAME product. The heptane phase containing the FAMEs was collected by a separatory funnel and rinsed twice with water, and finally, the organic solvent was recycled and water removed by a rotary evaporator under reduced pressure and at 65 °C for 1 h.
The C 16 FAME components containing the target compound were isolated from the methanolysis products via short-path distillation (Kugelrohr, Aldrich). For that, ~60 g, 350 g and 490 g of FAME products extracted from 0.3 kg, 1.8 kg and 2.5 kg seeds were distilled under high vacuum (Edwards Advanced Vacuum) and at temperatures ramping from 125 °C to 140 °C, with a gradient of 5 °C. The distillates collected at the same temperature were combined and checked by GC/MS before undergoing further processing.
To remove the saturated component from the short-path distillates, urea complexation was performed following the Hayes et al. protocol with some modifications 24 . On a scale of 300 ml methanol, 40 g of FAME distillate was added and the solution allowed to achieve homogeneity at 65 ± 1 °C. While solution temperature was maintained, 63 g (7 mol equivalent) of urea was added. Once the urea was totally dissolved, the solution was rapidly cooled down to 18-20 °C by shaking the reaction flask under cold tap water. The resultant slurry was then gravity filtrated, and the complexes were separated from the filtrate. Recovered complexes were washed with isooctane, a non-complexing solvent, to remove traces of the filtrate. Warm water (70 °C, 400 ml) and a small volume (10 ml) of n-heptane were applied to recover FAMEs from both filtrate in methanol and isolated complexes. The resulting FAME samples from top phase of filtrates and decomposed complexes were checked on GC/MS. Depending on the composition of saturated and unsaturated C 16 acids in the starting mixture, one or two rounds of complexation were performed to obtain the purified target acid in the filtrate.
Additional silver nitrate (AgNO 3 )-silica gel column chromatography was applied to the isolated FAME samples obtained from seeds of the 2016 field trial. After short-path distillation and urea complexation, the isolated samples containing mostly C 16 components were loaded on an AgNO 3-silica gel-based column (AgNO 3 5%) and eluted by n-heptane:diethyl ether:acetic acid (85:15:1 by volume). Fractions with the highest percentage of Z11-16 component were pooled (having 38% of Z11-16 methyl ester), and the majority of the pooled fractions was used in the synthesis of the final pheromones. A minor part of the pooled fraction was further separated on AgNO 3 -impregnated silica TLC plates (Silica 60, Merck) developed in n-heptane:diethyl ether:acetic acid (85:15:1 by volume). The monounsaturated band (containing 83% Z11-16 methyl ester) was eluted from the gel and used for synthesis of the final pheromones.
Chemical transformation of precursor acid to pheromone compounds. The purified precursor acid in the form of its methyl ester was converted into the corresponding alcohol using lithium aluminium hydride (LiAlH 4 ) (conversion rate 98%, yield 88%). The fatty alcohol product was oxidized into aldehyde using pyridinium chlorochromate according to Corey and Suggs 38 or acetylated into acetate pheromone by Fischer esterification with acetyl chloride nearly quantitatively.
Monitoring of the diamondback moth, Plutella xylostella. In the field trial (8 March to 3 May, 2017), four treatments including a blank control (N = 5), an optimal synthetic pheromone mixture of aldehyde, acetate and alcohol in a blend ratio of 30:70:0.1, and two seed oil-derived pheromone mixtures with different purities but with the AIs in the same amounts and proportions (as confirmed by GC analysis) were evaluated in an experimental field of Choy sum, Brassica rapa var. parachinensis in Guangzhou, China (22° 48′ 14.0″ N, 113° 26′ 26.0″ E, 54,000 m 2 ). The synthetic pheromone compounds with a purity of ≥99% were available from our laboratory collection. For each treatment, there were five replicates, and during the experimental period, traps (Csalomon, Budapest, Hungary) were replaced once with fresh pheromone baits by the end of week 6. The dispensers (red rubber septa, catalogue number 224100-020, Wheaton Science Products) were loaded with 100 µg AI per septum dissolved in 100 µl heptane and stored at −20 °C before use in the field. To install the bait in a trap, a pin was used to pierce the rubber septum through the sidewall of the trap.
Traps within a replicate were randomly placed in a row in the field on sticks at ~1 m height above ground level, separated from each other by 10-15 m. The traps were checked once per week. After each check, the traps were redistributed successively within a replicate row to eliminate any potential position effect.
Field experiment for trap shutdown. A trap shutdown experiment was performed in a field of common dry beans, Phaseolus vulgaris, in Fazenda Van den Broek, Parapanema, São Paulo, Brazil (23° 27′ 37.5″ S 48° 48′ 29.3″ W), targeting the noctuid species, H. armigera, a global agricultural pest that uses Z11-16:Ald as its major sex pheromone component, together with a minor component (Z)-9-hexadecenal (Z9-16:Ald) 28,39 . Trap catches were monitored in 30 m × 30 m plots using two Delta traps per plot placed 20 cm above the crop and baited with commercial H. armigera pheromone rubber septum lures (ISCAlure Armigera, ISCA Inc. USA). The experiment included three different treatments, each containing 1.8 g of the AI (Z11-16:Ald) per plot, two of them using plant-derived pheromone with a purity of 45% and 75%, respectively ( Supplementary Fig. 4) and the third using conventionally produced synthetic pheromone, with a purity of ≥98% (ISCA Inc. USA). Untreated plots receiving no pheromone mating disruption treatment served as a negative control. Number of replicates was four, each replicate with two monitoring traps checked weekly. Monitoring of the experimental plots before the installation of the mating disruption experiment indicated the presence of H. armigera at similar population densities across all fields.
Pheromone application to treatment plots in a SPLAT formulation was performed with a mechanical applicator. The splat portions were 0.5 g totalling 2,000 portions per ha. Thus, in each plot, 90 portions of 0.5 g were applied in six strips spaced 5 m apart. In each strip, 15 portions were applied 2 m apart. The different treatment plots were 50 m apart. Monitoring traps were placed in the middle of the plot between strips 2 and 3 and 4 and 5, respectively. Statistical analysis. The statistical analysis of trap catches (one-way analysis of variance (ANOVA) followed by a least significant difference (LSD) or Tukeyʼs test at P < 0.05 level) was performed with IBM SPSS Statistics 26.0.
Reporting summary. Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability
The data that support the findings of this study are available from the corresponding authors upon reasonable request.