# Production of Long Chain Fatty Alcohols Found in Bumblebee Pheromones by Yarrowia lipolytica

Jaroslav Hambalko, Peter Gajdoš, Jean-Marc Nicaud, Rodrigo

Ledesma-Amaro, Michal Tupec, Iva Pichová, Milan Čertík

► To cite this version:

Jaroslav Hambalko, Peter Gajdoš, Jean-Marc Nicaud, Rodrigo Ledesma-Amaro, Michal Tupec, et al. Production of Long Chain Fatty Alcohols Found in Bumblebee Pheromones by Yarrowia lipolytica. Frontiers in Bioengineering and Biotechnology, 2021, 8, pp. 1-12. 10.3389/fbioe.2020.593419. hal-03833046

HAL Id: hal-03833046

https://hal.inrae.fr/hal-03833046

Submitted on 28 Oct 2022

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
<figure><img src="image_2.png" /></figure>

Check for
updates

# Production of Long Chain Fatty Alcohols Found in Bumblebee Pheromones by Yarrowia lipolytica

Jaroslav Hambalko<fnref n="1" />, Peter Gajdoš<fnref n="1" />, Jean-Marc Nicaud<fnref n="2" />, Rodrigo Ledesma-Amaro³, Michal Tupec⁴, Iva Pichová⁴ and Milan Čertík<fnref n="1" />*

## OPEN ACCESS

**Edited by:**
M. Kalim Akhtar,
United Arab Emirates University,
United Arab Emirates

**Reviewed by:**
Matthew Wook Chang,
National University of Singapore,
Singapore
Marjorie A. Liénard,
Lund University, Sweden
Volkmar Passoth,
Swedish University of Agricultural
Sciences, Sweden

*Correspondence:
Milan Čertík,
milan-certik@stuba.sk

**Specialty section:**
This article was submitted to *Synthetic Biology*, a section of the journal *Frontiers in Bioengineering and Biotechnology*

**Received:** 10 August 2020
**Accepted:** 08 December 2020
**Published:** 08 January 2021

**Citation:**
Hambalko J, Gajdoš P, Nicaud J-M, Ledesma-Amaro R, Tupec M, Pichová I and Čertík M (2021) *Production of Long Chain Fatty Alcohols Found in Bumblebee Pheromones by Yarrowia lipolytica*. *Front. Bioeng. Biotechnol.* 8:593419. doi: 10.3389/fbioe.2020.593419

Fatty alcohols (FA-OH) are aliphatic unbranched primary alcohols with a chain of four or more carbon atoms. Besides potential industrial applications, fatty alcohols have important biological functions as well. In nature, fatty alcohols are produced as a part of a mixture of pheromones in several insect species, such as moths, termites, bees, wasps, etc. In addition, FA-OHs have a potential for agricultural applications, for example, they may be used as a suitable substitute for commercial insecticides. The insecticides have several drawbacks associated with their preparation, and they exert a negative impact on the environment. Currently, pheromone components are prepared mainly through the catalytic hydrogenation of plant oils and petrochemicals, which is an unsustainable, ecologically unfriendly, and highly expensive process. The biotechnological production of the pheromone components using engineered microbial strains and through the expression of the enzymes participating in the biosynthesis of these components is a promising approach that ensures ecological sustenance as well. The present study was aimed at evaluating the production of FA-OHs in the oleaginous yeast, *Yarrowia lipolytica*, with different lengths of fatty-acyl chains by expressing the fatty acyl-CoA reductase (FAR) *BlapFAR4* from *B. lapidarius*, producing C16:0-OH, C16:1Δ⁹-OH, and lower quantities of both C14:0-OH and C18:1Δ⁹-OH, and *BlucFAR1* from *B. lucorum*, producing FA-OHs with a chain length of 18–26 carbon atoms, in this yeast. Among the different novel *Y. lipolytica* strains used in the present study, the best results were obtained with JMY7086, which carried several lipid metabolism modifications and expressed the *BlucFAR1* gene under the control of a strong constitutive promoter *8UAS-pTEF*. JMY7086 produced only saturated fatty alcohols with chain lengths from 18 to 24 carbon atoms. The highest titer and accumulation achieved were 166.6 mg/L and 15.6 mg/g DCW of fatty alcohols, respectively. Unlike JMY7086, the *BlapFAR4*-expressing strain JMY7090 produced only 16 carbon atom-long FA-OHs with a titer of 14.6 mg/L.

**Keywords:** pheromone, fatty alcohol, reductase, metabolic engineering, *Bombus*, *Yarrowia lipolytica*

<footnote n="1">¹ Faculty of Chemical and Food Technology, Institute of Biotechnology, Slovak University of Technology, Bratislava, Slovakia,</footnote>

<footnote n="2">² French National Research Institute for Agriculture (INRAE), Food and Environment, AgroParisTech, Micalis Institute, Université Paris-Saclay, Jouy-en-Josas, France, ³ Department of Bioengineering and Imperial College Centre for Synthetic Biology, Faculty of Engineering, Imperial College London, London, United Kingdom, ⁴ Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Prague, Czechia</footnote>
INTRODUCTION

Fatty alcohols (FA-OH) are aliphatic unbranched primary alcohols with varying chain lengths ranging from 4 to 28 carbon atoms and containing either saturated or unsaturated carbon bonds (McNaught and Wilkinson, 1997). The properties and the potential applications of fatty alcohol depend on their molecular structure. Generally, fatty alcohols are used as fuels, solvents, detergents, cosmetics, lubricants, and pharmaceuticals, or may serve as precursors for other compounds such as waxes or polymers (Rutter and Rao, 2016; Wang G. et al., 2016; Wang W. et al., 2016; Borodina et al., 2018a; Cordova et al., 2019). In 2019, the global demand for fatty alcohols was estimated to be over two million tons, with an annual growth rate of 4.3%. Traditionally, these molecules are produced through the catalytic hydrogenation of petrochemicals or plant oils, which currently relies on fossil fuels or unsustainable palm farming and has significant environmental consequences such as deforestation or contribution to global warming (Rutter and Rao, 2016; Shah et al., 2016; Cordova et al., 2019). Therefore, there is an urgent requirement for a further efficient and ecologically-friendly process.

Fatty alcohols and their derivatives also have important biological functions. Insects have evolved an efficient mate-finding system that is based on volatile sex pheromones. In most species, sex pheromones are released either as a single component (Jurenka, 2004; Groot et al., 2016; Tupec et al., 2017) or as a specific blend of molecules in specific ratios, most of which are fatty acid (FA) derivatives, usually alcohols, aldehydes, or acetates. The sex pheromone blend of the bumblebees mainly comprises saturated, mono-unsaturated, and poly-unsaturated fatty alcohols with a chain length of 16-18 carbon atoms along with terpenoid compounds (Ayasse and Jarau, 2014). Sex pheromones are synthesized *de novo* in specialized cells known as pheromone glands, which have evolved from epidermal cells (Žáček et al., 2013; Tittiger and Blomquist, 2017). The pheromone biosynthesis process involves several key enzymes. In addition, the saturated FAs such as stearic and palmitic must undergo processing mediated by chain-shortening enzymes, desaturases, reductases, acetyltransferases, and oxidases, among others (Tillman et al., 1999; Jurenka, 2004; Matoušková et al., 2008; Matsumoto, 2010; Buček et al., 2013; Koutroumpa and Jacquin-Joly, 2014; Tupec et al., 2017).

The first pheromone to be identified and purified was a fatty alcohol named bombykol (10E,12Z-hexadeca-10,12-dien-1-ol), which was isolated from the silkmoth *Bombyx mori* (Butenandt et al., 1961). Enzymes involved in bombykol biosynthesis were described later (Moto et al., 2003, 2004). Since then, pheromones have been identified in several thousands of insect species and are known for their potential as attractants or repellents in agriculture (Koutroumpa and Jacquin-Joly, 2014; El-Sayed, 2019). Pheromones represent the most suitable substitute for insecticides. Insecticides have been in use in agriculture for over 50 years. However, the environmental damage caused by the insecticides and the development of insecticide resistance among insects and pests are emerging as serious threats. The most significant problems associated with the use of insecticides

