Transcriptome analysis of thermophilic methylotrophic Bacillus methanolicus MGA3 using RNA-sequencing provides detailed insights into its previously uncharted transcriptional landscape

Genome-Wide Transcription Start Sites Mapping in Methylorubrum Grown with Dichloromethane and Methanol

Dichloromethane (DCM, methylene chloride) is a toxic halogenated volatile organic compound massively used for industrial applications, and consequently often detected in the environment as a major pollutant. DCM biotransformation suggests a sustainable decontamination strategy of polluted sites. Among methylotrophic bacteria able to use DCM as a sole source of carbon and energy for growth, Methylorubrum extorquens DM4 is a longstanding reference strain. Here, the primary 5′-ends of transcripts were obtained using a differential RNA-seq (dRNA-seq) approach to provide the first transcription start site (TSS) genome-wide landscape of a methylotroph using DCM or methanol. In total, 7231 putative TSSs were annotated and classified with respect to their localization to coding sequences (CDSs). TSSs on the opposite strand of CDS (antisense TSS) account for 31% of all identified TSSs. One-third of the detected TSSs were located at a distance to the start codon inferior to 250 nt (average of 84 nt) with 7% of leaderless mRNA. Taken together, the global TSS map for bacterial growth using DCM or methanol will facilitate future studies in which transcriptional regulation is crucial, and efficient DCM removal at polluted sites is limited by regulatory processes.

Aerobic Utilization of Methanol for Microbial Growth and Production

Methanol is a reduced one-carbon (C1) compound. It supports growth of aerobic methylotrophs that gain ATP from reduced redox equivalents by respiratory phosphorylation in their electron transport chains. Notably, linear oxidation of methanol to carbon dioxide may yield three reduced redox equivalents if methanol oxidation is NAD-dependent as, e.g., in Bacillus methanolicus. Methanol has a higher degree of reduction per carbon than glucose (6 vs. 4), and thus, lends itself as an ideal carbon source for microbial production of reduced target compounds. However, C-C bond formation in the RuMP or serine cycle, a prerequisite for production of larger molecules, requires ATP and/or reduced redox equivalents. Moreover, heat dissipation and a high demand for oxygen during catabolic oxidation of methanol may pose challenges for fermentation processes. In this chapter, we summarize metabolic pathways for aerobic methanol utilization, aerobic methylotrophs as industrial production hosts, strain engineering, and methanol bioreactor processes. In addition, we provide technological and market outlooks.

Characterization of two 3-deoxy-d-Arabino-Heptulosonate 7-phosphate synthases from Bacillusmethanolicus

3-Deoxy-d-arabino-heptulosonate 7-phosphate (DAHP) synthase catalyzes the condensation of phosphoenolpyruvate (PEP) with d-erythrose 4-phosphate (E4P) and plays an important role in regulating carbon flux toward aromatic amino acid biosynthesis in bacteria and plants. Sequence analysis of the DAHP synthases AroG1 and AroG2 from Bacillus methanolicus MGA3 suggested this thermophilic, methylotrophic bacterium possesses two type Iβ DAHP synthases. This study describes production of AroG1 and AroG2 in Escherichia coli as hexa-histidine fused proteins, which were purified by affinity chromatography. Treatment with TEV protease afforded native proteins for characterization and kinetic analysis. AroG1 and AroG2 are, respectively, 30.1 kDa and 40.0 kDa proteins. Both enzymes have maximal activity over a pH range of 6.3-7.2. The apparent kinetic parameters at 50 °C and pH 7.2 for AroG1 are K_m^PEP 1100 ± 100 μM, K_m^E4P 530 ± 100 μM, and k_cat 10.3 ± 1.2 s^-1. The kinetic parameters for AroG2 are K_m^PEP 90 ± 20 μM, K_m^E4P 130 ± 40 μM, and k_cat 2.0 ± 0.2 s^-1. At 50 °C AroG2 retains 50% of its activity after 96 min whereas AroG1 retains less than 5% of its activity after 10 min. AroG2, which contains an N-terminal regulatory domain, is inhibited by chorismate and prephenate but not l-phenylalanine, l-tyrosine, or l-tryptophan. AroG1 is not inhibited by any of the molecules examined. Understanding DAHP synthase regulation in B. methanolicus is a first step toward generating biocatalysts that exploit the target-rich aromatic amino acid biosynthetic pathway for synthesis of chemicals from methanol.

Interrogating the Role of the Two Distinct Fructose-Bisphosphate Aldolases of Bacillus methanolicus by Site-Directed Mutagenesis of Key Amino Acids and Gene Repression by CRISPR Interference

The Gram-positive Bacillus methanolicus shows plasmid-dependent methylotrophy. This facultative ribulose monophosphate (RuMP) cycle methylotroph possesses two fructose bisphosphate aldolases (FBA) with distinct kinetic properties. The chromosomally encoded FBA^C is the major glycolytic aldolase. The gene for the major gluconeogenic aldolase FBA^P is found on the natural plasmid pBM19 and is induced during methylotrophic growth. The crystal structures of both enzymes were solved at 2.2 Å and 2.0 Å, respectively, and they suggested amino acid residue 51 to be crucial for binding fructose-1,6-bisphosphate (FBP) as substrate and amino acid residue 140 for active site zinc atom coordination. As FBA^C and FBA^P differed at these positions, site-directed mutagenesis (SDM) was performed to exchange one or both amino acid residues of the respective proteins. The aldol cleavage reaction was negatively affected by the amino acid exchanges that led to a complete loss of glycolytic activity of FBA^P. However, both FBA^C and FBA^P maintained gluconeogenic aldol condensation activity, and the amino acid exchanges improved the catalytic efficiency of the major glycolytic aldolase FBA^C in gluconeogenic direction at least 3-fold. These results confirmed the importance of the structural differences between FBA^C and FBA^P concerning their distinct enzymatic properties. In order to investigate the physiological roles of both aldolases, the expression of their genes was repressed individually by CRISPR interference (CRISPRi). The fba ^C RNA levels were reduced by CRISPRi, but concomitantly the fba ^P RNA levels were increased. Vice versa, a similar compensatory increase of the fba ^C RNA levels was observed when fba ^P was repressed by CRISPRi. In addition, targeting fba ^P decreased tkt ^P RNA levels since both genes are cotranscribed in a bicistronic operon. However, reduced tkt ^P RNA levels were not compensated for by increased RNA levels of the chromosomal transketolase gene tkt ^C.

