Effect of cerulenin on fatty acid composition and gene expression pattern of DHA-producing strain Colwellia psychrerythraea strain 34H
© Wan et al. 2016
Received: 25 October 2015
Accepted: 25 January 2016
Published: 6 February 2016
Colwellia psychrerythraea 34H is a psychrophilic bacterium able to produce docosahexaenoic acid (DHA). Polyketide synthase pathway is assumed to be responsible for DHA production in marine bacteria.
Five pfa genes from strain 34H were confirmed to be responsible for DHA formation by heterogeneous expression in Escherichia coli. The complexity of fatty acid profile of this strain was revealed by GC and GC–MS. Treatment of cells with cerulenin resulted in significantly reduced level of C16 monounsaturated fatty acid (C16:1Δ9t, C16:1Δ7). In contrast, the amount of saturated fatty acids (C10:0, C12:0, C14:0), hydroxyl fatty acids (3-OH C10:0 and 3-OH C12:0), as well as C20:4ω3, C20:5ω3 and C22:6ω3 were increased. RNA sequencing (RNA-Seq) revealed the altered gene expression pattern when C. psychrerythraea cells were treated with cerulenin. Genes involved in polyketide synthase pathway and fatty acid biosynthesis pathway were not obviously affected by cerulenin treatment. In contrast, several genes involved in fatty acid degradation or β-oxidation pathway were dramatically reduced at the transcriptional level.
Genes responsible for DHA formation in C. psychrerythraea was first cloned and characterized. We revealed the complexity of fatty acid profile in this DHA-producing strain. Cerulenin could substantially change the fatty acid composition by affecting the fatty acid degradation at transcriptional level. Acyl-CoA dehydrogenase gene family involved in the first step of β-oxidation pathway may be important to the selectivity of degraded fatty acids. In addition, inhibition of FabB protein by cerulenin may lead to the accumulation of malonyl-CoA, which is the substrate for DHA formation.
Docosahexaenoic acid (DHA, 22:6ω3) and eicosapentaenoic acid (EPA, 20:5) are the two main commercial omega-3 very long chain polyunsaturated fatty acids (VLCPUFAs) because of their health benefits. Their beneficial effects in human cardiovascular health, mental health, immune function and infant cognitive development are well documented . In addition, they are precursors for certain hormones and lipid-derived signaling molecules . Such VLCPUFAs can be synthesized in human with a low overall conversion rate using α-linolenic acid (ALA, 18:3) as substrate , rather than de novo synthesis. Alternatively, food or food additives rich in VLCPUFAs are a good choice to obtain VLCPUFAs directly.
Commercial VLCPUFAs are mainly produced from fish oils, algae oils or krill oils. Recently, it became more and more attractive to achieve commercial production of VLCPUFAs via metabolic engineering of microbes due to the limited natural resources. Microbes grow fast, are easy to be genetically manipulated and use versatile carbon sources. One of these successful examples, oleaginous yeast Yarrowia lipolytica has been genetically engineered to make EPA at 25 % dry cell weight and at 50 % of total lipids, which is much higher than those of other reported wild-type microbes . In general, VLCPUFAs could be synthesized through conventional fatty acid desaturase/elongase pathway in algae and animals. There are more than 20 fatty acid desaturases and fatty acid elongases involved in this process . Several key desaturases and elongases, such as Δ6 desaturase, Δ5 desaturase, Δ4 desaturase, Δ9 elongase, Δ8 desaturase, Δ5 elongase, have been transformed into microbes or plants to make them successfully produce VLCPUFAs [6–8]. Alternatively, VLCPUFAs can be synthesized via unconventional polyketide synthase (PKS) pathway as well. This pathway can only be found in narrow species of marine bacteria including Shewanella, Moritella, Vibrio, Photobacterium, Aureispira and Colwellia and some lower eukaryotes [9–12]. PKS system contains a set of essential domains for VLCPUFAs formation including malonyl-CoA:ACP acyltransferase (MAT), acyl carrier protein (ACP), 3-ketoacyl synthase (KS), 3-ketoacyl-ACP reductase (KR), acyltransferase (AT), chain length factor (CLF), enoyl reductase (ER), and 3-hydroxyacyl-ACP dehydrase/isomerase (HD/IS) . Compared to conventional desaturase/elongase pathway, there are only 3–5 genes (pfa gene cluster) responsible for de novo VLCPUFAs synthesis via PKS pathway . This cluster usually contains 1 MAT domain (pfaA), 1 KR domain (pfaA), 1 CLF domain (pfaC), 1 ER domain (pfaD), 2 KS domains (pfaA and pfaC), 3–6 ACP domains (pfaA), 1 AT (pfaB) and 2 DH/IS domains (pfaC) . In addition, PKS pathway was demonstrated to be anaerobic pathway and less NADPH was needed for the whole process compared with that of conventional pathway . Notably, PKS proteins typically show high sequence similarity to enzymes of fatty acid synthase (FAS) but not to any desaturase or elongase . Similarly to bacterial FAS, PKS make use of pfa gene cluster to accomplish the iterative cycle of condensation, reduction, dehydration and further reduction to produce VLCPUFAs. However, the details of synthetic route for VLCPUFAs via PKS and the mechanism of crosstalk between FAS and PKS are still unknown.
Colwellia psychrerythraea strain 34H, is a psychrotrophic gram-negative bacterium isolated from the sediment of Northeast Water Polynya, Greenland. The presence of a pfa gene cluster in the C. psychrerythraea genome suggested that this strain is able to produce DHA . Recently, C. psychrerythraea 34H was demonstrated to produce 2.3 % (w/w) of DHA and trace amount of EPA to the total fatty acids . Cerulenin [(2R,3S)-2,3-epoxy-4-oxo-7,10-trans,trans-dodecadienamide] is an FAS inhibitor originally developed as an antifungal antibiotic. It was reported to be an irreversible inhibitor of β–ketoacyl-ACP synthase I and II activities [17, 18] and it could increase the DHA content in Photobacterium, Shewanella and Moritella [19, 20]. However, the mechanism of this process is still unknown. In the present study we describe the DHA production by C. psychrerythraea at different temperatures and with the treatment of cerulenin. Five pfa genes were required for DHA production. Cerulenin treatment resulted in significant increased DHA production, along with the increased levels of medium chain, short chain saturated fatty acids (SFAs) and hydroxyl fatty acids (HFAs), at the expense of medium chain monounsaturated fatty acids (MUFAs). RNA deep sequencing revealed a certain number of genes either up regulated or down regulated, suggesting possible involvement of these genes in regulating of the fatty acid composition.
Chemicals, enzymes and kits
Fatty acid standards, 2-amino-2-methyl-1-propanol, 1 % trimethylchlorosilane in N,O-bis(trimethylsilyl)trifluoroacetamide (BSTFA) and cerulenin were purchased from Sigma-Aldrich (St. Louis, MO). Zwittergent, lysozyme and proteinase K were purchased from Merck, Sigma-Aldrich (St. Louis, MO) and Roche Diagnostic GmbH (Mannheim, Germany), respectively. RNA protect reagent, RNeasy Mini Kit and Qiaquick PCR spin columns were purchased from Qiagen (Qiagen, Venlo, The Netherlands). TURBO DNA-free Kit was purchased from Applied Biosystems (Foster City, CA). Ribo-Zero rRNA removal kit for Gram-negative bacteria was purchased from Epicentre. Superscript II Reverse Transcriptase, RNaseH and DNA polymerase I were purchased from Invitrogen (Waltham, MA). PrimeScript RT reagent kit with gDNA eraser and SYBR PremixEx Taq II kit were purchased from TaKaRa (Dalian, Japan).
Strains and growth conditions
The marine psychrophile bacterium C. psychrerythraea 34H was purchased from ATCC. This strain was cultured in marine broth 2216 (Difco Laboratories) at the temperature ranging from 0 to 15 °C without shaking. For solid media, agar was added at 15 g/L. Cerulenin dissolved in 50 % (v/v) ethanol was added to culture medium at various concentrations prior to cultivation. To investigate the effect of cerulenin on C. psychrerythraea, cells were cultivated at 10 °C with shaking at 200 rpm for 24 h after cerulenin added.
Fatty acid profile analysis
Colwellia psychrerythraea 34H cells were harvested and freeze-dried overnight. Fatty acid methyl esters (FAME) were prepared by direct trans-methylation method . The FAMEs were analyzed by GC-FID (7890A GC, Agilent Technologies, Palo Alto, CA, USA) equipped with a 30 m BPX70 column (0.25 mm inner diameter, 0.25 μm film thickness, SGE, Austin, Texas, USA). Peaks were integrated with Agilent Technologies ChemStation software. GC–MS analysis was essentially performed according to Zhou et al. . Dimethyloxazoline (DMOX) derivatives of FAMEs were prepared according to Fay and Richli  to determine the double bond positions. For hydroxyl fatty acids, the FAMEs were further derivatized by trimethylchlorosilane as described  to make trimethylsilyl (TMS) derivatives.
Heterogeneous expression of pfa gene cluster in Escherichia coli
The primers used in this section are listed in Additional file 1: Table S1. Long PCR was carried out with the high fidelity PrimeSTAR GXL DNA Polymerase (Takara Bio) according to manufacturer’s instruction. Specifically, the C. psychrerythraea pfaABCD sequence was amplified with primers mpfaA-F containing EcoRI restrction site and mpfaD-R containing SalI restrction site, respectively. The resulting DNA fragment of approximately 19 kbp was double digested with EcoRI and SalI, and then subcloned into the pColdI (Takara Bio) to generate pColdI-pfaA-D. pfaE was amplified with primers mpfaE-F containing KpnI restriction site and mpfaE-R containing BamHI restriction site, respectively. The resulting sequence of about 900 bp was double digested with KpnI and BamHI. Then it was inserted into pSTV28 (Takara Bio) to form pSTV28-pfaE. Two compatible vectors were co-transformed into E. coli DH5α by electroporation. Positive colonies were confirmed by PCR and sequencing. Transformants were then cultivated in Luria–Bertani (LB) medium with 50 mg/L ampicillin and 17 mg/L chloramphenicol. Half mL of cells pre-cultured at 30 °C was inoculated into 50 mL of fresh LB medium supplemented with 50 mg/L ampicillin and 17 mg/L chloramphenicol. When OD600 reached at about 0.6, 1 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) was added into the culture. Then cells were cultivated with shaking at 180 rpm at 15 °C for 24 h to allow gene expression. Cells were collected and dried.
RNA extraction, RNA enrichment and preparation of cDNA fragment libraries
Colwellia psychrerythraea 34H was grown in 2216 medium at 10 °C. Overnight grown cells were inoculated into culture flasks containing 3 mL of 2216 media each, under two treatments with three repeats, One was cultured at 10 °C (M10), while the other treatment also cultured at 10 °C but added cerulenin at a final concentration of 12 μg/mL (M10_C). Colwellia psychrerythraea cells were collected at exponentially growing phase (typical optical density at 600 nm~0.9), then fixed with two volumes of RNA protect reagent and harvested by centrifugation. The cellular pellet was then lysed with 280 μL of lysis buffer (10 % zwittergent, 15 mg/ml lysozime and 20 mg/ml proteinase K in TE buffer). Total RNA was isolated using the RNeasy Mini Kit according to manufacturer’s instructions. During this RNA purification, the residual DNA in the samples was removed using TURBO DNA-free Kit. NanoDrop ND-1000 spectrophotometer (NanoDrop Technology, Rockland, DE) and Experion Automated Electrophoresis were used to check the quantity and quality of RNA with the RNA StdSens Analysis Kit (Bio Rad). Genomic DNA contamination was tested by quantitative real-time PCR (qPCR) amplification of 16S rRNA using the iO5 iCycler (Bio-Rad). Total RNA samples were stored at −80 °C until further use.
Ribo-Zero rRNA Removal Kit was used to remove 16S and 23S rRNA from total RNA according to the manufacturer’s protocol with the exception that no more than 5 μg total RNA was treated per enrichment reaction. Each RNA sample was divided into multiple aliquots of ≤5 μg RNA and separate enrichment reactions were performed for each sample. Enriched mRNA samples were pooled and analyzed on the 2100 Bioanalzyer (Agilent) to asses the reduction of 16S and 23S rRNA before the preparation of cDNA fragment libraries. 100 ng of mRNA were treated with the Ambion RNA fragmentation reagents to generate 60-200 nucleotide RNA fragments. After fragmented RNA precipitation, the synthesis of the cDNA first strand was performed using random N6 primers and Superscript II Reverse Transcriptase. Second strand cDNA synthesis followed, using RNaseH and DNA Polymerase I. Double stranded cDNA was purified using Qiaquick PCR spin columns according to the manufacturer’s protocol.
Illumina Genome Analyzer RNA-Seq
This double-stranded cDNA was then processed for RNA-Seq using the Illumina Genomic DNA Sample Prep kit (Illumina, Inc., San Diego, CA). After loading the amplified material onto independent flow cells, sequencing was carried out by running 36 cycles on the Illumina Genome Analyzer IIx at the BGI (Beijing, China). Illumina RNA-Seq quality scores converted to phred format (http://www.phrap.com/phred/) analyzed the quality of the RNA-Seq reads.
Mapping to reference genome and annotated genes
Colwellia psychrerythraea 34H genome and gene information were available from NCBI (http://www.ncbi.nlm.nih.gov/genome/1067). After removing low quality reads and the reads containing adaptors, the remaining reads were aligned to the C. psychrerythraea 34H genome using software bowtie2 (version 2.0.2). Reads Per Kilobase per Million mapped reads (RPKM) were calculated for each transcript dividing number of mapped reads by the length of the transcripts and the number of total sequenced reads in this library.
Quantification of pfa transcription in C. psychrerythraea by Quantitative real-time reverse transcription PCR (qRT-PCR)
After RNA isolation as described in the previous section, first-strand cDNA was synthesized using a PrimeScript RT reagent kit with gDNA eraser. Two microliters of RT product was used as a template for qRT-PCR with SYBR PremixEx Taq II kit. The primers used in this study are listed in Additional file 2: Table S2.The qRT-PCR mixture (25 μL) contained 2 μL of cDNA, 0.4 μM each gene-specific primer. The qRT-PCR was performed as follows: 4 min at 94 °C followed by 40 cycles of 30 s at 94 °C, 30 s at 56 °C, and 25 s at 72 °C, followed by a melting cycle from 56 °C to 94 °C to check for amplification specificity. To assess for DNA contamination, RNA samples were run without reverse transcriptase. The constitutively expressed 16S ribosomal RNA gene was used as an internal control for relative quantification, using the comparative threshold method.
Optimal growth condition of C. psychrerythraea34H for DHA production
Effect of cerulenin (12 μg/mL) on the fatty acid composition of C. psychrerythraea 34H
Fatty acid composition (% of total peak area)
4.3 ± 0.2
0.4 ± 0.0
2.5 ± 0.1
0.3 ± 0.0
0.0 ± 0.0
0.2 ± 0.3
4.1 ± 0.0
0.3 ± 0.0
0.3 ± 0.0
0.0 ± 0.0
2.7 ± 0.0
2.9 ± 0.0
0.2 ± 0.1
28.3 ± 0.9
13.3 ± 0.2
9.0 ± 0.2
17.6 ± 0.3
0.4 ± 0.0
5.9 ± 0.1
0.1 ± 0.1
0.0 ± 0.0
0.2 ± 0.3
26.7 ± 0.4
0.3 ± 0.0
0.3 ± 0.0
0.1 ± 0.2
1.1 ± 0.0
0.3 ± 0.2
0.0 ± 0.0
19.7 ± 0.2
0.6 ± 0.0
3.0 ± 0.1
20.6 ± 0.3
0.7 ± 0.1
3.0 ± 0.1
0.7 ± 0.0
0.1 ± 0.1
2.6 ± 0.4
1.2 ± 0.2
0.1 ± 0.1
0.0 ± 0.0
0.0 ± 0.1
0.0 ± 0.0
2.3 ± 0.1
0.2 ± 0.0
8.0 ± 0.1
0.0 ± 0.0
3.7 ± 0.2
0.2 ± 0.2
0.0 ± 0.0
0.0 ± 0.0
0.0 ± 0.0
0.0 ± 0.0
0.3 ± 0.2
1.0 ± 0.1
0.1 ± 0.1
10.1 ± 0.1
0.3 ± 0.3
Five genes involved in PKS pathway are essential to DHA formation
Occurrence of multiple components of monounsaturated, polyunsaturated and hydroxyl fatty acids
In addition to DHA, trace amounts of other VLCPUFA including C20:4ω3 and C20:5ω3 were also found (Additional file 9: Figure S7). The identities of these components were confirmed by GC analysis comparing the FAME standards, or by GC–MS analysis of FAMEs and/or their DMOX derivatives.
Effect of cerulenin treatment on fatty acid composition of C. psychrerythraea 34H
Transcriptome analysis of C. psychrerythraea after cerulenin treatment revealed the differential response of pathways
Expression of most genes involved in PKS and fatty acid biosynthesis pathway was not affected by cerulenin
Differential expression of genes involved in fatty acid metabolism in C. psychrerythraea without or with 12 μg/mL cerulenin at 10 °C
Fold change (log2)
Gene involved in fatty acid biosynthesis
acetyl-CoA carboxylase subunit A
acetyl-CoA carboxylase subunit B
acetyl-CoA carboxylase subunit C
acetyl-CoA carboxylase subunit D
3-hydroxydecanoyl-ACP dehydratase, isomerase
3-ketoacyl-ACP synthase I
3-ketoacyl-ACP synthase II
3-ketoacyl-ACP synthase III
DNA-binding transcriptional repressor
Gene involved in fatty acid degradation
Long chain fatty acid-CoA ligase
Long chain fatty acid-CoA ligase
Fatty acid metabolism regulator
Gene involved in fatty acid modification or lipid metabolism
fatty acid cis/trans isomerase
Genes involved in PKS pathway
polyunsaturated fatty acid synthase
polyunsaturated fatty acid synthase
polyunsaturated fatty acid synthase
polyunsaturated fatty acid synthase
However, several genes involved in fatty acid degradation, such as fadD1 (CPS_RS05300) and fadB (CPS_RS01765), were down-regulated. Besides, the genes encoding for acyltransferases involved in polar lipid and neutral lipid formation, including glycerol-3-phosphate acyltransferase (CPS_RS10220, CPS_RS19415, CPS_RS00605), 1-acyl-sn-glycerol-3-phosphate acyltransferase (CPS_RS04525, CPS_RS21310) were not affected. Three out of six genes for other annotated acyltransferases (CPS_RS04435, CPS_RS08065, CPS_RS14560, CPS_RS21745, CPS_RS22505 and CPS_RS22520) had slight increased expression level. Their functions are to be further investigated.
Top up-regulated genes by cerulenin
Complexity of C. psychrerythraea fatty acid profile
Although the major fatty acids in C. psychrerythraea 34H cells were C16 fatty acids, the very complicated fatty acid profile was uncovered in this work. These FAs included short and medium chain SFAs, odd number carbon length FAs, short and medium chain MUFAs, multiple component of monounsaturated C16:1 or C18:1, branch fatty acids, PUFAs and 3-HFAs, although some of the components had very low levels. The existence of these components was confirmed by GC–MS analysis. The occurrence of trans monounsaturated C16:1 and PUFAs in 34H has also been recently reported . Under our experimental conditions, we detected more components that were not reported previously. These included trace amounts of C11:0, C13:0, three isomers of C14:1, 3-OH C10:0, C18:1Δ9trans, and branch fatty acids C16:0iso and C19:0iso. Such fatty acid profile could indicate the complexity of fatty acid synthesis and modification pathways in strain 34H.
We showed that the fatty acid profile was differentially affected by cerulenin treatment. Cerulenin, a specific inhibitor of condensing enzymes, had been reported to inhibit MUFA synthesis but not PUFA synthesis in deep-sea bacteria M. marina . This is consistent to our result that C16:1 and C18:1 MUFAs were significantly decreased, especially the trans fatty acids, while the VLCPUFA DHA increased. We also showed the enhanced levels of C10:0, C12:0, C14:0 and 3-hydroxy C12:0 after cerulenin treatment (Table 1). The condensing enzymes principally extend the fatty acid chain length by two carbons in each fatty acid synthesis cycle. Inhibition of the condensing enzymes by cerulenin could explain the accumulation of shorter chain length fatty acid, and decreased the longer chain length fatty acids. Indeed, C16:0 and C18:0 FAs were both reduced after cerulenin treatment. However, the significant increase of C14:0 levels might suggest the condensing step from C12:0 to C14:0 was not inhibited, same as the step from C8:0 to C10:0. It is possible that cerulenin might inhibit these steps differently, or these steps might be catalyzed by different condensing enzyme isomers. In addition, C16 and C18 MUFAs were decreased. Escherichia coli possesses at least three condensing enzymes: FabH selectively uses acetyl-CoA to initiate the fatty acid synthesis pathway, and FabB/FabF carry out the subsequent rounds of fatty acid elongation which condense malonyl-ACP with acyl-ACP to extend the acyl chain by two carbons . FabF is confirmed to be the sole condensing enzyme for subsequent elongation step in this strain . In contrast, Bacillus subtilis contains two FabH isozymes (FabHA and FabHB) which differ from E. coli FabH in that they could initiate both the straight and branched chain fatty acid synthesis. Similar to E. coli, only one initial condensing enzyme FabH (CPS_RS09575) and two other condensing enzymes, FabB (CPS_RS17030) and FabF (CPS_RS10240), were denoted in the genome of 34H.
No effect of PKS pathway at transcriptional level
In order to better understand the effect of cerulenin on the fatty acid level change, we performed the transcriptome analysis. To our surprise, there was no obvious change of PUFA PKS pathway at the transcription level after cerulenin treatment. The five pfa genes were confirmed to be essential to DHA production in strain 34H by heterologous co-expression in E. coli. However, all these five genes did not increase the expression level after cerulenin treatment, although the DHA level was significantly increased. This might suggest the increase of DHA content in C. psychrerythraea was not regulated by pfa gene cluster at the transcriptional level when cells were treated with cerulenin. It was reported that cultivation of P. profundum strain SS9 at elevated hydrostatic pressure or reduced temperature did not increase pfa gene expression despite the resulting increase in percentage composition of EPA under these conditions . We proposed that there could be some other genes or transcriptional factors related to the increase of DHA accumulation. Another explanation could be that the availability of malonyl-CoA for PKS pathway is higher after cerulenin treatment.
Besides, pfa genes are regulated by PfaR. PfaR is predicted to be an inactive transcriptase that does not seem to be affected for any effectors. Therefore, it has sense that cerulenin affecting fatty acid profile does not affect PfaR binding and thus does not affect PfaABCD transcription. A BLAST search against current databases using the sequence of PfaR revealed more than 100 putative PfaR homologs in marine bacteria, but not in M. marina (data not shown). In fact the PfaRABCD cluster is conserved in all these marine bacteria. Moreover, we carried out structural prediction and modeling of PfaR using the PHYRE. 263 residues (91 %) of PfaR could be modeled at >90 % accuracy using structure 4FCY (crystal structure of the bacteriophage mu transpososome) as a model, although the sequence identity between PfaR and 4FCY is only 14 % (unpublished data). MuA, like many other DNA transposases, and retroviral integrases share a conserved RNaseH-like or “DDE” catalytic domain, named for the 3 Mg2+ binding carboxylates in their active sites. However, these catalytic residues are not conserved in PfaR. Thus, we think PfaR is not an active transposase but rather a transcriptional regulator.
Similarly, the gene CPS_RS00390 coding for a cis/trans isomerase that could alter the ratio of cis- to trans-esterified fatty acids in phospholipids was not obviously changed at RNA level, although the trans monounsaturated fatty acids (C16:1Δ9t and C18:1Δ9t) were greatly reduced. Moreover, considering that the reduction level of C16:1Δ9t was much higher than that of C16:1Δ9c, a more complicated effect could exist rather than a simple effect of CPS_RS00390 expression level.
Down-regulation of some genes involved in fatty acids degradation
Bacterial FASII system contained basically the same set of enzymatic domains as PKS . Although FAS and PKS used the same substrates, acetyl-CoA and manoly-CoA, the final products were quite distinctive. In general, FAS was responsible for the essential short and medium chain fatty acids while PKS was only capable of production of VLCPUFAs . Until now, the mechanism of cross talk between FAS and PKS was still unknown. FAS pathway can be classified into fatty acids biosynthesis pathway and fatty acids degradation pathway. As shown in Table 2, most of the genes involved in the fatty acid biosynthesis pathway were not obviously affected in the presence of cerulenin. Notably, one of those genes, fabB encoding for 3-ketoacyl-ACP synthase, was responsible for the elongation of unsaturated fatty acids. FabB activity was demonstrated to be irreversibly inhibited by cerulenin . Although the content of MUFAs, such as C16:1Δ9t, C16:1Δ7c and C16:1Δ9c was indeed decreased (Table 1), the transcriptional expression of fabB was not down-regulated in the presence of cerulenin in this test, but slightly up-regulated (Table 2). It was also demonstrated that transcription of the operon containing the fabF gene from Bacillus subtilis increased eightfold in response to cerulenin treatment . We speculated that the minor increase of fabB or fabF might be to compensate the low activity of FabB/FabF protein.
In contrast, two genes (fadB and fadD1) involved in the fatty acid degradation pathway were dramatically down-regulated (Table 2). In the first step of long chain fatty acid degradation, fatty acids would be transported across the cell membrane via a transport and acyl-activation mechanism. Gene fadD encoding an inner membrane-bound acyl-CoA synthase was involved in this process. In general, there is only one acyl-CoA synthase in E. coli with broad substrate specificity . However, in C. psychrerythraea, two genes fadD1 and fadD2, were both denoted as acyl-CoA synthases. Only fadD1 was dramatically down-regulated while fadD2 was slightly up-regulated (Table 2). Therefore, two FadD enzymes might have different substrate specificities based on the changed fatty acids profile. This result awaits further investigation. The other down-regulated gene fadB was thought to encode for enoyl-CoA hydratase/3-hydroxyacyl-CoA dehydrogenase. This enzyme converted enoyl-CoA to 3-ketoayl-CoA via 3-hydroxylacyl-CoA through hydration and oxidation . Both fadD and fadB played critical roles in fatty acid degradation.
The final degradation of acyl-CoA is achieved by β-oxidation pathway. In E. coli, gene fadE encoding for acyl-CoA dehydrogenase, was demonstrated to be responsible for the conversion of acyl-CoA to enoyl-CoA. This gene was also thought to be the only acyl-CoA dehydrogenase gene in E. coli . According to the genome sequence of strain 34H, CPS_RS10085, denoted as fadE, codes for an acyl-CoA dehydrogenase. However, there were six other genes denoted as acyl-CoA dehydrogenase as well in this strain (Additional file 13: Figure S9). Notably, only CPS_RS02980 and CPS_RS10085 were up-regulated while the others including CPS_RS10075, CPS_RS05245, CPS_RS12125, CPS_RS20065 and CPS_RS02940 were down-regulated (Additional file 13: Figure S9). The expansion of acyl-CoA dehydrogenase gene family might be important to the selective oxidation of short (C4–C6), medium (C8–C14), long (C16–C18) and very long (>C20) chain fatty acids. This could be one of the reasons led to the higher accumulation of DHA after cerulenin treatment.
Other top regulated genes by cerulenin might contribute to the changes of fatty acid profile in strain 34H
Except for fatty acid metabolism related genes, other significantly regulated genes were also investigated. Three of the top 20 up-regulated genes are TetR family regulators (CPS_RS06200, CPS_RS07325 and CPS_RS07450, Additional file 12: Table S4). They belong to one-component systems in which a single polypeptide contains both a sensory domain and a DHA-binding domain . They are widely associated with many biological processes, such as antibiotic resistance and fatty acids synthesis [36, 37]. These proteins have several homologs in the Protein Data Bank (PDB). One of them is FadR. FadR binds one lauroyl (C12)-CoA molecule per FadR monomer, with its acyl chain moiety in the center of the FadR molecule, enclosed within a tunnel-like substrate-binding pocket surrounded by hydrophobic residues, and the CoA moiety interacting with basic residues on the protein surface . Recently, fadR was demonstrated to be an important transcriptional factor that was capable of positively control the fabA and fabB/F genes in E. coli . Over-expression of E. coli FadR could dramatically increase the fatty acid titer by 7.5-fold over wild-type strains. Moreover, the unsaturated fatty acids including 14:1, 16:1 and 18:1 were increased . Besides, FadR was demonstrated to dissociate with C16–C18 fatty acids but not with C10 fatty acids . Thus, cerulenin is not expected to affect FadR regulation.
The other most up-regulated gene was CPS_RS06205. The predicted protein structure of this gene showed high similarity to OYEs. OYEs have been found to be highly expressed in the EPA producing strain S. oneidensis . The expression of OYE in S. oneidensis was observed in the absence of cerulenin. Thus, CPS_RS06205 upregulation seems to be related to the metabolic changes that drive DHA accumulation and not to be the cerulenin effect itself.
CPS_RS04135 denoted as a CCR family gene, was dramatically up-regulated as well. It could reduce double bond in the accumulated MUFAs and other fatty acids to produce the extender units needed by PKS for DHA synthesis. Therefore, the accumulation of precursors might be one possible reason for DHA increase in cerulenin-treated cells. On the other side, CPS_RS07330 (14-fold) is homolog to the multidrug efflux protein AcrB. AcrA/B system is controlled by the TetR-like protein ArcR. This efflux system could be upregulated to cope with the toxicity of cerulenin. The ability to reduce the dye 2,3,5-triphenyltetrazolium chloride (TTC) has been shown to be directly associated with EPA production in marine Gram-negative bacteria in solid or liquid medium . The up-regulated oxidoreductases could be responsible of this reduction. Moreover, detoxification of H2O2 is other characteristic of PUFA producing deep-sea bacteria  that could also be due to the expression of different reductases.
Although C. psychrerythraea 34H cell produces DHA, it can only grow at relatively low temperature, and the transformation system for genetic modification has not been set up. However, this strain synthesizes DHA via PKS pathway. Compared to conventional desaturase/elongase pathway, PKS is an anaerobic pathway and less NADPH was needed for the whole process. More important, unlike desaturase/elongase pathway, the final products of PKS are only VLCPUFAs. Our current work is to seek and investigate good gene candidates from microbes harboring PKS pathway. These genes can be engineered into other expression systems such as oleaginous yeast Y. lipolytica  to produce high levels of VLCPUFAs.
In summary, the complexity of psychrophilic C. psychrerythraea fatty acid profile was revealed by GC or GC–MS. Cerulenin treatment of this strain resulted in enhanced accumulation of short and medium chain saturated fatty acids, 3-hydroxy fatty acids, and DHA, at the expense of C16 and C18 saturated and MUFAs. Five pfa genes involved in PKS pathway were demonstrated to be responsible for DHA production by heterogeneous expression in E. coli. However, genes involved in polyketide synthase pathway and fatty acid biosynthesis pathway were not obviously affected by cerulenin treatment. In contrast, several genes involved in fatty acid degradation or β-oxidation pathway were dramatically reduced at the transcriptional level. The results suggested that acyl-CoA dehydrogenase gene family involved in the first step of β-oxidation pathway may be important to the selectivity of degraded fatty acids. In addition, inhibition of FabB protein by cerulenin may lead to the accumulation of malonyl-CoA, which is the substrate for DHA formation.
XW, YG, FH and GM conceived the study, XW, YP and GM conducted experiments. XW, YP, GM and XZ analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.
We are grateful to Dr. Juan Maria Garcia-Lobo and Dr. Maria Cruz Rodriguez for RNA-Seq analysis performed in the massive sequencing service at the IBBTEC. This work was supported by the National Natural Science Foundation of China (3100078) and 948 project (2015-Z3) to X Wan and by Grants BIO2010-14809 and BFU2014-55534-C2-2-P from the Spanish Ministry of Economy and Competitiveness to GM.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Swanson D, Block R, Mousa SA. Omega-3 fatty acids EPA and DHA: health benefits throughout life. Adv Nutr. 2012;3:1–7.View ArticleGoogle Scholar
- Friedrichs W, Ruparel SB, Marciniak RA, deGraffenried L. Omega-3 fatty acid inhibition of prostate cancer progression to hormone independence is associated with suppression of mTOR signaling and androgen receptor expression. Nutr Cancer. 2011;63:771–7.View ArticleGoogle Scholar
- Valenzuela R, Bascunan K, Chamorro R, Barrera C, Sandoval J, Puigrredon C, Parraguez G, Orellana P, Gonzalez V, Valenzuela A. Modification of docosahexaenoic acid composition of milk from nursing women who received alpha linolenic acid from chia oil during gestation and nursing. Nutrients. 2015;7:6405–24.View ArticleGoogle Scholar
- Xue Z, Sharpe PL, Hong SP, Yadav NS, Xie D, Short DR, Damude HG, Rupert RA, Seip JE, Wang J, et al. Production of omega-3 eicosapentaenoic acid by metabolic engineering of Yarrowia lipolytica. Nat Biotechnol. 2013;31:734–40.View ArticleGoogle Scholar
- Gong Y, Wan X, Jiang M, Hu C, Hu H, Huang F. Metabolic engineering of microorganisms to produce omega-3 very long-chain polyunsaturated fatty acids. Prog Lipid Res. 2014;56:19–35.View ArticleGoogle Scholar
- Ruiz-Lopez N, Haslam RP, Usher S, Napier JA, Sayanova O. An alternative pathway for the effective production of the omega-3 long-chain polyunsaturates EPA and ETA in transgenic oilseeds. Plant Biotechnol J. 2015;. doi:10.1111/pbi.12328.Google Scholar
- Ruiz-Lopez N, Haslam RP, Napier JA, Sayanova O. Successful high-level accumulation of fish oil omega-3 long-chain polyunsaturated fatty acids in a transgenic oilseed crop. Plant J. 2014;77:198–208.View ArticleGoogle Scholar
- Guo B, Jiang M, Wan X, Gong Y, Liang Z, Hu C. Identification and heterologous expression of a delta4-fatty acid desaturase gene from Isochrysis sphaerica. J Microbiol Biotechnol. 2013;23:1413–21.View ArticleGoogle Scholar
- Ujihara T, Nagano M, Wada H, Mitsuhashi S. Identification of a novel type of polyunsaturated fatty acid synthase involved in arachidonic acid biosynthesis. FEBS Lett. 2014;588:4032–6.View ArticleGoogle Scholar
- Shulse CN, Allen EE. Widespread occurrence of secondary lipid biosynthesis potential in microbial lineages. PLoS One. 2011;6:e20146. doi:10.1371/journal.pone.0020146.View ArticleGoogle Scholar
- Shulse CN, Allen EE. Diversity and distribution of microbial long-chain fatty acid biosynthetic genes in the marine environment. Environ Microbiol. 2011;13:684–95.View ArticleGoogle Scholar
- Okuyama H, Orikasa Y, Nishida T, Watanabe K, Morita N. Bacterial genes responsible for the biosynthesis of eicosapentaenoic and docosahexaenoic acids and their heterologous expression. Appl Environ Microbiol. 2007;73:665–70.View ArticleGoogle Scholar
- Napier JA. Plumbing the depths of PUFA biosynthesis: a novel polyketide synthase-like pathway from marine organisms. Trends Plant Sci. 2002;7:51–4.View ArticleGoogle Scholar
- Metz JG, Roessler P, Facciotti D, Levering C, Dittrich F, Lassner M, Valentine R, Lardizabal K, Domergue F, Yamada A, et al. Production of polyunsaturated fatty acids by polyketide synthases in both prokaryotes and eukaryotes. Science. 2001;293:290–3.View ArticleGoogle Scholar
- Methe BA, Nelson KE, Deming JW, Momen B, Melamud E, Zhang X, Moult J, Madupu R, Nelson WC, Dodson RJ, et al. The psychrophilic lifestyle as revealed by the genome sequence of Colwellia psychrerythraea 34H through genomic and proteomic analyses. Proc Natl Acad Sci USA. 2005;102:10913–8.View ArticleGoogle Scholar
- Hashimoto M, Orikasa Y, Hayashi H, Watanabe K, Yoshida K, Okuyama H. Occurrence of trans monounsaturated and polyunsaturated fatty acids in Colwellia psychrerythraea strain 34H. J Basic Microbiol. 2015;55:838–45.View ArticleGoogle Scholar
- Trajtenberg F, Altabe S, Larrieux N, Ficarra F, de Mendoza D, Buschiazzo A, Schujman GE. Structural insights into bacterial resistance to cerulenin. FEBS J. 2014;281:2324–38.View ArticleGoogle Scholar
- Price AC, Choi KH, Heath RJ, Li Z, White SW, Rock CO. Inhibition of beta-ketoacyl-acyl carrier protein synthases by thiolactomycin and cerulenin. Struct Mech J Biol Chem. 2001;276:6551–9.Google Scholar
- Morita N, Nishida T, Tanaka M, Yano Y, Okuyama H. Enhancement of polyunsaturated fatty acid production by cerulenin treatment in polyunsaturated fatty acid-producing bacteria. Biotechnol Lett. 2005;27:389–93.View ArticleGoogle Scholar
- Allen EE, Facciotti D, Bartlett DH. Monounsaturated but not polyunsaturated fatty acids are required for growth of the deep-sea bacterium Photobacterium profundum SS9 at high pressure and low temperature. Appl Environ Microbiol. 1999;65:1710–20.Google Scholar
- Zhou XR, Singh SP, Green AG. Characterisation of the FAD2 gene family from Hiptage benghalensis. A ricinoleic acid accumulating plant. Phytochemistry. 2013;92:42–8.View ArticleGoogle Scholar
- Fay L, Richli U. Location of double bonds in polyunsaturated fatty acids by gas chromatography-mass spectrometry after 4,4-dimethyloxazoline derivatization. J Chromatogr A. 1991;541:89–98.View ArticleGoogle Scholar
- Orikasa Y, Tanaka M, Sugihara S, Hori R, Nishida T, Ueno A, Morita N, Yano Y, Yamamoto K, Shibahara A, et al. pfaB products determine the molecular species produced in bacterial polyunsaturated fatty acid biosynthesis. FEMS Microbiol Lett. 2009;295:170–6.View ArticleGoogle Scholar
- Jiang H, Zirkle R, Metz JG, Braun L, Richter L, Van Lanen SG, Shen B. The role of tandem acyl carrier protein domains in polyunsaturated fatty acid biosynthesis. J Am Chem Soc. 2008;130:6336–7.View ArticleGoogle Scholar
- Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40:D109–14.View ArticleGoogle Scholar
- Kelley LA, Sternberg MJ. Protein structure prediction on the web: a case study using the Phyre server. Nat Protoc. 2009;4:363–71.View ArticleGoogle Scholar
- Romano D, Contente ML, Molinari F, Eberini I, Ruvutuso E, Sensi C, Amaretti A, Rossi M, Raimondi S. Recombinant S. cerevisiae expressing old yellow enzymes from non-conventional yeasts: an easy system for selective reduction of activated alkenes. Microb Cell Fact. 2014;13:60. doi:10.1186/1475-2859-13-60.View ArticleGoogle Scholar
- Quade N, Huo L, Rachid S, Heinz DW, Muller R. Unusual carbon fixation gives rise to diverse polyketide extender units. Nat Chem Biol. 2012;8:117–24.View ArticleGoogle Scholar
- White SW, Zheng J, Zhang YM, Rock. The structural biology of type II fatty acid biosynthesis. Annu Rev Biochem. 2005;74:791–831.View ArticleGoogle Scholar
- Fujita Y, Matsuoka H, Hirooka K. Regulation of fatty acid metabolism in bacteria. Mol Microbiol. 2007;66:829–39.View ArticleGoogle Scholar
- Allen EE, Bartlett DH. Structure and regulation of the omega-3 polyunsaturated fatty acid synthase genes from the deep-sea bacterium Photobacterium profundum strain SS9. Microbiology. 2002;148:1903–13.View ArticleGoogle Scholar
- Hauvermale A, Kuner J, Rosenzweig B, Guerra D, Diltz S, Metz JG. Fatty acid production in Schizochytrium sp.: involvement of a polyunsaturated fatty acid synthase and a type I fatty acid synthase. Lipids. 2006;41:739–47.View ArticleGoogle Scholar
- Schujman GE, Choi KH, Altabe S, Rock CO, de Mendoza D. Response of Bacillus subtilis to cerulenin and acquisition of resistance. J Bacteriol. 2001;183:3032–40.View ArticleGoogle Scholar
- Campbell JW, Cronan JE Jr. The enigmatic Escherichia coli fadE gene is yafH. J Bacteriol. 2002;184:3759–64.View ArticleGoogle Scholar
- Cuthbertson L, Nodwell JR. The TetR family of regulators. Microbiol Mol Biol Rev. 2013;77:440–75.View ArticleGoogle Scholar
- Agari Y, Agari K, Sakamoto K, Kuramitsu S, Shinkai A. TetR-family transcriptional repressor thermus thermophilus FadR controls fatty acid degradation. Microbiology. 2011;157:1589–601.View ArticleGoogle Scholar
- Zhang F, Ouellet M, Batth TS, Adams PD, Petzold CJ, Mukhopadhyay A, Keasling JD. Enhancing fatty acid production by the expression of the regulatory transcription factor FadR. Metab Eng. 2012;14:653–60.View ArticleGoogle Scholar
- Zhang H, Zheng B, Gao R, Feng Y. Binding of Shewanella FadR to the fabA fatty acid biosynthetic gene: implications for contraction of the fad regulon. Protein Cell. 2015;6:667–79.View ArticleGoogle Scholar
- Brige A, Van den Hemel D, Carpentier W, De Smet L, Van Beeumen JJ. Comparative characterization and expression analysis of the four old yellow enzyme homologues from Shewanella oneidensis indicate differences in physiological function. Biochem J. 2006;394:335–44.View ArticleGoogle Scholar
- Ryan J, Farr H, Visnovsky S, Vyssotski M, Visnovsky G. A rapid method for the isolation of eicosapentaenoic acid-producing marine bacteria. J Microbiol Methods. 2010;82:49–53.View ArticleGoogle Scholar
- Tilay A, Annapure U. Novel simplified and rapid method for screening and isolation of polyunsaturated fatty acids producing marine bacteria. Biotechnol Res Int. 2012;2012:542721.View ArticleGoogle Scholar