include: (a) harmful effects on other organisms, including humans and the plants that rely on insects for pollination, (b) persistence of the insecticides in the biosphere, (c) worldwide spread, and (d) significant levels of pollution associated with the current methods of insecticide production (Hagström et al., 2013; Borodina et al., 2018b). In order to resolve some of these issues, synthetic pheromones were developed to control the pest insects in a species-specific manner and to maintain healthy agricultural practices (Hagström et al., 2013). For instance, tetradec-9-enyl acetate (C14:1Δ⁹-OAc) was reported to disrupt the mating efficiency of the fall armyworm when applied alone (Mitchell and McLaughlin, 1982). However, most agricultural applications of pheromones are limited by their high cost. Pheromones are expensive because their purity is paramount for eliciting a response in an insect, and the production of pheromones with such high levels of purity through chemical processes requires expensive and complicated methods while generating waste in huge quantities, which again requires disposal and increases costs (Borodina et al., 2018b).

The development of various genetic tools has allowed the characterization of heterologously produced fatty acyl-CoA reductases (FARs), which catalyze the reduction of fatty acyl-CoA precursors into the corresponding alcohols (Moto et al., 2003; Liénard et al., 2010; Antony et al., 2016; Ding et al., 2016; Tupec et al., 2019). The FAR genes are present in several species, including vertebrates, non-insect invertebrates, and fungi, with a particularly high number of FAR gene families reported in plant and insect genomes (Eirín-López et al., 2012; Buček et al., 2016; Tupec et al., 2019). Within the class Insecta, large quantities of long-chain alcohols have been identified in the pheromone mixtures of different bumblebees, including *Bombus lucorum*, *Bombus lapidarius*, and *Bombus terrestris*. Tupec et al. (2019) used heterologous expression in *S. cerevisiae* to demonstrate that bumblebees have evolved a specific FAR gene group that encodes reductases with unusual specificities and contributes to the biosynthesis of different fatty alcohols that form a part of bumble-specific pheromones.

The recent advancements in metabolic engineering and synthetic biology have enabled an environment-friendly production of FA-derived compounds, including FA-OHs and the biofuels from renewable feedstock using microbial biomass (Guo et al., 2016; Rutter and Rao, 2016). *Yarrowia lipolytica* is an oleaginous non-pathogenic yeast belonging to the Ascomycota phylum of kingdom Fungi (Abdel-Mawgoud et al., 2018), which could serve as a perfect cell factory for industrial applications (Groenewald et al., 2014). This yeast species is of great importance to researchers due to its high tolerance to a variety of organic substrates, higher salt levels in the environment, and a broad range of pH (Miller and Alper, 2019). Since the genome of *Y. lipolytica* was unraveled long ago and the tools for manipulating genomes and the knowledge of genetic engineering has also progressed dramatically, *Y. lipolytica* has become a suitable representative model organism for the production of natural biosynthetic products in the laboratory (Ledesma-Amaro and Nicaud, 2016). In this context, the present study was aimed to evaluate the ability of *Y. lipolytica* in the production of FA-OHs of different lengths, for which two
bumblebee FARs (Tupec et al., 2019) were expressed in this yeast. The *BlapFAR4* from *B. lapidarius* is capable of preferentially catalyzing the production of shorter FA-OHs (14–16 carbons), while *BlucFAR1* from *B. lucorum* prefers acyl chains containing 18–26 carbon atoms.

In the present study, we engineered multiple strains of *Y. lipolytica* for a redesigned lipid metabolism, with genes encoding *BlucFAR1* and *BlapFAR4*, to produce the FARs and obtain the corresponding fatty alcohols. The expression of *BlucFAR1* (JMY7086) stimulated the production of 18–24 carbon atom-long fatty alcohols, presenting the highest fatty alcohol production among all the strains (166.6 mg/L). Unlike JMY7086, the *BlapFAR4-expressing strain* (JMY7090) produced only 14.6 mg/L of fatty alcohols with a chain length of 16 carbons.

## MATERIALS AND METHODS

### Strains, Media Composition, and Culture Conditions

All the strains of *Escherichia coli* and *Y. lipolytica* used in the present study are listed in Table 1. The recombinant strains of *Y. lipolytica* were constructed by engineering the wild type strain W29 (ATCC 20 460). The *E. coli* strains were cultured in a lysogeny broth medium containing a suitable antibiotic (100 µg/mL of ampicillin or 50 µg/mL of kanamycin), according to the standard protocol described by Sambrook and Russell (2001). Minimal YNB, YNBURA, and YNBLeu media agar plates were used for the selection of *Y. lipolytica* transformants. The minimal YNB medium comprised 1.7 g/L yeast nitrogen base (without amino acids and ammonium sulfate; BD, Erembodegem, Belgium), 5 g/L NH4Cl, 50 mM phosphate buffer (pH 6.8), and 20 g/L glucose. The YNBURA and YNBLeu media were prepared by adding 0.1 g/L of uracil and leucine, respectively, to the YNB medium. The agar plates were prepared by adding 20 g/L agar to the respective medium. A rich YPD medium containing 10 g/L yeast extract (BD, Erembodegem, Belgium), 10 g/L peptone (BD, Erembodegem, Belgium), and 20 g/L glucose (Mikrochem, Pezinok, Slovakia) was prepared for obtaining the inoculum of *Y. lipolytica*. The medium for lipid production (MedA⁺) comprised 1.5 g/L yeast extract, 0.5 g/L NH₄Cl, 7 g/L KH₂PO₄, 5 g/L Na₂HPO₄·12H₂O, 0.1 g/L CaCl₂, 1.5 g/L MgSO₄·7H₂O, 10 mg/L ZnSO₄·7H₂O, 0.6 mg/L FeCl₃·6H₂O, 0.07 mg/L MnSO₄·H₂O, and 0.04 mg/L CuSO₄·5H₂O. The carbon source used was either glucose or crude glycerol (Mikrochem, Pezinok, Slovakia) in the concentration of 60 g/L. Owing to its high C/N ratio, this medium was suitable for the accumulation of lipids in yeasts. The MedA⁺ growth medium was prepared by modifying the MedA medium reported by Holdsworth et al. (1988). The yeast inoculum was prepared in 20 mL of the YPD medium in 100 mL flasks. Subsequently, 50 mL production medium in 250 mL baffled flasks was inoculated with a 24-h inoculum which had an optical density (OD₆₀₀) of 0.1. The cells were cultured at 28°C and 130 rpm inside an orbital shaker (Innova 40R, NB, Canada). In order to confirm fatty alcohol production, the strains of *Y. lipolytica* were cultured for 3 days, while the cultivation of the selected strain

<table><caption>**TABLE 1** | The *Escherichia coli* and *Yarrowia lipolytica* strains and plasmids used in the present study.</caption><thead><tr><th>Strain (host strain)</th><th>Plasmid/genotype</th><th>References</th></tr></thead><tbody><tr><td colspan="3"><strong>Escherichia coli</strong></td></tr><tr><td>JME1046</td><td>JMP62-pTEF-URA3ex</td><td>Lazar et al. (2013)</td></tr><tr><td>JME2563</td><td>JMP62-pTEF-LEU2ex</td><td>Dulermo et al. (2017)</td></tr><tr><td>JME2607</td><td>JMP62-8UAS-pTEF-RedStar2-LEU2ex</td><td>Dulermo et al. (2017)</td></tr><tr><td>JME3048</td><td>JMP62-8UAS-pTEF-URA3ex</td><td>Dulermo et al. (2017)</td></tr><tr><td>JME4149</td><td>JMP1046-BlucFAR1</td><td>This work</td></tr><tr><td>JME4151</td><td>JMP1046-BlapFAR4</td><td>This work</td></tr><tr><td>JME4303</td><td>JMP3048-BlucFAR1</td><td>This work</td></tr><tr><td>JME4305</td><td>JMP2607-BlapFAR4</td><td>This work</td></tr><tr><td colspan="3"><strong>Yarrowia lipolytica</strong></td></tr><tr><td>W29</td><td>MATA, wild type</td><td>Barth and Gaillardin (1996)</td></tr><tr><td>Po1d</td><td>MATA leu2-270 ura3-302 xpr2-322+pXPR2-SUC2</td><td>Barth and Gaillardin (1996)</td></tr><tr><td>JMY6697</td><td>Po1d, pTEF-BlucFAR1-URA3ex, LEU2</td><td>This work</td></tr><tr><td>JMY6698</td><td>Po1d, pTEF-BlapFAR4-URA3ex, LEU2</td><td>This work</td></tr><tr><td>JMY3501</td><td>W29 ura3-302 leu2-270 xpr2-322 Δpox1-6 Δtgl4+pXPR2-SUC2+pTEF-DGA2-LEU2ex+pTEF-GPD1-URA3ex</td><td>Lazar et al. (2014)</td></tr><tr><td>JMY3820</td><td>W29 ura3-302 leu2-270 xpr2-322 Δpox1-6 Δtgl4+pXPR2-SUC2+pTEF-DGA2+pTEF-GPD1</td><td>Lazar et al. (2014)</td></tr><tr><td>JMY7086</td><td>JMY3820, 8UAS-pTEF-BlucFAR1-URA3ex, LEU2</td><td>This work</td></tr><tr><td>JMY7090</td><td>JMY3820, 8UAS-pTEF-BlapFAR4-LEU2ex, URA3</td><td>This work</td></tr><tr><td>JMY7094</td><td>JMY3820, 8UAS-pTEF-BlucFAR1-URA3ex, 8UAS-pTEF-BlapFAR4-LEU2ex</td><td>This work</td></tr></tbody></table>