Developing a Riboswitch-Mediated Regulatory System for Metabolic Flux Control in Thermophilic Bacillus methanolicus

Genome-wide transcriptomic data obtained in RNA-seq experiments can serve as a reliable source for identification of novel regulatory elements such as riboswitches and promoters. Riboswitches are parts of the 5′ untranslated region of mRNA molecules that can specifically bind various metabolites and control gene expression. For that reason, they have become an attractive tool for engineering biological systems, especially for the regulation of metabolic fluxes in industrial microorganisms. Promoters in the genomes of prokaryotes are located upstream of transcription start sites and their sequences are easily identifiable based on the primary transcriptome data. Bacillus methanolicus MGA3 is a candidate for use as an industrial workhorse in methanol-based bioprocesses and its metabolism has been studied in systems biology approaches in recent years, including transcriptome characterization through RNA-seq. Here, we identify a putative lysine riboswitch in B. methanolicus, and test and characterize it. We also select and experimentally verify 10 putative B. methanolicus-derived promoters differing in their predicted strength and present their functionality in combination with the lysine riboswitch. We further explore the potential of a B. subtilis-derived purine riboswitch for regulation of gene expression in the thermophilic B. methanolicus, establishing a novel tool for inducible gene expression in this bacterium.

Extensive Reannotation of the Genome of the Model Streptomycete Streptomyces lividans TK24 Based on Transcriptome and Proteome Information

Streptomyces lividans TK24 is a relevant Gram-positive soil inhabiting bacterium and one of the model organisms of the genus Streptomyces. It is known for its potential to produce secondary metabolites, antibiotics, and other industrially relevant products. S. lividans TK24 is the plasmid-free derivative of S. lividans 66 and a close genetic relative of the strain Streptomyces coelicolor A3(2). In this study, we used transcriptome and proteome data to improve the annotation of the S. lividans TK24 genome. The RNA-seq data of primary 5'-ends of transcripts were used to determine transcription start sites (TSS) in the genome. We identified 5,424 TSS, of which 4,664 were assigned to annotated CDS and ncRNAs, 687 to antisense transcripts distributed between 606 CDS and their UTRs, 67 to tRNAs, and 108 to novel transcripts and CDS. Using the TSS data, the promoter regions and their motifs were analyzed in detail, revealing a conserved -10 (TAnnnT) and a weakly conserved -35 region (nTGACn). The analysis of the 5' untranslated region (UTRs) of S. lividans TK24 revealed 17% leaderless transcripts. Several cis-regulatory elements, like riboswitches or attenuator structures could be detected in the 5'-UTRs. The S. lividans TK24 transcriptome contains at least 929 operons. The genome harbors 27 secondary metabolite gene clusters of which 26 could be shown to be transcribed under at least one of the applied conditions. Comparison of the reannotated genome with that of the strain Streptomyces coelicolor A3(2) revealed a high degree of similarity. This study presents an extensive reannotation of the S. lividans TK24 genome based on transcriptome and proteome analyses. The analysis of TSS data revealed insights into the promoter structure, 5'-UTRs, cis-regulatory elements, attenuator structures and novel transcripts, like small RNAs. Finally, the repertoire of secondary metabolite gene clusters was examined. These data provide a basis for future studies regarding gene characterization, transcriptional regulatory networks, and usage as a secondary metabolite producing strain.

Charting the Metabolic Landscape of the Facultative Methylotroph Bacillus methanolicus

Methanol is inexpensive, is easy to transport, and can be produced both from renewable and from fossil resources without mobilizing arable lands. As such, it is regarded as a potential carbon source to transition toward a greener industrial chemistry. Metabolic engineering of bacteria and yeast able to efficiently consume methanol is expected to provide cell factories that will transform methanol into higher-value chemicals in the so-called methanol economy. Toward that goal, the study of natural methylotrophs such as Bacillus methanolicus is critical to understand the origin of their efficient methylotrophy. This knowledge will then be leveraged to transform such natural strains into new cell factories or to design methylotrophic capability in other strains already used by the industry.

Bacillus methanolicus MGA3 is a thermotolerant and relatively fast-growing methylotroph able to secrete large quantities of glutamate and lysine. These natural characteristics make B. methanolicus a good candidate to become a new industrial chassis organism, especially in a methanol-based economy. Intriguingly, the only substrates known to support B. methanolicus growth as sole sources of carbon and energy are methanol, mannitol, and, to a lesser extent, glucose and arabitol. Because fluxomics provides the most direct readout of the cellular phenotype, we hypothesized that comparing methylotrophic and nonmethylotrophic metabolic states at the flux level would yield new insights into MGA3 metabolism. In this study, we designed and performed a ¹³C metabolic flux analysis (¹³C-MFA) of the facultative methylotroph B. methanolicus MGA3 growing on methanol, mannitol, and arabitol to compare the associated metabolic states. On methanol, results showed a greater flux in the ribulose monophosphate (RuMP) pathway than in the tricarboxylic acid (TCA) cycle, thus validating previous findings on the methylotrophy of B. methanolicus. New insights related to the utilization of cyclic RuMP versus linear dissimilation pathways and between the RuMP variants were generated. Importantly, we demonstrated that the linear detoxification pathways and the malic enzyme shared with the pentose phosphate pathway have an important role in cofactor regeneration. Finally, we identified, for the first time, the metabolic pathway used to assimilate arabitol. Overall, those data provide a better understanding of this strain under various environmental conditions.

IMPORTANCE Methanol is inexpensive, is easy to transport, and can be produced both from renewable and from fossil resources without mobilizing arable lands. As such, it is regarded as a potential carbon source to transition toward a greener industrial chemistry. Metabolic engineering of bacteria and yeast able to efficiently consume methanol is expected to provide cell factories that will transform methanol into higher-value chemicals in the so-called methanol economy. Toward that goal, the study of natural methylotrophs such as Bacillus methanolicus is critical to understand the origin of their efficient methylotrophy. This knowledge will then be leveraged to transform such natural strains into new cell factories or to design methylotrophic capability in other strains already used by the industry.

