Transcription profiling of butanol producer Clostridium beijerinckii NRRL B-598 using RNA-Seq

Augusta: From RNA-Seq to gene regulatory networks and Boolean models

Computational models of gene regulations help to understand regulatory mechanisms and are extensively used in a wide range of areas, e.g., biotechnology or medicine, with significant benefits. Unfortunately, there are only a few computational gene regulatory models of whole genomes allowing static and dynamic analysis due to the lack of sophisticated tools for their reconstruction. Here, we describe Augusta, an open-source Python package for Gene Regulatory Network (GRN) and Boolean Network (BN) inference from the high-throughput gene expression data. Augusta can reconstruct genome-wide models suitable for static and dynamic analyses. Augusta uses a unique approach where the first estimation of a GRN inferred from expression data is further refined by predicting transcription factor binding motifs in promoters of regulated genes and by incorporating verified interactions obtained from databases. Moreover, a refined GRN is transformed into a draft BN by searching in the curated model database and setting logical rules to incoming edges of target genes, which can be further manually edited as the model is provided in the SBML file format. The approach is applicable even if information about the organism under study is not available in the databases, which is typically the case for non-model organisms including most microbes. Augusta can be operated from the command line and, thus, is easy to use for automated prediction of models for various genomes. The Augusta package is freely available at github.com/JanaMus/Augusta. Documentation and tutorials are available at augusta.readthedocs.io.

Transcriptomic studies of solventogenic clostridia, Clostridium acetobutylicum and Clostridium beijerinckii

Solventogenic clostridia are not a strictly defined group within the genus Clostridium but its representatives share some common features, i.e. they are anaerobic, non-pathogenic, non-toxinogenic and endospore forming bacteria. Their main metabolite is typically 1-butanol but depending on species and culture conditions, they can form other metabolites such as acetone, isopropanol, ethanol, butyric, lactic and acetic acids, and hydrogen. Although these organisms were previously used for the industrial production of solvents, they later fell into disuse, being replaced by more efficient chemical production. A return to a more biological production of solvents therefore requires a thorough understanding of clostridial metabolism. Transcriptome analysis, which reflects the involvement of individual genes in all cellular processes within a population, at any given (sampling) moment, is a valuable tool for gaining a deeper insight into clostridial life. In this review, we describe techniques to study transcription, summarize the evolution of these techniques and compare methods for data processing and visualization of solventogenic clostridia, particularly the species Clostridium acetobutylicum and Clostridium beijerinckii. Individual approaches for evaluating transcriptomic data are compared and their contributions to advancements in the field are assessed. Moreover, utilization of transcriptomic data for reconstruction of computational clostridial metabolic models is considered and particular models are described. Transcriptional changes in glucose transport, central carbon metabolism, the sporulation cycle, butanol and butyrate stress responses, the influence of lignocellulose-derived inhibitors on growth and solvent production, and other respective topics, are addressed and common trends are highlighted.

Genome-Wide Transcriptomic Analysis of n-Caproic Acid Production in Ruminococcaceae Bacterium CPB6 with Lactate Supplementation

n-Caproic acid (CA) is gaining increased attention due to its high value as a chemical feedstock. Ruminococcaceae bacterium strain CPB6 is an anaerobic mesophilic bacterium that is highly prolific in its ability to perform chain elongation of lactate to CA. However, little is known about the genome-wide transcriptional analysis of strain CPB6 for CA production triggered by the supplementation of exogenous lactate. In this study, cultivation of strain CPB6 was carried out in the absence and presence of lactate. Transcriptional profiles were analyzed using RNA-seq, and differentially expressed genes (DEGs) between the lactate-supplemented cells and control cells without lactate were analyzed. The results showed that lactate supplementation led to earlier CA p,roduction, and higher final CA titer and productivity. 295 genes were substrate and/or growth dependent, and these genes cover crucial functional categories. Specifically, 5 genes responsible for the reverse β-oxidation pathway, 11 genes encoding ATP-binding cassette (ABC) transporters, 6 genes encoding substrate-binding protein (SBP), and 4 genes encoding phosphotransferase system (PTS) transporters were strikingly upregulated in response to the addition of lactate. These genes would be candidates for future studies aiming at understanding the regulatory mechanism of lactate conversion into CA, as well as for the improvement of CA production in strain CPB6. The findings presented herein reveal unique insights into the biomolecular effect of lactate on CA production at the transcriptional level.