continued for 5 days, with 24-h interval sampling, to describe the kinetics of fatty alcohol formation. All the experiments were performed in three independent biological replicates.

### Plasmid and Strain Construction

The genes *BlucFar1* and *BlapFar4* were codon-optimized for *Yarrowia lipolytica* (Supplementary Figure 1). The synthetic fragments were digested using BamHI/AvrII, followed by insertion into the corresponding BamHI/AvrII sites of the already-available plasmids JME2607 and JME3048, which contained the 8UAS-pTEF promoter (Dulermo et al., 2017). The *JMP62-pTEF-LEU2ex* and *JMP62-pTEF-URA3ex* vectors were employed to complement the LEU and URA auxotrophy, respectively, in the final strain. The plasmids were digested with NotI prior to being used for *Y. lipolytica* transformation using the lithium acetate method (Le Dall et al., 1994). The transformants were selected on the YNBURA, YNBLeu, or YNB media, depending on their genotypes. Subsequently, the genomic DNA was derived from the yeast transformants as described by Querol et al. (1992). The positive transformants were confirmed using PCR. The PCR amplifications were performed in an Eppendorf 2720 thermal cycler using GoTaq DNA polymerases (Promega). The obtained PCR fragments were purified using a
QIAgen Purification Kit (Qiagen, Hilden, Germany), followed by verification through gel electrophoresis and sequencing. All the reactions were performed in accordance with the respective manufacturers' instructions.

## Analytical Methods

In order to isolate the biomass, the cell suspensions were centrifuged (2,880 × g, 5 min), washed twice with the saline solution (NaCl, 9 g/L), and then once with deionized water, and finally freeze-dried. The freeze-dried cells were subjected to lipid analysis and dry cell weight (DCW) determination. DCW was determined gravimetrically.

The residual glycerol amounts during growth profile analysis were determined by performing HPLC (Agilent Technologies, Santa Clara, CA, United States) using an Aminex HPX87H column (Bio-Rad, Hercules, CA, United States) coupled with an RI detector and a DAD detector. $H_2SO_4$ (5 mM) was used as the mobile phase with a flow rate of 0.6 mL/min, as described by Lazar et al. (2011).

The freeze-dried cells (approximately 10 mg) were added to a mixture of 1 mL $CH_2Cl_2$ (containing 0.1 mg of C13:0 as the internal standard) and 2 mL anhydrous methanolic HCl solution, and the suspension was incubated at 50°C for 3 h. After the incubation, 1 mL of water and 1 mL of hexane were added, and the whole suspension was vortexed vigorously. The organic layer containing fatty alcohols and fatty acid methyl esters (FAME) was separated through centrifugation (2,880 × g, 5 min) and analyzed using GC-6890 N (Agilent Technologies, Santa Clara, CA, United States). The samples (1 µL) were injected automatically into the DB-23 column (50% cyanopropyl-methylpolysiloxane, length 60 m, diameter 0.25 mm, film thickness 0.25 µm) and analyzed. The analysis conditions were: carrier gas-hydrogen, inlet (230°C; hydrogen flow: 37 mL/min; split-10:1), FID detector (250°C, hydrogen flow: 40 mL/min, air flow: 450 mL/min.), gradient (150°C-0 min; 150-170°C-5,0°C/min; 170-220°C-6,0°C/min; 220°C-6 min; 220-230°C-6°C/min; 230°C-1 min; 230-240°C-30°C/min; 240°C-6 min). The chromatograms were analyzed using the Agilent Open LAB CDS software. The fatty alcohols and fatty acids were quantified according to the individual peak area normalized with the internal standard (C13:0). Individual fatty acids were identified according to the C4-C24 FAME standard (Supelco, Bellefonte, PA, United States). The fatty alcohol standards were obtained from Nu-Chek Prep and Sigma-Aldrich. GC-MS (EI at 70 eV) was performed to confirm the identity of the obtained peaks according to their MS spectra.

## RESULTS

### Insertion of FAR Genes Into Wild Type

#### *Y. lipolytica*

The genes *BlucFAR1* from *B. lucorum* and *BlapFAR4* from *B. lapidarius* were overexpressed, under the control of the *pTEF* promoter, in the genetic background of Pold strain (**Supplementary Figure 2**). Pold was constructed from the wild type strain W29 in an earlier study (Barth and Gaillardin, 1996).

Both new strains JMY6697 (*BlucFAR1*) and JMY6698 (*BlapFAR4*) were cultured in YPD and two $MedA^+$ media, supplemented with either glycerol or glucose as the carbon source. The W29 strain was used as a control as it has the same genetic background as the host Pold. All the strains were quite similar in terms of growth and did not differ significantly regarding the amount of accumulated lipids. The lowest lipid accumulation was obtained with the YPD medium, which is consistent with the assumption that the YPD medium is not suitable for lipid overproduction as it does not have a high C/N ratio. In both the $MedA^+$ media, the lipid accumulation amounts were similar (**Supplementary Table 1**). However, trace amounts of FA-OHs were observed only under oleaginous conditions. The JMY6697 strain (*BlucFAR1*) produced saturated FA-OHs with a chain length of 18 or more carbon atoms, while the JMY6698 strain (*BlapFAR4*) produced saturated FA-OH with a 16 carbon atom-long chain only. No FA-OHs were secreted into the medium. Since insect FA-OHs were produced in low quantities in these strains, it was decided to construct *Y. lipolytica* strains with further improved lipid metabolism to achieve higher FA-OH production by the expressed FARs.

### Metabolic Redesigning of *Y. lipolytica* for Effective Fatty Alcohol Production

New strains capable of accumulating higher amounts of lipids were constructed. The JMY3820 strain, which had all the 6 *POX* genes and *TGL4* lipase deleted and the *DGA2* and *GPD1* overexpressed using the *pTEF* promoter, was selected as the host strain for FAR expression. The JMY3820 strain is an auxotrophic version of the JMY3501 strain. Both these strains were constructed from the JMY1233 (Beopoulos et al., 2008) strain in an earlier study (Lazar et al., 2014). In both these strains, the FARs were expressed under a stronger hybrid constitutive promoter named *8UAS-pTEF* (Dulermo et al., 2017). Therefore, in total, three new strains were constructed (**Supplementary Figure 3**): JMY7086 (*8UAS-pTEF-BlucFAR1*), JMY7090 (*8UAS-pTEF-BlapFAR4*), and JMY7094 (*8UAS-pTEF-BlucFAR1* and *8UAS-pTEF-BlapFAR4*). Changes in lipid metabolism are displayed in **Figure 1**.

### Biomass, Fatty Acid, and Fatty Alcohol Production

The yeast strains were cultured in two $MedA^+$ media with different carbon sources (glucose or glycerol) for 3 days. The C/N ratio of both the media was 80 and the concentration of carbon source was 60 g/L. The non-alcohol producing JMY3501 strain was cultured under the same conditions as the control for evaluating the influence of the alcohols on cell growth and lipid accumulation. Glycerol promoted higher biomass growth and total fatty acid including fatty alcohol accumulation (TFA) in all the strains, compared to glucose (**Figure 2**). Whether growing on glucose or glycerol, all the alcohol-producing strains produced quite similar amounts of lipid-free biomass, and the amount of lipid-free biomass produced by the control strain was lower compared to that of the alcohol-producing strains (5.8 g/L vs. approx. 6.8 g/L on glucose and 5.8 g/L vs. approx.).
<figure><img src="image_3.png" /></figure>