Establishment of a functional system for recombinant production of secreted proteins at 50 °C in the thermophilic Bacillus methanolicus

BACKGROUND

The suitability of bacteria as microbial cell factories is dependent on several factors such as price of feedstock, product range, production yield and ease of downstream processing. The facultative methylotroph Bacillus methanolicus is gaining interest as a thermophilic cell factory for production of value-added products from methanol. The aim of this study was to expand the capabilities of B. methanolicus as a microbial cell factory by establishing a system for secretion of recombinant proteins.

RESULTS

Native and heterologous signal peptides were tested for secretion of α-amylases and proteases, and we have established the use of the thermostable superfolder green fluorescent protein (sfGFP) as a valuable reporter protein in B. methanolicus. We demonstrated functional production and secretion of recombinant proteases, α-amylases and sfGFP in B. methanolicus MGA3 at 50 °C and showed that the choice of signal peptide for optimal secretion efficiency varies between proteins. In addition, we showed that heterologous production and secretion of α-amylase from Geobacillus stearothermophilus enables B. methanolicus to grow in minimal medium with starch as the sole carbon source. An in silico signal peptide library consisting of 169 predicted peptides from B. methanolicus was generated and will be useful for future studies, but was not experimentally investigated any further here.

CONCLUSION

A functional system for recombinant production of secreted proteins at 50 °C has been established in the thermophilic B. methanolicus. In addition, an in silico signal peptide library has been generated, that together with the tools and knowledge presented in this work will be useful for further development of B. methanolicus as a host for recombinant protein production and secretion at 50 °C.

Production of Value-Added Chemicals by Bacillus methanolicus Strains Cultivated on Mannitol and Extracts of Seaweed Saccharina latissima at 50°C

The facultative methylotroph Bacillus methanolicus MGA3 has previously been genetically engineered to overproduce the amino acids L-lysine and L-glutamate and their derivatives cadaverine and γ-aminobutyric acid (GABA) from methanol at 50°C. We here explored the potential of utilizing the sugar alcohol mannitol and seaweed extract (SWE) containing mannitol, as alternative feedstocks for production of chemicals by fermentation using B. methanolicus. Extracts of the brown algae Saccharina latissima harvested in the Trondheim Fjord in Norway were prepared and found to contain 12–13 g/l of mannitol, with conductivities corresponding to a salt content of ∼2% NaCl. Initially, 12 B. methanolicus wild type strains were tested for tolerance to various SWE concentrations, and some strains including MGA3 could grow on 50% SWE medium. Non-methylotrophic and methylotrophic growth of B. methanolicus rely on differences in regulation of metabolic pathways, and we compared production titers of GABA and cadaverine under such growth conditions. Shake flask experiments showed that recombinant MGA3 strains could produce similar and higher titers of cadaverine during growth on 50% SWE and mannitol, compared to on methanol. GABA production levels under these conditions were however low compared to growth on methanol. We present the first fed-batch mannitol fermentation of B. methanolicus and production of 6.3 g/l cadaverine. Finally, we constructed a recombinant MGA3 strain synthesizing the C30 terpenoids 4,4′-diaponeurosporene and 4,4′-diapolycopene, experimentally confirming that B. methanolicus has a functional methylerythritol phosphate (MEP) pathway. Together, our results contribute to extending the range of both the feedstocks for growth and products that can be synthesized by B. methanolicus.

Transaldolase in Bacillus methanolicus: biochemical characterization and biological role in ribulose monophosphate cycle

BACKGROUND

The Gram-positive facultative methylotrophic bacterium Bacillus methanolicus uses the sedoheptulose-1,7-bisphosphatase (SBPase) variant of the ribulose monophosphate (RuMP) cycle for growth on the C₁ carbon source methanol. Previous genome sequencing of the physiologically different B. methanolicus wild-type strains MGA3 and PB1 has unraveled all putative RuMP cycle genes and later, several of the RuMP cycle enzymes of MGA3 have been biochemically characterized. In this study, the focus was on the characterization of the transaldolase (Ta) and its possible role in the RuMP cycle in B. methanolicus.

RESULTS

The Ta genes of B. methanolicus MGA3 and PB1 were recombinantly expressed in Escherichia coli, and the gene products were purified and characterized. The PB1 Ta protein was found to be active as a homodimer with a molecular weight of 54 kDa and displayed K_M of 0.74 mM and V_max of 16.3 U/mg using Fructose-6 phosphate as the substrate. In contrast, the MGA3 Ta gene, which encodes a truncated Ta protein lacking 80 amino acids at the N-terminus, showed no Ta activity. Seven different mutant genes expressing various full-length MGA3 Ta proteins were constructed and all gene products displayed Ta activities. Moreover, MGA3 cells displayed Ta activities similar as PB1 cells in crude extracts.

CONCLUSIONS

While it is well established that B. methanolicus can use the SBPase variant of the RuMP cycle this study indicates that B. methanolicus possesses Ta activity and may also operate the Ta variant of the RuMP.

Structure of a Core Promoter in Bifidobacterium longum NCC2705

Bacterial promoters consist of core sequence motifs termed -35 and -10 boxes. The consensus motifs are TTGACA and TATAAT, respectively, which were identified from leading investigations on Escherichia coli However, the consensus sequences are not likely to fit genetically divergent bacteria. The sigma factor of the genus Bifidobacterium has a characteristic polar domain in the N terminus, suggesting the possibility of specific promoter recognition. We reevaluated the structure of Bifidobacterium longum NCC2705 promoters and compared them to other bacteria. Transcriptional start sites (TSSs) of the B. longum NCC2705 strain were identified using transcriptome sequencing (RNA-Seq) analysis to extract promoter regions. Conserved motifs of a bifidobacterial promoter were determined using regions upstream of TSSs and a hidden Markov model. As a result, consensus motifs of the -35 and -10 boxes were TTGTGC and TACAAT, respectively. To assess each base of both motifs, we constructed 37 plasmids based on pKO403-TPCTcon, including the hup promoter connected with a chloramphenicol acetyltransferase as a reporter gene. This reporter assay showed two optimal motifs of the -35 and -10 boxes, namely, TTGNNN and TANNNT, respectively. We further analyzed spacer lengths between the -35 and -10 boxes via a bioinformatics approach. The spacer lengths predominant in bacteria have been generally reported to be approximately 17 bp. In contrast, the predominant spacer lengths in the genus Bifidobacterium and related species were 11 bp, in addition to 17 bp. A reporter assay to assess the spacer lengths indicated that the 11-bp spacer length produced unusually high activity.IMPORTANCE The structures of sigma factors vary among bacterial strains, indicating that recognition rules may also vary. Therefore, we investigated the promoter structure of Bifidobacterium longum NCC2705 using a bioinformatics approach and wet analyses. The most frequent and optimal motifs were similar to other bacterial consensus motifs. The optimal spacer length between the two boxes was reported to be 17 bp. It is widely applied to a bioinformatics approach for other bacteria. Unexpectedly, conserved spacer lengths were 11 bp as well as 17 bp in the genus Bifidobacterium Moreover, the sigma factor of the genus Bifidobacterium has a characteristic domain in the N terminus which may contribute to the additional functions. Hence, it would be valuable to reevaluate the promoter in other organisms.

