Metabolomic analysis of the salt-sensitive mutants reveals changes in amino acid and fatty acid composition important to long-term salt stress in Synechocystis sp. PCC 6803

Mechanisms on salt tolerant of Paenibacillus polymyxa SC2 and its growth-promoting effects on maize seedlings under saline conditions

Soil salinity negatively affects microbial communities, soil fertility, and agricultural productivity and has become a major agricultural problem worldwide. Plant growth-promoting rhizobacteria (PGPR) with salt tolerance can benefit plant growth under saline conditions and diminish the negative effects of salt stress on plants. In this study, we aimed to understand the salt-tolerance mechanism of Paenibacillus polymyxa at the genetic and metabolic levels and elucidate the mechanism of strain SC2 in promoting maize growth under saline conditions. Under salt stress, we found that strain SC2 promoted maize seedling growth, which was accompanied by a significant upregulation of genes encoding for the biosynthesis of peptidoglycan, polysaccharide, and fatty acid, the metabolism of purine and pyrimidine, and the transport of osmoprotectants such as trehalose, glycine betaine, and K+ in strain SC2. To further enhance the salt resistance of strain SC2, three mutants (SC2-11, SC2-13, and SC2-14) with higher capacities for salt resistance and exopolysaccharide synthesis were obtained via atmospheric and room-temperature plasma mutagenesis. In saline-alkaline soil, the mutants showed better promoting effect on maize seedlings than wild-type SC2. The fresh weight of maize seedlings was increased by 68.10% after treatment with SC2-11 compared with that of the control group. The transcriptome analysis of maize roots demonstrated that SC2 and SC2-11 could induce the upregulation of genes related to the plant hormone signal transduction, starch and sucrose metabolism, reactive oxygen species scavenging, and auxin and ethylene signaling under saline-alkaline stress. In addition, various transcription factors, such as zinc finger proteins, ethylene-responsive-element-binding protein, WRKY, myeloblastosis proteins, basic helix-loop-helix proteins, and NAC proteins, were up-regulated in response to abiotic stress. Moreover, the microbial community composition of maize rhizosphere soil after inoculating with strain SC2 was varied from the one after inoculating with mutant SC2-11. Our results provide new insights into the various genes involved in the salt resistance of strain SC2 and a theoretical basis for utilizing P. polymyxa in saline-alkaline environments.

Identification of adult resistant genes to stripe rust in wheat from southwestern China based on GWAS and WGCNA analysis

KEY MESSAGE

In this study, genome-wide association studies combined with transcriptome data analysis were utilized to reveal potential candidate genes for stripe rust resistance in wheat, providing a basis for screening wheat varieties for stripe rust resistance. Wheat stripe rust, which is caused by the wheat stripe rust fungus (Puccinia striiformis f. sp. tritici, Pst) is one of the world's most devastating diseases of wheat. Genetic resistance is the most effective strategy for controlling diseases. Although wheat stripe rust resistance genes have been identified to date, only a few of them confer strong and broad-spectrum resistance. Here, the resistance of 335 wheat germplasm resources (mainly wheat landraces) from southwestern China to wheat stripe rust was evaluated at the adult stage. Combined genome-wide association study (GWAS) and weighted gene co-expression network analysis (WGCNA) based on RNA sequencing from stripe rust resistant accession Y0337 and susceptible accession Y0402, five candidate resistance genes to wheat stripe rust (TraesCS1B02G170200, TraesCS2D02G181000, TraesCS4B02G117200, TraesCS6A02G189300, and TraesCS3A02G122300) were identified. The transcription level analyses showed that these five genes were significantly differentially expressed between resistant and susceptible accessions post inoculation with Pst at different times. These candidate genes could be experimentally transformed to validate and manipulate fungal resistance, which is beneficial for the development of the wheat cultivars resistant to stripe rust.

Integrated Comparative Transcriptome and Weighted Gene Co-Expression Network Analysis Provide Valuable Insights into the Mechanisms of Pinhead Initiation in Chinese Caterpillar Mushroom Ophiocordyceps sinensis (Ascomycota)

The initiation and formation of the "pinhead" is the key node in growth process of Ophiocordyceps sinensis (Chinese Cordyceps). The research on the mechanism of changes in this growth stage is the basis for realizing the industrialization of its artificial cultivation. Clarifying the mechanisms of pinhead initiation is essential for its further application. Here, we performed a comprehensive transcriptome analysis of pinhead initiation process in O. sinensis. Comparative transcriptome analysis revealed remarkable variation in gene expression and enriched pathways at different pinhead initiation stages. Gene co-expression network analysis by WGCNA identified 4 modules highly relevant to different pinhead initiation stages, and 23 hub genes. The biological function analysis and hub gene annotation of these identified modules demonstrated that transmembrane transport and nucleotide excision repair were the topmost enriched in pre-pinhead initiation stage, carbohydrate metabolism and protein glycosylation were specially enriched in pinhead initiation stage, nucleotide binding and DNA metabolic process were over-represented after pinhead stage. These key regulators are mainly involved in carbohydrate metabolism, synthesis of proteins and nucleic acids. This work excavated the candidate pathways and hub genes related to the pinhead initiation stage, which will serve as a reference for realizing the industrialization of artificial cultivation in O. sinensis.