**FIGURE 1** | Lipid metabolism of manipulated *Y. lipolytica* strains producing fatty alcohols. Overexpressed genes are shown in green. Deleted genes are shown in red. Dashed lines indicate multiple steps. DHAP-dihydroxycetone phosphate, G3P-glycerol-3-phosphate, GA3P-glyceraldehyde-3-phosphate, DAG-diacylglycerol, TAG-triacylglycerol, FFA-free fatty acid, LB-lipid body, P-peroxisome. Enzymes: GPD1-glycerol-3-phosphate dehydrogenase, GUT1-glycerol kinase, GUT2-glycerol-3-phosphate dehydrogenase, ACL1-ATP-citrate lyase 1, ACL2-ATP citrate lyase 2, FAS1-fatty acid synthase 1, FAS2-fatty acid synthase, LRO1-phospholipid:diacylglycerol acyltransferase, DGA1 and DGA2-acyl-CoA:diacylglycerol acyltransferases 1, TGL3 and TGL4-intracellular lipases, FAA1-fatty acyl-CoA synthetase, FAO1-fatty alcohol oxidase, FADH-fatty alcohol oxidase, ACS-acyl-CoA synthase, POX1-6-acyl-CoA oxidases, MFE1-multifunctional enzyme, POT1-thiolase, and FAR1 and FAR4-fatty acyl-CoA reductases.

7.5 g/L on glycerol). The most remarkable difference was observed in the FAs accumulation in yeasts, which was inversely proportional to the FA-OH accumulation. The control strain JMY3501 accumulated approximately 60% of TFA/DCW, while the alcohol-producing strains accumulated only 17-45% of it. Even a higher production of fatty alcohols did not promote the secretion of fatty alcohols into the medium.

After 3 days of cultivation, all the yeast strains grown on glycerol accumulated more FA-OHs compared to the yeast growing on glucose (9.77 mg/g in JMY7086 on glycerol vs. 6.17 mg/g in JMY7086 on glucose) (Table 2). The titer of FA-OH was also higher on glycerol (104.78 mg/L on glycerol vs. 49.97 mg/L on glucose in JMY7086). The rates of the biosynthesis of individual FA-OHs facilitated by *BlucFAR1* reductase were different, although the final amounts of FA-OH produced were quite similar, as evidenced by the growth profile (Table 3). The JMY7086 strain (*BlucFAR1*) could produce only saturated fatty alcohols with chain length ranging from C18 to C24. All the FA-OHs were produced in similar final amounts, with C22:0-OH being the most abundant one. When the *BlapFAR4* expression was driven by a stronger promoter *8UAS-pTEF* in the JMY7090 and JMY7094 strains accumulating higher amounts of lipids, an additional unsaturated FA-OH C16:1Δ⁹-OH was identified compared to the previous experiment in which only C16:0-OH was produced. The best FA-OH-producing strain

was the JMY7086 strain expressing the *BlucFAR1* gene and yielding 9.77 mg/g DCW of FA-OH (104.78 mg/L). JMY7090 (*BlapFAR4*) produced the least amount of FA-OH at 1.04 mg/g DCW (14.6 mg/L). The expected higher FA-OH production in the JMY7094 strain expressing both the reductases could not be achieved as this strain produced only 5.68 mg/g DCW and 65.04 mg/L yield. On the basis of these results, only glycerol was selected for use as the carbon source, while JMY7086 was selected as the FA-OH-producing strain for subsequent experiments.

## Fatty Acid Profile

The FA profiles of the new strains and a control strain were compared (Figure 3). In all strains, oleic acid was determined as the major FA in the intracellular lipids. The FA analysis clearly demonstrated that the expression of FARs influenced the FA profile. The FAR-expressing yeast strains that produced higher amounts of FA-OHs presented a greater change in the FA profiles. The biggest change was observed for the JMY7086 strain, which accumulated the highest levels of FA-OH among all the strains. With the increasing FA-OH amount, the stearic acid content increased dramatically (8.9% in JMY3501 vs. 27.7% in JMY7086), while the oleic acid content decreased (54.3% in JMY3501 vs. 35.4% in JMY7086). In addition, the palmitoleic acid level decreased and reached up to 1.9% in JMY7086 vs. 6.7% in JMY3501. When comparing JMY3501 and JMY7086, the ratio
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Biomass and lipid accumulation of JMY3501 (control), JMY7086 (8UAS-pTEF-BlucFAR1), JMY7090 (8UAS-pTEF-BlapFAR4), and JMY7094 (8UAS-pTEF-BlucFAR1, 8UAS-pTEF-BlapFAR4) cultured in two different media. The yeast cells were cultured in a medium with a C/N ratio of 80 and 60 g/L of the carbon substrate. The carbon source was either glucose or glycerol. ■ Lipid free biomass (g/L) and □ total fatty acids, including fatty alcohols (g/L). The number inside the TFA is lipid accumulation, expressed as the TFA to DCW ratio. Each value is an average of the values obtained from three independent experiments.</figcaption></figure>

<table><caption>**TABLE 2** | Fatty alcohol production in JMY7086 (8UAS-pTEF-BlucFAR1), JMY7090 (8UAS-pTEF-BlapFAR4), and JMY7094 (8UAS-pTEF-BlucFAR1, 8UAS-pTEF-BlapFAR4) cultivated for 72 h using different carbon sources.</caption><thead><tr><th></th><th>JMY7086</th><th>JMY7090</th><th>JMY7094</th><th>JMY7086</th><th>JMY7090</th><th>JMY7094</th></tr><tr><th></th><th colspan="3">Glucose</th><th colspan="3">Glycerol</th></tr></thead><tbody><tr><td colspan="7"><strong>FA-OH (mg/g DCW)</strong></td></tr><tr><td>C16:0-OH</td><td>-</td><td>0.13</td><td>0.53</td><td>-</td><td>0.68</td><td>1.65</td></tr><tr><td>C16:1Δ<sup>9</sup>-OH</td><td>-</td><td>0.09</td><td>0.18</td><td>-</td><td>0.36</td><td>0.49</td></tr><tr><td>C18:0-OH</td><td>0.82</td><td>-</td><td>0.44</td><td>2.58</td><td>-</td><td>0.93</td></tr><tr><td>C20:0-OH</td><td>0.82</td><td>-</td><td>0.66</td><td>2.09</td><td>-</td><td>0.78</td></tr><tr><td>C22:0-OH</td><td>2.02</td><td>-</td><td>1.56</td><td>2.93</td><td>-</td><td>1.13</td></tr><tr><td>C24:0-OH</td><td>2.51</td><td>-</td><td>0.93</td><td>2.17</td><td>-</td><td>0.70</td></tr><tr><td>Total</td><td>6.17</td><td>0.22</td><td>4.29</td><td>9.77</td><td>1.04</td><td>5.68</td></tr><tr><td colspan="7"><strong>FA-OH (mg/L)</strong></td></tr><tr><td>C16:0-OH</td><td>-</td><td>1.66</td><td>5.09</td><td>-</td><td>9.61</td><td>18.92</td></tr><tr><td>C16:1Δ<sup>9</sup>-OH</td><td>-</td><td>1.07</td><td>1.74</td><td>-</td><td>4.99</td><td>5.64</td></tr><tr><td>C18:0-OH</td><td>6.66</td><td>-</td><td>4.27</td><td>27.69</td><td>-</td><td>10.62</td></tr><tr><td>C20:0-OH</td><td>6.61</td><td>-</td><td>6.35</td><td>22.38</td><td>-</td><td>8.96</td></tr><tr><td>C22:0-OH</td><td>16.33</td><td>-</td><td>15.05</td><td>31.44</td><td>-</td><td>12.94</td></tr><tr><td>C24:0-OH</td><td>20.36</td><td>-</td><td>9.02</td><td>23.26</td><td>-</td><td>7.96</td></tr><tr><td>Total</td><td>49.97</td><td>2.73</td><td>41.52</td><td>104.78</td><td>14.60</td><td>65.04</td></tr></tbody></table>

of FAs with chains longer than 20 carbon atoms rose by 2.6-fold for C20:0, 3.4-fold for C22:0, and 2-fold for C24:0. All these changes could be observed in all the FA-OH-producing strains, depending on the amount of FA-OH. The decrease in the palmitic acid level and increase in the linoleic acid level were obvious in JMY7086 and JMY7094 (both expressing *BlucFAR1*), while no such dependence was observed in JMY7090 (expressing only *BlapFAR4*).

## Daily Production Levels of the Individual Fatty Acids and Fatty Alcohols in the Strain JMY7086

The alcohol-producing JMY7086 strain and the control strain JMY3501 were cultured in the *MedA<sup>+</sup>* medium containing