Improving sustainable hydrogen production from green waste: [FeFe]-hydrogenases quantitative gene expression RT-qPCR analysis in presence of autochthonous consortia

BACKGROUND

Bio-hydrogen production via dark fermentation of low-value waste is a potent and simple mean of recovering energy, maximising the harvesting of reducing equivalents to produce the cleanest fuel amongst renewables. Following several position papers from companies and public bodies, the hydrogen economy is regaining interest, especially in combination with circular economy and the environmental benefits of short local supply chains, aiming at zero net emission of greenhouse gases (GHG). The biomasses attracting the largest interest are agricultural and urban green wastes (pruning of trees, collected leaves, grass clippings from public parks and boulevards), which are usually employed in compost production, with some concerns over the GHG emission during the process. Here, an alternative application of green wastes, low-value compost and intermediate products (partially composted but unsuitable for completing the process) is studied, pointing at the autochthonous microbial consortium as an already selected source of implementation for biomass degradation and hydrogen production. The biocatalysts investigated as mainly relevant for hydrogen production were the [FeFe]-hydrogenases expressed in Clostridia, given their very high turnover rates.

RESULTS

Bio-hydrogen accumulation was related to the modulation of gene expression of multiple [FeFe]-hydrogenases from two strains (Clostridium beijerinckii AM2 and Clostridium tyrobutyricum AM6) isolated from the same waste. Reverse Transcriptase quantitative PCR (RT-qPCR) was applied over a period of 288 h and the RT-qPCR results showed that C. beijerinckii AM2 prevailed over C. tyrobutyricum AM6 and a high expression modulation of the 6 different [FeFe]-hydrogenase genes of C. beijerinckii in the first 23 h was observed, sustaining cumulative hydrogen production of 0.6 to 1.2 ml H₂/g VS (volatile solids). These results are promising in terms of hydrogen yields, given that no pre-treatment was applied, and suggested a complex cellular regulation, linking the performance of dark fermentation with key functional genes involved in bio-H₂ production in presence of the autochthonous consortium, with different roles, time, and mode of expression of the involved hydrogenases.

CONCLUSIONS

An applicative outcome of the hydrogenases genes quantitative expression analysis can be foreseen in optimising (on the basis of the acquired functional data) hydrogen production from a nutrient-poor green waste and/or low added value compost, in a perspective of circular bioeconomy.

Sporulation in solventogenic and acetogenic clostridia

The Clostridium genus harbors compelling organisms for biotechnological production processes; while acetogenic clostridia can fix C1-compounds to produce acetate and ethanol, solventogenic clostridia can utilize a wide range of carbon sources to produce commercially valuable carboxylic acids, alcohols, and ketones by fermentation. Despite their potential, the conversion by these bacteria of carbohydrates or C1 compounds to alcohols is not cost-effective enough to result in economically viable processes. Engineering solventogenic clostridia by impairing sporulation is one of the investigated approaches to improve solvent productivity. Sporulation is a cell differentiation process triggered in bacteria in response to exposure to environmental stressors. The generated spores are metabolically inactive but resistant to harsh conditions (UV, chemicals, heat, oxygen). In Firmicutes, sporulation has been mainly studied in bacilli and pathogenic clostridia, and our knowledge of sporulation in solvent-producing or acetogenic clostridia is limited. Still, sporulation is an integral part of the cellular physiology of clostridia; thus, understanding the regulation of sporulation and its connection to solvent production may give clues to improve the performance of solventogenic clostridia. This review aims to provide an overview of the triggers, characteristics, and regulatory mechanism of sporulation in solventogenic clostridia. Those are further compared to the current knowledge on sporulation in the industrially relevant acetogenic clostridia. Finally, the potential applications of spores for process improvement are discussed.Key Points• The regulatory network governing sporulation initiation varies in solventogenic clostridia.• Media composition and cell density are the main triggers of sporulation.• Spores can be used to improve the fermentation process.

Identification and Validation of Reference Genes in Clostridium beijerinckii NRRL B-598 for RT-qPCR Using RNA-Seq Data