Essentiality of the Maltase AmlE in Maltose Utilization and Its Transcriptional Regulation by the Repressor AmlR in the Acarbose-Producing Bacterium Actinoplanes sp. SE50/110

Actinoplanes sp. SE50/110 is the wild type of industrial production strains of the fine-chemical acarbose (acarviosyl-maltose), which is used as α-glucosidase inhibitor in the treatment of type II diabetes. Although maltose is an important building block of acarbose, the maltose/maltodextrin metabolism has not been studied in Actinoplanes sp. SE50/110 yet. Bioinformatic analysis located a putative maltase gene amlE (ACSP50_2474, previously named malL; Wendler et al., 2015a), in an operon with an upstream PurR/LacI-type transcriptional regulator gene, named amlR (ACSP50_2475), and a gene downstream (ACSP50_2473) encoding a GGDEF-EAL-domain-containing protein putatively involved in c-di-GMP signaling. Targeted gene deletion mutants of amlE and amlR were constructed by use of the CRISPR/Cas9 technology. By growth experiments and functional assays of ΔamlE, we could show that AmlE is essential for the maltose utilization in Actinoplanes sp. SE50/110. Neither a gene encoding a maltose phosphorylase (MalP) nor MalP enzyme activity were detected in the wild type. By this, the maltose/maltodextrin system appears to be fundamentally different from other described prokaryotic systems. By sequence similarity analysis and functional assays from the species Streptomyces lividans TK23, S. coelicolor A3(2) and S. glaucescens GLA.O, first hints for a widespread lack of MalP and presence of AmlE in the class Actinobacteria were given. Transcription of the aml operon is significantly repressed in the wild type when growing on glucose and repression is absent in an ΔamlR deletion mutant. Although AmlR apparently is a local transcriptional regulator of the aml operon, the ΔamlR strain shows severe growth inhibitions on glucose and – concomitantly – differential transcription of several genes of various functional classes. We ascribe these effects to ACSP50_2473, which is localized downstream of amlE and presumably involved in the metabolism of the second messenger c-di-GMP. It can be assumed, that maltose does not only represent the most important carbon source of Actinoplanes sp. SE50/110, but that its metabolism is coupled to the nucleotide messenger system of c-di-GMP.

The Draft Genome Sequence of Methylophilus sp. D22, Capable of Growing Under High Concentration of Methanol

In this study, a wild type Methylophilus sp. strain D22 belonging to the family Methylophilus was isolated and characterized, which shows high tolerance towards methanol, as it can grow under 50 g/L of methanol. Methylophilus sp. strain D22 was isolated from the lake sludge in Nanjing Tech University, China. The assembled draft genome contains one circular chromosome with 3,004,398 bp, 49.7% of GC content, and 2107 predicted encoding proteins. Sequence-based genomic analysis demonstrates that the assimilation pathway of ribulose monophosphate (RuMP) pathway and dissimilation pathway of tetrahydromethanopterin (H₄MPT) pathway are co-existing and contribute to the high methanol utilization efficiency.

Proteomic Response to Rising Temperature in the Marine Cyanobacterium Synechococcus Grown in Different Nitrogen Sources

Synechococcus is one of the most important contributors to global primary productivity, and ocean warming is predicted to increase abundance and distribution of Synechococcus in the ocean. Here, we investigated molecular response of an oceanic Synechococcus strain WH8102 grown in two nitrogen sources (nitrate and urea) under present (25°C) and predicted future (28°C) temperature conditions using an isobaric tag (IBT)-based quantitative proteomic approach. Rising temperature decreased growth rate, contents of chlorophyll a, protein and sugar in the nitrate-grown cells, but only decreased protein content and significantly increased zeaxanthin content of the urea-grown cells. Expressions of CsoS2 protein involved in carboxysome formation and ribosomal subunits in both nitrate- and urea-grown cells were significantly decreased in rising temperature, whereas carbohydrate selective porin and sucrose-phosphate synthase (SPS) were remarkably up-regulated, and carbohydrate degradation associated proteins, i.e., glycogen phosphorylase kinase, fructokinase and glucose-6-phosphate dehydrogenase, were down-regulated in the urea-grown cells. Rising temperature also increased expressions of three redox-sensitive enzymes (peroxiredoxin, thioredoxin, and CP12) in both nitrate- and urea-grown cells. Our results indicated that rising temperature did not enhance cell growth of Synechococcus; on the contrary, it impaired cell functions, and this might influence cell abundance and distribution of Synechococcus in a future ocean.

Characterization of D-Arabitol as Newly Discovered Carbon Source of Bacillus methanolicus