Weighted gene co-expression network analysis identifies genes related to HG Type 0 resistance and verification of hub gene GmHg1

INTRODUCTION

The soybean cyst nematode (SCN) is a major disease in soybean production thatseriously affects soybean yield. At present, there are no studies on weighted geneco-expression network analysis (WGCNA) related to SCN resistance.

METHODS

Here, transcriptome data from 36 soybean roots under SCN HG Type 0 (race 3) stresswere used in WGCNA to identify significant modules.

RESULTS AND DISCUSSION

A total of 10,000 differentially expressed genes and 21 modules were identified, of which the module most related to SCN was turquoise. In addition, the hub gene GmHg1 with high connectivity was selected, and its function was verified. GmHg1 encodes serine/threonine protein kinase (PK), and the expression of GmHg1 in SCN-resistant cultivars ('Dongnong L-204') and SCN-susceptible cultivars ('Heinong 37') increased significantly after HG Type 0 stress. Soybean plants transformed with GmHg1-OX had significantly increased SCN resistance. In contrast, the GmHg1-RNAi transgenic soybean plants significantly reduced SCN resistance. In transgenic materials, the expression patterns of 11 genes with the same expression trend as the GmHg1 gene in the 'turquoise module' were analyzed. Analysis showed that 11genes were co-expressed with GmHg1, which may be involved in the process of soybean resistance to SCN. Our work provides a new direction for studying the Molecular mechanism of soybean resistance to SCN.

Comparative transcriptome and WGCNA reveal key genes involved in lignocellulose degradation in Sarcomyxa edulis

The developmental transcriptomes of Sarcomyxa edulis were assessed to explore the molecular mechanisms underlying lignocellulose degradation. Six stages were analyzed, spanning the entire developmental process: growth of mycelium until occupying half the bag (B1), mycelium under low-temperature stimulation after occupying the entire bag (B2), appearance of mycelium in primordia (B3), primordia (B4), mycelium at the harvest stage (B5), and mature fruiting body (B6). Samples from all six developmental stages were used for transcriptome sequencing, with three biological replicates for all experiments. A co-expression network of weighted genes associated with extracellular enzyme physiological traits was constructed using weighted gene co-expression network analysis (WGCNA). We obtained 19 gene co-expression modules significantly associated with lignocellulose degradation. In addition, 12 key genes and 8 kinds of TF families involved in lignocellulose degradation pathways were discovered from the four modules that exhibited the highest correlation with the target traits. These results provide new insights that advance our understanding of the molecular genetic mechanisms of lignocellulose degradation in S. edulis to facilitate its utilization by the edible mushroom industry.

WGCNA Analysis Identifies the Hub Genes Related to Heat Stress in Seedling of Rice (Oryza sativa L.)

Frequent high temperature weather affects the growth and development of rice, resulting in the decline of seed–setting rate, deterioration of rice quality and reduction of yield. Although some high temperature tolerance genes have been cloned, there is still little success in solving the effects of high temperature stress in rice (Oryza sativa L.). Based on the transcriptional data of seven time points, the weighted correlation network analysis (WGCNA) method was used to construct a co–expression network of differentially expressed genes (DEGs) between the rice genotypes IR64 (tolerant to heat stress) and Koshihikari (susceptible to heat stress). There were four modules in both genotypes that were highly correlated with the time points after heat stress in the seedling. We further identified candidate hub genes through clustering and analysis of protein interaction network with known–core genes. The results showed that the ribosome and protein processing in the endoplasmic reticulum were the common pathways in response to heat stress between the two genotypes. The changes of starch and sucrose metabolism and the biosynthesis of secondary metabolites pathways are possible reasons for the sensitivity to heat stress for Koshihikari. Our findings provide an important reference for the understanding of high temperature response mechanisms and the cultivation of high temperature resistant materials.

Metabolic Responses of a Model Green Microalga Euglena gracilis to Different Environmental Stresses

Euglena gracilis, a green microalga known as a potential candidate for jet fuel producers and new functional food resources, is highly tolerant to antibiotics, heavy metals, and other environmental stresses. Its cells contain many high-value products, including vitamins, amino acids, pigments, unsaturated fatty acids, and carbohydrate paramylon as metabolites, which change contents in response to various extracellular environments. However, mechanism insights into the cellular metabolic response of Euglena to different toxic chemicals and adverse environmental stresses were very limited. We extensively investigated the changes of cell biomass, pigments, lipids, and paramylon of E. gracilis under several environmental stresses, such as heavy metal CdCl₂, antibiotics paromomycin, and nutrient deprivation. In addition, global metabolomics by Ultra-high-performance liquid chromatography tandem mass spectrometry (UHPLC-MS/MS) was applied to study other metabolites and potential regulatory mechanisms behind the differential accumulation of major high-valued metabolites. This study collects a comprehensive update on the biology of E. gracilis for various metabolic responses to stress conditions, and it will be of great value for Euglena cultivation and high-value [154mm][10mm]Q7metabolite production.

