Plant-based bioproduction of insect sex pheromones has been proposed as an innovative strategy to increase the sustainability of pest control in agriculture. Here, we describe the engineering of transgenic plants producing (Z)-11-hexadecenol (Z11-16OH) and (Z)-11-hexadecenyl acetate (Z11-16OAc), two main volatile components in many Lepidoptera sex pheromone blends. We assembled multigene DNA constructs encoding the pheromone biosynthetic pathway and stably transformed them into Nicotiana benthamiana plants. The constructs contained the Amyelois transitella AtrΔ11 desaturase gene, the Helicoverpa armigera fatty acyl reductase HarFAR gene, and the Euonymus alatus diacylglycerol acetyltransferase EaDAct gene in different configurations. All the pheromone-producing plants showed dwarf phenotypes, the severity of which correlated with pheromone levels. All but one of the recovered lines produced high levels of Z11-16OH, but very low levels of Z11-16OAc, probably as a result of recurrent truncations at the level of the EaDAct gene. Only one plant line (SxPv1.2) was recovered that harboured an intact pheromone pathway and which produced moderate levels of Z11-16OAc (11.8 μg g-1 FW) and high levels of Z11-16OH (111.4 μg g-1). Z11-16OAc production was accompanied in SxPv1.2 by a partial recovery of the dwarf phenotype. SxPv1.2 was used to estimate the rates of volatile pheromone release, which resulted in 8.48 ng g-1 FW per day for Z11-16OH and 9.44 ng g-1 FW per day for Z11-16OAc. Our results suggest that pheromone release acts as a limiting factor in pheromone biodispenser strategies and establish a roadmap for biotechnological improvements.

Plant-based bio-production of insect sex pheromones has been proposed as an innovative strategy to increase the sustainability of pest control in agriculture. Here we describe the engineering of transgenic plants producing (Z)-11-hexadecen-1-ol (Z11-16OH) and (Z)-11-hexadecenyl acetate (Z11-16OAc), two main volatile components in many Lepidoptera sex pheromone blends. We assembled multiple multigene DNA constructs encoding the pheromone biosynthetic pathway and stably transformed them in Nicotiana benthamiana plants. The constructs comprised the Amyelois transitella AtrΔ11 desaturase gene, the Helicoverpa armigera farnesyl reductase HarFAR gene, and the Euonymus alatus diacylglycerol acetyltransferase EaDAct gene in different configurations. All the pheromone-producing plants showed dwarf phenotypes, whose severity correlated with pheromone levels. All but one of the recovered lines produced high levels of Z11-16OH but very low levels of Z11-16OAc, probably as a result of recurrent truncations at the level of the EaDAct gene. Only one plant line (SxPv1.2) was recovered harbouring an intact pheromone pathway and producing moderate levels of Z11-16OAc (11.8 μg g-1 FW), next to high levels of Z11-16OH (111.4 μg g-1). Z11-16OAc production was accompanied in SxPv1.2 by a partial recovery of the dwarf phenotype. SxPv1.2 was used to estimate the rates of volatile pheromone release, which resulted in 8.48 ng g-1 FW per day for Z11-16OH and 9.44 ng g-1 FW per day for Z11-16OAc. Our results suggest that pheromone release acts as a limiting factor in pheromone bio-dispenser strategies and establish a roadmap for biotechnological improvements.

We provide integrated protein sequence-based predictions via https://bio2byte.be/b2btools/. The aim of our predictions is to identify the biophysical behaviour or features of proteins that are not readily captured by structural biology and/or molecular dynamics approaches. Upload of a FASTA file or text input of a sequence provides integrated predictions from DynaMine backbone and side-chain dynamics, conformational propensities, and derived EFoldMine early folding, DisoMine disorder, and Agmata β-sheet aggregation. These predictions, several of which were previously not available online, capture ‘emergent’ properties of proteins, i.e. the inherent biophysical propensities encoded in their sequence, rather than context-dependent behaviour (e.g. final folded state). In addition, upload of a multiple sequence alignment (MSA) in a variety of formats enables exploration of the biophysical variation observed in homologous proteins. The associated plots indicate the biophysical limits of functionally relevant protein behaviour, with unusual residues flagged by a Gaussian mixture model analysis. The prediction results are available as JSON or CSV files and directly accessible via an API. Online visualisation is available as interactive plots, with brief explanations and tutorial pages included. The server and API employ an email-free token-based system that can be used to anonymously access previously generated results.