Gene expression analysis through reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR) depends on correct data normalization by reference genes with stable expression. Although Clostridium beijerinckii NRRL B-598 is a promising Gram-positive bacterium for the industrial production of biobutanol, validated reference genes have not yet been reported. In this study, we selected 160 genes with stable expression based on an RNA sequencing (RNA-Seq) data analysis, and among them, seven genes (zmp, rpoB1, rsmB, greA, rpoB2, topB2, and rimO) were selected for experimental validation by RT-qPCR and gene ontology (GO) enrichment analysis. According to statistical analyses, zmp and greA were the most stable and suitable reference genes for RT-qPCR normalization. Furthermore, our methodology can be useful for selection of the reference genes in other strains of C. beijerinckii and it also suggests that the RNA-Seq data can be used for the initial selection of novel reference genes, however, their validation is required.

Changes in efflux pump activity of Clostridium beijerinckii throughout ABE fermentation

Pumping toxic substances through a cytoplasmic membrane by protein transporters known as efflux pumps represents one bacterial mechanism involved in the stress response to the presence of toxic compounds. The active efflux might also take part in exporting low-molecular-weight alcohols produced by intrinsic cell metabolism; in the case of solventogenic clostridia, predominantly acetone, butanol and ethanol (ABE). However, little is known about this active efflux, even though some evidence exists that membrane pumps might be involved in solvent tolerance. In this study, we investigated changes in overall active efflux during ABE fermentation, employing a flow cytometric protocol adjusted for Clostridia and using ethidium bromide (EB) as a fluorescence marker for quantification of direct efflux. A fluctuation in efflux during the course of standard ABE fermentation was observed, with a maximum reached during late acidogenesis, a high efflux rate during early and mid-solventogenesis and an apparent decrease in EB efflux rate in late solventogenesis. The fluctuation in efflux activity was in accordance with transcriptomic data obtained for various membrane exporters in a former study. Surprisingly, under altered cultivation conditions, when solvent production was attenuated, and extended acidogenesis was promoted, stable low efflux activity was reached after an initial peak that appeared in the stage comparable to standard ABE fermentation. This study confirmed that efflux pump activity is not constant during ABE fermentation and suggests that undisturbed solvent production might be a trigger for activation of pumps involved in solvent efflux. KEY POINTS: • Flow cytometric assay for efflux quantification in Clostridia was established. • Efflux rate peaked in late acidogenesis and in early solventogenesis. • Impaired solventogenesis led to an overall decrease in efflux.

Deeper below the surface-transcriptional changes in selected genes of Clostridium beijerinckii in response to butanol shock

The main bottleneck in the return of industrial butanol production from renewable feedstock through acetone–butanol–ethanol (ABE) fermentation by clostridia, such as Clostridium beijerinckii, is the low final butanol concentration. The problem is caused by the high toxicity of butanol to the production cells, and therefore, understanding the mechanisms by which clostridia react to butanol shock is of key importance. Detailed analyses of transcriptome data that were obtained after butanol shock and their comparison with data from standard ABE fermentation have resulted in new findings, while confirmed expected population responses. Although butanol shock resulted in upregulation of heat shock protein genes, their regulation is different than was assumed based on standard ABE fermentation transcriptome data. While glucose uptake, glycolysis, and acidogenesis genes were downregulated after butanol shock, solventogenesis genes were upregulated. Cyclopropanation of fatty acids and formation of plasmalogens seem to be significant processes involved in cell membrane stabilization in the presence of butanol. Surprisingly, one of the three identified Agr quorum‐sensing system genes was upregulated. Upregulation of several putative butanol efflux pumps was described after butanol addition and a large putative polyketide gene cluster was found, the transcription of which seemed to depend on the concentration of butanol.

Phenotypic and Genomic Analysis of Clostridium beijerinckii NRRL B-598 Mutants With Increased Butanol Tolerance