glycerol as a carbon source for 5 days, with sample retrieval every 24 h (Figure 4). After the first 24 h, no alcohol was detected in JMY7086, and the cells behaved similar to those in the control. Both the strains produced 7 g/L of DCW, of which a slightly lower FA content was accumulated in JMY7086 (15% in JMY7086 vs. 19% in JMY3501). The residual glycerol in the medium was 46 g/L. After 48 h, alcohol production had begun and differences among the strains could be observed. JMY7086 consumed less glycerol (22 g/L) compared to JMY3501 (25 g/L), the lipid-free biomass was approximately the same as that in the control (JMY3501 7.4 g/L vs. JMY7086 7.8 g/L), the DCW value was lower (JMY3501 12 g/L vs. JMY7086 10 g/L), and the accumulated TFA was less than that of the control strain (37% in JMY3501 vs. 20% in JMY7086). After 72 h, both the strains reached the respective stationary phases, and
<table><caption>TABLE 3 | Kinetics of the individual fatty acids and fatty alcohols in *Y. lipolytica*.</caption><thead><tr><th rowspan="2">Strain</th><th colspan="5">JMY7086</th><th colspan="5">JMY3501</th></tr><tr><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th>Time (days)</th><th colspan="10">FA (mg/g DCW)</th></tr></thead><tbody><tr><td>FA</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>C14:0</td><td>0.33</td><td>0.32</td><td>0.42</td><td>0.44</td><td>0.45</td><td>0.36</td><td>0.83</td><td>1.10</td><td>1.30</td><td>1.35</td></tr><tr><td>C15:0</td><td>0.20</td><td>0.18</td><td>0.54</td><td>0.69</td><td>0.83</td><td>0.22</td><td>0.37</td><td>0.80</td><td>0.95</td><td>1.09</td></tr><tr><td>C16:0</td><td>26.48</td><td>25.90</td><td>42.65</td><td>45.51</td><td>48.83</td><td>32.49</td><td>81.38</td><td>123.88</td><td>136.32</td><td>139.45</td></tr><tr><td>C16:1Δ<sup>9</sup></td><td>5.64</td><td>4.09</td><td>5.19</td><td>5.29</td><td>5.30</td><td>7.46</td><td>23.92</td><td>40.11</td><td>44.46</td><td>44.74</td></tr><tr><td>C16:0-OH</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td></tr><tr><td>C17:0</td><td>0.51</td><td>0.54</td><td>0.62</td><td>0.76</td><td>0.85</td><td>0.67</td><td>1.03</td><td>3.26</td><td>3.65</td><td>3.77</td></tr><tr><td>C16:1Δ<sup>9</sup>-OH</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td></tr><tr><td>C18:0</td><td>28.06</td><td>50.83</td><td>69.34</td><td>80.88</td><td>81.26</td><td>30.13</td><td>48.59</td><td>53.24</td><td>58.55</td><td>56.07</td></tr><tr><td>C18:1Δ<sup>9</sup></td><td>57.62</td><td>66.54</td><td>93.81</td><td>97.81</td><td>96.68</td><td>65.84</td><td>180.29</td><td>323.85</td><td>347.97</td><td>341.47</td></tr><tr><td>C18:1Δ<sup>11</sup></td><td>0.63</td><td>0.68</td><td>0.80</td><td>0.89</td><td>0.94</td><td>0.75</td><td>1.63</td><td>2.70</td><td>2.96</td><td>3.08</td></tr><tr><td>C18:2Δ<sup>9,12</sup></td><td>19.18</td><td>19.31</td><td>24.04</td><td>26.01</td><td>30.72</td><td>18.01</td><td>18.14</td><td>27.88</td><td>29.69</td><td>31.93</td></tr><tr><td>C18:0-OH</td><td>0.00</td><td>0.54</td><td>2.12</td><td>4.72</td><td>1.59</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td></tr><tr><td>C20:0</td><td>1.89</td><td>4.29</td><td>5.67</td><td>7.07</td><td>6.81</td><td>1.79</td><td>3.15</td><td>5.10</td><td>5.35</td><td>5.05</td></tr><tr><td>C20:0-OH</td><td>0.00</td><td>0.66</td><td>1.92</td><td>3.27</td><td>1.06</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td></tr><tr><td>C22:0</td><td>2.11</td><td>5.25</td><td>4.60</td><td>5.95</td><td>5.04</td><td>1.83</td><td>2.91</td><td>3.81</td><td>4.13</td><td>3.86</td></tr><tr><td>C22:0-OH</td><td>0.00</td><td>1.57</td><td>3.05</td><td>4.41</td><td>2.71</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td></tr><tr><td>C24:0</td><td>5.89</td><td>13.00</td><td>8.55</td><td>9.92</td><td>8.52</td><td>6.10</td><td>10.25</td><td>10.46</td><td>11.65</td><td>10.99</td></tr><tr><td>C24:0-OH</td><td>0.00</td><td>2.17</td><td>2.31</td><td>3.17</td><td>0.79</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td></tr><tr><td>TFA</td><td>148.55</td><td>195.87</td><td>265.63</td><td>296.78</td><td>292.38</td><td>165.66</td><td>372.50</td><td>596.21</td><td>646.98</td><td>642.86</td></tr></tbody></table>

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

the amount of DCW reached its highest value and remained nearly the same for the next few days. Interestingly, the peak for the lipid-free biomass was observed on the second day in

JMY3501, which decreased on the third day and then remained constant in the following days (approx. 5.6 g/L). Meanwhile, the amount of lipid-free biomass in JMY7086 increased to
<figure><img src="image_6.png" /></figure>

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

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

**FIGURE 4** | Kinetics of alcohol production. (A) JMY3501 (control) and (B,C) JMY7086 (8UAS-ρTEF-BlucFAR1). Both the strains were cultivated for 5 days in a glycerol medium with an initial glycerol concentration of 60 g/L. Gly, kinetics of substrate consumption (◊); CW (cell weight), DCW biomass growth (●), TFA free DCW (○), TFA, total fatty acid accumulation, including that of fatty alcohols (◊), FA-OH, fatty alcohols, mg/g DCW (▲), mg/L (Δ).

8.2 g/L on the third day. At this time, glycerol was completely consumed by JMY3501, while 6 g/L of residual glycerol remained in the medium with JMY7086. The TFA accumulation reached its maximum value 24 h later, i.e., on the fourth day, with the FA-OH-producing strain accumulating less than half of the TFA amount accumulated by the control strain (65 vs. 29%). After the fourth day, no significant change in the lipid content was observed. The peak of FA-OH production was achieved on the 4th day. The yeast cells accumulated 5% FA-OH of TFA, which represented 15.6 mg/g DCW and a 166.6 mg/L titer of alcohols. On the 5th day, the FA-OH content decreased, which occurred most probably due to the depletion of the carbon substrate in the medium. In order to survive, the cells attempted to derive energy from the fat stored inside the lipid bodies. The deletion of the genes encoding lipase (TFL4) and acyl-CoA oxidases (POX1-6) prevented the cells from consuming lipids, although the presence of functional fatty alcohol oxidase (FAO1) and fatty alcohol dehydrogenase (HFD4) continued converting the fatty alcohols into fatty acids.

The effect of fatty alcohols on fatty acid accumulation is best described in **Table 3**. On the first day of cultivation, no huge difference was observed between the two strains. While JMY7086 accumulated slightly fewer FAs, its FA profile was quite similar to the FA profile of JMY3501. On the second day and beyond, the differences between the strains began deepening. JMY3501 began accumulating higher amounts of palmitic acid, palmitoleic acid, and oleic acid, while the increased amounts of palmitic and oleic acids accumulated in JMY7086 were much lower than those in the control. In addition, from the second day onward, JMY7086 did not accumulate palmitoleic acid at all. On the contrary, JMY7086 exhibited an increased amount of accumulated stearic acid, which was higher than that of the control. Although the relative contents of other FAs were different, the absolute amounts did not differ significantly.

DISCUSSION

