Genome-wide determination of transcription start sites reveals new insights into promoter structures in the actinomycete Corynebacterium glutamicum

Structure-Guided Protein Engineering of Glyceraldehyde-3-phosphate Dehydrogenase from Corynebacterium glutamicum for Dual NAD/NADP Cofactor Specificity

Since the discovery of l-glutamate-producing Corynebacterium glutamicum, it has evolved to be an industrial workhorse. For biobased chemical production, suppling sufficient amounts of the NADPH cofactor is crucial. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), a glycolytic enzyme that converts glyceraldehyde-3-phosphate (G3P) to 1,3-bisphosphoglycerate and produces NADH, is a major prospective solution for the cofactor imbalance issue. In this study, we determined the crystal structure of GAPDH from C. glutamicum ATCC13032 (CgGAPDH). Based on the structural information, we generated six CgGAPDH variants, CgGAPDHL36S, CgGAPDHL36S/T37K, CgGAPDHL36S/T37K/P192S, CgGAPDHL36S/T37K/F100V/P192S, CgGAPDHL36S/T37K/F100L/P192S, and CgGAPDHL36S/T37K/F100I/P192S, that can produce both NADH and NAPDH. The final CgGAPDHL36S/T37K/F100V/P192S variant showed a 212-fold increase in enzyme activity for NADP as well as 200% and 30% increased activity for the G3P substrate under NAD and NADP cofactor conditions, respectively. In addition, crystal structures of CgGAPDH variants in complex with NAD(P) permit the elucidation of differences between wild-type CgGAPDH and variants in relation to cofactor stabilization.

Deep annotation of long noncoding RNAs by assembling RNA-seq and small RNA-seq data

Long noncoding RNAs (lncRNAs) are increasingly being recognized as modulators in various biological processes. However, due to their low expression, their systematic characterization is difficult to determine. Here, we performed transcript annotation by a newly developed computational pipeline, termed RNA-seq and small RNA-seq combined strategy (RSCS), in a wide variety of cellular contexts. Thousands of high-confidence potential novel transcripts were identified by the RSCS, and the reliability of the transcriptome was verified by analysis of transcript structure, base composition, and sequence complexity. Evidenced by the length comparison, the frequency of the core promoter and the polyadenylation signal motifs, and the locations of transcription start and end sites, the transcripts appear to be full length. Furthermore, taking advantage of our strategy, we identified a large number of endogenous retrovirus-associated lncRNAs, and a novel endogenous retrovirus-lncRNA that was functionally involved in control of Yap1 expression and essential for early embryogenesis was identified. In summary, the RSCS can generate a more complete and precise transcriptome, and our findings greatly expanded the transcriptome annotation for the mammalian community.

Overlapping SigH and SigE sigma factor regulons in Corynebacterium glutamicum

The sigma H (σΗ) and sigma E (σE) subunits of Corynebacterium glutamicum RNA polymerase belong to Group 4 of sigma factors, also called extracytoplasmic function (ECF) sigma factors. Genes of the C. glutamicum σΗ regulon that are involved in heat and oxidative stress response have already been defined, whereas the genes of the σE regulon, which is involved in cell surface stress response, have not been explored until now. Using the C. glutamicum RES167 strain and its derivative C. glutamicum ΔcseE with a deletion in the anti-σΕ gene, differential gene expression was analyzed by RNA sequencing. We found 296 upregulated and 398 downregulated genes in C. glutamicum ΔcseE compared to C. glutamicum RES167. To confirm the functional link between σΕ and the corresponding promoters, we tested selected promoters using the in vivo two-plasmid system with gfpuv as a reporter gene and by in vitro transcription. Analyses with RNAP+σΗ and RNAP+σΕ, which were previously shown to recognize similar promoters, proved that the σΗ and σE regulons significantly overlap. The σE-controlled genes were found to be involved for example in protein quality control (dnaK, dnaJ2, clpB, and clpC), the regulation of Clp proteases (clgR), and membrane integrity maintenance. The single-promoter analyses with σΗ and σΕ revealed that there are two groups of promoters: those which are exclusively σΗ-specific, and the other group of promoters, which are σΗ/σE-dependent. No exclusively σE-dependent promoter was detected. We defined the consensus sequences of exclusively σΗ-regulated promotors to be -35 GGAAt and - 10 GTT and σΗ/σE-regulated promoters to be -35 GGAAC and - 10 cGTT. Fifteen genes were found to belong to the σΗ/σΕ regulon. Homology modeling showed that there is a specific interaction between Met170 in σΗ and the nucleotides -31 and - 30 within the non-coding strand (AT or CT) of the σΗ-dependent promoters. In σE, Arg185 was found to interact with the nucleotides GA at the same positions in the σE-dependent promoters.