N-Butanol, a valuable solvent and potential fuel extender, can be produced via acetone-butanol-ethanol (ABE) fermentation. One of the main drawbacks of ABE fermentation is the high toxicity of butanol to producing cells, leading to cell membrane disruption, low culture viability and, consequently, low produced concentrations of butanol. The goal of this study was to obtain mutant strains of Clostridium beijerinckii NRRL B-598 with improved butanol tolerance using random chemical mutagenesis, describe changes in their phenotypes compared to the wild-type strain and reveal changes in the genome that explain improved tolerance or other phenotypic changes. Nine mutant strains with stable improved features were obtained by three different approaches and, for two of them, ethidium bromide (EB), a known substrate of efflux pumps, was used for either selection or as a mutagenic agent. It is the first utilization of this approach for the development of butanol-tolerant mutants of solventogenic clostridia, for which generally there is a lack of knowledge about butanol efflux or efflux mechanisms and their regulation. Mutant strains exhibited increase in butanol tolerance from 36% up to 127% and the greatest improvement was achieved for the strains for which EB was used as a mutagenic agent. Additionally, increased tolerance to other substrates of efflux pumps, EB and ethanol, was observed in all mutants and higher antibiotic tolerance in some of the strains. The complete genomes of mutant strains were sequenced and revealed that improved butanol tolerance can be attributed to mutations in genes encoding typical stress responses (chemotaxis, autolysis or changes in cell membrane structure), but, also, to mutations in genes X276_07980 and X276_24400, encoding efflux pump regulators. The latter observation confirms the importance of efflux in butanol stress response of the strain and offers new targets for rational strain engineering.

How to outwit nature: Omics insight into butanol tolerance

The energy crisis, depletion of oil reserves, and global climate changes are pressing problems of developed societies. One possibility to counteract that is microbial production of butanol, a promising new fuel and alternative to many petrochemical reagents. However, the high butanol toxicity to all known microbial species is the main obstacle to its industrial implementation. The present state of the art review aims to expound the recent advances in modern omics approaches to resolving this insurmountable to date problem of low butanol tolerance. Genomics, transcriptomics, and proteomics show that butanol tolerance is a complex phenomenon affecting multiple genes and their expression. Efflux pumps, stress and multidrug response, membrane transport, and redox-related genes are indicated as being most important during butanol challenge, in addition to fine-tuning of global regulators of transcription (Spo0A, GntR), which may further improve tolerance. Lipidomics shows that the alterations in membrane composition (saturated lipids and plasmalogen increase) are very much species-specific and butanol-related. Glycomics discloses the pleiotropic effect of CcpA, the role of alternative sugar transport, and the production of exopolysaccharides as alternative routes to overcoming butanol stress. Unfortunately, the strain that simultaneously syntheses and tolerates butanol in concentrations that allow its commercialization has not yet been discovered or produced. Omics insight will allow the purposeful increase of butanol tolerance in natural and engineered producers and the effective heterologous expression of synthetic butanol pathways in strains hereditary butanol-resistant up to 3.2 - 4.9% (w/v). Future breakthrough can be achieved by a detailed study of the membrane proteome, of which 21% are proteins with unknown functions.

Synthetic Biology Tools for Genome and Transcriptome Engineering of Solventogenic Clostridium

Strains of Clostridium genus are used for production of various value-added products including fuels and chemicals. Development of any commercially viable production process requires a combination of both strain and fermentation process development strategies. The strain development in Clostridium sp. could be achieved by random mutagenesis, and targeted gene alteration methods. However, strain improvement in Clostridium sp. by targeted gene alteration method was challenging due to the lack of efficient tools for genome and transcriptome engineering in this organism. Recently, various synthetic biology tools have been developed to facilitate the strain engineering of solventogenic Clostridium. In this review, we consolidated the recent advancements in toolbox development for genome and transcriptome engineering in solventogenic Clostridium. Here we reviewed the genome-engineering tools employing mobile group II intron, pyrE alleles exchange, and CRISPR/Cas9 with their application for strain development of Clostridium sp. Next, transcriptome engineering tools such as untranslated region (UTR) engineering and synthetic sRNA techniques were also discussed in context of Clostridium strain engineering. Application of any of these discussed techniques will facilitate the metabolic engineering of clostridia for development of improved strains with respect to requisite functional attributes. This might lead to the development of an economically viable butanol production process with improved titer, yield and productivity.

Transcriptional analysis of amino acid, metal ion, vitamin and carbohydrate uptake in butanol-producing Clostridium beijerinckii NRRL B-598

