- Open Access
Engineering the fatty acid synthesis pathway in Synechococcus elongatus PCC 7942 improves omega-3 fatty acid production
Biotechnology for Biofuelsvolume 11, Article number: 239 (2018)
The microbial production of fatty acids has received great attention in the last few years as feedstock for the production of renewable energy. The main advantage of using cyanobacteria over other organisms is their ability to capture energy from sunlight and to transform CO2 into products of interest by photosynthesis, such as fatty acids. Fatty acid synthesis is a ubiquitous and well-characterized pathway in most bacteria. However, the activity of the enzymes involved in this pathway in cyanobacteria remains poorly explored.
To characterize the function of some enzymes involved in the saturated fatty acid synthesis in cyanobacteria, we genetically engineered Synechococcus elongatus PCC 7942 by overexpressing or deleting genes encoding enzymes of the fatty acid synthase system and tested the lipid profile of the mutants. These modifications were in turn used to improve alpha-linolenic acid production in this cyanobacterium. The mutant resulting from fabF overexpression and fadD deletion, combined with the overexpression of desA and desB desaturase genes from Synechococcus sp. PCC 7002, produced the highest levels of this omega-3 fatty acid.
The fatty acid composition of S. elongatus PCC 7942 can be significantly modified by genetically engineering the expression of genes coding for the enzymes involved in the first reactions of fatty acid synthesis pathway. Variations in fatty acid composition of S. elongatus PCC 7942 mutants did not follow the pattern observed in Escherichia coli derivatives. Some of these modifications can be used to improve omega-3 fatty acid production. This work provides new insights into the saturated fatty acid synthesis pathway and new strategies that might be used to manipulate the fatty acid content of cyanobacteria.
There is an ample interest in turning cyanobacteria into biosynthetic platforms [1,2,3]. The main advantage of photosynthetic microorganisms is their capacity to economically transform sunlight and CO2 into value-added products without requiring organic compounds as raw materials. Therefore, the cyanobacterial production of fatty acids (FAs) has received great attention in the last few years as a sustainable feedstock for the production of high-energy–density biofuels and other FA-derived oleochemical compounds [4, 5].
The main FA synthesis pathway is highly conserved throughout bacterial species, and even in all kingdoms of life. Its importance comes from the fact that it is crucial for the cells, as it constitutes the first step in the formation of membrane lipids. The lipid composition of bacterial membranes is modified in response to environmental cues, allowing bacteria to survive under unfavorable conditions [6, 7]. In cyanobacteria, FAs are doubly important, since they are constituents of glycoglycerolipids, which form the thylakoid membranes, where photosynthesis takes place . FA composition varies among cyanobacteria. This led Murata et al.  to classify them into four groups, a scheme recently updated to five . Members of Group 1, such as Synechococcus elongatus PCC 7942 (Se7942), encode the Δ9 desaturase gene desC and contain saturated and mono-unsaturated FAs while lack polyunsaturated FAs. A richer FA composition can be found in other cyanobacteria, such as Group 3α Synechococcus sp. PCC 7002 (Ss7002), characterized by the production of trienoic alpha-linolenic acid (ALA, C18:3Δ9,12,15) by Δ9 (DesC), Δ12 (DesA) and Δ15 (DesB) desaturases.
FA synthesis has been widely studied in Escherichia coli, which serves as a general model for FA production in other bacterial species [11, 12]. FA synthesis requires two enzyme complexes, acetyl-CoA carboxylase (AccABCD), which catalyzes the committed step at the beginning of the FA pathway (conversion of acetyl-CoA into malonyl-CoA), and the FA synthase (FAS) system that catalyzes the remaining reactions . Based on the organization of their catalytic units, there are two groups of FAS. Independently of the architecture of the FAS system, FA synthesis involves two steps, initiation and elongation  (Fig. 1a). Type I FAS carries out all steps of FA biosynthesis as a multimeric protein complex and is present in fungi and yeast genomes . Type II FA synthases catalyze these reactions in most bacteria and consist of a series of enzymes, encoded by individual genes  (Fig. 1b). In many bacteria, such as E. coli, genes encoding core enzymes of the FAS system are organized in the fab (FA biosynthesis) cluster  (Fig. 1b, upper panel). On the contrary, in the case of Se7942, fab genes are non-clustered and rather scattered along the chromosome (Fig. 1b, lower panel). Moreover, the organization and distribution of the fab genes along the chromosome are different in both microorganisms. The composition of E. coli and Se7942 fab genes is similar except for the absence of the fabA and fabB in Se7942.
The enzymes involved in the FA synthesis pathway display high amino acid similarity between cyanobacterial and E. coli homologs [17, 18]. In fact, individual subunits of the Ss7002 and E. coli FAS complexes were qualitatively interchangeable in vitro, using purified protein components . Nevertheless, the turnover rate of Ss7002 and E. coli FAS differed in the respective reconstituted systems. The metabolic flux of the Ss7002 FAS system was exclusively limited by FabH, a β-ketoacyl-ACP synthase that initiates FA synthesis by condensing malonyl-ACP with acetyl-CoA to form acetoacetyl-ACP (Fig. 1a) . In contrast, the reconstituted E. coli system was limited by its dehydratase (FabZ) and enoyl reductase (FabI) . Congruently, when the fabZ, fabI and fabG reductase genes were overexpressed in vivo, the resulting E. coli strain produced 50% more FAs . Overproduction of FAs of different chain lengths from C6 to C16 was also obtained in E. coli, based on in silico predictive models of metabolic flux, when the fabD strain was engineered to overproduce the fabZ and acyl-ACP thioesterase genes . Finally, the abundance of FA synthesis proteins was compared between a FA high-yield E. coli strain (fadE deficient) and the wild-type (wt) strain by mass spectrometry . Higher levels of acetyl-CoA carboxylase, thioesterase and FabZ were detected in the FA high-yield strain, implying that they are key enzymes for improving FA production in E. coli.
In this work, we genetically modified the expression of several enzymes involved in the FA synthesis pathway in the model cyanobacterium Se7942 and evaluated their impact on the lipid profile. Our results showed remarkable differences in the functionality of some of these enzymes between Se7942 and E. coli. This analysis allowed the selection of genetic modifications that, combined with the introduction of Ss7002 genes coding for DesA and DesB desaturases, rendered an engineered Se7942 strain able to produce ALA up to 22.6% of its total FA profile. These modifications might be implemented in other cyanobacteria, especially in omega-3 natural producers, to increase their yield.
Overproduction of enzymes involved in the initiation step of fatty acid synthesis pathway changes the lipid profile of Se7942
FabD, the first enzyme of the fab cluster involved in saturated FA synthesis, is a malonyl-CoA ACP transacylase (MCAT). It catalyzes the malonyl-CoA transfer to ACP, generating malonyl-ACP and free CoA  (Figs. 1a, 2a). fabD overexpression in E. coli produced an altered FA composition, increasing the proportion of cis-vaccenic acid (C18:1n7) while decreasing palmitoleic acid (C16:1) . To test if FabD overproduction in Se7942 had the same effect, strain MSM24, containing a Ptrc::fabD cassette in neutral site 1 (NS1), was constructed. As shown in Fig. 2b, the FA profile exhibited small changes in C14:0 and C18:1, but only in the case of C14:0 this change was statistically significant. In E. coli, changes in C14:0 were not observed, while an increase in C18:1 was detected as a consequence of a concomitant decrease in C16:1 . These results demonstrated that fabD overexpression did not produce the same effect in both bacteria.
The next target in the FA synthesis pathway to be modified in this study was the β-ketoacyl-ACP synthase III (KAS III), which catalyzes the formation of 3-ketoacyl-ACP by condensation of acetyl-CoA with malonyl-ACP (Fig. 3a). This enzyme is encoded by the fabH gene and produces the precursor needed for the elongation cycle  (Fig. 1a). Subsequent elongation steps are performed by FabB and FabF, described in more detail below.
In E. coli, fabH deletion improved the production of C18 species, while reducing C16 . Overexpression of fabH had an opposite effect: C14 and C16 increased at the expense of C18:1 . In Ss7002, Kuo and Khosla determined that FabH was the rate-limiting enzyme of FAS using in vitro reconstituted systems and suggested that, either overexpressing the endogenous fabH gene or by replacement with its E. coli ortholog, the FA flux should be increased . Attempts to delete fabH in Ss7002 proved impossible, rendering always non-pure deficient mutants and thus suggesting that fabH is essential in this cyanobacterium .
Based on this information, two different modifications were carried out in Se7942: a mutant that overexpressed fabH under the Ptrc promoter (MSM28), and a merodiploid mutant in which some chromosomal copies of the fabH gene were deleted (MSM23). To generate the latter, we used the gene replacement system designed by Matsuoka et al. . A kanamycin resistance cassette and a wt copy of the rps12 gene of Synechocystis sp. PCC 6803 (Ss6803) were inserted into the Se7942 fabH open reading frame (Fig. 3b, upper panel). After purification of the transformant colonies by successive streakings, the fabH mutant was checked by PCR using primers 111 and 112 (Additional file 3: Table S3). Two amplicon bands were observed in the electrophoresis gel. The largest band (~ 2.6 kb) belonged to the chromosomes with the fabH deletion and the smallest one (~ 2.1 kb) to the wt chromosome (Fig. 3b, lower panel). Using a transposon mutagenesis approach, fabH was described as an essential gene in Se7942 , a probable reason for obtaining here only merodiploid cells that contained both inactivated and fabH wt copies. On the other hand, we were able to delete this gene in Se7942 only when an extra copy of the fabH gene was previously inserted under the control of a strong constitutive promoter, Ptrc (strain MSM34, Additional file 1: Figures S1 and S2B). Lower expression of the extra copy of fabH, driven by a weaker promoter PnrsB with a theophylline-dependent riboswitch  in the absence of inducers, did not allow complete deletion of the wt gene (strain MSM35, Additional file 1: Figure S1 and S2C).
The fabH mutants were used to compare their FA composition with the wt strain (Fig. 3c). Lauric acid (C12:0) was detected only in MSM28, the strain that overexpresses fabH. FabH overproduction resulted also in an increase in medium-chain FAs (MCFAs) (C12:0 and C14:0) and a decrease in long-chain FAs (LCFAs) (C18:0 and C18:1), as happened in E. coli . In contrast, the merodiploid mutant, MSM23, showed higher production of LCFAs (C18 and C18:1) and lower yield of C16:1, also similar to the FA profile described for the E. coli fabH derivative .
FabF overproduction increases transformation of C14 and C16 fatty acids into C18
FabB and FabF, both independently catalyze the first reaction in the second elongation cycle (Fig. 4a). FabB is a KAS I, which catalyzes the condensation of cis-3-decenoyl-ACP, cis-5-dodecenoyl-ACP, and cis-7-tetradodecenoyl-ACP with malonyl-ACP . In E. coli, FabB showed activity with saturated C4–C14 FAs as substrates and was involved in the control of unsaturation . Furthermore, fabB overexpression increased the amount of C18:1 . The fabB gene is present in Ss7002 but no homolog is found in Se7942. Kuo and Khosla determined that fabB was not essential in Ss7002 and speculated that it could play an additional role in fast-growing Ss7002 . Based on these results, we tried to identify the function of this gene by analyzing the effects of its overexpression in Se7942.
To overproduce FabB in Se7942, the Ss7002 fabB gene was cloned under PnrsB and Ptrc promoters, generating strains MSM19 and MSM27, respectively (see Additional file 2 for details). As shown in Fig. 4b, FabB overproduction did not lead to major changes in the FA profile. The most significant variation was the increase in C14:0. This fact suggests that FabB has more affinity for MCFAs, as happens in E. coli . On the other hand, in MSM27 the amount of C16:0 slightly diminished (Fig. 4b, purple bar). This effect was not observed in MSM19, which carries a weaker promoter to produce FabB than MSMS27, and consequently, it could be related with the amount of FabB protein obtained during overexpression (Additional file 1: Figures S2A and S2B, respectively).
FabF also catalyzes the reaction shown in Fig. 4a. FabF is a KAS II and catalyzes the same condensation reaction as KAS I enzymes. In E. coli, KAS II is predominantly responsible for the elongation of palmitoleic acid (C16:1) required for the cis-vaccenic acid (C18:1n-7) synthesis . In E. coli, fabF overexpression is lethal, as high levels of FabF block the access of FabB to the product generated by FabD . Due to the lack of fabB in Se7942, this reaction is only catalyzed by FabF, and thus the blocking effect observed in E. coli should not take place. To verify this end, we constructed Se7942 mutant MSM29, which overexpresses fabF. This mutant showed a significant decrease in MCFAs (C14:0, C16:0 and C16:1) linked to a significant increase in LCFAs (C18:0 and C18:1) (Fig. 5). These results indicate that high FabF levels improved the transformation of C14 and C16 into C18.
FadD deletion increases transformation of C16:1 into C18
The metabolic pathway to degrade MCFAs and LCFAs is called β-oxidation. Despite being highly conserved, it is not present in cyanobacteria, according to the literature and genome annotation. The first step in this pathway is catalyzed by FadD, a long-chain fatty acyl-CoA ligase (acyl-CoA synthetase) responsible for the activation of exogenous LCFAs into acyl-CoAs  (Fig. 6a). FadD is the single gene of this pathway identified in cyanobacteria . In E. coli, fadD disruption caused the accumulation of free FAs in the cytosol [35, 37]. In Se7942, fadD deletion produced a change in the FA profile of membrane lipids, particularly an increase in the degree of saturation .
Based on this information, we disrupted the fadD gene in Se7942, obtaining strain MSM42 (Additional file 1: Figure S3A), which is a pure fadD mutant (Additional file 1: Figure S3B). The lipid profiles of the engineered and wt strains were analyzed and compared. As shown in Fig. 6b, MSM42 exhibited a significant increase in C18 FAs (C18:0 and C18:1) at the expense of C16:1 content decrease.
Overexpression of desA and desB desaturase genes leads to increase alpha-linolenic acid production at low temperature
Three desaturases are essential to produce ALA in cyanobacteria, DesC, DesA and DesB (Fig. 7a). DesC introduces the first cis double bond in the Δ9 position of stearic acid (C18:0), producing oleic acid (C18:1). The next double bond is introduced at the Δ12 position of C18:1 by DesA, producing linoleic acid (C18:2). Finally, DesB introduces a third double bond at the Δ15 (ω3) position of C18:2, rendering ALA. DesC and DesB may also use C16 FAs as substrates [39, 40]. Se7942 encodes desC but lacks both desA and desB desaturase genes , whose role in ALA production is well documented in Ss7002 . The desA mRNA accumulates faster in Ss7002  than in other cyanobacteria, such as Ss6803 . Thus, desA and desB from Ss7002 were chosen as candidate genes to be introduced in Se7942 for producing ALA in this organism.
Two Se7942 derivative strains were constructed by integrating Ss7002 desA and desB genes in the NS1 site: MSM17 (PnrsB::desA-desB) and MSM26 (Ptrc::desA-desB). Optimal conditions to express these genes were previously defined in Ss7002. Both desA and desB transcripts were more abundant at low temperature, specifically at 22 °C [44, 45]. Therefore, ALA production in Se7942 was tested at this temperature. Both mutant strains produced ALA (Fig. 7b). In the case of MSM17, which expresses desAB genes from the PnrsB promoter, ALA production (1.14% of total lipids) was significantly lower than that of Ss7002. Replacement of PnrsB by a stronger promoter, Ptrc (strain MSM26), led to an increase of approximately 6.6 times in ALA production yield, reaching levels as high as those obtained from Ss7002.
Based on data recently published by Yang et al.  in the Cyanomics database, a condition that might improve the expression of both desaturases, i.e., high light intensity, was identified. This database comprises genomic, transcriptomic and proteomic datasets of Ss7002. This condition had been already tested by Sakamoto et al., who found an increase of 3.8 times in ALA production by Ss7002 under high light intensity (250 µmol photons m−2 s−1) . The influence of light intensity on ALA production was thus investigated in the engineered Se7942 strain MSM26. No statistically significant difference was observed under low light intensity (60 µmol photons m−2 s−1), 7.61 ± 1.45%, versus 8.78 ± 1.79% under high light intensity (250 µmol photons m−2 s−1).
Overproduction of FabF and DesAB in a Se7942 fadD deletion mutant improves alpha-linolenic acid yields
By construction of a series of Se7942 mutants, we observed that both fabF overexpression and fadD deletion improved the ALA production, whereas fabH partial deletion did not (data not shown). The combined effect of both modifications was checked in the mutant strain MSM45, which lacked fadD and overexpressed fabF and both desaturase genes (desA and desB). We compared the FA composition of MSM45 versus wt Se7942 and several differences were observed (Fig. 8). First, a significant reduction in C14:0, C16:1 and C18:1 was produced in MSM45, as well as a significant increase in C16:0 and C18:0. Strain MSM45 produced ALA to as high levels as 22.6%. This amount was higher than the ALA production by Ss7002 under high light condition (19%) .
Se7942 is a model cyanobacterium because of its small genome size and simple metabolism, making it ideal for the study of lipid biosynthesis in photoautotrophic bacteria. To gain a detailed understanding of the FA synthesis pathway in Se7942, a deeper knowledge on the activity of the enzymes involved was needed. To fulfill this aim, we engineered several modifications in the production of these enzymes in Se7942. Our results showed important differences in the activity of some enzymes involved in FA synthesis in Se7942 regarding those in E. coli, despite their amino acid sequence homology (ranging from 42 to 52% identity). Whereas fabD overexpression altered the FA composition of E. coli , this modification did not have a remarkable effect in the Se7942 FA profile. On the other hand, the E. coli fabH mutant is viable and exhibits a small-colony phenotype , suggesting the implication of another condensing enzyme performing less efficiently this reaction in the initiation cycle of the FA biosynthesis. In Se7942, the fabH deletion was lethal, when assayed by random barcode transposon site sequencing (RB-TnSeq) . Neither was this gene removed from all chromosomal copies of Ss7002 . Our attempts to delete fabH in Se7942 also resulted in merodiploid mutants, an indication of the indispensability of this gene, and that no other Se7942 KAS enzyme (e.g., FabF) can act as a functional replacement for FabH. Overexpression of fabH led to increase in MCFA production, specifically lauric (C12:0) and myristic (C14:0) acids. Notably, C12:0 is not naturally present in Se7942. Increased levels of MCFAs were reported in an engineered Ss7002 strain in which fabH was replaced by Chaetoceros GSL56 KASIII gene . MCFAs are appreciated in the fuel industry, which makes fabH overexpression in cyanobacteria potentially useful. These results suggest that in cyanobacteria FabH may have more affinity toward short-chain FAs, a feature described previously in other bacteria .
Regarding the two other KAS, FabB and FabF, relevant differences to E. coli were found. In E. coli, FabB plays an important role in controlling the unsaturation degree of FAs . It is noteworthy that FabB is absent in Se7942 and most cyanobacteria, but present in Ss7002, where it is dispensable in unsaturated FA synthesis . In cyanobacteria, unsaturated FA synthesis is controlled downstream of the FA synthases, because their synthesis is performed by desaturases, which act on FAS products . The growth temperature of cyanobacteria is responsible for the up- or downregulation of the desaturase gene expression . Overproduction of FabB from Ss7002 in Se7942 did not increase C18:1 yield, contrary to the E. coli case , suggesting that FabB did not play a role in the cyanobacterial unsaturated FA synthesis. Kuo and Khosla speculated that FabB could have an additional role, still unknown . Overexpression of fabF was lethal in E. coli, likely due to the fact that high levels of FabF block the access of FabB to the product generated by FabD . This lethality was not observed when the endogenous fabF gene was overexpressed in Se7942, ruled out any FabB blockage by FabF since FabB is absent in Se7942.
Taking all these results together, we suggest that the enzymes involved in FA synthesis in Se7942 play a different role from that observed in E. coli. Moreover, it also suggests that these enzymes could be regulated in a different way in both bacteria. Even though the biochemistry of FA synthesis in E. coli is well known and documented, its genetic regulation still remains unclear. In E. coli, the key regulators are the long-chain acyl-ACP end products, which exert a negative regulatory feedback on key enzymes of the synthesis pathway . Moreover, this pathway is also controlled by transcription factors, mainly FabR (TetR family) and FadR (GntR family), among others . FabR and FadR also control the expression of FA degradation genes involved in the β-oxidation cycle . As mentioned above, it appears that cyanobacteria lack the major FA degradative metabolic pathway β-oxidation, generally thought to be universal. In several Gram-positive bacteria, there are other families of regulators that control the expression of fab genes. In Bacillus subtilis, FapR is a global negative regulator of genes involved in FA and phospholipid biosynthesis . This regulator, a member of the DeoR family, is common to other species of the Bacillus genus, as well as Listeria, Clostridium and Staphylococcus  and senses malonyl-CoA, which releases its transcriptional repression [52, 53]. In Streptococcus pneumoniae, a different FA synthesis regulator was identified, FabT, which belongs to the MarR family . FabT homologs were also located in other groups of Gram‐positive bacteria, such as Enterococcus, Clostridium and Lactococcus .
Little is known about the regulation of the gene expression of FA synthesis in cyanobacteria. In most bacteria, fab genes are organized into operons, and this organization enables the identification of cognate transcription factors based on shared synteny . In cyanobacteria , fab genes are scattered throughout the genome. Moreover, due to the fact that cyanobacteria do not code for a FA degradation pathway, some authors suggested that the regulatory of mechanisms for FA metabolism may be largely different from those of heterotrophic bacteria [56, 57]. There is scarce information on the transcriptional regulation of FA biosynthesis in cyanobacteria. Deletion of a cyanobacteria-specific transcriptional regulator gene, cyAbrB2, was found to enhance the production of free FAs in engineered strains of Ss6803 . Besides, the LexA transcription factor was described as a regulator of FA synthesis also in Ss6803, and more specifically as a repressor of genes involved in the initiation step of FA synthesis (fabD, fabH and fabF) and the first reductive step in the elongation cycle (fabG) . Gene lexA deletion largely increased the production of FAs in an engineered Ss6803 strain . No LexA homolog was found in Se7942 . The information related to FA synthesis regulation in Se7942 focuses on the signal transduction protein PII. The PII mutant strain of Se7942 showed a decreased intracellular acetyl-CoA content and enhanced activity of the acetyl-CoA–carboxylase complex through transcriptional activation of the accABCD genes . However, PII does not regulate the enzymes involved in the FAS system. These and our results may suggest that the regulation of FA synthesis enzymes in cyanobacteria is independent of FA degradation by the β-oxidation cycle and could involve a global regulator, which has yet to be found.
Omega-3 FAs are naturally produced in some cyanobacteria, members of the groups 3α and 4 , which encode Δ12 and Δ15 desaturase genes. Such is the case of Ss6803 , Synechococcus sp. NKBG 15041c  and Ss7002 . The first two produced high levels of ALA (23 and 50% of total FA, respectively) when engineered to overexpress, respectively, Δ6 + Δ15 and Δ6 FA desaturase genes [63, 64]. On the other hand, Ss7002 produced the highest ALA levels (19% of total FA) in a cyanobacterial wt strain when grown at high light intensity either at 22 °C or at 38 °C and shifted to 22 °C for 12 h . The introduction of the Δ12 desaturase gene desA from Ss6803 into Se7942 caused a modification of its FA composition: dienoic FA C16:2Δ9,12 and C18:2Δ9,12 emerged at significant levels at the expense of C16:1Δ9 and C18:1Δ9 . This modification led to an enhancing in chilling tolerance  and protection against photoinhibition . When transformed with Ss6803 Δ12 and Δ15 desaturase genes, desA and desB, Se7942 was able to desaturase FA at Δ9, Δ12 and Δ15 positions . ALA production by the transformant strain at 34 °C was low (1% of total lipids) but significantly increased (up to 5%) when incubated at 22 °C. Taking advantage of the high similarity in codon usage between Se7942 and Ss7002, the Δ12 and Δ15 desaturase genes of the latter, desA and desB, were chosen to produce transgenically ALA in Se7942, which lacks these two essential enzymes of the desaturation pathway. The desB gene of Ss7002 is transcribed only at low temperatures (22 °C) [44, 45]. Thus, an inducible promoter was chosen to drive the desAB expression at 22 °C once enough biomass was achieved at 30 °C, because Se7942 growth is impaired at low temperature . The resulting strain, MSM26, produced ALA at the same level of Ss7002 (Fig. 7b), at the expense of C18 and C16:1 FAs. Se7942 mutants either overexpressing fabF or lacking fadD exhibited a FA profile enriched in C18 and thus these mutations were evaluated in combination with desAB expression. This strategy resulted in a significant improvement in the yield of this omega-3 FA in relation to the available data for Ss7002 [43, 45, 68, 69], and comparable to the engineered Se6803 .
In this study, the activity of the enzymes involved in the first reactions of FA synthesis pathway in Se7942 was examined by deleting or overexpressing their respective coding genes and analyzing the resulting FA profiles. In general, the mutant strains showed changes in their FA composition, regarding both the wt Se7942 strain and E. coli mutants. In addition, the modifications that improved the synthesis of C16:0 and C18:0 (overexpression of fabF and deletion of fadD) used in combination with the expression of Ss7002 desAB desaturase genes was revealed as a feasible strategy to engineer Se7942 for ALA production. To the best of our knowledge, this is the first report on increased ALA production in cyanobacteria using modifications in the expression of the genes belonging to the FAS system. The combination of fadD deletion and fabF overexpression is a strategy that may be evaluated in omega-3 natural producers to raise ALA yield.
Strains and culture conditions
All strains used in this study are listed in Table 1. Ss7002 was cultured in A + medium , while Se7942 and GRPS1 were cultured in BG11 medium . All cyanobacterial strains were grown at 30 °C by bubbling 1% CO2 with continuous light at 60 µmol photons m−2 s−1. Antibiotics used for selecting Se7942 were neomycin at 5 or 25 µg/ml (Neo5 or Neo25), spectinomycin at 10 or 20 µg/ml (Sp10 or Sp20) and chloramphenicol at 5 or 10 µg/ml (Cm5 or Cm10).
Escherichia coli strain DH5α was used for transformation with recombinant plasmids. It was grown in Luria–Bertani (LB) medium at 37 °C under shaking (150 rpm) and, when required, was supplemented with antibiotics: kanamycin, 50 µg/ml (Km50), chloramphenicol, 25 µg/ml (Cm25), spectinomycin, 100 µg/ml (Sp100) and ampicillin, 100 µg/ml (Ap100).
Natural transformation of Se7942
Transformation of Se7942 was performed as previously described . Briefly, a Se7942 culture sample equivalent to 10 μg of chlorophyll (around 4 × 109 cyanobacterial cells) was mixed with 500 ng plasmid DNA for 24 h in the dark at 30 °C. Transformation mixtures were deposited onto nitrocellulose filters (Millipore) and incubated for 24 h on BG11 plates at 30 °C with continuous light. Transformant colonies were selected in BG11 supplemented with the corresponding antibiotic. Mutant segregation was achieved by repeatedly transferring individual transformant colonies to fresh selective plates. Mutant genotypes were confirmed by PCR followed by DNA sequencing using specific primers (see Additional file 3: Tables S1 and S3).
Cultivation of cyanobacteria and engineered strains for fatty acid production
Se7942 was cultured in BG11 medium in a Multi-cultivator MC 1000-OD (Photon Systems Instruments, Drasov, Czech Republic) with cool white light (60 µmol photons m−2 s−1) at 30 °C and 3% CO2. From solid BG11 agar medium, individual colonies were inoculated into 60 ml of BG11 and grown for 5–7 days to be used as inocula. Cultures were inoculated at initial OD720 = 0.05 in 60 ml of fresh BG11 medium (without antibiotics). At OD720 0.5–0.6, the temperature was decreased from 30 °C to 22 °C, IPTG (Ptrc promoter) or NiSO4 (PnrsB promoter) was added to 1 mM or 5 µM final concentrations, respectively. Samples were taken 24 h post-induction.
Chromosomal gene deletion system based on a dominant streptomycin-sensitive rps12 gene
Gene deletion was carried out in the Se7942 derivative strain GRPS1 , using the rps12-mediated gene replacement method . In this approach, a transformation is performed using a plasmid that contains rps12 and nptI genes flanked by USDT (upstream sequences of the deletion target) and DSDT (downstream sequences of the deletion target). Transformant colonies were selected using Neo25. Individual colonies were checked for the Neor/Sms phenotype and the absence of the deleted region was assayed by PCR with primers that hybridized into the flanking regions of the deletion target (30 cycles of 94 °C during 60 s, 50 °C 30 s, and 72 °C 30 s).
Promoter activity measurements: determination of luciferase activity
Promoters used to overexpress genes in Se7942 were tested by placing them upstream promoter-less reporter genes luxAB from Photorhabdus luminescens . The different constructions were introduced into NS1 of Se7942 to test for luciferase activity. Bioluminescence measurements were carried out essentially as described previously . 1.0 ml aliquots were taken at selected time points from shaking cultures and the final OD720 was adjusted to 0.5. 100 µl samples were transferred in triplicates to white 96-well microtiter plates (Thermo Fisher Scientific). 100 mM decanal stock solution was prepared in methanol and freshly diluted with BG11 for 2 mM decanal ready-to-use solution prior measuring. 100 µl of 2 mM decanal ready-to-use solution was added to the samples, and the plate was immediately placed into a plate reader (Wallac Victor 2 1420 multilabel counter, PerkinElmer Life Sciences). Bioluminescence was measured for 30 min at 25 °C. The maximum light emission (generally around 10 min after start measuring) was used as the bioluminescence value and normalized to the OD720. Results were presented in relative bioluminescence units. Experiments were carried out in triplicates.
Fatty acid analysis
Total FA content of cyanobacterial cultures was determined by gas chromatography flame ionization detection (GC-FID) analysis. Cultured cells (60 ml) were harvested by centrifugation 24 h post-induction. Preparation of FA methyl esters (FAMEs) was done according to . Briefly, saponification was conducted by adding to each sample 1 ml of saponification reagent (45 g NaOH, 150 ml methanol, 150 ml deionized H20). Samples were vortexed for 10 s, heated 5 min at 100 °C, vortexed 10 s, and heated again, for 25 min. Tubes were then rapidly cooled to room temperature. Methylation was accomplished by adding 2 ml of methylation reagent (325 ml 6 N HCl, 275 ml methanol), vortexing for 10 s, and heating for 10 min at 80 °C. After cooling to room temperature, FAMEs were removed from the acidic aqueous phase and transferred to an organic phase by adding 1.25 ml extraction solvent (hexane:methyl tert-butyl ether [MTBE] 1:1), mixing end-over-end in a laboratory rotator for 10 min, and removing the lower aqueous phase. The organic phase was washed in 3 ml base wash (10.8 g NaOH, 900 ml distilled water) for 5 min with end-over-end mixing. Two-thirds of the upper solvent phase were removed for FAME analysis. Samples were stored at − 20 °C. FAME analysis by GC–FID employed an Agilent 7890A gas chromatograph equipped with a capillary column DB-23 (Agilent Technologies, Santa Clara, CA, USA). The initial oven temperature was set at 130 °C and programmed up to 215 °C at a rate of 2.75 °C/min, and up to 230 °C at a rate of 40 °C/min. The injector temperature was set at 270 °C, and the FID detector at 280 °C. FAMEs were identified by comparing the retention time with Supelco® 37 component FAME standard mix (Sigma-Aldrich, St. Louis, MO, USA).
All measurements were carried out at least in triplicates. Data were expressed as the mean ± SD or mean + SD of n separate experiments, depending on data representation. The significance of differences between groups was evaluated by one-way ANOVA followed by Dunnett’s multiple comparison test or by an unpaired Student’s t test. Statistical analyses were carried out using GraphPad Prism software (GraphPad Software Inc., San Diego, CA). Differences were considered statistically significant at p < 0.05.
Synechococcus elongatus PCC 7942
Synechococcus sp. PCC 7002
fatty acid synthase
- fab :
fatty acid biosynthesis
malonyl-CoA ACP transacylase
neutral site 1
Synechocystis sp. PCC 6803
medium-chain fatty acid
long-chain fatty acid
random barcode transposon site sequencing
upstream sequences of the deletion target
downstream sequences of the deletion target
gas chromatography–flame ionization detection
fatty acid methyl ester
analysis of variance
neutral site 2
Knoot CJ, Ungerer J, Wangikar PP, Pakrasi HB. Cyanobacteria: promising biocatalysts for sustainable chemical production. J Biol Chem. 2018;293(14):5044–52.
Lau NS, Matsui M, Abdullah AA. Cyanobacteria: photoautotrophic microbial factories for the sustainable synthesis of industrial products. Biomed Res Int. 2015;2015:754934.
Nozzi NE, Oliver JW, Atsumi S. Cyanobacteria as a platform for biofuel production. Front Bioeng Biotechnol. 2013;1:7.
Ducat DC, Way JC, Silver PA. Engineering cyanobacteria to generate high-value products. Trends Biotechnol. 2011;29(2):95–103.
Wijffels RH, Kruse O, Hellingwerf KJ. Potential of industrial biotechnology with cyanobacteria and eukaryotic microalgae. Curr Opin Biotechnol. 2013;24(3):405–13.
Koga Y. Thermal adaptation of the archaeal and bacterial lipid membranes. Archaea. 2012;2012:789652.
Roy H, Dare K, Ibba M. Adaptation of the bacterial membrane to changing environments using aminoacylated phospholipids. Mol Microbiol. 2009;71(3):547–50.
Holzl G, Dormann P. Structure and function of glycoglycerolipids in plants and bacteria. Prog Lipid Res. 2007;46(5):225–43.
Murata N, Wada H, Gombos Z. Modes of fatty-acid desaturation in cyanobacteria. Plant Cell Physiol. 1992;33(7):933–41.
Los DA, Mironov KS. Modes of fatty acid desaturation in cyanobacteria: an update. Life (Basel). 2015;5(1):554–67.
Janssen HJ, Steinbuchel A. Fatty acid synthesis in Escherichia coli and its applications towards the production of fatty acid based biofuels. Biotechnol Biofuels. 2014;7(1):7.
Handke P, Lynch SA, Gill RT. Application and engineering of fatty acid biosynthesis in Escherichia coli for advanced fuels and chemicals. Metab Eng. 2011;13(1):28–37.
Chan DI, Vogel HJ. Current understanding of fatty acid biosynthesis and the acyl carrier protein. Biochem J. 2010;430(1):1–19.
Schweizer E, Hofmann J. Microbial type I fatty acid synthases (FAS): major players in a network of cellular FAS systems. Microbiol Mol Biol Rev. 2004;68(3):501–17.
White SW, Zheng J, Zhang YM, Rock. The structural biology of type II fatty acid biosynthesis. Annu Rev Biochem. 2005;74:791–831.
Zhang Y, Cronan JE Jr. Polar allele duplication for transcriptional analysis of consecutive essential genes: application to a cluster of Escherichia coli fatty acid biosynthetic genes. J Bacteriol. 1996;178(12):3614–20.
Kuo J, Khosla C. The initiation ketosynthase (FabH) is the sole rate-limiting enzyme of the fatty acid synthase of Synechococcus sp. PCC 7002. Metab Eng. 2014;22:53–9.
Liu Y, Zhang Y, Cao X, Xue S. Cloning, purification, crystallization and preliminary X-ray crystallographic analysis of MCAT from Synechocystis sp. PCC 6803. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2013;69(Pt 11):1256–9.
Yu X, Liu T, Zhu F, Khosla C. In vitro reconstitution and steady-state analysis of the fatty acid synthase from Escherichia coli. Proc Natl Acad Sci USA. 2011;108(46):18643–8.
Ranganathan S, Tee TW, Chowdhury A, Zomorrodi AR, Yoon JM, Fu Y, Shanks JV, Maranas CD. An integrated computational and experimental study for overproducing fatty acids in Escherichia coli. Metab Eng. 2012;14(6):687–704.
Tao H, Zhang Y, Cao X, Deng Z, Liu T. Absolute quantification of proteins in the fatty acid biosynthetic pathway using protein standard absolute quantification. Synth Syst Biotechnol. 2016;1(3):150–7.
Molnos J, Gardiner R, Dale GE, Lange R. A continuous coupled enzyme assay for bacterial malonyl-CoA:acyl carrier protein transacylase (FabD). Anal Biochem. 2003;319(1):171–6.
Magnuson K, Oh W, Larson TJ, Cronan JE Jr. Cloning and nucleotide sequence of the fabD gene encoding malonyl coenzyme A-acyl carrier protein transacylase of Escherichia coli. FEBS Lett. 1992;299(3):262–6.
Gajiwala KS, Margosiak S, Lu J, Cortez J, Su Y, Nie Z, Appelt K. Crystal structures of bacterial FabH suggest a molecular basis for the substrate specificity of the enzyme. FEBS Lett. 2009;583(17):2939–46.
Yao Z, Davis RM, Kishony R, Kahne D, Ruiz N. Regulation of cell size in response to nutrient availability by fatty acid biosynthesis in Escherichia coli. Proc Natl Acad Sci USA. 2012;109(38):E2561–8.
Tsay JT, Oh W, Larson TJ, Jackowski S, Rock CO. Isolation and characterization of the beta-ketoacyl-acyl carrier protein synthase III gene (fabH) from Escherichia coli K-12. J Biol Chem. 1992;267(10):6807–14.
Gu H, Jinkerson RE, Davies FK, Sisson LA, Schneider PE, Posewitz MC. Modulation of medium-chain fatty acid synthesis in Synechococcus sp. PCC 7002 by replacing FabH with a chaetoceros ketoacyl-ACP synthase. Front Plant Sci. 2016;7:690.
Matsuoka M, Takahama K, Ogawa T. Gene replacement in cyanobacteria mediated by a dominant streptomycin-sensitive rps12 gene that allows selection of mutants free from drug resistance markers. Microbiology. 2001;147(Pt 8):2077–87.
Rubin BE, Wetmore KM, Price MN, Diamond S, Shultzaberger RK, Lowe LC, Curtin G, Arkin AP, Deutschbauer A, Golden SS. The essential gene set of a photosynthetic organism. Proc Natl Acad Sci USA. 2015;112(48):E6634–43.
Nakahira Y, Ogawa A, Asano H, Oyama T, Tozawa Y. Theophylline-dependent riboswitch as a novel genetic tool for strict regulation of protein expression in Cyanobacterium Synechococcus elongatus PCC 7942. Plant Cell Physiol. 2013;54(10):1724–35.
Feng Y, Cronan JE. Escherichia coli unsaturated fatty acid synthesis: complex transcription of the fabA gene and in vivo identification of the essential reaction catalyzed by FabB. J Biol Chem. 2009;284(43):29526–35.
Edwards P, Nelsen JS, Metz JG, Dehesh K. Cloning of the fabF gene in an expression vector and in vitro characterization of recombinant fabF and fabB encoded enzymes from Escherichia coli. FEBS Lett. 1997;402(1):62–6.
Clark DP, DeMendoza D, Polacco ML, Cronan JE. β-Hydroxydecanoyl thioester dehydrase does not catalyze a rate-limiting step in Escherichia coli unsaturated fatty acid synthesis. Biochemistry. 1983;22(25):5897–902.
Subrahmanyam S, Cronan JE Jr. Overproduction of a functional fatty acid biosynthetic enzyme blocks fatty acid synthesis in Escherichia coli. J Bacteriol. 1998;180(17):4596–602.
Pech-Canul A, Nogales J, Miranda-Molina A, Alvarez L, Geiger O, Soto MJ, Lopez-Lara IM. FadD is required for utilization of endogenous fatty acids released from membrane lipids. J Bacteriol. 2011;193(22):6295–304.
Kaczmarzyk D, Fulda M. Fatty acid activation in cyanobacteria mediated by acyl-acyl carrier protein synthetase enables fatty acid recycling. Plant Physiol. 2010;152(3):1598–610.
Xu P, Gu Q, Wang W, Wong L, Bower AG, Collins CH, Koffas MA. Modular optimization of multi-gene pathways for fatty acids production in E. coli. Nat Commun. 2013;4:1409.
Ruffing AM, Jones HD. Physiological effects of free fatty acid production in genetically engineered Synechococcus elongatus PCC 7942. Biotechnol Bioeng. 2012;109(9):2190–9.
Los DA, Murata N. Membrane fluidity and its roles in the perception of environmental signals. Biochim Biophys Acta. 2004;1666(1–2):142–57.
Tocher DR, Leaver MJ, Hodgson PA. Recent advances in the biochemistry and molecular biology of fatty acyl desaturases. Prog Lipid Res. 1998;37(2–3):73–117.
Wada H, Murata N. Temperature-induced changes in the fatty acid composition of the Cyanobacterium, Synechocystis PCC6803. Plant Physiol. 1990;92(4):1062–9.
Ludwig M, Bryant DA. Transcription Profiling of the Model Cyanobacterium Synechococcus sp. Strain PCC 7002 by Next-Gen (SOLiD) Sequencing of cDNA. Front Microbiol. 2011;2:41.
Sakamoto T, Higashi S, Wada H, Murata N, Bryant DA. Low-temperature-induced desaturation of fatty acids and expression of desaturase genes in the cyanobacterium Synechococcus sp. PCC 7002. FEMS Microbiol Lett. 1997;152(2):313–20.
Ludwig M, Bryant DA. Synechococcus sp. strain PCC 7002 transcriptome: acclimation to temperature, salinity, oxidative stress, and mixotrophic growth conditions. Front Microbiol. 2012;3:354.
Sakamoto T, Bryant DA. Temperature-regulated mRNA accumulation and stabilization for fatty acid desaturase genes in the cyanobacterium Synechococcus sp. strain PCC 7002. Mol Microbiol. 1997;23(6):1281–92.
Yang Y, Feng J, Li T, Ge F, Zhao J. CyanOmics: an integrated database of omics for the model cyanobacterium Synechococcus sp. PCC 7002. Database (Oxford). 2015. https://doi.org/10.1093/database/bau127.
Khandekar SS, Gentry DR, Van Aller GS, Warren P, Xiang H, Silverman C, Doyle ML, Chambers PA, Konstantinidis AK, Brandt M, et al. Identification, substrate specificity, and inhibition of the Streptococcus pneumoniae beta-ketoacyl-acyl carrier protein synthase III (FabH). J Biol Chem. 2001;276(32):30024–30.
Murata N, Wada H. Acyl-lipid desaturases and their importance in the tolerance and acclimatization to cold of cyanobacteria. Biochem J. 1995;308(Pt 1):1–8.
My L, Rekoske B, Lemke JJ, Viala JP, Gourse RL, Bouveret E. Transcription of the Escherichia coli fatty acid synthesis operon fabHDG is directly activated by FadR and inhibited by ppGpp. J Bacteriol. 2013;195(16):3784–95.
Zhu K, Zhang YM, Rock CO. Transcriptional regulation of membrane lipid homeostasis in Escherichia coli. J Biol Chem. 2009;284(50):34880–8.
Schujman GE, Paoletti L, Grossman AD, de Mendoza D. FapR, a bacterial transcription factor involved in global regulation of membrane lipid biosynthesis. Dev Cell. 2003;4(5):663–72.
Albanesi D, Reh G, Guerin ME, Schaeffer F, Debarbouille M, Buschiazzo A, Schujman GE, de Mendoza D, Alzari PM. Structural basis for feed-forward transcriptional regulation of membrane lipid homeostasis in Staphylococcus aureus. PLoS Pathog. 2013;9(1):e1003108.
Schujman GE, Guerin M, Buschiazzo A, Schaeffer F, Llarrull LI, Reh G, Vila AJ, Alzari PM, de Mendoza D. Structural basis of lipid biosynthesis regulation in Gram-positive bacteria. EMBO J. 2006;25(17):4074–83.
Lu YJ, Rock CO. Transcriptional regulation of fatty acid biosynthesis in Streptococcus pneumoniae. Mol Microbiol. 2006;59(2):551–66.
Eckhardt TH, Skotnicka D, Kok J, Kuipers OP. Transcriptional regulation of fatty acid biosynthesis in Lactococcus lactis. J Bacteriol. 2013;195(5):1081–9.
Kizawa A, Kawahara A, Takashima K, Takimura Y, Nishiyama Y, Hihara Y. The LexA transcription factor regulates fatty acid biosynthetic genes in the cyanobacterium Synechocystis sp. PCC 6803. Plant J. 2017;92(2):189–98.
Ruffing AM. Metabolic engineering and systems biology for free fatty acid production in Cyanobacteria. In: Los DA, editor. Cyanobacteria: omics and manipulation. Norfolk: Caister Academic Press; 2017. p. 161–86.
Kawahara A, Sato Y, Saito Y, Kaneko Y, Takimura Y, Hagihara H, Hihara Y. Free fatty acid production in the cyanobacterium Synechocystis sp. PCC 6803 is enhanced by deletion of the cyAbrB2 transcriptional regulator. J Biotechnol. 2016;220:1–11.
Cassier-Chauvat C, Veaudor T, Chauvat F. Comparative genomics of DNA recombination and repair in Cyanobacteria: biotechnological implications. Front Microbiol. 2016;7:1809.
Verma E, Chakraborty S, Tiwari B, Mishra AK. Transcriptional regulation of acetyl CoA and lipid synthesis by PII protein in Synechococcus PCC 7942. J Basic Microbiol. 2018;58(2):187–97.
Wada H, Murata N. Synechocystis PCC6803 mutants defective in desaturation of fatty acids. Plant Cell Physiol. 1989;30(7):971–8.
Yu R, Yamada A, Watanabe K, Yazawa K, Takeyama H, Matsunaga T, Kurane R. Production of eicosapentaenoic acid by a recombinant marine cyanobacterium, Synechococcus sp. Lipids. 2000;35(10):1061–4.
Chen G, Qu S, Wang Q, Bian F, Peng Z, Zhang Y, Ge H, Yu J, Xuan N, Bi Y, He Q. Transgenic expression of delta-6 and delta-15 fatty acid desaturases enhances omega-3 polyunsaturated fatty acid accumulation in Synechocystis sp. PCC6803. Biotechnol Biofuels. 2014;7(1):32.
Yoshino T, Kakunaka N, Liang Y, Ito Y, Maeda Y, Nomaguchi T, Matsunaga T, Tanaka T. Production of omega3 fatty acids in marine cyanobacterium Synechococcus sp. strain NKBG 15041c via genetic engineering. Appl Microbiol Biotechnol. 2017;101(18):6899–905.
Wada H, Gombos Z, Murata N. Enhancement of chilling tolerance of a cyanobacterium by genetic manipulation of fatty acid desaturation. Nature. 1990;347(6289):200–3.
Gombos Z, Kanervo E, Tsvetkova N, Sakamoto T, Aro EM, Murata N. Genetic enhancement of the ability to tolerate photoinhibition by introduction of unsaturated bonds into membrane glycerolipids. Plant Physiol. 1997;115(2):551–9.
Sakamoto T, Los DA, Higashi S, Wada H, Nishida I, Ohmori M, Murata N. Cloning of omega 3 desaturase from cyanobacteria and its use in altering the degree of membrane-lipid unsaturation. Plant Mol Biol. 1994;26(1):249–63.
Sakamoto T, Bryant DA. Synergistic effect of high-light and low temperature on cell growth of the Delta12 fatty acid desaturase mutant in Synechococcus sp. PCC 7002. Photosynth Res. 2002;72(3):231–42.
Sakamoto T, Shen G, Higashi S, Murata N, Bryant DA. Alteration of low-temperature susceptibility of the cyanobacterium Synechococcus sp. PCC 7002 by genetic manipulation of membrane lipid unsaturation. Arch Microbiol. 1998;169(1):20–8.
Stevens SE, Patterson COP, Myers J. The production of hydrogen peroxide by blue-green algae: a survey. J Phycol. 1973;9(4):427–30.
Rippka R, Deruelles J, Waterbury JB, Herdman M, Stanier RY. Generic assignments, strain histories and properties of pure cultures of Cyanobacteria. Microbiology. 1979;111(1):1–61.
Golden SS, Brusslan J, Haselkorn R. Genetic engineering of the cyanobacterial chromosome. Methods Enzymol. 1987;153:215–31.
Takahama K, Matsuoka M, Nagahama K, Ogawa T. High-frequency gene replacement in cyanobacteria using a heterologous rps12 gene. Plant Cell Physiol. 2004;45(3):333–9.
Meighen EA. Molecular biology of bacterial bioluminescence. Microbiol Rev. 1991;55(1):123–42.
Eisenhut M, Georg J, Klahn S, Sakurai I, Mustila H, Zhang P, Hess WR, Aro EM. The antisense RNA As1_flv4 in the Cyanobacterium Synechocystis sp. PCC 6803 prevents premature expression of the flv4-2 operon upon shift in inorganic carbon supply. J Biol Chem. 2012;287(40):33153–62.
Miller L, Berger T. Bacteria identification by gas chromatography of whole cell fatty acids. Hewlett-Packard application note. Hewlett-Packard: Avondale; 1985. p. 228–41.
MSM carried out all the experimental procedures of the study, participated in its design and drafted the manuscript. MPGB helped with some of the experiments and drafted the manuscript. FDLC conceived the study, participated in its design and coordination and drafted the manuscript. All authors read and approved the final manuscript.
We thank the Pasteur Culture Collection (Institute Pasteur, Paris, France) for supplying Se7942 and Ss7002, A. Contreras and K. Takahama for kindly sharing pUAGC280 and pAM1580, and GRPS1, respectively. We are also grateful to A. Fernández and F. Romero for GC–FID analysis performed at Biomar Microbial Technologies, León, Spain, and R. Fernández-López, M. Robledo, R. Ruiz, S. Redondo and C. Palencia-Gándara for helpful discussions.
The authors declare that they have no competing interests.
Availability of supporting data
All data generated or analyzed during this study are included in this published article and its Additional files.
Consent for publication
All authors agree to publication.
Ethics approval and consent to participate
Work in the FDLC laboratory was financed by the Spanish Ministry of Economy and Competitivity (MINECO) Grant BFU2014-55534-C2-1-P. MSM. was recipient of a Ph.D. fellowship (BES-2012-057387) from MINECO.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.