Functional evaluation of LTR-derived lncRNAs in porcine oocytes and zygotes with RNA-seq and small RNA-seq

Long noncoding RNAs (lncRNAs) are increasingly being recognized as modulators of early embryonic development in mammals. However, they are seldom investigated in pigs. Here, to annotate full-length RNA transcripts, we performed annotation using a newly developed computational pipeline—an RNA-seq and small RNA-seq combined strategy—using our previously obtained RNA-seq and small RNA-seq data from porcine oocytes and zygotes. As evidenced by the length comparison, the frequency of the core promoter, and the polyadenylation signal motifs, the transcripts appear to be full-length. Furthermore, our strategy allowed the identification of a large number of endogenous retrovirus-associated lncRNAs (ERV-lncRNAs) and found that some of them were highly expressed in porcine zygotes, as compared to oocytes. Through the knockdown strategy, two ERV-lncRNAs (TCONS_00035465 and TCONS_00031520) were identified as playing potential roles in the early embryo development of pigs, laying a foundation for future research.

Identification of Rhodococcus erythropolis Promoters Controlled by Alternative Sigma Factors Using In Vivo and In Vitro Systems and Heterologous RNA Polymerase

Rhodococcus erythropolis CCM2595 is a bacterial strain, which has been studied for its capability to degrade phenol and other toxic aromatic compounds. Its cell wall contains mycolic acids, which are also an attribute of other bacteria of the Mycolata group, such as Corynebacterium and Mycobacterium species. We suppose that many genes upregulated by phenol stress in R. erythropolis are controlled by the alternative sigma factors of RNA polymerase, which are active in response to the cell envelope or oxidative stress. We developed in vitro and in vivo assays to examine the connection between the stress sigma factors and genes activated by various extreme conditions, e.g., heat, cell surface, and oxidative stress. These assays are based on the procedures of such tests carried out in the related species, Corynebacterium glutamicum. We showed that the R. erythropolis CCM2595 genes frmB1 and frmB2, which encode S-formylglutathione hydrolases (named corynomycolyl transferases in C. glutamicum), are controlled by SigD, just like the homologous genes cmt1 and cmt2 in C. glutamicum. The new protocol of the in vivo and in vitro assays will enable us to classify R. erythropolis promoters according to their connection to sigma factors and to assign the genes to the corresponding sigma regulons. The complex stress responses, such as that induced by phenol, could, thus, be analyzed with respect to the gene regulation by sigma factors.

The ldhA Gene Encoding Fermentative l-Lactate Dehydrogenase in Corynebacterium Glutamicum Is Positively Regulated by the Global Regulator GlxR

