- Open Access
Integrated OMICS guided engineering of biofuel butanol-tolerance in photosynthetic Synechocystissp. PCC 6803
Biotechnology for Biofuelsvolume 6, Article number: 106 (2013)
Photosynthetic cyanobacteria have been recently proposed as a ‘microbial factory’ to produce butanol due to their capability to utilize solar energy and CO2 as the sole energy and carbon sources, respectively. However, to improve the productivity, one key issue needed to be addressed is the low tolerance of the photosynthetic hosts to butanol.
In this study, we first applied a quantitative transcriptomics approach with a next-generation RNA sequencing technology to identify gene targets relevant to butanol tolerance in a model cyanobacterium Synechocystis sp. PCC 6803. The results showed that 278 genes were induced by the butanol exposure at all three sampling points through the growth time course. Genes encoding heat-shock proteins, oxidative stress related proteins, transporters and proteins involved in common stress responses, were induced by butanol exposure. We then applied GC-MS based metabolomics analysis to determine the metabolic changes associated with the butanol exposure. The results showed that 46 out of 73 chemically classified metabolites were differentially regulated by butanol treatment. Notably, 3-phosphoglycerate, glycine, serine and urea related to general stress responses were elevated in butanol-treated cells. To validate the potential targets, we constructed gene knockout mutants for three selected gene targets. The comparative phenotypic analysis confirmed that these genes were involved in the butanol tolerance.
The integrated OMICS analysis provided a comprehensive view of the complicated molecular mechanisms employed by Synechocystis sp. PCC 6803 against butanol stress, and allowed identification of a series of potential gene candidates for tolerance engineering in cyanobacterium Synechocystis sp. PCC 6803.
Due to its high energy content and superior chemical properties such as low volatility and corrosiveness, and its compatibility with the existing fuel storage and distribution infrastructure, butanol has been proposed as a good candidate for next-generation transportation biofuel [1, 2]. Traditionally, bio-butanol can be produced by anaerobic Gram-positive bacteria, such as Clostridium acetobutylicum through a so-called acetone-butanol-ethanol (ABE) fermentation process [3, 4]. Although significant improvements have been made in the past decades to increase efficiency of the ABE process through a combination of strain screening, genetic engineering and process optimization [5–8], butanol production from the fermentation processes is still not competitive economically. As one of the alternatives, photosynthetic cyanobacteria have recently attracted significant attention as a ‘microbial factory’ to produce biofuels and chemicals due to their capability to utilize solar energy and CO2 as the sole energy and carbon sources, respectively [9, 10]. Recent synthetic biology efforts have led to successful production of n-butanol, isobutyraldehyde and isobutanol in cyanobacterium Synechococcus elongatus PCC 7942 [11, 12], demonstrating the potentials of using engineered photosynthetic microbes for large-scale production of butanol or other biofuel products in the future.
Currently, the butanol production by the synthetic cyanbacterial systems is at a level of a few dozen or hundred milligrams per liter , much lower than the native Clostridium or even synthetic Escherichia coli systems [13–15]. To improve productivity, one of the key issues needed to be addressed is the low tolerance of the photosynthetic hosts to butanol [16, 17]. The tolerance mechanism of native Clostridium strains to butanol has been well-studied [16–19]. For example, analysis of butanol tolerant transposon-insertion mutants of Clostridium beijerinckii NCIMB 8052 have led to the discovery that butanol-tolerance is associated with reduced activity of the enzyme, glycerol dehydrogenase . Recently a functionally unknown protein (encoded by SMB_G1518) with a hypothetical alcohol interacting domain was also found negatively related to butanol tolerance . In E. coli, a global transcription factor cyclic AMP receptor protein (CRP) was also engineered for increasing butanol tolerance . However, currently information related to biofuel tolerance in cyanobacteria is very limited.
Recently various genome-wide approaches, such as genomic library enrichment and whole-genome sequencing of tolerant mutants were also employed to identify genes conferring enhanced tolerance to n-butanol in E. coli[23, 24]. The results showed that microbes tend to employ multiple and synergistic resistance mechanisms in dealing with a single stress , and to fully interpret the complicated and synergistic tolerance mechanism, genome-wide based analytical approaches are necessary . In a previous study, we investigated responses of Synechocystis sp. PCC 6803 (hereafter Synechocystis) to butanol using an iTRAQ - LC-MS/MS based proteomics, the results identified 303 proteins differentially regulated by butanol . To further decipher responses at transcript and metabolite levels, and to identify gene targets relevant to butanol tolerance, in this study, we applied an integrated approach coupling quantitative RNA-seq transcriptomics approach, quantitative reverse-transcript PCR (qRT-PCR) and GC-MS based metabolomics to analyze cellular responses of Synechocystis to butanol exposure. The transcriptomic result revealed very similar response patterns as those identified by the previous proteomic analysis that multiple resistance mechanisms may be utilized in coping with butanol stress in Synechocystis; and the metabolomic analysis showed that 46 chemically classified metabolites were differentially regulated by butanol treatment, including 3-phosphoglycerate, glycine and urea which were elevated in butanol-treated cells. The integrated analysis led to the identification of a series of potential gene targets and pathways for tolerance engineering, we then constructed gene knockout mutants for three selected butanol-induced genes, sll0690, slr0947 and slr1295, and comparative phenotype analyses showed that their disruptions led to increased sensitivity to butanol, suggesting the gene targets identified can be used for engineering butanol tolerance in Synechocystis.
Results and discussion
Overview of RNA-Seq transcriptomics analysis
To make the transcriptomics data comparable with previous proteomics data, we used the identical sampling conditions for transcriptomics as our previous proteomic analysis . As described previously, Synechocystis was grown in BG11 supplemented with 0.20% (v/v) butanol and cell samples of both control and butanol treatment were collected by centrifugation (8,000 × g for 10 min at 4°C) at 24 h, 48 h and 72 h, corresponded to middle-exponential, exponential-stationary transition and stationary phases of the cell growth, respectively.
A total of 79.5-million raw sequencing reads was obtained from the RNA-seq transcriptomics analysis of six samples, with average reads of 13.2-million reads. After a two-step standard data filtering process, first to eliminate reads with low-quality bases (such as multiple N) and reads shorter than 20 bp, and then to eliminate sequence reads mapped to non-coding RNA of Synechocystis, a total of 27.5-million qualified mRNA-based sequence reads were identified (Table 1). The qualified sequence reads have an average genome mapping ratio of 66.4%. To assess the analytical reproducibility between biological replicates, we collected two biological replicates for butanol treated samples at 72 h, and plotted them using the normalized Reads Per Kilobase of Gene per Million Mapped Reads (RPKM) values, the result showed a correlation coefficient around 0.991 (Figure 1), indicating the overall good quality of RNA-sequencing based transcriptomics technology. The sequence reads matched to all 3189 coding genes in Synechocystis genome (data not shown), suggesting excellent sequencing depth and overall transcript coverage.
Using a strict criterion of 1.5-fold change at all three time points (i.e., 24, 48 and 72 h), we determined that 278 genes were induced upon butanol exposure, out of which 70 important genes with known functional categories were listed in Table 2. Functional category analysis of the induced genes showed that the most affected functional categories were “hypothetical proteins”, representing a total of more than 40% of all the up-regulated genes, consistent with the fact that nearly half of the genes in the Synechocystis genome are still annotated as hypothetical up to now . Based on their expression level and regulation patterns by butanol, a subset of 10 genes was randomly selected for quantitative RT-PCR validation. Comparative RT-PCR analysis was performed for the genes between the butanol-treated sample and control at 48 h. The results showed very similar trends between qRT-PCR and RNA-Seq transcriptomics data (Table 3), suggesting a good quality of RNA-seq data.
Potential gene targets related to butanol tolerance
Our previous proteomic analysis found that the Synechocystis cells employed a combination of approaches to cope with butanol stress, and the responses included an induced common stress response, modifications of cell envelope, and induction of multiple transporters and signal transduction proteins against butanol stress . Transcriptomic analysis showed very similar responses:
Heat-shock and general stress proteins: early analysis of butanol tolerance in both native and unnatural producing microorganisms showed that heat-shock proteins were relevant to tolerance [7, 17]. Our quantitative proteomics found that DnaJ1 (Slr0093) was significantly induced at 48 h after butanol treatment . At transcriptional level, we found that four genes involved in heat shock and general stress responses were induced (i.e., slr0093, sll1988, sll1388 and slr1854). In addition, slr1204 encoding a putative serine protease (HtrA) and slr0835 encoding a MoxR protein homolog were also up-regulated significantly by butanol (Table 2). HtrA-type serine proteases participate in folding and degradation of aberrant proteins and in processing and maturation of native proteins, and htrA mutation often conferred a pleiotropic phenotype that can include high sensitivity to various stress . The MoxR family AAA+ proteins are ubiquitous proteins that employ the energy obtained from ATP hydrolysis to remodel proteins, DNA or RNA. Early studies have showed that some members of this protein group can potentially function as molecular chaperones involved in the assembly of protein complexes , and be involved in stress resistance and virulence in Francisella tularensis .
Oxidative stress response: early studies showed that solvent like ethanol or butanol can challenge cells by causing increased production of highly reactive oxygen species (ROS) . Transcriptomic analysis found that butanol induced expression of slr1828 and sll0248 genes encoding a petF-like ferredoxin and a flavodoxin protein (Table 2), respectively, consistent with the up-regulation of these two proteins in proteomics dataset. In addition, transcriptomic analysis showed that other genes involved in oxidative stress response, such as ssl2250 encoding a bacterioferritin-associated ferredoxin, slr1846 encoding a putative monothiol glutaredoxin and slr1795 encoding a peptide methionine sulfoxide reductase were also up-regulated (Table 2). Recent study showed that bacterioferritin comigratory proteins, along with glutathione peroxidase-reductase, were responsible for detoxification of bentazone-derived peroxide in a S. elongatus PCC7942 mutant Mu2 . Monothiol glutaredoxins was found with roles in actin cytoskeleton remodeling and cellular defenses against oxidative stress caused by ROS accumulation in Saccharomyces cerevisiae and Schizosaccharomyces pombe [33, 34]. In addition, monothiol glutaredoxin (Slr1846) was found up-regulated by ethanol in Synechocystis . Specific modifications of certain amino acid side chains are common during oxidative stress. Cysteine and methionine both contain a sulfur atom in their side chains and are among the most easily oxidized amino acids. Methionine sulfoxides can be reduced back to the methionines by peptide methionine sulfoxide reductase (MSR), providing cells with a mechanism to repair proteins damaged by reactive oxygen species rather than having them degraded and then re-synthesizing them de novo . Induction of the methionine sulfoxide reductase by oxidative stress has been found in anaerobic Desulfovibrio vulgaris, E. coli, S. cerevisiae and Synechocystis [35–38].
Transporters: transcriptomics analysis identified 19 membrane transporters were up-regulated. Among them only two genes, sll0689 and slr1512 which were in the same operon with butanol-induced slr1515, were identified in the previous proteomics analysis . Interestingly, the up-regulated transporters involved a wide range of putative substrates, including iron, Na+/H+, nitrate/nitrite, phosphate, sodium, potassium, urea, bicarbonate and sulfate (Table 2). Moreover, many of these transporters were induced at significantly high fold changes, such as slr2131 encoding a RND multidrug efflux transporter up-regulated 12.87 fold at 72 h, and sll1428 encoding a probable sodium-dependent transporter up-regulated 32.0 folds. Other up-regulated genes included sll1697 which encodes a well-studied multidrug efflux pump NorA . Exact functions of these transporters in butanol tolerance may worth further investigation.
Protein translocation: Bacteria have two major protein translocation systems, one of which is catalyzed by the Sec-dependent protein translocation system, and another is the Twin-arginine (Tat) protein translocation system [40, 41]. Our proteomic analysis showed that SecE protein (Ssl3335) of Sec-dependent translocation system and Tha4 protein (Slr1047) of the Tat translocation system were up-regulated by butanol. Trnascriptomic analysis showed that ssr3307 encoding a preprotein translocase SecG subunit, sll0616 encoding a preprotein translocase SecA subunit and slr1046 encoding a putative TatA protein, were up-regulated by butanol. Genes slr1046 and slr1047 were organized in the same operon. The results confirmed that enhanced protein translocation systems may be an important mechanism against butanol stress.
Cell envelope: Cell envelope is the important barrier in protecting cells. Consistent with proteomic results, our transcriptomic analysis also found that many genes involved in cell envelope function were up-regulated upon butanol exposure, such as sll2010 encoding UDP-N-acetylmuramoylalanine--D-glutamate ligase, slr0528 encoding UDP-N-acetylmuramoylalanyl-D-glutamate-2, 6-diaminopimelate ligase and sll088 encoding UDP-N-acetylglucosamine--peptide n-acetylglucosaminyltransferase (Table 3). Their up-regulation was supposed to strengthen cell wall structure against butanol stress.
Regulatory genes: Previous proteomic analysis showed that several signal transduction proteins involved in cell mobility (i.e. Che type) and nitrate induction, and repression of genes encoding nitrate respiration enzymes (i.e. NarL subfamily) were up-regulated by butanol . Transcriptomics analysis identified 11 butanol-induced signal transduction genes. The induced genes included two Che type response regulators (i.e. slr1042, slr1037) and one putative phototaxis histidine kinase (sll0043) involved in cell mobility, and one gene (ssl0707) involved in nitrogen metabolism. Gene ssl0707 encodes a nitrogen regulatory protein P-II belonging to the NtcA regulon in cyanobacteria . Although the transcriptomic results confirmed that regulation of cell mobility and nitrogen responses are important in combating butanol stress, none of regulatory genes/proteins was identified in both transciptomic and proteomic datasets, suggesting the complicity of signal transduction in Synechocystis, and also the insufficiency to use any single ‘omics’ approach to characterize the complexity of biological systems . To compare the proteomic and transcriptomic datasets quantitatively, 11 common genes/proteins up-regulated in both transcriptomics and proteomics datasets were listed in Table 4. The results also showed the very similar trends of up-regulation.
In our previous proteomic analysis, using a cutoff of 1.5-fold change and a p-value less than 0.05, we determined that 63 and 79 proteins were up-regulated between control and butanol treatments conditions at 24 h and 48 h, respectively; among which 35 proteins were up-regulated at both time points . Comparison of proteomic and transcriptomic datasets showed that among the 278 genes up-regulated by butanol, 17 induced genes also had their corresponding proteins up-regulated (Table 4), 10 genes had their corresponding proteins down-regulated, and 251 induced genes have their protein levels unchanged. The finding that a relatively low number of genes and proteins shared the same up-regulation patterns, was probably due to the fact we used highly strict criteria in determining induced genes (i.e., up-regulated at all three time points in this study). In spite of low correlation between the two datasets, the patterns of metabolic changes key to the butanol tolerance seemed similar, as described above for each of the functional categories.
One goal of the integrated OMICS analysis is to achieve a complete coverage of cellular molecules by using complementary techniques targeting different levels of information (i.e., RNA, protein or metabolites) . In this study, our transcriptomic analysis also revealed new cellular responses which were not observed in the previous proteomic analysis : i) Enhanced production of storage compounds: Polyhydroxyalkanoates (PHAs) are common carbon storage compounds that are accumulated during unbalanced growth conditions . Two genes involved in PHA biosynthesis, slr1994 encoding a PHA-specific acetoacetyl-CoA reductase and slr1993 encoding a PHA-specific beta-ketothiolase were found up-regulated by butanol (Table 2). Cyanophycin is a non-ribosomally synthesized peptide, composed of arginine and aspartic acid, accumulates when cells are grown under all unbalanced nutrient conditions except nitrogen starvation, and has been considered as a primary nitrogen reserve compound in cyanobacteria . Transcriptomic analysis showed that the key gene involved in cyanophycin synthesis, slr2002 encoding cyanophycin synthetase was up-regulated by butanol (Table 2). Although PHA and cyanophycin accumulation has been reported for many natural stress conditions, it may worth further investigation how these pathways respond to butanol stress; ii) Enhanced carotenoid biosynthesis: three genes involved in carotenoid biosynthesis were up-regulated: slr1254 encoding phytoene desaturase, slr0940 encoding zeta-carotene desaturase and slr0899 encoding cyanate lyase (Table 2). The results were consistent with the increased photosynthetic activity of Synechocystis upon butanol stress . Carotenoid biosynthesis has been found up-regulated by strong light in Synechococcus PCC7942 , and in stress-tolerant mutants of Haematococcus pluvialis. The results provided further evidences that the integrated OMICS approach could be advantageous in revealing global cellular responses.
Metabolomic signatures related to butanol response
GC-MS based metabolomic analysis was used to characterize the time-series metabolic responses of Synechocystis to butanol exposure, with unperturbed cultures as controls. Cell samples used for metabolomic analysis were collected at 24, 48 and 72 h, respectively, the identical time points of sampling for transcriptomic analysis. Three biological replicates were collected for each time point and treatment, thereby yielding a total of 18 samples. The analysis showed that a total of 73 metabolites were chemically identified with great confidence. Although more metabolites were detected in butanol-treated samples (70.4 ± 2.74) than the control samples (64.12 ± 4.01), the number of metabolites identified varied only slightly within control or treatment bins, implying an overall good analytical quality. To further assess the reproducibility of GC-MS metabolomics, we analyzed three technical replicates of one selected sample, and the results showed that most of the metabolites were identified in technical replicates (Date not shown).
The score plot of principal component analysis (PCA) was applied to evaluate the similarities and differences between the 18 metabolomic profiles (Figure 2). The score plot revealed the following features: i) the samples with or without butanol treatment at different time points were distinctly separated, suggesting significant metabolic differences between samples; ii) for the control samples, metabolic changes along the time courses were relatively small, as showed by the clustering patterns of 9 samples; and iii) when compared with controls, significant metabolic changes were observed for butanol-treated samples, especially for samples with 48 and 72 h butanol treatments. One of the butanol-treated biological replicates was slightly different from other two biological replicates at 48 h and 72 h, probably due to the fact the long-term butanol treatment has caused significant cell aggregation , which increased the sample heterogeneity. Nevertheless, the overall similar response patterns can still be observed in these replicate samples according to their position in the score plot (Figure 2). Using a cutoff ratio of 1.5 fold between butanol-treated and control samples, and change in at least 5 out of 9 replicate ratios in any time point, we determined 46 metabolites were differentially regulated, in which 35, 41 and 38 metabolites were detected in 24, 48 and 72 h, respectively (Table 5). Pattern analysis showed the 48 metabolites can be divided into at least 6 clusters according to their changes along the treatment time courses. For example, Cluster I included 7 metabolites up-regulated in all three time points, while Cluster II included 7 metabolites up-regulated only in 48 and 72 h after butanol exposure (Table 5).
Metabolomic analysis has identified several metabolites induced by butanol treatment, including 3-phosphoglycerate (3-PG) and glycerol 1-phosphate induced significantly in all three time points, serine induced at 24 and 48 h, and glycine induced at 48 and 72 h after butanol exposure, respectively. The findings were consistent with early studies which showed 3-phosphoglycerate is increasingly withdrawn from the Calvin cycle in S. elongatus PCC 7942 under iron limitation stress . In addition, phosphoglycerate kinase that catalyzes the production of 3-phosphoglycerate from 1,3-bisphosphoglycerate was also found induced in Anabaena sp. PCC7120 under arsenic stress . Moreover, early study has shown that the intracellular levels of organic acids (glyceric, glycolic and glyoxylic acids) and amino acids (glycine and serine) were elevated in salt-treated Anabaena sp. PCC 7120 as compared to those in the control cells . The results suggested that these metabolites could be important part of metabolic responses to both butanol and general environmental stresses.
Previous proteomic study found that a common stress response of Synechocystis under various environmental perturbations, irrespective of amplitude and duration, is the activation of atypical pathways for the acquisition of carbon and nitrogen from urea and arginine, as evidenced by the significant up-regulation of urease that converts urea into CO2 and ammonia, under most conditions . Our metabolomic analysis showed that urea was induced by butanol, especially at 48 and 72 h. Previous proteomic analysis showed that cyanophycinase, involved in the breakdown of cyanophycin, a storage molecule for excess carbon and nitrogen, into arginine and aspartic acid, was moderately up-regulated under several conditions . Arginine and aspartic acid can be further converted to glutamate and succinate, respectively . Metabolomic analysis showed that aspartic acid was significantly induced at all three time points, and succinic acid and L-glutamic acid were both induced at 48 and 72 h by butanol treatment. These results implied that a similar up-regulated degradation of cyanophycin may also occur under butanol stress.
Integrated transcriptomic and metabolomic analysis has been proposed as a powerful tool to build the relationship between information elements (i.e., genes/transcripts) and functional elements (i.e., metabolites) in cells [25, 52, 53]. In one recent study, integrated transcriptomic and metabolomic approach was used to determine the infection mechanism of Rhodococcus fascians into Arabidopsis thaliana. The transcriptomic analysis showed a significant impact of infection on the primary metabolism of the host, which was then confirmed by subsequent metabolite analysis, for example, invertase transcripts and activities strongly enhanced upon infection, may related to the increase in the hexose:sucrose ratio . In another study to compare the aerobic and anaerobic fermentations of Zymomonas mobilis, researchers found that greater amounts of end products such as acetate, lactate and acetoin were detected under aerobic conditions, while no change in terms of gene expression was found between aerobic and anaerobic conditions in the early exponential growth phase , implying the importance to applying integrated technology in uncovering related molecular mechanism. In this study, although only small number of metabolites can be chemically classified in Synechocystis, the metabolomic analysis found increased abundances of aspartic acid and serine, which was consistent with the induction of slr0550 encoding dihydrodipicolinate synthase involved in aspartate pathway, and sll0455 encoding homoserine dehydrogenase involved in serine pathway, respectively (Table 5). In addition, increased abundance of glutamic acid inside the cells was correlated with up-regulation of sll1883 encoding bifunctional ornithine acetyltransferase/N-acetylglutamate synthase protein, sll0461 encoding gamma-glutamyl phosphate reductase, slr0288 encoding glutamate--ammonia ligase, and slr1898 encoding acetylglutamate kinase that are involved in metabolism of glutamate family amino acids (Table 5). Moreover, metabolomic analysis showed the increased abundances of intermediates in the glycolysis pathway, such as glucose-6-P and 3-PG, consistent with the induction of two key genes, slr0752 encoding phosphopyruvate hydratase and sll0745 encoding 6-phosphofructokinase in the glycolysis pathway. Consistent with this result, up-regulation of glycolysis has been reported for various microbes under stress condition [56, 57]. In a recent 13C-based flux analysis, a thermophilic ethanol-tolerant Geobacillus thermoglucosidasius M10EXG was found to prefer glycolysis, the pentose phosphate pathway and the TCA cycle for glucose metabolism . On the other hand, for some of differentially regulated metabolites identified, such as urea and cyanophycin, no change was observed for their functionally-related genes in the transcriptomic datasets, which may be due to multiple factors, such as the snapshot nature of the analysis and the different stability of RNA molecules . Nevertheless, the results further demonstrated that transcriptomic and metabolomic technologies could be complementary to each other, allowing better decipherment of cellular responses of Synechocystis under butanol stress.
Validation of potential tolerance targets
Three genes, sll0690, slr0947 and slr1295 which were found induced by butanol exposure at all three time points (i.e. 24, 48 and 72 h) (Table 2), were selected for construction of knockout mutants and for validation of their involvement in butanol resistance. sll0690 encoding a probable transcription regulator, was up-regulated 5–6 folds, slr0947 encoding an OmpR-type DNA-binding response regulator, was up-regulated 2.4-5.5 folds, and slr1295 encoding an iron transport system substrate-binding protein was up-regulated 1.8-5.9 folds by butanol, respectively. Two corresponding proteins of the genes, Slr0947 and Slr1295, were identified in our previous proteomic analysis, in which they were also slightly up-regulated 1.16-1.55 and 1.16-1.57 folds after butanol treatment for 48 h, respectively . After confirmed by PCR and sequencing, the mutants were grown in parallel with wild type Synechocystis in both normal BG11 medium and the BG11 medium supplemented with 0.25% (v/v) butanol. Comparative analysis showed that although there is no visible difference in terms of growth patterns between the wild type and all three mutants in the regular BG11 medium (Figure 3A), gene disruption of sll0690, slr0947 and slr1295 led to increased butanol sensitivity, suggesting they were involved in butanol resistance (Figure 3B). Currently little is known how these genes are involved in butanol tolerance, although early studies have found that the slr1295 gene product, a periplasm-located component of an iron transporter, has a function in protecting photosystem (PS) II  and was induced under salt-stress condition ; and the slr0947 gene was involved in the regulation of the coupling of phycobilisomes to photosynthetic reaction centers, and reduction of the copy number of slr0947 resulted in decreased efficiency of energy transfer from phycobilisomes to photosystem II relative to photosystem I .
RNA-Seq based transcriptomics coupled with RT-PCR and GC-MS metabolomics were used to determine gene targets related to butanol tolerance in Synechocystis. Although the overall cellular responses revealed by transcriptomics and metabolomics were very similar to those revealed by our previous proteomic analysis, the genes/proteins involved in each type of responses were not always identical, consistent with recent conclusions that only a weak correlation exists between large-scale transcriptomic and proteomic datasets so that an integrative analysis of multiple levels of gene expression would be necessary and valuable . A comprehensive transcriptomic and metabolomic analysis with proteomic analysis led to identification of putative gene targets which may be involved in butanol tolerance. By constructing KO mutants and analyzing their butanol resistance, we validated three potential gene targets identified by the integrated OMICS approaches. In the future, once further functional characterization of these candidate genes completed, it is possible they can serve as target genes to engineer more robust butanol-tolerant cyanobacterial hosts.
Materials and methods
Bacterial growth conditions and butanol treatment
Synechocystis sp. PCC 6803 was grown in BG11 medium (pH 7.5) as described previously [26, 27]. Butanol of 0.20% (v/v) was added at the beginning of cultivation. Cells were collected by centrifugation at 8,000 × g for 10 min at 4°C.
RNA preparation and cDNA synthesis
Approximately 10 mg of cell pellets were frozen by liquid nitrogen immediately after centrifugation and cell walls were broken with mechanical cracking at low temperature. Cell pellets were then resuspended in Trizol reagent (Ambion, Austin, TX) and mixed well by vortex. Total RNA extraction was achieved using a miRNeasy Mini Kit (Qiagen, Valencia, CA). Contaminating DNA in RNA samples was removed with DNase I according to the instruction in the miRNeasy Mini Manual (Qiagen, Valencia, CA). The RNA quality and quantity were determined using Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA) and subjected to cDNA synthesis. The RNA integrity number (RIN) of every RNA sample used for sequencing was more than 8.0. For each sample, 500 ng total RNA were subjected to cDNA synthesis using a NuGEN Ovation® Prokaryotic RNA-Seq System according to manufacturer’s protocol (NuGEN, San Carlos, CA). The resulting double-stranded cDNA was purified using the MinElute Reaction Cleanup Kit (Qiagen, Valencia, CA).
RNA-seq library preparation
The double-stranded cDNA obtained was subjected to library preparation using the Illumina TruSeqTM RNA Sample Preparation Kit (Illumina, San Diego, CA), through a four-step protocol of end repairing, adenylate 3’ ends adding, adapter ligation, and cDNA template enrichment. To determine the quality of libraries, a Qubit® 2.0 Fluorometer and Qubit™ dsDNA HS (Invitrogen, Grand Island, NY) were first used to determine the DNA concentration of the libraries, and then FlashGel DNA Cassette (Lonza, USA) or Agilent Technologies 2100 Bioanalyzer (Agilent, Santa Clara, CA) was used to determine the product size of the libraries, with good libraries typically around 300 bp.
RNA 2 × 100 bp paired-end sequencing was performed using Illumina’s Solexa Genome Analyzer II using the standard protocol. The cDNA library of each sample was loaded to a single lane of an Illumina flow cell. The image deconvolution and calculation of quality value were performed using Goat module (Firecrest v.1.4.0 and Bustard v.1.4.0 programs) of Illumina pipeline v.1.4.
Transcriptomics data analysis
Sequence reads were pre-processed using FASTX Toolkit (v. 0.0.13) to remove low-quality bases, and reads shorter than 20 bp. The qualified sequence reads were then mapped to non-coding RNA (ncRNA) sequences using Bowtie (v. 2.0.0) with default settings. Genome sequences (including ncRNA sequences) and annotation information of Synechocystis sp. PCC 6803 were downloaded from NCBI (Downloaded on April 22, 2012) . Reads that mapped to ncRNA sequences were excluded from further analysis in this study. For paired-end Illumina reads, both pairs were removed if either pair mapped to rRNA. Remaining reads were mapped to the Synechocystis genome using Bowtie (v. 2.0.0) with the default parameters. For gene expression determination, we performed a standard calculation of Reads Per Kilobase of Gene per Million Mapped Reads (RPKM) . We performed comparative transcriptome analysis for all three time points (i.e., 24, 48 and 72 h). To identify the reliable gene targets related to butanol tolerance, only the genes with 1.5-fold induction by butanol at all three time points were regarded as up-regulated genes.
Quantitative real-time RT-PCR analysis
The identical RNA samples used for transcriptomics analysis as described above were used for RT-qPCR analysis. cDNAs were synthesized using RevertAidTM Reverse Transcriptase (Fermentas, Glen Burnie, MD). The qPCR reaction was carried out in 20 μl reactions containing 10 μl of SYBR® Green PCR Master Mix (Applied Biosystems, Foster City, CA), and 2 μl of each PCR primer at 2 mM, employing the StepOnePlus™ Real-Time PCR System (Applied Biosystems, Foster City, CA), under the following condition: 50°C for 2 min and 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. Quantification of gene expression was determined according to standard process of RT-PCR which used serial dilutions of known concentration of chromosome DNA as template to make a standard curve. The rnpB gene (6803 s01) encoding RNase P subunit B was used as an internal control according to the previous publication . Three technical replicates were performed for each gene. Data analysis was carried out using the StepOnePlus analytical software (Applied Biosystems, Foster City, CA). Data was presented as ratios of the amount of normalized transcript in the treatment to that from the control. The gene ID and their related primer sequences used for real-time RT-PCR analysis were listed in Additional file 1: Table S1.
GC-MS based metabolomics analysis
All chemicals used for metabolome isolation and GC/MS analysis were obtained from Sigma-Aldrich (Taufkirchen, Germany). Cells were collected from control and butanol-treated (0.2% v/v) cultures at 24, 48 and 72 h, respectively. Three biological replicates were established for each sample, and every sample was analyzed three times. For each sample, cells from 5 to 20 mL culture, equivalent to 108 cells mL-1, were collected by centrifugation at 8000 × g for 10 min at 4°C (Eppendorf 5430R, Hamburg, Germany). The cell pellets were frozen in liquid nitrogen and then stored at −80°C before use. i) Metabolome extraction: cells were re-suspended in 1 mL cold 10:3:1 (v/v/v) methanol: chloroform: H2O solution (MCW), and frozen in liquid nitrogen and thawed for five times. Supernatants were collected by centrifugation at 14,000 × g for 3 min at 4°C (Eppendorf 5430R, Hamburg, Germany). To normalize variations across samples, an internal standard (IS) solution (100 μg/mL U- 13C-sorbitol, 5 μL) was added to 100 μL supernatant in a 1.5-mL microtube before it was dried by vacuum centrifugation for 2–3 h (4°C). ii) Sample derivatization: derivatization was conducted according to the two-stage technique by Roessner et al. (2001) . The samples were dissolved in 10 μL methoxyamine hydrochloride (40 mg/mL in pyridine) and shaken at 30°C for 90 min, then were added with 90 μL N-methyl- N -(trimethylsilyl) trifluoroacetamide (MSTFA) and incubated at 37°C for 30 min to trimethylsilylate the polar functional groups. The derivate samples were collected by centrifugation at 14,000 × g for 3 min before GC/MS analysis. iii) GC-MS analysis: sample analysis was performed on a GC-MS system-GC 7890 coupled to an MSD 5975 (Agilent Technologies, Inc., Santa Clara, CA, USA) equipped with a HP-5MS capillary column (30 m × 250 mm id). 2 μL derivatized sample was injected in splitless mode at 230°C injector temperature. The GC was operated at constant flow of 1 mL/min helium. The temperature program started isocratic at 45°C for 2 min, followed by temperature ramping of 5°C/ min to a final temperature of 280°C, and then held constant for additional 2 min. The range of mass scan was m/z 38–650. iv) Data processing and statistical analysis: The mass fragmentation spectrum was analyzed using the Automated Mass Spectral Deconvolution and Identification System (AMDIS)  to identify the compounds by matching the data with Fiehn Library  and the mass spectral library of the National Institute of Standards and Technology (NIST). Peak areas of all identified metabolites were normalized against the internal standard and the acquired relative abundances for each identified metabolite were used for future data analysis.
All metabolomic profile data was first normalized by internal control and cell numbers, and then subjected to Principal Component Analysis using software SIMCA-P 11.5 . Differentially regulated metabolites were determined using a threshold of fold change greater than 1.5 between butanol-treated samples and controls. For each time point, three biological replicates of butanol-treated samples were compared with three biological replicates of control, generating 9 ratios. For each ratio, r > 1.5 was assigned as “ + 1”, r < −1.5 as “-1”, and −1.5 < r < 1.5 as “0”. The sums of the nine ratios for each metabolite at any time point were provided in Table 5.
Construction and analysis of knockout mutants
A fusion PCR based method was employed for the construction of gene knockout fragments . Briefly, for the gene target selected, three sets of primers were designed to amplify a linear DNA fragment containing the chloramphenicol resistance cassette (amplified from a plasmid pACYC184) with two flanking arms of DNA upstream and downstream of the targeted gene. The linear fused PCR amplicon was used directly for transformation into Synechocystis by natural transformation. The chloramphenicol-resistant transformants were obtained and passed several times on fresh BG11 plates supplemented with 10 μg/ml chloramphenicol to achieve complete chromosome segregation. Three genes, sll0690, slr0947 and slr1295 that have been found differentially regulated by butanol exposure, were selected for construction of gene knockout mutants. PCR primers for mutant construction and validation were listed in Additional file 1: Table S1. Full segregation for sll0690 and slr1295 genes was confirmed by PCR. For Δslr0947 mutant, we found that it contained trace amount of original wild-type band in the DNA gels even after more than ten passages, it may worth further investigation whether slr0947 is a lethal gene for the condition. Comparative growth analysis of the wild type 6803 and the mutants were performed in 100-mL flasks each with 10 mL BG11 medium with or without 0.25% (v/v) butanol. Cultivation conditions are the same as described above. Growth analysis was performed in biological triplicates.
Gas chromatography–Mass spectrometry
Isobaric tag for relative and absolute quantitation
Liquid chromatography-tandem mass spectrometry
Principal component analysis
Reads per kilobase of gene per million mapped reads
Reactive oxygen species
Durre P: Biobutanol: an attractive biofuel. Biotechnol J. 2007, 2: 1525-1534. 10.1002/biot.200700168.
Durre P: Fermentative butanol production: bulk chemical and biofuel. Ann NY Acad Sci. 2008, 1125: 353-362. 10.1196/annals.1419.009.
Durre P, Fischer RJ, Kuhn A, Lorenz K, Schreiber W, Sturzenhofecker B, Ullmann S, Winzer K, Sauer U: Solventogenic enzymes of Clostridium acetobutylicum: catalytic properties, genetic organization, and transcriptional regulation. FEMS Microbiol Rev. 1995, 17: 251-262.
Lee SY, Park JH, Jang SH, Nielsen LK, Kim J, Jung KS: Fermentative butanol production by Clostridia. Biotechnol Bioeng. 2008, 101: 209-228. 10.1002/bit.22003.
Ezeji TC, Qureshi N, Blaschek H: Acetone–butanol–ethanol production from concentrated substrate: reduction in substrate inhibition by fed-batch technique and product inhibition by gas stripping. Appl Micriobiol Biotechnol. 2004, 63: 653-658. 10.1007/s00253-003-1400-x.
Ezeji TC, Qureshi N, Blaschek HP: Butanol fermentation research: upstream and downstream manipulations. Chem Rec. 2004, 4: 305-314. 10.1002/tcr.20023.
Zheng YN, Li LZ, Xian M, Ma YJ, Yang JM, Xu X, He DZ: Problems with the microbial production of butanol. J Ind Microbiol Biotechnol. 2009, 36: 1127-1138. 10.1007/s10295-009-0609-9.
Lutke-Eversloh T, Bahl H: Metabolic engineering of Clostridium acetobutylicum: recent advances to improve butanol production. Curr Opin Biotechnol. 2011, 22: 634-647. 10.1016/j.copbio.2011.01.011.
Zhu XG, Long SP, Ort DR: Improving photosynthetic efficiency for greater yield. Annu Rev Plant Biol. 2010, 61: 235-261. 10.1146/annurev-arplant-042809-112206.
Machado IM, Atsumi S: Cyanobacterial biofuel production. J Biotechnol. 2012, 162: 50-56. 10.1016/j.jbiotec.2012.03.005.
Atsumi S, Higashide W, Liao J: Direct photosynthetic recycling of carbon dioxide to isobutyraldehyde. Nat Biotechnol. 2009, 27: 1177-1180. 10.1038/nbt.1586.
Lan E, Liao J: Metabolic engineering of cyanobacteria for 1-butanol production from carbon dioxide. Metab Eng. 2011, 13: 353-363. 10.1016/j.ymben.2011.04.004.
Atsumi S, Cann AF, Connor MR, Shen CR, Smith KM, Brynildsen MP, Chou KJ, Hanai T, Liao JC: Metabolic engineering of Escherichia coli for 1-butanol production. Metab Eng. 2008, 10: 305-311. 10.1016/j.ymben.2007.08.003.
Inui M, Suda M, Kimura S, Yasuda K, Suzuki H, Toda H, Yamamoto S, Okino S, Suzuki N, Yukawa H: Expression of Clostridium acetobutylicum butanol synthetic genes in Escherichia coli. Appl Microbiol Biotechnol. 2008, 77: 1305-1316. 10.1007/s00253-007-1257-5.
Shen CR, Liao JC: Metabolic engineering of Escherichia coli for 1-butanol and 1-propanol production via the keto-acid pathways. Metab Eng. 2008, 10: 312-320. 10.1016/j.ymben.2008.08.001.
Nicolaou SA, Gaida SM, Papoutsakis ET: A comparative view of metabolite and substrate stress and tolerance in microbial bioprocessing: From biofuels and chemicals, to biocatalysis and bioremediation. Metab Eng. 2010, 12: 307-331. 10.1016/j.ymben.2010.03.004.
Dunlop MJ: Engineering microbes for tolerance to next-generation biofuels. Biotechnol Biofuels. 2011, 4: 32-10.1186/1754-6834-4-32.
Papoutsakis ET: Engineering solventogenic clostridia. Curr Opin Biotechnol. 2008, 19: 420-429. 10.1016/j.copbio.2008.08.003.
Ezeji T, Milne C, Price ND, Blaschek HP: Achievements and perspectives to overcome the poor solvent resistance in acetone and butanol-producing microorganisms. Appl Microbiol Biotechnol. 2010, 85: 1697-1712. 10.1007/s00253-009-2390-0.
Liyanage H, Young M, Kashket ER: Butanol tolerance of Clostridium beijerinckii NCIMB 8052 associated with down-regulation of gldA by antisense RNA. J Mol Microbiol Biotechnol. 2000, 2: 87-93.
Jia K, Zhang Y, Li Y: Identification and Characterization of two functionally unknown genes involved in butanol tolerance of Clostridium acetobutylicum. PLoS One. 2012, 7: e38815-10.1371/journal.pone.0038815.
Zhang H, Chong H, Ching CB, Song H, Jiang R: Engineering global transcription factor cyclic AMP receptor protein of Escherichia coli for improved 1-butanol tolerance. Appl Microbiol Biotechnol. 2012, 94: 1107-1117. 10.1007/s00253-012-4012-5.
Atsumi S, Wu TY, Machado IMP, Huang WC, Chen PY, Pellegrini M, Liao JC: Evolution, genomic analysis, and reconstruction of isobutanol tolerance in Escherichia coli. Mol Syst Biol. 2010, 6: 449-
Reyes LH, Almario MP, Kao KC: Genomic library screens for genes involved in n-butanol tolerance in Escherichia coli. PLoS One. 2011, 6 (3): e17678-10.1371/journal.pone.0017678.
Zhang W, Li F, Nie L: Integrating multiple ‘omics’ analysis for microbial biology: application and methodologies. Microbiology. 2010, 156: 287-301. 10.1099/mic.0.034793-0.
Tian XX, Chen L, Wang JX, Qiao J, Zhang WW: Quantitative proteomics reveals dynamic responses of Synechocystis sp. PCC 6803 to next-generation biofuel butanol. J Proteomics. 2012, 78: 326-345.
Kaneko T, Nakamura Y, Sasamoto S, Watanabe A, Kohara M, Matsumoto M, Shimpo S, Yamada M, Tabata S: Structural analysis of four large plasmids harboring in a unicellular cyanobacterium, Synechocystis sp. PCC 6803. DNA Res. 2003, 10: 221-228. 10.1093/dnares/10.5.221.
Noone D, Howell A, Collery R, Devine KM: YkdA and YvtA, HtrA-like serine proteases in Bacillus subtilis, engage in negative autoregulation and reciprocal cross-regulation of ykdA and yvtA gene expression. J Bacteriol. 2001, 183: 654-663. 10.1128/JB.183.2.654-663.2001.
Snider J, Gutsche I, Lin M, Baby S, Cox B, Butland G, Greenblatt J, Emili A, Houry WA: Formation of a distinctive complex between the inducible bacterial lysine decarboxylase and a novel AAA + ATPase. J Biol Chem. 2006, 281: 1532-1546.
Dieppedale J, Sobral D, Dupuis M, Dubail I, Klimentova J, Stulik J, Postic G, Frapy E, Meibom KL, Barel M, Charbit A: Identification of a putative chaperone involved in stress resistance and virulence in Francisella tularensis. Infect Immun. 2011, 79: 1428-1439. 10.1128/IAI.01012-10.
Rutherford BJ, Dahl RH, Price RE, Szmidt HL, Benke PI, Mukhopadhyay A, Keasling JD: Functional genomic study of exogenous n-butanol stress in Escherichia coli. Appl Environ Microbiol. 2010, 76: 1935-1945. 10.1128/AEM.02323-09.
Das PK, Bagchi SN: Role of bacterioferritin comigratory protein and glutathione peroxidase-reductase system in promoting bentazone tolerance in a mutant of Synechococcus elongatus PCC7942. Protoplasma. 2012, 249: 65-74. 10.1007/s00709-011-0262-9.
Kim HG, Kim BC, Park EH, Lim CJ: Stress-dependent regulation of a monothiol glutaredoxin gene from Schizosaccharomyces pombe. Can J Microbiol. 2005, 51: 613-620. 10.1139/w05-034.
Pujol-Carrion N, De la Torre-Ruiz MA A: Glutaredoxins Grx4 and Grx3 of Saccharomyces cerevisiae play a role in actin dynamics through their Trx domains, which contributes to oxidative stress resistance. Appl Environ Microbiol. 2010, 76: 7826-7835. 10.1128/AEM.01755-10.
Qiao JJ, Wang JX, Chen L, Tian XX, Zhang WW: Quantitative iTRAQ LC-MS/MS proteomics reveals metabolic response to biofuel ethanol in cyanobacterial Synechocystis sp. PCC 6803. J Proteome Res. 2012, 11 (11): 5286-5300. 10.1021/pr300504w.
Zhang W, Culley DE, Hogan M, Vitirit L, Brockman FJ: Oxidative stress and heat-shock responses in Desulfovibrio vulgaris by genome-wide transcriptomic analysis. Antonie Van Leeuwenhoek. 2006, 90: 41-55. 10.1007/s10482-006-9059-9.
Moskovitz J, Berlett B, Poston JM, Stadtman ER: The yeast peptide-methionine sulfoxide reductase functions as an antioxidant in vivo. Proc Natl Acad Sci USA. 1997, 94: 9585-9589. 10.1073/pnas.94.18.9585.
Moskovitz J, Lescher E, Berlett BS, Azare J, Poston JM, Stadtman ER: Over-expression of peptide-methionine sulfoxide reductase in Saccharomyces cerevisiae and human T cells provides them with high resistance to oxidative stress. Proc Natl Acad Sci USA. 1998, 95: 14071-14075. 10.1073/pnas.95.24.14071.
Truong-Bolduc QC, Ding Y, Hooper DC: Posttranslational modification influences the effects of MgrA on norA expression in Staphylococcus aureus. J Bacteriol. 2008, 190: 7375-7381. 10.1128/JB.01068-08.
Economou A: Bacterial protein translocase: a unique molecular machine with an army of substrates. FEBS Lett. 2000, 476: 18-21. 10.1016/S0014-5793(00)01662-8.
Cline K, McCaffery M: Evidence for a dynamic and transient pathway through the TAT protein transport machinery. EMBO J. 2007, 26: 3039-3049. 10.1038/sj.emboj.7601759.
von Wobeser E, Ibelings BW, Bok J, Krasikov V, Huisman J, Matthijs HCP: Concerted changes in gene expression and Cell physiology of the cyanobacterium Synechocystis sp. Strain PCC 6803 during transitions between nitrogen and light-limited growth. Plant Physiol. 2011, 155: 1445-1457. 10.1104/pp.110.165837.
Kadouri D, Jurkevitch E, Okon Y, Castro-Sowinski S: Ecological and agricultural significance of bacterial polyhydroxyalkanoates. Crit Rev Microbiol. 2005, 31: 55-67. 10.1080/10408410590899228.
Kolodny NH, Bauer D, Bryce K, Klucevsek K, Lane A, Medeiros L, Mercer W, Moin S, Park D, Petersen J, Wright J, Yuen C, Wolfson AJ, Allen MM: Effect of nitrogen source on cyanophycin synthesis in Synechocystis sp. strain PCC 6308. J Bacteriol. 2006, 188: 934-940. 10.1128/JB.188.3.934-940.2006.
Schäfer L, Sandmann M, Woitsch S, Sandmann G: Coordinate up-regulation of carotenoid biosynthesis as a response to light stress in Synechococcus PCC7942. Plant Cell Environ. 2006, 29: 1349-1356. 10.1111/j.1365-3040.2006.01515.x.
Sandesh Kamath B, Vidhyavathi R, Sarada R, Ravishankar GA: Enhancement of carotenoids by mutation and stress induced carotenogenic genes in Haematococcus pluvialis mutants. Bioresour Technol. 2008, 99: 8667-8673. 10.1016/j.biortech.2008.04.013.
Nodop A, Pietsch D, Höcker R, Becker A, Pistorius EK, Forchhammer K, Michel KP: Transcript profiling reveals new insights into the acclimation of the mesophilic fresh-water cyanobacterium Synechococcus elongatus PCC 7942 to iron starvation. Plant Physiol. 2008, 147: 747-763. 10.1104/pp.107.114058.
Pandey S, Rai R, Rai LC: Proteomics combines morphological, physiological and biochemical attributes to unravel the survival strategy of Anabaena sp. PCC7120 under arsenic stress. J Proteomics. 2012, 75: 921-937. 10.1016/j.jprot.2011.10.011.
Srivastava AK, Alexova R, Jeon YJ, Kohli GS, Neilan BA: Assessment of salinity-induced photorespiratory glycolate metabolism in Anabaena sp. PCC 7120. Microbiology. 2011, 157: 911-917. 10.1099/mic.0.045682-0.
Wegener KM, Singh AK, Jacobs JM, Elvitigala T, Welsh EA, Keren N, Gritsenko MA, Ghosh BK, Camp DG, Smith RD, Pakrasi HB: Global proteomics reveal an atypical strategy for carbon/nitrogen assimilation by a cyanobacterium under diverse environmental perturbations. Mol Cell Proteomics. 2010, 9: 2678-2689. 10.1074/mcp.M110.000109.
Quintero MJ, Muro-Pastor AM, Herrero A, Flores E: Arginine catabolism in the cyanobacterium Synechocystis sp. Strain PCC 6803 involves the urea cycle and arginase pathway. J Bacteriol. 2000, 182: 1008-1015. 10.1128/JB.182.4.1008-1015.2000.
Urbanczyk-Wochniak E, Luedemann A, Kopka J, Selbig J, Roessner-Tunali U, Willmitzer L, Fernie AR: Parallel analysis of transcript and metabolic profiles: a new approach in systems biology. EMBO Rep. 2003, 4: 989-993. 10.1038/sj.embor.embor944.
Urbanczyk-Wochniak E, Willmitzer L, Fernie AR: Integrating profiling data: using linear correlation to reveal coregulation of transcript and metabolites. Methods Mol Biol. 2007, 358: 77-85. 10.1007/978-1-59745-244-1_5.
Depuydt S, Trenkamp S, Fernie AR, Elftieh S, Renou JP, Vuylsteke M, Holsters M, Vereecke D: An integrated genomics approach to define niche establishment by Rhodococcus fascians. Plant Physiol. 2009, 149: 1366-1386. 10.1104/pp.108.131805.
Yang S, Tschaplinski TJ, Engle NL, Carroll SL, Martin SL, Davison BH, Palumbo AV, Rodriguez M, Brown SD: Transcriptomic and metabolomic profiling of Zymomonas mobilis during aerobic and anaerobic fermentations. BMC Genomics. 2009, 10: 34-10.1186/1471-2164-10-34.
Chaussee MA, McDowell EJ, Rieck LD, Callegari EA, Chaussee MS: Proteomic analysis of a penicillin-tolerant rgg mutant strain of Streptococcus pyogenes. J Antimicrob Chemother. 2006, 58: 752-759. 10.1093/jac/dkl319.
Guo W, Hao H, Dai M, Wang Y, Huang L, Peng D, Wang X, Wang H, Yao M, Sun Y, Liu Z, Yuan Z: Development of quinoxaline 1, 4-dioxides resistance in Escherichia coli and molecular change under resistance selection. PLoS One. 2012, 7: e43322-10.1371/journal.pone.0043322.
Tang YJ, Sapra R, Joyner D, Hazen TC, Myers S, Reichmuth D, Blanch H, Keasling JD: Analysis of metabolic pathways and fluxes in a newly discovered thermophilic and ethanol-tolerant Geobacillus strain. Biotechnol Bioeng. 2009, 102: 1377-1386. 10.1002/bit.22181.
Tölle J, Michel KP, Kruip J, Kahmann U, Preisfeld A, Pistorius EK: Localization and function of the IdiA homologue Slr1295 in the cyanobacterium Synechocystis sp. strain PCC 6803. Microbiology. 2002, 148: 3293-3305.
Huang F, Fulda S, Hagemann M, Norling B: Proteomic screening of salt-stress-induced changes in plasma membranes of Synechocystis sp. strain PCC 6803. Proteomics. 2006, 6: 910-920. 10.1002/pmic.200500114.
Ashby MK, Mullineaux CW: Cyanobacterial ycf27 gene products regulate energy transfer from phycobilisomes to photosystems I and II. FEMS Microbiol Lett. 1999, 181: 253-260. 10.1111/j.1574-6968.1999.tb08852.x.
Nie L, Wu G, Culley DE, Scholten JC, Zhang W: Integrative analysis of transcriptomic and proteomic data: challenges, solutions and applications. Crit Rev Biotechnol. 2007, 27: 63-75. 10.1080/07388550701334212.
Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-10.1186/gb-2010-11-10-r106.
Kloft N, Rasch G, Forchhammer K: Protein phosphatase PphA from Synechocystis sp. PCC 6803: the physiological framework of PII-P dephosphorylation. Microbiology. 2005, 151: 1275-1283. 10.1099/mic.0.27771-0.
Roessner U, Luedemann A, Brust D, Fiehn O, Linke T, Willmitzer L, Fernie A: Metabolic profiling allows comprehensive phenotyping of genetically or environmentally modified plant systems. Plant Cell. 2001, 13: 11-29.
Stein SE: An integrated method for spectrum extraction and compound identification from GC/MS data. Mass Spectrom. 1999, 10 (8): 770-781.
Fiehn O: Metabolomics-the link between genotypes and phenotypes. Plant Mol Biol. 2002, 48: 155-171. 10.1023/A:1013713905833.
Laiakis EC, Morris GA, Fornace AJ, Howie SR: Metabolomic analysis in severe childhood pneumonia in the Gambia, West Africa: findings from a pilot study. PLoS One. 2010, 5 (9): e12655-10.1371/journal.pone.0012655.
Wang HL, Postier BL, Burnap RL: Optimization of fusion PCR for in vitro construction of gene knockout fragments. Biotechniques. 2002, 33 (26): 28-30.
The research was supported by grants from National Basic Research Program of China (“973” program, project No. 2011CBA00803 and No. 2012CB721101) and National High-tech R&D Program (“863” program, project No. 2012AA02A707). The authors would also like to thank Tianjin University and the “985 Project” of Ministry of Education for their generous supports in establishing the research laboratory.
The authors declare no competing interests.
XR, ZS, LC and JW carried out cultivation and transcriptomics analysis. LC, XR, ZS and JL carried out the RT-PCR, mutant construction and phenotypic analysis. HZ, LC, JW, JQ and WZ finished the statistical analysis for transcripomic data. MS, XT, LC and WZ carried out the metabolomics experiments and analysis. LC, JW and WZ conceived of the study, participated in its design and coordination. HZ, JW, LC and WZ drafted the manuscript. All authors read and approved the final manuscript.