New insights into the function of the proteins IsiC and IsiD from Synechocystis sp. PCC 6803 under iron limitation

Iron is a common cofactor in biological processes such as respiration, photosynthesis, and nitrogen fixation. The genes isiC and isiD encode unknown proteins, and the growth of ΔisiC and ΔisiD mutants is inhibited under iron-deficient conditions. To study the regulatory mechanisms of IsiC and IsiD during iron starvation, we carried out transcriptome and metabolome sequencing. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the photosynthesis, nitrogen metabolism, and ABC transporter pathways play a vital role in regulating iron deficiency. Upon iron repletion, IsiC and IsiD also have a regulatory effect on these pathways. Additionally, KEGG analysis of the differential metabolites of wild type (WT) and mutants showed that they were all enriched in starch and sucrose metabolism after iron limitation. Weighted gene co-expression network analysis (WGCNA) constructed a co-expression network of differentially expressed genes with phenotypes and metabolites, and finally identified five modules. The turquoise module was positively correlated with iron deficiency. In contrast, the WT and blue module exhibited a negative correlation, and the mutants ΔisiC and ΔisiD were positively correlated with the gray and brown modules, respectively. WGCNA also analyzed the relationship between metabolites and phenotypes, and the green module was related to iron starvation. The co-expression network determined the hub genes and metabolites of each module. This study lays a foundation for a better understanding of cyanobacteria in response to iron deficiency. KEY POINTS: • Nitrogen metabolism and ABC transporters are involved in iron regulation. • Starch and sucrose metabolism is related to the regulation of iron deficiency. • WGCNA analyzes the correlation between genes and metabolites.

Salt-Tolerant Synechococcus elongatus UTEX 2973 Obtained via Engineering of Heterologous Synthesis of Compatible Solute Glucosylglycerol

The recently isolated cyanobacterium Synechococcus elongatus UTEX 2973 (Syn2973) is characterized by a faster growth rate and greater tolerance to high temperature and high light, making it a good candidate chassis for autotrophic photosynthetic microbial cell factories. However, Syn2973 is sensitive to salt stress, making it urgently important to improve the salt tolerance of Syn2973 for future biotechnological applications. Glucosylglycerol, a compatible solute, plays an important role in resisting salt stress in moderate and marine halotolerant cyanobacteria. In this study, the salt tolerance of Syn2973 was successfully improved by introducing the glucosylglycerol (GG) biosynthetic pathway (OD₇₅₀ improved by 24% at 60 h). In addition, the salt tolerance of Syn2973 was further enhanced by overexpressing the rate-limiting step of glycerol-3-phosphate dehydrogenase and downregulating the gene rfbA, which encodes UDP glucose pyrophosphorylase. Taken together, these results indicate that the growth of the end-point strain M-2522-GgpPS-drfbA was improved by 62% compared with the control strain M-pSI-pSII at 60 h under treatment with 0.5 M NaCl. Finally, a comparative metabolomic analysis between strains M-pSI-pSII and M-2522-GgpPS-drfbA was performed to characterize the carbon flux in the engineered M-2522-GgpPS-drfbA strain, and the results showed that more carbon flux was redirected from ADP-GLC to GG synthesis. This study provides important engineering strategies to improve salt tolerance and GG production in Syn2973 in the future.

Dynamic Inventory of Intermediate Metabolites of Cyanobacteria in a Diurnal Cycle

Cyanobacteria are gaining importance both as hosts for photoautotrophic production of chemicals and as model systems for studies of diurnal lifestyle. The proteome and transcriptome of cyanobacteria have been closely examined under diurnal growth, whereas the downstream effects on the intermediary metabolism have not received sufficient attention. The present study focuses on identifying the cellular metabolites whose inventories undergo dramatic changes in a fast-growing cyanobacterium, Synechococcus elongatus PCC 11801. We identified and quantified 67 polar metabolites, whose inventory changes significantly during diurnal growth, with some metabolites changing by 100-fold. The Calvin-Benson-Bassham cycle intermediates peak at midday to support fast growth. The hitherto unexplored γ-glutamyl peptides act as reservoirs of amino acids. Interestingly, several storage molecules or their precursors accumulate during the dark phase, dispelling the notion that all biosynthetic activity takes place in the light phase. Our results will guide metabolic modeling and strain engineering of cyanobacteria.

•

We identify and quantify 67 polar intermediate metabolites in cyanobacteria via LC-MS

•