Bacterial metabolism shifts from aerobic respiration to fermentation at the transition from exponential to stationary growth phases in response to limited oxygen availability. Corynebacterium glutamicum, a Gram-positive, facultative aerobic bacterium used for industrial amino acid production, excretes l-lactate, acetate, and succinate as fermentation products. The ldhA gene encoding l-lactate dehydrogenase is solely responsible for l-lactate production. Its expression is repressed at the exponential phase and prominently induced at the transition phase. ldhA is transcriptionally repressed by the sugar-phosphate-responsive regulator SugR and l-lactate-responsive regulator LldR. Although ldhA expression is derepressed even at the exponential phase in the sugR and lldR double deletion mutant, a further increase in its expression is still observed at the stationary phase, implicating the action of additional transcription regulators. In this study, involvement of the cAMP receptor protein-type global regulator GlxR in the regulation of ldhA expression was investigated. The GlxR-binding site found in the ldhA promoter was modified to inhibit or enhance binding of GlxR. The ldhA promoter activity and expression of ldhA were altered in proportion to the binding affinity of GlxR. Similarly, l-lactate production was also affected by the binding site modification. Thus, GlxR was demonstrated to act as a transcriptional activator of ldhA.

Physiology and Transcriptional Analysis of (p)ppGpp-Related Regulatory Effects in Corynebacterium glutamicum

The alarmone species ppGpp and pppGpp are elementary components of bacterial physiology as they both coordinate the bacterial stress response and serve as fine-tuners of general metabolism during conditions of balanced growth. Since the regulation of (p)ppGpp metabolism and the effects of (p)ppGpp on cellular processes are highly complex and show massive differences between bacterial species, the underlying molecular mechanisms have so far only been insufficiently investigated for numerous microorganisms. In this study, (p)ppGpp physiology in the actinobacterial model organism Corynebacterium glutamicum was analyzed by phenotypic characterization and RNAseq-based transcriptome analysis. Total nutrient starvation was identified as the most effective method to induce alarmone production, whereas traditional induction methods such as the addition of serine hydroxamate (SHX) or mupirocin did not show a strong accumulation of (p)ppGpp. The predominant alarmone in C. glutamicum represents guanosine tetraphosphate, whose stress-associated production depends on the presence of the bifunctional RSH enzyme Rel. Interestingly, in addition to ppGpp, another substance yet not identified accumulated strongly under inducing conditions. A C. glutamicum triple mutant (Δrel,ΔrelS,ΔrelH) unable to produce alarmones [(p)ppGpp⁰ strain] exhibited unstable growth characteristics and interesting features such as an influence of illumination on its physiology, production of amino acids as well as differences in vitamin and carotenoid production. Differential transcriptome analysis using RNAseq provided numerous indications for the molecular basis of the observed phenotype. An evaluation of the (p)ppGpp-dependent transcriptional regulation under total nutrient starvation revealed a complex interplay with the involvement of ribosome-mediated transcriptional attenuation, the stress-responsive sigma factors σ^B and σ^H and transcription factors such as McbR, the master regulator of sulfur metabolism. In addition to the differential regulation of genes connected with various cell functions, the transcriptome analysis revealed conserved motifs within the promoter regions of (p)ppGpp-dependently and independently regulated genes. In particular, the representatives of translation-associated genes are both (p)ppGpp-dependent transcriptionally downregulated and show a highly conserved and so far unknown TTTTG motif in the -35 region, which is also present in other actinobacterial genera.

Overlap of Promoter Recognition Specificity of Stress Response Sigma Factors SigD and SigH in Corynebacterium glutamicum ATCC 13032