Fatty alcohol production using microorganisms could reduce
the production costs for insect pheromones and render the
<table><caption>**TABLE 4** | Comparison of pheromone production between *Y. lipolytica* and *S. cerevisiae* expressing the same genes.</caption><thead><tr><th rowspan="3"></th><th colspan="2">BlucFAR1</th><th colspan="2">BlapFAR4</th></tr><tr><th>JMY7086</th><th>S. cerevisiae</th><th>JMY7090</th><th>S. cerevisiae</th></tr><tr><th colspan="4">Fatty alcohols (mg/L)</th></tr></thead><tbody><tr><td>C16:0-OH</td><td>-</td><td>0.1</td><td>9.61</td><td>35</td></tr><tr><td>C16:1-OH</td><td>-</td><td>-</td><td>4.99</td><td>44</td></tr><tr><td>C18:0-OH</td><td>50.5</td><td>3.3</td><td>-</td><td>-</td></tr><tr><td>C20:0-OH</td><td>35.0</td><td>1</td><td>-</td><td>-</td></tr><tr><td>C22:0-OH</td><td>47.2</td><td>1.3</td><td>-</td><td>-</td></tr><tr><td>C24:0-OH</td><td>33.9</td><td>0.7</td><td>-</td><td>-</td></tr><tr><td>C26:0-OH</td><td>-</td><td>0.5</td><td>-</td><td>-</td></tr></tbody></table>

production process environmentally safer (Borodina et al., 2018b; Petkevicius et al., 2020). Although the literature has reported several approaches for producing fatty alcohols using *E. coli* and *S. cerevisiae*, one promising tool is the use of oleaginous microorganisms. So far, among all the oleaginous yeasts, *Y. lipolytica* is the most researched one (Borodina et al., 2018a). Host systems, either bacterial (*E. coli* or cyanobacteria) or yeast (*S. cerevisiae*, *R. toruloides*, or *Y. lipolytica*) that express non-insect enzyme sequences are reported to have reached several grams of fatty alcohols per liter of cultivation medium after optimization (Tupec et al., 2017). A fatty alcohol titer of 6.33 g/L was achieved in fed-batch fermentation in *E. coli* via the deletion of all fatty acyl-CoA thioesterases and *ldhA*, *pta*, and *ackA* genes, to starve cells of fatty acids and remove competing pathways (Liu et al., 2016). The expression of *Tyto alba* FAR (TaFAR1) in *Y. lipolytica* resulted in the production of 690 mg/L of hexadecan-1-ol (Wang G. et al., 2016), while the *Y. lipolytica* transformed with the *MhFAR* gene encoding the reductase from *Marinobacter hydrocarbonoclasticus* produced 6 g/L of alcohols representing an accumulation of 36 mg alcohols/g of cells (Cordova et al., 2019). The same reductase stimulated the production of 770 mg/L of fatty alcohols in *L. starkeyi* (Wang W. et al., 2016) and 8 g/L of fatty alcohols in *Rhodosporidium toruloides* (Fillet et al., 2015). Since the pheromone blends often contain unusual components, it is an advantage to use insect fatty acyl biosynthetic enzymes for their production.

In the present work, *Y. lipolytica* strains were engineered to have different genetic backgrounds with genes encoding reductases from *Bombus lucorum* (*BlucFAR1*) and *Bombus lapidarius* (*BlapFAR4*). When these genes were introduced into the Po1d strain under the control of the *pTEF* promoter, alcohols were detected only in trace amounts (in strains JMY6697 and JMY6698). Therefore, the second strategy was adopted, which involved redesigning the strain JMY3820, which had all the six *POX* genes and *TGL4* lipase deleted and the *DGA2* and *GPD1* overexpressed using the *pTEF* promoter, with the incorporation of FAR genes driven by a stronger promoter *8UAS-pTEF*. This genotype is particularly effective in accumulating high amounts of lipids (Dulermo et al., 2015). Deletions of *POX* genes and *TGL4* gene prevents degradation of storage lipids and oxidation of fatty acids. On the other hand,

overexpression of the *DGA2* and *GPD1* drives the formation of the TAGs and incorporation of fatty acids into the TAGs. The premise was to implement the push and pull strategy (Tai and Stephanopoulos, 2013), with the higher lipid production as the driving force for the higher fatty acid synthesis. With the higher fatty acyl-CoA pool, more fatty acyl-CoAs can be redirected to the reduction by FARs. As the cells become more saturated with TAGs, more fatty acyl-CoA becomes available as a substrate for reductases. This strategy produced three new strains: JMY7086 (*8UAS-pTEF-BlucFAR1*), JMY7090 (*8UAS-pTEF-BlapFAR4*), and JMY7094 (*8UAS-pTEF-BlucFAR1*, *8UAS-pTEF-BlapFAR4*). A co-expression of both the reductases (JMY7094) did not meet the expectations. It might have affected the FA profiles or reduced the activity of the elongases (Figure 3). This could have led to limited availability of the substrate and ultimately a lower yield of alcohols. For instance, co-expression of *BlucFAR1* and *BlapFAR4* resulted in a higher level of C16:0-OH and C16:1Δ⁹-OH and a decreased level of longer FA-OHs (Table 2). Another explanation could be as a consequence of an increased degree of FA-OH degradation. Table 3 presents the degradation of long-chain fatty alcohol between Day 4 and Day 5, which might have been induced by the presence of C16 fatty alcohols, which would also explain the lower level of total FA-OHs in JMY7090. It may also be interpreted by the level of *BlucFAR1* expression, possibly due to the integration site of *BlucFAR1* in the double mutant or titration of the *pTEF* promoter. All of the results did not indicate alcohol toxicity on the cells and rather indicated the inhibition of FA formation and lipid accumulation by the alcohols or the by-products of their metabolism. It is also possible that in addition to being affected by the total FA-OH amount, the FA profile is affected by the FA-OH composition as well. For example, Figure 3 clearly shows that the JMY7086 strain produces more long chain fatty acids (>C18). Indirectly, this tells us that the production of long-chain fatty alcohols may well be driving the elongation of fatty acids - the so-called metabolic driving force. This also corresponds to the push and pull strategy described by Tai and Stephanopoulos (2013). If FAR is taking substrates (long-chain fatty acyl-CoAs) from lipid metabolism, the cells are overcompensating by producing more long-chain fatty acyl-CoAs. Although the metabolic engineering of the selected *Y. lipolytica* strains enabled the production of the desired FA-OHs, there is scope for further improvement. A combination of lipid metabolism manipulation in *Y. lipolytica* and engineering of the FA-OH metabolism could produce higher amounts of specific FA-OHs synthesized by a particular FAR in low concentrations. It is well recognized that in certain cases, *Y. lipolytica* secretes fatty acids and fatty alcohols into the medium, which increases the production while simplifying the purification. Most often, this occurs spontaneously when the cells accumulate high amounts of toxic products that they are unable to degrade. The excretion of fatty alcohols may also be achieved by optimizing the fermentation process. The deletion of *FAO1* (fatty alcohol oxidase) and *HFD4* (fatty alcohol dehydrogenase) could prevent product degradation by the cells (Dahlin et al., 2019; Holkenbrink et al., 2020). Another approach is to manipulate the *FAA1* gene
(fatty acyl-CoA synthetase) or the lipases to increase the substrate availability for the functioning of FARs, or to use dodecane-mediated extraction fermentation for the secretion of alcohols into the medium (Cordova et al., 2019).

In a previous study, the FARs used in the present work were expressed in *Saccharomyces cerevisiae* (Tupec et al., 2019), and several discrepancies were observed between the composition of FA-OH synthesized by the FARs in labial gland of bumblebee and in the *S. cerevisiae* expression system. The reductase *BlapFAR4* produced poly-unsaturated fatty alcohols in the yeast cells, which are not present in the labial gland of male *B. lapidarius*. It is possible that a pool of fatty acyl CoA precursors in the yeasts differs from those present in the labial gland. The concentrations of C18:1Δ⁹-CoA, C18:2Δ⁹,¹²-CoA, and C18:3Δ⁹,¹²,¹⁵-CoA in the labial gland of *B. lapidarius* might be extremely low, to the extent that detectable accumulation of the corresponding alcohols is not possible. This implies that the ratio of fatty alcohols in the pheromones is dependent not only on the FAR substrate specificity, rather also on substrate availability, which is affected greatly by the other mechanisms such as the conversion of fatty acyls into fatty acyl-CoAs, fatty acyl-CoAs hydrolysis, etc. The differences in the quantity and ratio of fatty alcohols between *Y. lipolytica* and *S. cerevisiae* (Table 4) most likely reflects the differences in metabolism of these species. While *Y. lipolytica* is considered as oleaginous yeast, *S. cerevisiae* do not accumulate significant quantities of lipids. FA composition of lipids of these two yeast species is also different. While *Y. lipolytica* naturally prefers to accumulate the FAs with a chain length of 18 carbon atoms, *S. cerevisiae* prefers 16 carbon atom-long FAs. It could, therefore, be assumed that this difference in the accumulation of FAs was reflected in the fatty acyl-CoA pool as well. Since the fatty acyl-CoA pool is the source of substrates for FAR, a higher quantity of either 16 or 18 carbon long FA could influence the quantity and quality of FA-OHs. While the enzyme *BlucFAR1* produced the highest FA-OH content in *Y. lipolytica* (166.6 mg/L), the same enzyme produced only a small amount of FA-OH in *S. cerevisiae* (6.9 mg/L). On the contrary, *BlapFAR4* worked much better in *S. cerevisiae* (79 mg/L) than in *Y. lipolytica* (14.6 mg/L). Therefore, variations in yields could be attributed to the difference in the substrate availability between *Y. lipolytica* and *S. cerevisiae*. However, the expression level of the enzymes and in particular their activities also might have an effect on the synthesis of FA-OH.