Bacillus methanolicus is a Gram-positive, thermophilic, methanol-utilizing bacterium. As a facultative methylotroph, B. methanolicus is also known to utilize D-mannitol, D-glucose and, as recently discovered, sugar alcohol D-arabitol. While metabolic pathways for utilization of methanol, mannitol and glucose are known, catabolism of arabitol has not yet been characterized in B. methanolicus. In this work we present the elucidation of this hitherto uncharted pathway. In order to confirm our predictions regarding genes coding for arabitol utilization, we performed differential gene expression analysis of B. methanolicus MGA3 cells grown on arabitol as compared to mannitol via transcriptome sequencing (RNA-seq). We identified a gene cluster comprising eight genes that was up-regulated during growth with arabitol as a sole carbon source. The RNA-seq results were subsequently confirmed via qRT-PCR experiments. The transcriptional organization of the gene cluster identified via RNA-seq was analyzed and it was shown that the arabitol utilization genes are co-transcribed in an operon that spans from BMMGA3_RS07325 to BMMGA3_RS07365. Since gene deletion studies are currently not possible in B. methanolicus, two complementation experiments were performed in an arabitol negative Corynebacterium glutamicum strain using the four genes discovered via RNA-seq analysis as coding for a putative PTS for arabitol uptake (BMMGA3_RS07330, BMMGA3_RS07335, and BMMGA3_RS07340 renamed to atlABC) and a putative arabitol phosphate dehydrogenase (BMMGA3_RS07345 renamed to atlD). C. glutamicum is a natural D-arabitol utilizer that requires arabitol dehydrogenase MtlD for arabitol catabolism. The C. glutamicum mtlD deletion mutant was chosen for complementation experiments. Heterologous expression of atlABCD as well as the arabitol phosphate dehydrogenase gene atlD from B. methanolicus alone restored growth of the C. glutamicum ΔmtlD mutant with arabitol. Furthermore, D-arabitol phosphate dehydrogenase activities could be detected in crude extracts of B. methanolicus and these were higher in arabitol-grown cells than in methanol- or mannitol-grown cells. Thus, B. methanolicus possesses an arabitol inducible operon encoding, amongst others, a putative PTS system and an arabitol phosphate dehydrogenase for uptake and activation of arabitol as growth substrate.

Function of L-Pipecolic Acid as Compatible Solute in Corynebacterium glutamicum as Basis for Its Production Under Hyperosmolar Conditions

Pipecolic acid or L-PA is a cyclic amino acid derived from L-lysine which has gained interest in the recent years within the pharmaceutical and chemical industries. L-PA can be produced efficiently using recombinant Corynebacterium glutamicum strains by expanding the natural L-lysine biosynthetic pathway. L-PA is a six-membered ring homolog of the five-membered ring amino acid L-proline, which serves as compatible solute in C. glutamicum.

Here, we show that de novo synthesized or externally added L-PA partially is beneficial for growth under hyper-osmotic stress conditions. C. glutamicum cells accumulated L-PA under elevated osmotic pressure and released it after an osmotic down shock. In the absence of the mechanosensitive channel YggB intracellular L-PA concentrations increased and its release after osmotic down shock was slower. The proline permease ProP was identified as a candidate L-PA uptake system since RNAseq analysis revealed increased proP RNA levels upon L-PA production. Under hyper-osmotic conditions, a ΔproP strain showed similar growth behavior than the parent strain when L-proline was added externally. By contrast, the growth impairment of the ΔproP strain under hyper-osmotic conditions could not be alleviated by addition of L-PA unless proP was expressed from a plasmid. This is commensurate with the view that L-proline can be imported into the C. glutamicum cell by ProP and other transporters such as EctP and PutP, while ProP appears of major importance for L-PA uptake under hyper-osmotic stress conditions.

Metabolic and proteomic alteration in phytohormone-producing endophytic Bacillus amyloliquefaciens RWL-1 during methanol utilization

INTRODUCTION

Methanol utilization by bacteria is important for various industrial processes. Methylotrophic bacteria are taxonomically diverse and some species promote plant growth and induce stress tolerance. However, methylotrophic potential of bacterial endophytes is poorly understood.

OBJECTIVE

The current study aimed to evaluate the metabolomic and proteomic changes in endophytic Bacillus amyloliquefaciens RWL-1 caused by its methanol utilization and the resultant influence on its phytohormone production.

METHODS

B. amyloliquefaciens RWL-1 was grown in LB medium with different concentrations [0 (control), 0.5, 1, 1.5, 2, 2.5, 3, 3.5, and 4%) of methanol to examine its methylotrophic potential. SDS-PAGE analysis was carried out for bacterial protein confirmation. Moreover, the phytohormones (indole 3 acetic acid (IAA), gibberellins (GAs), abscisic acid (ABA)) produced by RWL-1 in methanol supplemented medium were quantified by GC-MS/SIM (6890N Network GC system, and 5973 Network Mass Selective Detector; Agilent Technologies, Santa Clara, CA, USA), while the antioxidants were estimated spectrophotometrically (T60 UV-VIS spectrophotometer, Leicester, UK). The amino acid quantification was carried out by amino acid analyzer (HITACHI L-8900, Japan). Furthermore, Nano-liquid chromatography (LC)-MS/MS analysis was performed with an Agilent system (Wilmington, DE, USA) for proteomic analysis while mascot algorithm (Matrix science, USA) was used to identify peptide sequences present in the protein sequence database.

RESULTS

RWL-1 showed significant growth in media supplemented with 2 and 3.5% methanol, when compared with other concentrations. Mass spectroscopy analysis revealed that RWL-1 utilizes methanol efficiently as a carbon source. In the presence of methanol, RWL-1 produced significantly higher levels of IAA but lower levels of ABA, when compared with the control. Further, enzymatic antioxidants and functional amino acids were significantly up-regulated, with predominant expression of glutamic acid and alanine. Nano-liquid chromatography, quadrupole time-of-flight analysis, and quantitative analysis of methanol-treated bacterial cells showed expression of eight different types of proteins, including detoxification proteins, unrecognized and unclassified enzymes with antioxidant properties, proteases, metabolism enzymes, ribosomal proteins, antioxidant proteins, chaperones, and heat shock proteins.

CONCLUSION

Results demonstrate that RWL-1 can significantly enhance its growth by utilizing methanol, and could produce phytohormones when growing in methanol-supplemented media, with increased expression of specific proteins and different biochemicals. These results will be useful in devising strategies for utilizing methylotrophic bacterial endophytes as alternative promoters of plant growth. Understanding RWL-1 ability to utilize methanol. The survival and phytohormones production by Bacillus amyloliquefaciens RWL-1 in methanol supplemented media whistle inducing metabolic and proteomic changes.

DNA sequence encodes the position of DNA supercoils