In-depth knowledge of cell metabolism and nutrient uptake mechanisms can lead to the development of a tool for improving acetone-butanol-ethanol (ABE) fermentation performance and help to overcome bottlenecks in the process, such as the high cost of substrates and low production rates. Over 300 genes potentially encoding transport of amino acids, metal ions, vitamins and carbohydrates were identified in the genome of the butanol-producing strain Clostridium beijerinckii NRRL B-598, based on similarity searches in protein function databases. Transcriptomic data of the genes were obtained during ABE fermentation by RNA-Seq experiments and covered acidogenesis, solventogenesis and sporulation. The physiological roles of the selected 81 actively expressed transport genes were established on the basis of their expression profiles at particular stages of ABE fermentation. This article describes how genes encoding the uptake of glucose, iron, riboflavin, glutamine, methionine and other nutrients take part in growth, production and stress responses of C. beijerinckii NRRL B-598. These data increase our knowledge of transport mechanisms in solventogenic Clostridium and may be used in the selection of individual genes for further research.

A transcriptional response of Clostridium beijerinckii NRRL B-598 to a butanol shock

BACKGROUND

One of the main obstacles preventing solventogenic clostridia from achieving higher yields in biofuel production is the toxicity of produced solvents. Unfortunately, regulatory mechanisms responsible for the shock response are poorly described on the transcriptomic level. Although the strain Clostridium beijerinckii NRRL B-598, a promising butanol producer, has been studied under different conditions in the past, its transcriptional response to a shock caused by butanol in the cultivation medium remains unknown.

RESULTS

In this paper, we present a transcriptional response of the strain during a butanol challenge, caused by the addition of butanol to the cultivation medium at the very end of the acidogenic phase, using RNA-Seq. We resequenced and reassembled the genome sequence of the strain and prepared novel genome and gene ontology annotation to provide the most accurate results. When compared to samples under standard cultivation conditions, samples gathered during butanol shock represented a well-distinguished group. Using reference samples gathered directly before the addition of butanol, we identified genes that were differentially expressed in butanol challenge samples. We determined clusters of 293 down-regulated and 301 up-regulated genes whose expression was affected by the cultivation conditions. Enriched term "RNA binding" among down-regulated genes corresponded to the downturn of translation and the cluster contained a group of small acid-soluble spore proteins. This explained phenotype of the culture that had not sporulated. On the other hand, up-regulated genes were characterized by the term "protein binding" which corresponded to activation of heat-shock proteins that were identified within this cluster.

CONCLUSIONS

We provided an overall transcriptional response of the strain C. beijerinckii NRRL B-598 to butanol shock, supplemented by auxiliary technologies, including high-pressure liquid chromatography and flow cytometry, to capture the corresponding phenotypic response. We identified genes whose regulation was affected by the addition of butanol to the cultivation medium and inferred related molecular functions that were significantly influenced. Additionally, using high-quality genome assembly and custom-made gene ontology annotation, we demonstrated that this settled terminology, widely used for the analysis of model organisms, could also be applied to non-model organisms and for research in the field of biofuels.

Acidogenesis, solventogenesis, metabolic stress response and life cycle changes in Clostridium beijerinckii NRRL B-598 at the transcriptomic level

Clostridium beijerinckii NRRL B-598 is a sporulating, butanol and hydrogen producing strain that utilizes carbohydrates by the acetone-butanol-ethanol (ABE) fermentative pathway. The pathway consists of two metabolic phases, acidogenesis and solventogenesis, from which the latter one can be coupled with sporulation. Thorough transcriptomic profiling during a complete life cycle and both metabolic phases completed with flow cytometry, microscopy and a metabolites analysis helped to find out key genes involved in particular cellular events. The description of genes/operons that are closely involved in metabolism or the cell cycle is a necessary condition for metabolic engineering of the strain and will be valuable for all C. beijerinckii strains and other Clostridial species. The study focused on glucose transport and catabolism, hydrogen formation, metabolic stress response, binary fission, motility/chemotaxis and sporulation, which resulted in the composition of the unique image reflecting clostridial population changes. Surprisingly, the main change in expression of individual genes was coupled with the sporulation start and not with the transition from acidogenic to solventogenic metabolism. As expected, solvents formation started at pH decrease and the accumulation of butyric and acetic acids in the cultivation medium.