## REFERENCES

* Abdel-Mawgoud, A. M., Markham, K. A., Palmer, C. M., Liu, N., Stephanopoulos, G., and Alper, H. S. (2018). Metabolic engineering in the host *Yarrowia lipolytica*. *Metab. Eng.* 50, 192–208. doi: 10.1016/j.ymben.2018.07.016

* Antony, B., Ding, B. J., Moto, K., Aldosari, S. A., and Aldawood, A. S. (2016). Two fatty acyl reductases involved in moth pheromone biosynthesis. *Sci. Rep.* 6:29927. doi: 10.1038/srep29927

* Ayasse, M., and Jarau, S. (2014). Chemical ecology of bumble bees. *Ann. Rev. Entomol.* 59, 299–319. doi: 10.1146/annurev-ento-011613-161949

Taken together, the results of the present work demonstrated that bioengineered *Y. lipolytica* strain with deleted *POX* and *TGL4* lipase genes and overexpressed *DGA2* and *GPD1* genes is suitable for the efficient production of FA-OHs with different chain lengths. To our knowledge, this is the first study showing the ability of *Y. lipolytica* to produce FA-OHs of different lengths and the first study achieving production of FA-OHs longer than 20 carbons by *Y. lipolytica*.

## DATA AVAILABILITY STATEMENT

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

## AUTHOR CONTRIBUTIONS

IP, RL-A, J-MN, and MC defined the concept of the study. MT isolated the genes. IP provided the gene sequences. RL-A designed and prepared the constructs. PG and JH designed and performed the cultivation experiments and performed the GC MS analyses. JH prepared the manuscript. IP, MT, RL-A, J-MN, MC, and PG revised the manuscript. All authors contributed to the article and approved the submitted version.

## FUNDING

The work was supported by grant (APVV-17-0262) from the Slovak Research and Development Agency, Slovak Republic.

## ACKNOWLEDGMENTS

We thank Tatiana Klempová for valuable comments and Aleš Buček for discussion.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2020.593419/full#supplementary-material

* Barth, G., and Gaillardin, C. (1996). "Yarrowia lipolytica," in *Nonconventional Yeasts in Biotechnology*, ed. K. Wolf, (Berlin: Springer), 313–388.

* Beopoulos, A., Mrozova, Z., Thevenieau, F., Le Dall, M. Z., Hapala, I., Papanikolaou, S., et al. (2008). Control of lipid accumulation in the yeast *Yarrowia lipolytica*. *Appl. Environ. Microbiol.* 74:24. doi: 10.1128/AEM.01412-08

* Borodina, I., Holkenbrink, C., Dam, M. I., and Löfstedt, C. (2018a). *Methods for Producing Fatty Alcohols and Derivatives Thereof in Yeast*. Patent No WO2018109163. Geneva: World Intellectual Property Organization.

* Borodina, I., Holkenbrink, C., Dam, M. I., Löfstedt, C., Ding, B., and Wang, H. (2018b). *Production of Desaturated Fatty Alcohols and Desaturated Fatty*
*Acyl Acetates in Yeast*. Patent No WO2018109167. Geneva: World Intellectual Property Organization.

Buček, A., Brabcová, J., Vogel, H., Prchalová, D., Kindl, J., Valterová, I., et al. (2016). Exploring complex pheromone biosynthetic processes in the bumblebee male labial gland by RNA sequencing. *Insect Mol. Biol.* 25, 295–314. doi: 10.1111/imb.12221

Buček, A., Vogel, H., Matoušková, P., Prchalová, D., Žádek, P., Vrkoslav, V., et al. (2013). The role of desaturases in the biosynthesis of marking pheromones in bumblebee males. *Insect biochemistry and molecular biology* 43, 724–731. doi: 10.1016/j.ibmb.2013.05.003

Butenandt, A., Beckmann, R., and Hecker, E. (1961). Über den sexuallockstoff des seidenspinners, I. Der biologische test und die isolierung des reinen sexuallockstoffes bombykol. *Biol. Chem.* 324, 71–83. doi: 10.1515/bchm2.1961.324.1.71

Cordova, L. T., Butler, J., and Alper, H. S. (2019). Direct production of fatty alcohols from glucose using engineered strains of *Yarrowia lipolytica*. *Metab. Eng. Commun.* 10:e00105. doi: 10.1016/j.mec.2019.e00105

Dahlin, J., Holkenbrink, C., Marella, E. R., Wang, G., Liebal, U., Lieven, C., et al. (2019). Multi-omics analysis of fatty alcohol production in engineered yeasts *Saccharomyces cerevisiae* and *Yarrowia lipolytica*. *Front. Genet.* 10:747. doi: 10.3389/fgene.2019.00747

Ding, B. J., Lager, I., Bansal, S., Durrett, T. P., Stymne, S., and Löfstedt, C. (2016). The yeast ATF1 acetyltransferase efficiently acetylates insect pheromone alcohols: implications for the biological production of moth pheromones. *Lipids* 51, 469–475. doi: 10.1007/s11745-016-4122-4

Dulermo, R., Brunel, F., Dulermo, T., Ledesma-Amaro, R., Vion, J., Trassaert, M., et al. (2017). Using a vector pool containing variable-strength promoters to optimize protein production in *Yarrowia lipolytica*. *Microb. Cell Fact.* 16:31.

Dulermo, T., Lazar, Z., Dulermo, R., Rakicka, M., Haddouche, R., and Nicaud, J.-M. (2015). Analysis of ATP-citrate lyase and malic enzyme mutants of *Yarrowia lipolytica* points out the importance of mannitol metabolism in fatty acid synthesis. *Biochim. Biophys. Acta* 1851, 1107–1117. doi: 10.1016/j.jbbalip.2015.04.007

Eirín-López, J. M., Rebordinos, L., Rooney, A. P., and Rozas, J. (2012). The birth-and-death evolution of multigene families revisited. *Genome Dynamics* 7, 170–196. doi: 10.1159/000337119

El-Sayed, A. M. (2019). *The Pherobase: Database of Pheromones and Semiochemicals*. Available online at: http://www.pherobase.com (accessed October 20, 2020).

Fillet, S., Gibert, J., Suárez, B., Lara, A., Ronchel, C., and Adrio, J. L. (2015). Fatty alcohols production by oleaginous yeast. *J. Indust. Microbiol. Biotechnol.* 42, 1463–1472. doi: 10.1007/s10295-015-1674-x

Groenewald, M., Boekhout, T., Neuvéglise, C., Gaillardin, C., van Dijck, P. W., and Wyss, M. (2014). *Yarrowia lipolytica*: safety assessment of an oleaginous yeast with a great industrial potential. *Crit. Rev. Microbiol.* 40, 187–206. doi: 10.3109/1040841X.2013.770386

Groot, A. T., Dekker, T., and Heckel, D. G. (2016). The genetic basis of pheromone evolution in moths. *Annu. Rev. Entomol.* 61, 99–117. doi: 10.1146/annurev-ento-010715-023638

Guo, W., Sheng, J., Zhao, H., and Feng, X. (2016). Metabolic engineering of *Saccharomyces cerevisiae* to produce 1-hexadecanol from xylose. *Microb. Cell Fact.* 15:24. doi: 10.1186/s12934-016-0423-9

Hagström, ÅK., Wang, H. L., Liénard, M. A., Lassance, J. M., Johansson, T., and Löfstedt, C. (2013). A moth pheromone brewery: production of (Z)-11-hexadecenol by heterologous co-expression of two biosynthetic genes from a noctuid moth in a yeast cell factory. *Microb. Cell Fact.* 12:125. doi: 10.1186/1475-2859-12-125