A number of metabolites show large variations during the diurnal cycle

•

Intermediates of the CBB cycle peak at midday, coinciding with peak in growth rate

•

Gamma-glutamyl dipeptides identified as new storage compounds that peak at dawn

Engineering salt tolerance of photosynthetic cyanobacteria for seawater utilization

Photosynthetic cyanobacteria are capable of utilizing sunlight and CO₂ as sole energy and carbon sources, respectively. With genetically modified cyanobacteria being used as a promising chassis to produce various biofuels and chemicals in recent years, future large-scale cultivation of cyanobacteria would have to be performed in seawater, since freshwater supplies of the earth are very limiting. However, high concentration of salt is known to inhibit the growth of cyanobacteria. This review aims at comparing the mechanisms that different cyanobacteria respond to salt stress, and then summarizing various strategies of developing salt-tolerant cyanobacteria for seawater cultivation, including the utilization of halotolerant cyanobacteria and the engineering of salt-tolerant freshwater cyanobacteria. In addition, the challenges and potential strategies related to further improving salt tolerance in cyanobacteria are also discussed.

Dynamic Changes in Metabolite Accumulation and the Transcriptome during Leaf Growth and Development in Eucommia ulmoides

Eucommia ulmoides Oliver is widely distributed in China. This species has been used mainly in medicine due to the high concentration of chlorogenic acid (CGA), flavonoids, lignans, and other compounds in the leaves and barks. However, the categories of metabolites, dynamic changes in metabolite accumulation and overall molecular mechanisms involved in metabolite biosynthesis during E. ulmoides leaf growth and development remain unknown. Here, a total of 515 analytes, including 127 flavonoids, 46 organic acids, 44 amino acid derivatives, 9 phenolamides, and 16 vitamins, were identified from four E. ulmoides samples using ultraperformance liquid chromatography-mass spectrometry (UPLC-MS) (for widely targeted metabolites). The accumulation of most flavonoids peaked in growing leaves, followed by old leaves. UPLC-MS analysis indicated that CGA accumulation increased steadily to a high concentration during leaf growth and development, and rutin showed a high accumulation level in leaf buds and growing leaves. Based on single-molecule long-read sequencing technology, 69,020 transcripts and 2880 novel loci were identified in E. ulmoides. Expression analysis indicated that isoforms in the flavonoid biosynthetic pathway and flavonoid metabolic pathway were highly expressed in growing leaves and old leaves. Co-expression network analysis suggested a potential direct link between the flavonoid and phenylpropanoid biosynthetic pathways via the regulation of transcription factors, including MYB (v-myb avian myeloblastosis viral oncogene homolog) and bHLH (basic/helix-loop-helix). Our study predicts dynamic metabolic models during leaf growth and development and will support further molecular biological studies of metabolite biosynthesis in E. ulmoides. In addition, our results significantly improve the annotation of the E. ulmoides genome.

Comparison of gene co-networks analysis provide a systems view of rice (Oryza sativa L.) response to Tilletia horrida infection

The biotrophic soil-borne fungus Tilletia horrida causes rice kernel smut, an important disease affecting the production of rice male sterile lines in most hybrid rice growing regions of the world. There are no successful ways of controlling this disease and there has been little study of mechanisms of resistance to T. horrida. Based on transcriptional data of different infection time points, we found 23, 782 and 23, 718 differentially expressed genes (fragments per kilobase of transcript sequence per million, FPKM >1) in Jiangcheng 3A (resistant to T. horrida) and 9311A (susceptible to T. horrida), respectively. In order to illuminate the differential responses of the two rice male sterile lines to T. horrida, we identified gene co-expression modules using the method of weighted gene co-expression network analysis (WGCNA) and compared the different biological functions of gene co-expression networks in key modules at different infection time points. The results indicated that gene co-expression networks in the two rice genotypes were different and that genes contained in some modules of the two groups may play important roles in resistance to T. horrida, such as DTH8 and OsHop/Sti1a. Furthermore, these results provide a global view of the responses of two different phenotypes to T. horrida, and assist our understanding of the regulation of expression changes after T. horrida infection.

The drnf1 Gene from the Drought-Adapted Cyanobacterium Nostoc flagelliforme Improved Salt Tolerance in Transgenic Synechocystis and Arabidopsis Plant