Corynebacterium glutamicum ATCC 13032 harbors five sigma subunits of RNA polymerase belonging to Group IV, also called extracytoplasmic function (ECF) σ factors. These factors σ^C, σ^D, σ^E, σ^H, and σ^M are mostly involved in stress responses. The role of σ^D consists in the control of cell wall integrity. The σ^D regulon is involved in the synthesis of components of the mycomembrane which is part of the cell wall in C. glutamicum. RNA sequencing of the transcriptome from a strain overexpressing the sigD gene provided 29 potential σ^D-controlled genes and enabled us to precisely localize their transcriptional start sites. Analysis of the respective promoters by both in vitro transcription and the in vivo two-plasmid assay confirmed that transcription of 11 of the tested genes is directly σ^D-dependent. The key sequence elements of all these promoters were found to be identical or closely similar to the motifs -35 GTAAC^A/_G and -10 GAT. Surprisingly, nearly all of these σ^D-dependent promoters were also active to a much lower extent with σ^Hin vivo and one (Pcg0607) also in vitro, although the known highly conserved consensus sequence of the σ^H-dependent promoters is different (-35 GGAA^T/_C and -10 GTT). In addition to the activity of σ^H at the σ^D-controlled promoters, we discovered separated or overlapping σ^A- or σ^B-regulated or σ^H-regulated promoters within the upstream region of 8 genes of the σ^D-regulon. We found that phenol in the cultivation medium acts as a stress factor inducing expression of some σ^D-dependent genes. Computer modeling revealed that σ^H binds to the promoter DNA in a similar manner as σ^D to the analogous promoter elements. The homology models together with mutational analysis showed that the key amino acids, Ala 60 in σ^D and Lys 53 in σ^H, bind to the second nucleotide within the respective -10 promoter elements (GAT and GTT, respectively). The presented data obtained by integrating in vivo, in vitro and in silico approaches demonstrate that most of the σ^D-controlled genes also belong to the σ^H-regulon and are also transcribed from the overlapping or closely located housekeeping (σ^A-regulated) and/or general stress (σ^B-regulated) promoters.

Leaderless mRNAs in the Spotlight: Ancient but Not Outdated!

Previously, leaderless mRNAs (lmRNAs) were perceived to make up only a minor fraction of the transcriptome in bacteria. However, advancements in RNA sequencing technology are uncovering vast numbers of lmRNAs, particularly in archaea, Actinobacteria , and extremophiles and thus underline their significance in cellular physiology and regulation. Due to the absence of conventional ribosome binding signals, lmRNA translation initiation is distinct from canonical mRNAs and can therefore be differentially regulated. The ribosome’s inherent ability to bind a 5′-terminal AUG can stabilize and protect the lmRNA from degradation or allow ribosomal loading for downstream initiation events. As a result, lmRNAs remain translationally competent during a variety of physiological conditions, allowing them to contribute to multiple regulatory mechanisms. Furthermore, the abundance of lmRNAs can increase during adverse conditions through the upregulation of lmRNA transcription from alternative promoters or by the generation of lmRNAs from canonical mRNAs cleaved by an endonucleolytic toxin. In these ways, lmRNA translation can continue during stress and contribute to regulation, illustrating their importance in the cell. Due to their presence in all domains of life and their ability to be translated by heterologous hosts, lmRNAs appear further to represent ancestral transcripts that might allow us to study the evolution of the ribosome and the translational process.

Iron and Zinc Regulate Expression of a Putative ABC Metal Transporter in Corynebacterium diphtheriae

Corynebacterium diphtheriae, a Gram-positive, aerobic bacterium, is the causative agent of diphtheria and cutaneous infections. While mechanisms required for heme iron acquisition are well known in C. diphtheriae, systems involved in the acquisition of other metals such as zinc and manganese remain poorly characterized. In this study, we identified a genetic region that encodes an ABC-type transporter (iutBCD) and that is flanked by two genes (iutA and iutE) encoding putative substrate binding proteins of the cluster 9 family, a related group of transporters associated primarily with the import of Mn and Zn. We showed that IutA and IutE are both membrane proteins with comparable Mn and Zn binding abilities. We demonstrated that the iutABCD genes are cotranscribed and repressed in response to iron by the iron-responsive repressor DtxR. Transcription of iutE was positively regulated in response to iron availability in a DtxR-dependent manner and was repressed in response to Zn by the Zn-dependent repressor Zur. Electrophoretic mobility shift assays showed that DtxR does not bind to the iutE upstream region, which indicates that DtxR regulation of iutE is indirect and that other regulatory factors controlled by DtxR are likely responsible for the iron-responsive regulation. Analysis of the iutE promoter region identified a 50-bp sequence at the 3' end of the iutD gene that is required for the DtxR-dependent and iron-responsive activation of the iutE gene. These findings indicate that transcription of iutE is controlled by a complex mechanism that involves multiple regulatory factors whose activity is impacted by both Zn and Fe.IMPORTANCE Vaccination against diphtheria prevents toxin-related symptoms but does not inhibit bacterial colonization of the human host by the bacterium. Thus, Corynebacterium diphtheriae remains an important human pathogen that poses a significant health risk to unvaccinated individuals. The ability to acquire iron, zinc, and manganese is critical to the pathogenesis of many disease-causing organisms. Here, we describe a gene cluster in C. diphtheriae that encodes a metal importer that is homologous to broadly distributed metal transport systems, some with important roles in virulence in other bacterial pathogens. Two metal binding components of the gene cluster encode surface exposed proteins, and studies of such proteins may guide the development of second-generation vaccines for C. diphtheriae.

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.