The three-dimensional organization of DNA is increasingly understood to play a decisive role in vital cellular processes. Many studies focus on the role of DNA-packaging proteins, crowding, and confinement in arranging chromatin, but structural information might also be directly encoded in bare DNA itself. Here, we visualize plectonemes (extended intertwined DNA structures formed upon supercoiling) on individual DNA molecules. Remarkably, our experiments show that the DNA sequence directly encodes the structure of supercoiled DNA by pinning plectonemes at specific sequences. We develop a physical model that predicts that sequence-dependent intrinsic curvature is the key determinant of pinning strength and demonstrate this simple model provides very good agreement with the data. Analysis of several prokaryotic genomes indicates that plectonemes localize directly upstream of promoters, which we experimentally confirm for selected promotor sequences. Our findings reveal a hidden code in the genome that helps to spatially organize the chromosomal DNA.

Genetic Tools and Techniques for Recombinant Expression in Thermophilic Bacillaceae

Although Escherichia coli and Bacillus subtilis are the most prominent bacterial hosts for recombinant protein production by far, additional species are being explored as alternatives for production of difficult-to-express proteins. In particular, for thermostable proteins, there is a need for hosts able to properly synthesize, fold, and excrete these in high yields, and thermophilic Bacillaceae represent one potentially interesting group of microorganisms for such purposes. A number of thermophilic Bacillaceae including B.methanolicus, B.coagulans, B.smithii, B.licheniformis, Geobacillus thermoglucosidasius, G. kaustophilus, and G. stearothermophilus are investigated concerning physiology, genomics, genetic tools, and technologies, altogether paving the way for their utilization as hosts for recombinant production of thermostable and other difficult-to-express proteins. Moreover, recent successful deployments of CRISPR/Cas9 in several of these species have accelerated the progress in their metabolic engineering, which should increase their attractiveness for future industrial-scale production of proteins. This review describes the biology of thermophilic Bacillaceae and in particular focuses on genetic tools and methods enabling use of these organisms as hosts for recombinant protein production.

Transcriptome sequencing of the human pathogen Corynebacterium diphtheriae NCTC 13129 provides detailed insights into its transcriptional landscape and into DtxR-mediated transcriptional regulation

Background

The human pathogen Corynebacterium diphtheriae is the causative agent of diphtheria. In the 1990s a large diphtheria outbreak in Eastern Europe was caused by the strain C. diphtheriae NCTC 13129. Although the genome was sequenced more than a decade ago, not much is known about its transcriptome. Our aim was to use transcriptome sequencing (RNA-Seq) to close this knowledge gap and gain insights into the transcriptional landscape of a C. diphtheriae tox⁺ strain.

Results

We applied two different RNA-Seq techniques, one to retrieve 5′-ends of primary transcripts and the other to characterize the whole transcriptional landscape in order to gain insights into various features of the C. diphtheriae NCTC 13129 transcriptome. By examining the data we identified 1656 transcription start sites (TSS), of which 1202 were assigned to genes and 454 to putative novel transcripts. By using the TSS data promoter regions recognized by the housekeeping sigma factor σ^A and its motifs were analyzed in detail, revealing a well conserved −10 but an only weakly conserved −35 motif, respectively. Furthermore, with the TSS data 5’-UTR lengths were explored. The observed 5’-UTRs range from zero length (leaderless transcripts), which make up 20% of all genes, up to over 450 nt long leaders, which may harbor regulatory functions. The C. diphtheriae transcriptome consists of 471 operons which are further divided into 167 sub-operon structures. In a differential expression analysis approach, we discovered that genetic disruption of the iron-sensing transcription regulator DtxR, which controls expression of diphtheria toxin (DT), causes a strong influence on general gene expression. Nearly 15% of the genome is differentially transcribed, indicating that DtxR might have other regulatory functions in addition to regulation of iron metabolism and DT. Furthermore, our findings shed light on the transcriptional landscape of the DT encoding gene tox and present evidence for two tox antisense RNAs, which point to a new way of transcriptional regulation of toxin production.

Conclusions

This study presents extensive insights into the transcriptome of C. diphtheriae and provides a basis for future studies regarding gene characterization, transcriptional regulatory networks, and regulation of the tox gene in particular.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4481-8) contains supplementary material, which is available to authorized users.

RNAseq analysis of α-proteobacterium Gluconobacter oxydans 621H

Background

The acetic acid bacterium Gluconobacter oxydans 621H is characterized by its exceptional ability to incompletely oxidize a great variety of carbohydrates in the periplasm. The metabolism of this α-proteobacterium has been characterized to some extent, yet little is known about its transcriptomes and related data. In this study, we applied two different RNAseq approaches. Primary transcriptomes enriched for 5′-ends of transcripts were sequenced to detect transcription start sites, which allow subsequent analysis of promoter motifs, ribosome binding sites, and 5´-UTRs. Whole transcriptomes were sequenced to identify expressed genes and operon structures.

Results

Sequencing of primary transcriptomes of G. oxydans revealed 2449 TSSs, which were classified according to their genomic context followed by identification of promoter and ribosome binding site motifs, analysis of 5´-UTRs including validation of predicted cis-regulatory elements and correction of start codons. 1144 (41%) of all genes were found to be expressed monocistronically, whereas 1634 genes were organized in 571 operons. Together, TSSs and whole transcriptome data were also used to identify novel intergenic (18), intragenic (328), and antisense transcripts (313).

Conclusions

This study provides deep insights into the transcriptional landscapes of G. oxydans. The comprehensive transcriptome data, which we made publicly available, facilitate further analysis of promoters and other regulatory elements. This will support future approaches for rational strain development and targeted gene expression in G. oxydans. The corrections of start codons further improve the high quality genome reference and support future proteome analysis.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4415-x) contains supplementary material, which is available to authorized users.

Development of a formaldehyde biosensor with application to synthetic methylotrophy