Environmental abiotic stresses are limiting factors for less tolerant organisms, including soil plants. Abiotic stress tolerance-associated genes from prokaryotic organisms are supposed to have a bright prospect for transgenic application. The drought-adapted cyanobacterium Nostoc flagelliforme is arising as a valuable prokaryotic biotic resource for gene excavation. In this study, we evaluated the salt-tolerant function and application potential of a candidate gene drnf1 from N. flagelliforme, which contains a P-loop NTPase (nucleoside-triphosphatase) domain, through heterologous expression in two model organisms Synechocystis sp. PCC 6803 and Arabidopsis thaliana. It was found that DRNF1 could confer significant salt tolerance in both transgenic organisms. In salt-stressed transgenic Synechocystis, DRNF1 could enhance the respiration rate; slow-down the accumulation of exopolysaccharides; up-regulate the expression of salt tolerance-related genes at a higher level, such as those related to glucosylglycerol synthesis, Na⁺/H⁺ antiport, and sugar metabolism; and maintain a better K⁺/Na⁺ homeostasis, as compared to the wild-type strain. These results imply that DRNF1 could facilitate salt tolerance by affecting the respiration metabolism and indirectly regulating the expression of important salt-tolerant genes. Arabidopsis was employed to evaluate the salt tolerance-conferring potential of DRNF1 in plants. The results show that it could enhance the seed germination and shoot growth of transgenic plants under saline conditions. In general, a novel prokaryotic salt-tolerant gene from N. flagelliforme was identified and characterized in this study, enriching the candidate gene pool for genetic engineering in plants.

Comparison of gene co-networks reveals the molecular mechanisms of the rice (Oryza sativa L.) response to Rhizoctonia solani AG1 IA infection

Rhizoctonia solani causes rice sheath blight, an important disease affecting the growth of rice (Oryza sativa L.). Attempts to control the disease have met with little success. Based on transcriptional profiling, we previously identified more than 11,947 common differentially expressed genes (TPM > 10) between the rice genotypes TeQing and Lemont. In the current study, we extended these findings by focusing on an analysis of gene co-expression in response to R. solani AG1 IA and identified gene modules within the networks through weighted gene co-expression network analysis (WGCNA). We compared the different genes assigned to each module and the biological interpretations of gene co-expression networks at early and later modules in the two rice genotypes to reveal differential responses to AG1 IA. Our results show that different changes occurred in the two rice genotypes and that the modules in the two groups contain a number of candidate genes possibly involved in pathogenesis, such as the VQ protein. Furthermore, these gene co-expression networks provide comprehensive transcriptional information regarding gene expression in rice in response to AG1 IA. The co-expression networks derived from our data offer ideas for follow-up experimentation that will help advance our understanding of the translational regulation of rice gene expression changes in response to AG1 IA.

Homeostasis of Second Messenger Cyclic-di-AMP Is Critical for Cyanobacterial Fitness and Acclimation to Abiotic Stress

Second messengers are intracellular molecules regulated by external stimuli known as first messengers that are used for rapid organismal responses to dynamic environmental changes. Cyclic di-AMP (c-di-AMP) is a relatively newly discovered second messenger implicated in cell wall homeostasis in many pathogenic bacteria. C-di-AMP is synthesized from ATP by diadenylyl cyclases (DAC) and degraded by specific c-di-AMP phosphodiesterases (PDE). C-di-AMP DACs and PDEs are present in all sequenced cyanobacteria, suggesting roles for c-di-AMP in the physiology and/or development of these organisms. Despite conservation of these genes across numerous cyanobacteria, the functional roles of c-di-AMP in cyanobacteria have not been well-investigated. In a unique feature of cyanobacteria, phylogenetic analysis indicated that the broadly conserved DAC, related to CdaA/DacA, is always co-associated in an operon with genes critical for controlling cell wall synthesis. To investigate phenotypes regulated by c-di-AMP in cyanobacteria, we overexpressed native DAC (sll0505) and c-di-AMP PDE (slr0104) genes in the cyanobacterium Synechocystis sp. PCC 6803 (hereafter Synechocystis) to increase and decrease intracellular c-di-AMP levels, respectively. DAC- and PDE-overexpression strains, showed abnormal aggregation phenotypes, suggesting functional roles for regulating c-di-AMP homeostasis in vivo. As c-di-AMP may be implicated in osmotic responses in cyanobacteria, we tested whether sorbitol and NaCl stresses impacted expression of sll0505 and slr0104 or intracellular c-di-AMP levels in Synechocystis. Additionally, to determine the range of cyanobacteria in which c-di-AMP may function, we assessed c-di-AMP levels in two unicellular cyanobacteria, i.e., Synechocystis and Synechococcus elongatus PCC 7942, and two filamentous cyanobacteria, i.e., Fremyella diplosiphon and Anabaena sp. PCC 7120. C-di-AMP levels responded differently to abiotic stress signals in distinct cyanobacteria strains, whereas salt stress uniformly impacted another second messenger cyclic di-GMP in cyanobacteria. Together, these results suggest regulation of c-di-AMP homeostasis in cyanobacteria and implicate a role for the second messenger in maintaining cellular fitness in response to abiotic stress.

Exploring stress tolerance mechanism of evolved freshwater strain Chlorella sp. S30 under 30 g/L salt