Holdsworth, J. E., Veenhuis, M., and Ratledge, C. (1988). Enzyme activities in oleaginous yeasts accumulating and utilizing exogenous or endogenous lipids. *J. Gen. Microbiol.* 1988, 2907–2915. doi: 10.1099/00221287-134-11-2907

Holkenbrink, C., Ding, B.-J., Wang, H.-L., Dam, M. I., Petkevicius, K., Kildegaard, K., et al. (2020). Production of moth sex pheromones for pest control by yeast fermentation. *Metab. Eng.* 62, 312–321. doi: 10.1101/2020.07.15.205047

Jurenka, R. (2004). Insect pheromone biosynthesis. *Top. Curr. Chem.* 239, 97–132. doi: 10.1007/b95450

Koutroumpa, F. A., and Jacquin-Joly, E. (2014). Sex in the night: fatty acid-derived sex pheromones and corresponding membrane pheromone receptors in insects. *Biochimie* 107, 15–21. doi: 10.1016/j.biochi.2014.07.018

Lazar, Z., Dulermo, T., Neuvéglise, C., Crutz-Le Coq, A. M., and Nicaud, J. M. (2014). Hexokinase—A limiting factor in lipid production from fructose in *Yarrowia lipolytica*. *Metab. Eng.* 26, 89–99. doi: 10.1016/j.ymben.2014.09.008

Lazar, Z., Rossignol, T., Verbeke, J., Crutz-Le Coq, A. M., Nicaud, J.-M., Robak, M. (2013). Optimized invertase expression and secretion cassette for improving *Yarrowia lipolytica* growth on sucrose for industrial applications. *J. Ind. Microbiol. Biotechnol.* 40, 1273–1283. doi: 10.1007/s10295-013-1323-1

Lazar, Z., Walczak, E., and Robak, M. (2011). Simultaneous production of citric acid and invertase by *Yarrowia lipolytica* SUC+ transformants. *Bioresour. Technol.* 102, 6982–6989. doi: 10.1016/j.biortech.2011.04.032

Le Dall, M.-T., Nicaud, J.-M., Gaillardin, C. (1994). Multi-copy integration in the yeast *Yarrowia lipolytica*. *Curr. Genet.* 26, 38–44. doi: 10.1007/BF00326302

Ledesma-Amaro, R., and Nicaud, J.-M. (2016). *Yarrowia lipolytica* as a biotechnological chassis to produce usual and unusual fatty acids. *Prog. Lipid Res.* 61, 40–50. doi: 10.1016/j.plipres.2015.12.001

Liénard, M. A., Hagström, ÅK., Lassance, J.-M., and Löfstedt, C. (2010). Evolution of multi-component pheromone signals in small ermine moths involves a single fatty-acyl reductase gene. *Proc. Natl. Acad. Sci. U.S.A.* 107, 10955–10960. doi: 10.1073/pnas.1000823107

Liu, Y., Chen, S., Chen, J., Zhou, J., Wang, Y., Yang, M., et al. (2016). High production of fatty alcohols in *Escherichia coli* with fatty acid starvation. *Microb. Cell Fact.* 15:129. doi: 10.1186/s12934-016-0524-5

Matoušková, P., Luxová, A., Matoušková, J., Jiroš, P., Svatoš, A., Valterová, I., et al. (2008). A delta(9) desaturase from *Bombus lucorum* males: investigation of the biosynthetic pathway of marking pheromones. *Chembiochem* 9, 2534–2541. doi: 10.1002/cbic.200800374

Matsumoto, S. (2010). Molecular mechanisms underlying sex pheromone production in moths. *Biosci. Biotechnol. Biochem.* 74, 223–231. doi: 10.1271/bbb.90756

McNaught, A. D., and Wilkinson, A. (1997). *IUPAC. Compendium of Chemical Terminology*, 2 Edn. Oxford: Blackwell Scientific Publications. doi: 10.1351/goldbook

Miller, K. K., and Alper, H. S. (2019). *Yarrowia lipolytica*: more than an oleaginous workhorse. *Appl. Microbiol. Biotechnol.* 103, 9251–9262. doi: 10.1007/s00253-019-10200-x

Mitchell, E. R., and McLaughlin, J. R. (1982). Suppression of mating and oviposition by fall armyworm and mating by corn earworm in corn, using the air permeation technique. *J. Econom. Entomol.* 75, 270–274. doi: 10.1093/jee/75.2.270

Moto, K., Suzuki, M. G., Hull, J. J., Kurata, R., Takahashi, S., Yamamoto, M., et al. (2004). Involvement of a bifunctional fatty-acyl desaturase in the biosynthesis of the silkmoth, *Bombyx mori*, sex pheromone. *Proc. Natl. Acad. Sci. U.S.A.* 101, 8631–8636. doi: 10.1073/PNAS.0402056101

Moto, K., Yoshiga, T., Yamamoto, M., Takahashi, S., Okano, K., Ando, T., et al. (2003). Pheromone gland-specific fatty-acyl reductase of the silkmoth, *Bombyx mori*. *Proc. Natl. Acad. Sci. U.S.A.* 100, 9156–9161. doi: 10.1073/pnas.1531993100

Petkevicius, K., Löfstedt, C. H., and Borodina, I. (2020). Insect sex pheromone production in yeasts and plants. *Curr. Onin. Biotechnol.* 65, 259–267. doi: 10. 1016/j.copbio.2020.07.011

Querol, A., Barrio, E., Ramon, D. (1992). A comparative study of different methods of yeast strain characterization. *Syst. Appl. Microbiol.* 15, 439–446. doi: 10.1016/ S0723-2020(11)02192-5

Rutter, C. D., and Rao, C. V. (2016). Production of 1-decanol by metabolically engineered *Yarrowia lipolytica*. *Metab. Eng.* 38, 139–147. doi: 10.1016/j.ymben. 2016.07.011

Sambrook, J., and Russell, D. W. (2001). *Molecular Cloning: A Laboratorymanual*, 3 Edn. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press.

Shah, J., Arslan, E., Cirucci, J., O'Brien, J., and Moss, D. (2016). Comparison of oleo- vs petro-sourcing of fatty alcohols via cradle-to-gate life cycle assessment. *J. Surfact. Deterg.* 19, 1333–1351. doi: 10.1007/s11743-016-1867-y

Tai, M., and Stephanopoulos, G. (2013). Engineering the push and pull of lipid biosynthesis in oleaginous yeast *Yarrowia lipolytica* for biofuel production. *Metab. Eng.* 15, 1–9. doi: 10.1016/j.ymben.2012.08.007
Tillman, J. A., Seybold, S. J., Jurenka, R. A., and Blomquist, G. J. (1999). Insect pheromones—an overview of biosynthesis and endocrine regulation. *Insect Biochem. Mol. Biol.* 29, 481–514. doi: 10.1016/s0965-1748(99)00016-8

Tittiger, C., and Blomquist, G. J. (2017). Pheromone biosynthesis in bark beetles. *Curr. Opin. Insect Sci.* 24, 68–74. doi: 10.1016/j.cois.2017.09.005

Tupec, M., Buček, A., Janoušek, V., Vogel, H., Prchalová, D., Kindl, J., et al. (2019). Expansion of the fatty acyl reductase gene family shaped pheromone communication in Hymenoptera. *eLife* 8:e39231. doi: 10.7554/eLife.39231

Tupec, M., Buček, A., Valterová, I., and Pichová, I. (2017). Biotechnological potential of insect fatty acid-modifying enzymes. *Zeitschr. Naturfors.* C 72, 387–403. doi: 10.1515/znc-2017-0031

Wang, G., Xiong, X., Ghogare, R., Wang, P., Meng, Y., and Chen, S. (2016). Exploring fatty alcohol-producing capability of Yarrowia lipolytica. *Biotechnol. Biof.* 9:107. doi: 10.1186/s13068-016-0512-3

Wang, W., Wei, H., Knoshaug, E., Van Wychen, S., Xu, Q., Himmel, M. E., et al. (2016). Fatty alcohol production in *Lipomyces starkeyi*

and *Yarrowia lipolytica*. *Biotechnol. Biof.* 9:227. doi: 10.1186/s13068-016-0647-2

Žáček, P., Prchalová-Horòáková, D., Tykva, R., Kindl, J., Vogel, H., Svatoš, A., et al. (2013). De novo biosynthesis of sexual pheromone in the labial gland of bumblebee males. *Chembiochem* 14, 361–371. doi: 10.1002/cbic.201200684

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Hambalko, Gajdoš, Nicaud, Ledesma-Amaro, Tupec, Pichová and Čertík. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.