Formaldehyde is a prevalent environmental toxin and a key intermediate in single carbon metabolism. The ability to monitor formaldehyde concentration is, therefore, of interest for both environmental monitoring and for metabolic engineering of native and synthetic methylotrophs, but current methods suffer from low sensitivity, complex workflows, or require expensive analytical equipment. Here we develop a formaldehyde biosensor based on the FrmR repressor protein and cognate promoter of Escherichia coli. Optimization of the native repressor binding site and regulatory architecture enabled detection at levels as low as 1 µM. We then used the sensor to benchmark the in vivo activity of several NAD-dependent methanol dehydrogenase (Mdh) variants, the rate-limiting enzyme that catalyzes the first step of methanol assimilation. In order to use this biosensor to distinguish individuals in a mixed population of Mdh variants, we developed a strategy to prevent cross-talk by using glutathione as a formaldehyde sink to minimize intercellular formaldehyde diffusion. Finally, we applied this biosensor to balance expression of mdh and the formaldehyde assimilation enzymes hps and phi in an engineered E. coli strain to minimize formaldehyde build-up while also reducing the burden of heterologous expression. This biosensor offers a quick and simple method for sensitively detecting formaldehyde, and has the potential to be used as the basis for directed evolution of Mdh and dynamic formaldehyde control strategies for establishing synthetic methylotrophy.

Detailed transcriptome analysis of the plant growth promoting Paenibacillus riograndensis SBR5 by using RNA-seq technology

BACKGROUND

The plant growth promoting rhizobacterium Paenibacillus riograndensis SBR5 is a promising candidate to serve as crop inoculant. Despite its potential in providing environmental and economic benefits, the species P. riograndensis is poorly characterized. Here, we performed for the first time a detailed transcriptome analysis of P. riograndensis SBR5 using RNA-seq technology.

RESULTS

RNA was isolated from P. riograndensis SBR5 cultivated under 15 different growth conditions and combined together in order to analyze an RNA pool representing a large set of expressed genes. The resultant total RNA was used to generate 2 different libraries, one enriched in 5'-ends of the primary transcripts and the other representing the whole transcriptome. Both libraries were sequenced and analyzed to identify the conserved sequences of ribosome biding sites and translation start motifs, and to elucidate operon structures present in the transcriptome of P. riograndensis. Sequence analysis of the library enriched in 5'-ends of the primary transcripts was used to identify 1082 transcription start sites (TSS) belonging to novel transcripts and allowed us to determine a promoter consensus sequence and regulatory sequences in 5' untranslated regions including riboswitches. A putative thiamine pyrophosphate dependent riboswitch upstream of the thiamine biosynthesis gene thiC was characterized by translational fusion to a fluorescent reporter gene and shown to function in P. riograndensis SBR5.

CONCLUSIONS

Our RNA-seq analysis provides insight into the P. riograndensis SBR5 transcriptome at the systems level and will be a valuable basis for differential RNA-seq analysis of this bacterium.

Time-resolved transcriptome analysis and lipid pathway reconstruction of the oleaginous green microalga Monoraphidium neglectum reveal a model for triacylglycerol and lipid hyperaccumulation

BACKGROUND

Oleaginous microalgae are promising production hosts for the sustainable generation of lipid-based bioproducts and as bioenergy carriers such as biodiesel. Transcriptomics of the lipid accumulation phase, triggered efficiently by nitrogen starvation, is a valuable approach for the identification of gene targets for metabolic engineering.

RESULTS

An explorative analysis of the detailed transcriptional response to different stages of nitrogen availability was performed in the oleaginous green alga Monoraphidium neglectum. Transcript data were correlated with metabolic data for cellular contents of starch and of different lipid fractions. A pronounced transcriptional down-regulation of photosynthesis became apparent in response to nitrogen starvation, whereas glucose catabolism was found to be up-regulated. An in-depth reconstruction and analysis of the pathways for glycerolipid, central carbon, and starch metabolism revealed that distinct transcriptional changes were generally found only for specific steps within a metabolic pathway. In addition to pathway analyses, the transcript data were also used to refine the current genome annotation. The transcriptome data were integrated into a database and complemented with data for other microalgae which were also subjected to nitrogen starvation. It is available at https://tdbmn.cebitec.uni-bielefeld.de.

CONCLUSIONS

Based on the transcriptional responses to different stages of nitrogen availability, a model for triacylglycerol and lipid hyperaccumulation is proposed, which involves transcriptional induction of thioesterases, differential regulation of lipases, and a re-routing of the central carbon metabolism. Over-expression of distinct thioesterases was identified to be a potential strategy to increase the oleaginous phenotype of M. neglectum, and furthermore specific lipases were identified as potential targets for future metabolic engineering approaches.

Genome improvement of the acarbose producer Actinoplanes sp. SE50/110 and annotation refinement based on RNA-seq analysis

Actinoplanes sp. SE50/110 is the natural producer of acarbose, which is used in the treatment of diabetes mellitus type II. However, until now the transcriptional organization and regulation of the acarbose biosynthesis are only understood rudimentarily. The genome sequence of Actinoplanes sp. SE50/110 was known before, but was resequenced in this study to remove assembly artifacts and incorrect base callings. The annotation of the genome was refined in a multi-step approach, including modern bioinformatic pipelines, transcriptome and proteome data. A whole transcriptome RNA-seq library as well as an RNA-seq library enriched for primary 5'-ends were used for the detection of transcription start sites, to correct tRNA predictions, to identify novel transcripts like small RNAs and to improve the annotation through the correction of falsely annotated translation start sites. The transcriptome data sets were also applied to identify 31 cis-regulatory RNA structures, such as riboswitches or RNA thermometers as well as three leaderless transcribed short peptides found in putative attenuators upstream of genes for amino acid biosynthesis. The transcriptional organization of the acarbose biosynthetic gene cluster was elucidated in detail and fourteen novel biosynthetic gene clusters were suggested. The accurate genome sequence and precise annotation of the Actinoplanes sp. SE50/110 genome will be the foundation for future genetic engineering and systems biology studies.

l-lysine production by Bacillus methanolicus: Genome-based mutational analysis and l-lysine secretion engineering