Enhancement of stress tolerance to high concentration of salt and CO₂ is beneficial for CO₂ capture by microalgae. Adaptive evolution was performed for improving the tolerance of a freshwater strain, Chlorella sp. AE10, to 30 g/L salt. A resulting strain denoted as Chlorella sp. S30 was obtained after 46 cycles (138 days). The stress tolerance mechanism was analyzed by comparative transcriptomic analysis. Although the evolved strain could tolerate 30 g/L salt, high salinity caused loss to photosynthesis, oxidative phosphorylation, fatty acid biosynthesis and tyrosine metabolism. The related genes of antioxidant enzymes, CO₂ fixation, amino acid biosynthesis, central carbon metabolism and ABC transporter proteins were up-regulated. Besides the up-regulation of several genes in Calvin-Benson cycle, they were also identified in C4 photosynthetic pathway and crassulacean acid metabolism pathway. They were essential for the survival and CO₂ fixation of Chlorella sp. S30 under 30 g/L salt and 10% CO₂.

Combining metabolomics and network analysis to improve tacrolimus production in Streptomyces tsukubaensis using different exogenous feedings

Tacrolimus is widely used as an immunosuppressant in the treatment of various autoimmune diseases. However, the low fermentation yield of tacrolimus has thus far restricted its industrial applications. To solve this problem, the time-series response mechanisms of the intracellular metabolism that were highly correlated with tacrolimus biosynthesis were investigated using different exogenous feeding strategies in S. tsukubaensis. The metabolomic datasets, which contained 93 metabolites, were subjected to weighted correlation network analysis (WGCNA), and eight distinct metabolic modules and seven hub metabolites were identified to be specifically associated with tacrolimus biosynthesis. The analysis of metabolites within each metabolic module suggested that the pentose phosphate pathway (PPP), shikimate and aspartate pathway might be the main limiting factors in the rapid synthesis phase of tacrolimus accumulation. Subsequently, all possible key-limiting steps in the above metabolic pathways were further screened using a genome-scale metabolic network model (GSMM) of S. tsukubaensis. Based on the prediction results, two newly identified targets (aroC and dapA) were overexpressed experimentally, and both of the engineered strains showed higher tacrolimus production. Moreover, the best strain, HT-aroC/dapA, that was engineered to simultaneously enhanced chorismate and lysine biosynthesis was able to produce 128.19 mg/L tacrolimus, 1.64-fold higher than control (78.26 mg/L). These findings represent a valuable addition to our understanding of tacrolimus accumulation in S. tsukubaensis, and pave the way to further production improvements.

Using scale and feather traits for module construction provides a functional approach to chicken epidermal development

Gene co-expression network analysis has been a research method widely used in systematically exploring gene function and interaction. Using the Weighted Gene Co-expression Network Analysis (WGCNA) approach to construct a gene co-expression network using data from a customized 44K microarray transcriptome of chicken epidermal embryogenesis, we have identified two distinct modules that are highly correlated with scale or feather development traits. Signaling pathways related to feather development were enriched in the traditional KEGG pathway analysis and functional terms relating specifically to embryonic epidermal development were also enriched in the Gene Ontology analysis. Significant enrichment annotations were discovered from customized enrichment tools such as Modular Single-Set Enrichment Test (MSET) and Medical Subject Headings (MeSH). Hub genes in both trait-correlated modules showed strong specific functional enrichment toward epidermal development. Also, regulatory elements, such as transcription factors and miRNAs, were targeted in the significant enrichment result. This work highlights the advantage of this methodology for functional prediction of genes not previously associated with scale- and feather trait-related modules.

A partial proteome reference map of Tetragenococcus halophilus and comparative proteomic and physiological analysis under salt stress

Tetragenococcus halophilus, a moderately halophilic Gram-positive lactic acid bacteria, was widely existed in many food fermentation systems, where salt stress is an environmental condition commonly encountered.

MALDI-TOF Mass Spectrometry Enables a Comprehensive and Fast Analysis of Dynamics and Qualities of Stress Responses of Lactobacillus paracasei subsp. paracasei F19

Lactic acid bacteria (LAB) are widely used as starter cultures in the manufacture of foods. Upon preparation, these cultures undergo various stresses resulting in losses of survival and fitness. In order to find conditions for the subsequent identification of proteomic biomarkers and their exploitation for preconditioning of strains, we subjected Lactobacillus (Lb.) paracasei subsp. paracasei TMW 1.1434 (F19) to different stress qualities (osmotic stress, oxidative stress, temperature stress, pH stress and starvation stress). We analysed the dynamics of its stress responses based on the expression of stress proteins using MALDI-TOF mass spectrometry (MS), which has so far been used for species identification. Exploiting the methodology of accumulating protein expression profiles by MALDI-TOF MS followed by the statistical evaluation with cluster analysis and discriminant analysis of principle components (DAPC), it was possible to monitor the expression of low molecular weight stress proteins, identify a specific time point when the expression of stress proteins reached its maximum, and statistically differentiate types of adaptive responses into groups. Above the specific result for F19 and its stress response, these results demonstrate the discriminatory power of MALDI-TOF MS to characterize even dynamics of stress responses of bacteria and enable a knowledge-based focus on the laborious identification of biomarkers and stress proteins. To our knowledge, the implementation of MALDI-TOF MS protein profiling for the fast and comprehensive analysis of various stress responses is new to the field of bacterial stress responses. Consequently, we generally propose MALDI-TOF MS as an easy and quick method to characterize responses of microbes to different environmental conditions, to focus efforts of more elaborate approaches on time points and dynamics of stress responses.