Native promoters of Corynebacterium glutamicum and its application in L-lysine production

OBJECTIVE

To identify useful native promoters of Corynebacterium glutamicum for fine-tuning of gene expression in metabolic engineering.

RESULTS

Sixteen native promoters of C. glutamicum were characterized. These promoters covered a strength range of 31-fold with small increments and exhibited relatively stable activity during the whole growth phase using β-galactosidase as the reporter. The mRNA level and enzymatic activity of the lacZ reporter gene exhibited high correlation (R ² = 0.96) under the control of these promoters. Sequence analysis found that strong promoters had high similarity of the -10 hexamer to the consensus sequence and preference of the AT-rich UP element upstream the -35 region. To test the utility of the promoter library, the characterized native promoters were applied to modulate the sucCD-encoded succinyl-CoA synthetase expression for L-lysine overproduction.

CONCLUSIONS

The native promoters with various strengths realize the efficient and precise regulation of gene expression in metabolic engineering of C. glutamicum.

Insight of Genus Corynebacterium: Ascertaining the Role of Pathogenic and Non-pathogenic Species

This review gathers recent information about genomic and transcriptomic studies in the Corynebacterium genus, exploring, for example, prediction of pathogenicity islands and stress response in different pathogenic and non-pathogenic species. In addition, is described several phylogeny studies to Corynebacterium, exploring since the identification of species until biological speciation in one species belonging to the genus Corynebacterium. Important concepts associated with virulence highlighting the role of Pld protein and Tox gene. The adhesion, characteristic of virulence factor, was described using the sortase mechanism that is associated to anchorage to the cell wall. In addition, survival inside the host cell and some diseases, were too addressed for pathogenic corynebacteria, while important biochemical pathways and biotechnological applications retain the focus of this review for non-pathogenic corynebacteria. Concluding, this review broadly explores characteristics in genus Corynebacterium showing to have strong relevance inside the medical, veterinary, and biotechnology field.

Construction of an expression vector that uses the aph promoter for protein expression in Corynebacterium glutamicum

Corynebacterium glutamicum is an attractive host for the production of heterologous proteins despite its traditional use in fermentative production of amino acids. To enhance the expression levels of target genes, the development of useful promoters is required in the construction of expression systems. Here, we developed a new promoter, the aph promoter from aminoglycoside-3'-phosphotransferase gene, and used it to construct monocistronic and bicistronic expression systems that host different ribosome binding site (RBS) sequences. First, the expression level of the reporter protein, enhanced green fluorescent protein (EGFP), varied with changes in the RBS sequences in the constructed vectors. The results showed that the fluorescence intensities of the bicistronic group were higher than those of the monocistronic group and that RM3E showed the highest fluorescence intensity, which was 42-fold higher than the lowest (RA2E') among these groups. Next, taking advantage of the optimized aph promoter, we successfully employed this aph promoter for α-amylase and V_HH (camelid antibody fragment) expression. The secretion of α-amylase improved 1.5-fold after promoter mutation. This promoter will be useful for heterologous protein production in C. glutamicum cells.