Bacillus methanolicus is a methylotrophic bacterium with an increasing interest in academic research and for biotechnological applications. This bacterium was previously applied for methanol-based production of l-glutamate, l-lysine and the five-carbon diamine cadaverine by wild type, classical mutant and recombinant strains. The genomes of two different l-lysine secreting B. methanolicus classical mutant strains, NOA2#13A52-8A66 and M168-20, were sequenced. We focused on mutational mapping in genes present in l-lysine and other relevant amino acid biosynthetic pathways, as well as in the primary cell metabolism important for precursor supply. In addition to mutations in the aspartate pathway genes dapG, lysA and hom-1, new mutational target genes like alr, proA, proB1, leuC, odhA and pdhD were identified. Surprisingly, no mutations were found in the putative l-lysine transporter gene lysE^MGA3. Inspection of the wild type B. methanolicus strain PB1 genome sequence identified two homologous putative l-lysine transporter genes, lysE^PB1 and lysE2^PB1. The biological role of these putative l-lysine transporter genes, together with the heterologous l-lysine exporter gene lysE^Cg from Corynebacterium glutamicum, were therefore investigated. Our results demonstrated that the titer of secreted l-lysine in B. methanolicus was significantly increased by overexpression of lysE^Cg while overexpression of lysE^MGA3, lysE^PB1 and lysE2^PB1 had no measurable effect.

Genome-Based Genetic Tool Development for Bacillus methanolicus: Theta- and Rolling Circle-Replicating Plasmids for Inducible Gene Expression and Application to Methanol-Based Cadaverine Production

Bacillus methanolicus is a thermophilic methylotroph able to overproduce amino acids from methanol, a substrate not used for human or animal nutrition. Based on our previous RNA-seq analysis a mannitol inducible promoter and a putative mannitol activator gene mtlR were identified. The mannitol inducible promoter was applied for controlled gene expression using fluorescent reporter proteins and a flow cytometry analysis, and improved by changing the -35 promoter region and by co-expression of the mtlR regulator gene. For independent complementary gene expression control, the heterologous xylose-inducible system from B. megaterium was employed and a two-plasmid gene expression system was developed. Four different replicons for expression vectors were compared with respect to their copy number and stability. As an application example, methanol-based production of cadaverine was shown to be improved from 6.5 to 10.2 g/L when a heterologous lysine decarboxylase gene cadA was expressed from a theta-replicating rather than a rolling-circle replicating vector. The current work on inducible promoter systems and compatible theta- or rolling circle-replicating vectors is an important extension of the poorly developed B. methanolicus genetic toolbox, valuable for genetic engineering and further exploration of this bacterium.

Transcriptome Analysis Reveals the Genetic Basis of the Resveratrol Biosynthesis Pathway in an Endophytic Fungus (Alternaria sp. MG1) Isolated from Vitis vinifera

Alternaria sp. MG1, an endophytic fungus previously isolated from Merlot grape, produces resveratrol from glucose, showing similar metabolic flux to the phenylpropanoid biosynthesis pathway, currently found solely in plants. In order to identify the resveratrol biosynthesis pathway in this strain at the gene level, de novo transcriptome sequencing was conducted using Illumina paired-end sequencing. A total of 22,954,434 high-quality reads were assembled into contigs and 18,570 unigenes were identified. Among these unigenes, 14,153 were annotated in the NCBI non-redundant protein database and 5341 were annotated in the Swiss-Prot database. After KEGG mapping, 2701 unigenes were mapped onto 115 pathways. Eighty-four unigenes were annotated in major pathways from glucose to resveratrol, coding 20 enzymes for glycolysis, 10 for phenylalanine biosynthesis, 4 for phenylpropanoid biosynthesis, and 4 for stilbenoid biosynthesis. Chalcone synthase was identified for resveratrol biosynthesis in this strain, due to the absence of stilbene synthase. All the identified enzymes indicated a reasonable biosynthesis pathway from glucose to resveratrol via glycolysis, phenylalanine biosynthesis, phenylpropanoid biosynthesis, and stilbenoid pathways. These results provide essential evidence for the occurrence of resveratrol biosynthesis in Alternaria sp. MG1 at the gene level, facilitating further elucidation of the molecular mechanisms involved in this strain's secondary metabolism.

Quantitative metabolomics of the thermophilic methylotroph Bacillus methanolicus

Background

The gram-positive bacterium Bacillus methanolicus MGA3 is a promising candidate for methanol-based biotechnologies. Accurate determination of intracellular metabolites is crucial for engineering this bacteria into an efficient microbial cell factory. Due to the diversity of chemical and cell properties, an experimental protocol validated on B. methanolicus is needed. Here a systematic evaluation of different techniques for establishing a reliable basis for metabolome investigations is presented.

Results

Metabolome analysis was focused on metabolites closely linked with B. methanolicus central methanol metabolism. As an alternative to cold solvent based procedures, a solvent-free quenching strategy using stainless steel beads cooled to −20 °C was assessed. The precision, the consistency of the measurements, and the extent of metabolite leakage from quenched cells were evaluated in procedures with and without cell separation. The most accurate and reliable performance was provided by the method without cell separation, as significant metabolite leakage occurred in the procedures based on fast filtration. As a biological test case, the best protocol was used to assess the metabolome of B. methanolicus grown in chemostat on methanol at two different growth rates and its validity was demonstrated.

Conclusion

The presented protocol is a first and helpful step towards developing reliable metabolomics data for thermophilic methylotroph B. methanolicus. This will definitely help for designing an efficient methylotrophic cell factory.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-016-0483-x) contains supplementary material, which is available to authorized users.

Methanol-based cadaverine production by genetically engineered Bacillus methanolicus strains

Methanol is regarded as an attractive substrate for biotechnological production of value-added bulk products, such as amino acids and polyamines. In the present study, the methylotrophic and thermophilic bacterium Bacillus methanolicus was engineered into a microbial cell factory for the production of the platform chemical 1,5-diaminopentane (cadaverine) from methanol. This was achieved by the heterologous expression of the Escherichia coli genes cadA and ldcC encoding two different lysine decarboxylase enzymes, and by increasing the overall L-lysine production levels in this host. Both CadA and LdcC were functional in B. methanolicus cultivated at 50°C and expression of cadA resulted in cadaverine production levels up to 500 mg l⁻¹ during shake flask conditions. A volume-corrected concentration of 11.3 g l⁻¹ of cadaverine was obtained by high-cell density fed-batch methanol fermentation. Our results demonstrated that efficient conversion of L-lysine into cadaverine presumably has severe effects on feedback regulation of the L-lysine biosynthetic pathway in B. methanolicus. By also investigating the cadaverine tolerance level, B. methanolicus proved to be an exciting alternative host and comparable to the well-known bacterial hosts E. coli and Corynebacterium glutamicum. This study represents the first demonstration of microbial production of cadaverine from methanol.