Seawater cultivation of freshwater cyanobacterium Synechocystis sp. PCC 6803 drastically alters amino acid composition and glycogen metabolism

Water use assessment is important for bioproduction using cyanobacteria. For eco-friendly reasons, seawater should preferably be used for cyanobacteria cultivation instead of freshwater. In this study, we demonstrated that the freshwater unicellular cyanobacterium Synechocystis sp. PCC 6803 could be grown in a medium based on seawater. The Synechocystis wild-type strain grew well in an artificial seawater (ASW) medium supplemented with nitrogen and phosphorus sources. The addition of HEPES buffer improved cell growth overall, although the growth in ASW medium was inferior to that in the synthetic BG-11 medium. The levels of proteins involved in sugar metabolism changed depending on the culture conditions. The biosynthesis of several amino acids including aspartate, glutamine, glycine, proline, ornithine, and lysine, was highly up-regulated by cultivation in ASW. Two types of natural seawater (NSW) were also made available for the cultivation of Synechocystis cells, with supplementation of both nitrogen and phosphorus sources. These results revealed the potential use of seawater for the cultivation of freshwater cyanobacteria, which would help to reduce freshwater consumption during biorefinery using cyanobacteria.

Combinatory strategy for characterizing and understanding the ethanol synthesis pathway in cyanobacteria cell factories

BACKGROUND

Photosynthetic production of chemicals and fuels by recycling CO2 in cyanobacteria is a promising solution facing energy shortage and resource declination. Ethanol is an attractive and demonstrative biofuel product, and ethanol synthesis in cyanobacteria has been achieved by assembling of a pathway consisting of pyruvate decarboxylase (PDCzm) and alcohol dehydrogenase II (slr1192). For enabling more powerful ethanol photosynthetic production, an optimized and balanced catalyzing route was required. In this work, we provided a paradigm for systematically characterizing and optimizing the PDCzm-slr1192 pathway from engineered cyanobacteria strains, combining in vitro reconstitution, genetic engineering and feeding-cultivation.

RESULTS

We reconstituted the PDCzm-slr1192 pathway in vitro and performed specific titration assays for enzymes, substrates, cofactors, and metal ions. In the in vitro system, K 50 of PDCzm was 0.326 μM, with a V max of 2.074 μM/s; while for slr1192, the values were 0.109 μM and 1.722 μM/s, respectively. Titration response discrepancy indicated that PDCzm rather than slr1192 was the rate-limiting factor for ethanol synthesis. In addition, a 4:6 concentration ratio of PDCzm-slr1192 would endow the reaction with a maximal specific catalytic activity. Titration assays for other components were also performed. K m values for NADPH, pyruvate, TPP, Mg(2+) and acetaldehyde were 0.136, 6.496, 0.011, 0.104, and 0.393 mM, respectively. We further constructed Synechocystis mutant strains with diverse PDCzm-slr1192 concentrations and ratios, and compared the growth and ethanol synthesis performances. The results revealed that activities of PDCzm indeed held control over the ethanol generation capacities. We performed pyruvate-feeding treatment with the newly developed Syn-YQ4 strain, and confirmed that improvement of pyruvate supply would direct more carbon flow to ethanol formation.

CONCLUSIONS

We systematically characterized and optimized the PDCzm-slr1192 pathway in engineered cyanobacteria for ethanol production. Information gained from in vitro monitoring and genetic engineering revealed that for further enhancing ethanol synthesis capacities, PDCzm activities needed enhancement, and the PDCzm-slr1192 ratio should be improved and held to about 1:1.5. Considering actual metabolites concentrations of cyanobacteria cells, enhancing pyruvate supply was also a promising strategy for further updating the current ethanol photosynthetic cell factories.

Metabolomic analysis reveals mechanism of antioxidant butylated hydroxyanisole on lipid accumulation in Crypthecodinium cohnii

The heterotrophic dinoflagellate alga Crypthecodinium cohnii is known to accumulate lipids with a high fraction of docosahexaenoic acid (DHA). In this study, we first evaluated two antioxidant compounds, butylated hydroxyanisole (BHA) and propyl gallate (PG), for their effects on lipid accumulation in C. cohnii. The results showed that antioxidant BHA could increase lipid accumulation in C. cohnii by 8.80% at a final concentration of 30 μM, while PG had no obvious effect on lipid accumulation at the tested concentrations. To decipher the molecular mechanism responsible for the increased lipid accumulation by BHA, we employed an integrated GC-MS and LC-MS metabolomic approach to determine the time-series metabolic profiles with or without BHA, and then subjected the metabolomic data to a principal component analysis (PCA) and a weighted gene coexpression network analysis (WGCNA) network analyses to identify the key metabolic modules and metabolites possibly relevant to the increased lipid accumulation. LC-MS analysis showed that several metabolites, including NADPH, could be important for the stimulation role of BHA on lipid accumulation. Meanwhile GC-MS and network analyses allowed identification of eight metabolic modules and nine hub metabolites possibly relevant to the stimulation role of BHA in C. cohnii. The study provided a metabolomics view of the BHA mode of action on lipid accumulation in C. cohnii, and the information could be valuable for a better understanding of antioxidant effects on lipid accumulation in other microalgae as well.

Protein Network Signatures Associated with Exogenous Biofuels Treatments in Cyanobacterium Synechocystis sp. PCC 6803

Although recognized as a promising microbial cell factory for producing biofuels, current productivity in cyanobacterial systems is low. To make the processes economically feasible, one of the hurdles, which need to be overcome is the low tolerance of hosts to toxic biofuels. Meanwhile, little information is available regarding the cellular responses to biofuels stress in cyanobacteria, which makes it challenging for tolerance engineering. Using large proteomic datasets of Synechocystis under various biofuels stress and environmental perturbation, a protein co-expression network was first constructed and then combined with the experimentally determined protein–protein interaction network. Proteins with statistically higher topological overlap in the integrated network were identified as common responsive proteins to both biofuels stress and environmental perturbations. In addition, a weighted gene co-expression network analysis was performed to distinguish unique responses to biofuels from those to environmental perturbations and to uncover metabolic modules and proteins uniquely associated with biofuels stress. The results showed that biofuel-specific proteins and modules were enriched in several functional categories, including photosynthesis, carbon fixation, and amino acid metabolism, which may represent potential key signatures for biofuels stress responses in Synechocystis. Network-based analysis allowed determination of the responses specifically related to biofuels stress, and the results constituted an important knowledge foundation for tolerance engineering against biofuels in Synechocystis.

Metabolomic basis of laboratory evolution of butanol tolerance in photosynthetic Synechocystis sp. PCC 6803

BACKGROUND

Recent efforts demonstrated the potential application of cyanobacteria as a "microbial cell factory" to produce butanol directly from CO2. However, cyanobacteria have very low tolerance to the toxic butanol, which limits the economic viability of this renewable system.

RESULTS

Through a long-term experimental evolution process, we achieved a 150% increase of the butanol tolerance in a model cyanobacterium Synechocystis sp. PCC 6803 after a continuous 94 passages for 395 days in BG11 media amended with gradually increased butanol concentration from 0.2% to 0.5% (v/v). To decipher the molecular mechanism responsible for the tolerance increase, we employed an integrated GC-MS and LC-MS approach to determine metabolomic profiles of the butanol-tolerant Synechocystis strains isolated from several stages of the evolution, and then applied PCA and WGCNA network analyses to identify the key metabolites and metabolic modules related to the increased tolerance. The results showed that unstable metabolites of 3-phosphoglyceric acid (3PG), D-fructose 6-phosphate (F6P), D-glucose 6-phosphate (G6P), NADPH, phosphoenolpyruvic acid (PEP), D-ribose 5-phosphate (R5P), and stable metabolites of glycerol, L-serine and stearic acid were differentially regulated during the evolution process, which could be related to tolerance increase to butanol in Synechocystis.

CONCLUSIONS

The study provided the first time-series description of the metabolomic changes related to the gradual increase of butanol tolerance, and revealed a metabolomic basis important for rational tolerance engineering in Synechocystis.

Metabolomic and network analysis of astaxanthin-producing Haematococcus pluvialis under various stress conditions

Various combinations of acetate (Ac), Fe(2+) and high light (HL) stress conditions were evaluated to maximize astaxanthin accumulation and biomass production in Haematococcus pluvialis, and then GC-MS and LC-MS based metabolomics were applied to determine molecular mechanisms responsible for enhancing astaxanthin accumulation under the stress conditions. With the optimized analytical protocols, the GC-MS and LC-MS analyses allowed identification of 93 stable and 24 unstable intracellular metabolites from H. pluvialis, respectively. In addition, a metabolic network was constructed based on GC-MS metabolomic datasets using a weighted correlation network analysis (WGCNA) approach. The network analysis uncovered 2, 1 and 1 distinguished metabolic modules highly associated with HL, Fe(2+) & HL, and Ac & Fe(2+) & HL conditions, respectively. Finally, LC-MS analysis found that AKG, Glu and R5P may be metabolites associated with the Fe(2+) & HL condition. The study provided the first metabolomic view of cell growth and astaxanthin accumulation in H. pluvialis.