Abstract
Heme is an essential cofactor for almost all living cells by acting as prosthetic group for various proteins or serving as alternative iron source. However, elevated levels are highly toxic for cells. Several corynebacterial species employ two paralogous, heme-responsive two-component systems (TCS), ChrSA and HrrSA, to cope with heme stress and to maintain intracellular heme homeostasis. Significant cross-talk at the level of phosphorylation between these systems was previously demonstrated. In this study, we have performed a laboratory evolution experiment to adapt Corynebacterium glutamicum to increasing heme levels. Isolated strains showed a highly increased tolerance to heme growing at concentrations of up to 100 μM. The strain featuring the highest heme tolerance harbored a frameshift mutation in the catalytical and ATPase-domain (CA-domain) of the chrS gene, converting it into a catalytically-inactive pseudokinase (ChrS_CA-fs). Reintroduction of the respective mutation in the parental C. glutamicum strain confirmed high heme tolerance and showed a drastic upregulation of hrtBA encoding a heme export system, conserved in Firmicutes and Actinobacteria. The strain encoding the ChrS pseudokinase variant showed significantly higher heme tolerance than a strain lacking chrS. Mutational analysis revealed that induction of hrtBA in the evolved strain is solely mediated via the cross-phosphorylation of the response regulator (RR) ChrA by the kinase HrrS and BACTH assays revealed the formation of heterodimers between HrrS and ChrS. Overall, our results emphasize an important role of the ChrS pseudokinase in high heme tolerance of the evolved C. glutamicum and demonstrate the promiscuity in heme-dependent signaling of the paralogous two-component systems facilitating fast adaptation to changing environmental conditions.
Introduction
Heme constitutes 95% of functional iron in the human body and is a key molecule for almost all living cells (; ) acting as cofactor for many important proteins, including cytochromes, hydroxylases, catalases, peroxidases (; ), and serving as alternative iron source (). Nevertheless, elevated levels of this iron-bound protoporphyrin are highly cytotoxic. While this toxicity partially originates from the redox-active iron, causing the formation of reactive oxygen species (), a not yet unraveled porphyrin-related toxicity is furthermore suggested (). Consequently, organisms have evolved sophisticated mechanisms ensuring heme homeostasis (; ). Among those, systems enhancing heme tolerance play an important role in both pathogenic and non-pathogenic prokaryotes. Known strategies include mechanisms of (i) heme sequestration (e.g., HemS of Yersinia enteroliticia), (ii) heme degradation (e.g., IsdG of Bacillus anthracis), and (iii) heme export by the HrtBA system (; ; ). The heme-dedicated ATP-binding cassette efflux pump HrtBA is a highly conserved system and predominantly found in Firmicutes and Actinobacteria (). Recent structural studies shed light on the mechanism HrtBA employs to sequester and extrude heme from the cytoplasmic membrane ().
In Gram-positive bacteria, two-component systems (TCS) play a predominant role in the regulation of heme homeostasis (; ; ; ; ). The prototypical TCS consists of a membrane-bound histidine kinase (HK), which undergoes autophosphorlyation at a conserved histidine residue upon stimulus perception. Subsequently, the phosphoryl group is transferred to a conserved aspartate residue of a cytoplasmic response regulator (RR) resulting in an appropriate output, e.g., altering gene expression (; ; ). HKs may be composed of multiple domains with a significant architectural diversity, but typically consist of an N-terminal transmembrane domain and a C-terminal transmitter domain. The transmitter domain can be split up in the dimerization and histidine phosphotransfer (DHp) domain and the catalytical and ATPase (CA) domain (; ). The CA-domain comprises four sequence motifs, including N, G1, F, and G2 boxes, which bind ATP in a pocket using an ATP lid and are consequently necessary for the autophosphorylation reaction (; ). The DHp-domain possesses the H box motif harboring the conserved histidine residue which is phosphorylated upon stimulus perception (), as well as the X box, which is required for dimerization. The DHp- and CA-domain are connected via a flexible linker, which probably also supports keeping the RR in place during the phosphotransfer reaction (). Many HKs are bifunctional possessing also a phosphatase motif and subsequently acting both as kinase and phosphatase for the RR (; ; ). Furthermore, also catalytically inactive variants of kinases have been identified, referred to as pseudokinases, that can act as important signaling modulators by various mechanisms (; ; ; ).
Gene duplication events facilitate the evolution of TCS signaling enabling the integration of new input signals and diversification of the gene regulatory network. Members of the Corynebacteriaceae family, including the Gram-positive soil bacterium Corynebacterium glutamicum, represent an interesting example of a recent gene duplication event, encoding two paralogous two-component systems that both respond to the multifaceted molecule heme (Figure 1). After sensing heme availability via intramembrane interaction (; ), the TCS ChrSA acts as an activator of the hrtBA operon encoding the heme export system (). In contrast, the paralogous TCS HrrSA is a global regulator of heme homeostasis controlling more than 200 genomic targets including inter alia genes involved in heme biosynthesis, respiration as well as hmuO, encoding heme oxygenase (). Strikingly, a high level of cross-phosphorylation between the systems was observed (), while phosphatase activity of these HKs remains specific to their cognate RR ().
FIGURE 1
In this study, we addressed the question how this underlying signaling cascade consisting of two paralogous TCSs facilitates fast adaptation to high heme levels, such as encountered by pathogenic species in the mammalian host, using the non-pathogenic C. glutamicum as a model. Understanding the mechanisms underlying microbial heme tolerance are not only important for the control of bacterial infections but also of biotechnological relevance for the engineering of a microbial production host demanding high product tolerance (
Materials and methods
Bacterial strains and growth conditions
Bacterial strains used in this study are listed in Supplementary Table 1. For standard cultivation, C. glutamicum cells ATCC 13032 (wild type) and derivatives were streaked on agar plates (17 g/l) containing brain heart infusion (BHI) (Difco, BD, Heidelberg, Germany) (37 g/l) and inoculated at 30°C overnight. One single colony was picked and incubated for approximately 8 h at 30°C in 5 ml BHI in reaction tubes (for cultivation in shake flasks) or in 1 ml BHI in deep-well plates (VWR International, PA, United States) (for microtiter cultivation). This first pre-culture was used to inoculate the second pre-culture 1:10 in CGXII minimal medium (
For the ALE experiment, the main culture was grown in deep-well plates for 1–3 days and then freshly transferred at an OD600 of 1 for the next batch. After the 13th inoculation, glycerol stocks of each population were frozen at −80°C. This allowed a restreaking of each potentially heterogeneous population on BHI-agar plates and picking of single evolved clones. Online monitoring of bacterial growth was performed using the BioLector® microtiter cultivation system of Beckman Coulter GmbH (Baesweiler) (
Escherichia coli strains including DH5α and BTH101 were cultivated in Lysogeny Broth (10 g/l tryptone, 5 g/l yeast extract, 10 g/l NaCl) medium at 37°C in a rotary shaker and if needed for selection, 50 μg/ml kanamycin or 100 μg/ml ampicillin was added to the medium.
Recombinant DNA work
Standard molecular methods were performed according to
Plasmids were constructed by amplifying DNA fragments using the respective oligonucleotides (Supplementary Tables 2, 3) and enzymatically ligated into a pre-cut vector backbone using Gibson assembly (
For the deletion of genes in the genome of C. glutamicum, the suicide vector pK19-mobsacB was used (
Whole genome sequencing
Whole genome resequencing of C. glutamicum strains isolated during the ALE experiment was performed using next generation sequencing (NGS). Genomic DNA was prepared using the NucleoSpin microbial DNA kit (Macherey-Nagel, Düren, Germany) according to manufacturer’s instructions. Concentrations of the purified genomic DNA were measured using Qubit 2.0 fluorometer (Invitrogen, Carlsbad, CA, United States) according to manufacturer’s instructions. The purified genomic DNA was used for the preparation for genome sequencing using NEBNext Ultra II DNA Library Kit for Illumina (New England BioLabs, Frankfurt am Main) and MiSeq Reagent Kit v2 (300-cycles) (Illumina, San Diego, CA, United States), according to manufacturer’s instructions. A MiSeq system (Illumina, San Diego, CA, United States) was used for sequencing. Data analysis and base calling were accomplished with the Illumina instrument software. FASTQ output files were analyzed for single nucleotide polymorphisms using PathoSystems Resource Integration Center (PATRIC) 3.6.12 (
Gradient plates
For heme gradient plates, the different mutant strains were cultivated in triplicates as described above using deep-well plates for the first and second pre-culture. Subsequently, cultures were harvested and resuspended to an OD600 of 1 in 0.9% NaCl. For each spot, 2 μl of the respective samples were spotted on the gradient plates. The gradient plates were always prepared freshly. Therefore, 30 ml of iron-free CGXII with 2% (w/v) glucose, and 17 g/l Bacto Agar (Difco, BD, Heidelberg, Germany) was poured into a squared agar plate, which was in inclined position. Then, after drying of the first layer, the incline was removed and further 30 ml of iron-free CGXII containing 2% glucose, and 15 μM heme were poured in the plates so that a heme gradient results.
Bacterial two-hybrid assays
Bacterial two-hybrid plate assays for the qualitative assessment of protein-protein interactions
Bacterial two-hybrid assays were performed based on the BACTH kit according to manufacturer’s instructions (Euromedex, Souffelweyersheim, France). This method is based on the two fragments T25 and T18 of the catalytical domain of the adenylate cyclase from Bordetella pertussis, which is only active when these two fragments are physically in close contact. Therefore, T25 and T18 were each fused once to ChrS, ChrS-Ala245fs, and HrrS. If the HKs interact with each other, this allows a functional complementation of T25 and T18, leading to cAMP synthesis, which binds to the catabolic activator protein (CAP). cAMP/CAP complexes are pleiotropic regulators of gene transcription in E. coli and therefore turn on the expression of e.g., the lac operon.
Therefore, E. coli BTH101, which lack adenylate cyclase activity, were transformed with two plasmids of heterologous proteins fused once to T25 and once to T18. This approach was directly diluted and spotted as 100, 10–1, and 10–2 dilutions on LB plates containing 40 μg/ml X-Gal, 50 μg/ml kanamycin, 100 μg/ml ampicillin, and 0.5 mM IPTG and incubated approximately 24 h at 30°C. Bacteria producing interacting proteins will form blue colonies. Otherwise, the colonies remain white. This allowed to check also for heterogeneity in expression. Additionally, the approach was also plated on LB plates only with antibiotics, to allow picking for further biological replicates for the β-galactosidase assay.
Bacterial two-hybrid β-galactosidase measurements
Transformed E. coli strains (compare 1.5.1) were re-cultivated in deep-well plates as triplicates in LB media with 50 μg/ml kanamycin, 100 μg/ml ampicillin and 0.5 mM IPTG. OD600 of the overnight cultures was measured in a Tecan Reader (Thermo Fisher Scientific, Massachusetts, United States). The β-galactosidase assay was adapted according to a previous study for 96-well plates (
With A420 being the absorbance of the yellow o-nitrophenol, A550 the scatter from the cell debris, 1.75 is the factor which needs to be multiplied with A550 to approximate scatter observed at 420 nm, t is time in minutes, v is the volume of the culture employed in the plate and A600 for the cell density. The value of A600 was calibrated to proper OD600.
DNA microarrays
For the analysis of the transcriptome, C. glutamicum wild type and C. glutamicum ChrS-Ala245fs were cultivated in triplicates as described above in 50 ml CGXII medium containing 2% glucose, no FeSO4 and 4 μM heme in shake flasks. Cells were harvested after 6 h, when the wild type reached an OD600 around 2.5 and the mutated strain around 5. The cell suspension was centrifuged at 4,250 × g, 10 min, 4°C in falcons filled with ice. The resulting pellets were frozen in liquid nitrogen and stored at −80°C. The following RNA preparation, cDNA synthesis, microchip hybridization, scanning, and overall evaluation was performed as described in previous studies (
Results
Adaptive laboratory evolution of Corynebacterium glutamicum toward tolerance of high heme levels depends on the heme exporter HrtBA
Previous studies reported a crucial role of the two-component system ChrSA for heme tolerance of C. glutamicum (
FIGURE 2

Adaptive laboratory evolution of C. glutamicum to high heme levels. The C. glutamicum wild type (WT) (A) as well as the ΔhrtBA strain (B) were inoculated at a starting-OD600 of 1 in CGXII medium containing 2% glucose and either 36 μM iron (blue) or increasing concentrations of heme (2.5–15 μM, shades of orange). Data represent the average of three biological replicates including standard deviations depicted as error bars. (C) Schematic representation of the adaptive laboratory evolution (ALE) experiment. The heme concentration was increased from 10 to 60 μM heme for the WT in overall 13 repetitive batch cultures (depicted in the bar graph). For the ΔhrtBA strain, concentrations > 15 μM remained toxic throughout the experiment and therefore were not further increased. (D) Growth of each three single clones derived from four different evolved C. glutamicum wild type populations on 100 μM heme. (E) Growth of each three single clones of the three evolved ΔhrtBA populations on 10 μM heme.
To elucidate mechanisms promoting high heme tolerance, we performed an ALE experiment applying increasing concentrations of heme to C. glutamicum wild type and the ΔhrtBA strain (Figure 2C). The ALE was accomplished in CGXII minimal medium, with 13 repetitive batch cultivations and started from each four independent single colonies, yielding four evolving populations. Batch cultures were started on 10 μM heme and were stepwise increased to finally reach 60 μM heme in the case of C. glutamicum wild type. From C. glutamicum wild type populations, four single clones were isolated from agar plates after the 13th batch cultivation and then further characterized in liquid culture (Figure 2D). Further analysis on earlier inoculation steps was also performed (Supplementary Figure 1). Remarkably, all isolates from the 13th inoculation were able to grow in the re-cultivation on medium containing up to 100 μM heme where growth of the parental strain was completely inhibited. By contrast, heme levels above 15 μM heme remained toxic to the ΔhrtBA strain and could not be increased throughout the ALE experiment without killing the cells. Only three of the four starting populations survived 13 inoculations and for these, no significant adaptation was observed throughout the ALE experiment (Figure 2E). These results already underlined that HrtBA represents a key factor for the adaptation of C. glutamicum to high heme levels.
Mutations in the catalytic- and ATPase domain of ChrS lead to significantly improved growth on heme
Whole genome sequencing of the four isolated clones from the ALE experiment revealed that all of them possess a mutation in the gene encoding the HK ChrS (Table 1). Remarkably, in all four cases, the catalytic- and ATPase (CA) domain was affected. For clone 2, a stop codon was inserted directly at the beginning of the CA-domain, while clone 3 possesses an amino acid exchange inside and close to the start of this domain. Strikingly, clone 1 and 4 (evolved in two independent cell lines) showed the exact same frameshift mutation for alanine at position 245, immediately after the dimerization and histidine phosphotransfer (DHp) domain of ChrS. This mutation also affects the linker between the CA- and DHp-domain. Since the strains carrying this mutation showed the highest heme tolerance, evolved clone number 1, in the following referred to as 1.fs, was further analyzed within this study.
TABLE 1
| Clone | Type | Variance nucleotides | Variance amino acids | Locus |
| 1.fs, 4.fs | Frameshift | 733delG | Ala245fs | cg2201 Sensor histidine kinase ChrS |
| 2.* | Non-sense | 862C > T | Gln288* | |
| 3. < | Missense | 916A > C | Thr306Pro |
Key mutations identified in C. glutamicum strains featuring increased heme tolerance.
The * is the symbol for the insertion of a stop-codon. Similar, the > indicates an amino acid exchange.
Reintegration of this mutation into the wild type parental strain confirmed that this frameshift mutation in chrS led to highly improved growth on heme (Figure 3A and Supplementary Figure 2). Protein structure prediction via AlphaFold2 shows the truncated CA-domain of ChrS-Ala245fs with additional 60 amino acids caused by the frameshift resulting in a presumably catalytically-inactive pseudokinase variant (Figure 3B;
FIGURE 3

The ChrS-Ala245fs pseudokinase promotes heme tolerance. Data represent the average of three biological replicates including standard deviations depicted as error bars. Cells were inoculated at a starting-OD600 of 1 in CGXII medium containing 2% glucose and the indicated amount of heme or iron. (A) Growth of C. glutamicum carrying the reintegrated ChrS-Ala245fs allele (gray) compared to the evolved clone 1.fs (orange) and the wild type (WT) (blue). (B) Predicted protein structures of ChrS (blue) and the truncated ChrS-Ala245fs variant (orange, with 60 additional amino acids caused by the frameshift, shown in red). Prediction was performed using AlphaFold2 (
The 1.fs strain showed normal growth under standard conditions. However, it outcompeted the wild type in the presence of heme at all concentrations tested, but not under iron starved conditions (Figures 3C,D and Supplementary Figure 2). A heme gradient plate experiment further confirmed the improved heme tolerance of the evolved clone 1.fs, and showed that this is also the case for the reintegrated pseudokinase variant ChrS-Ala245fs (Figure 3E).
The ChrS pseudokinase promotes increased hrtBA expression and is crucial for the improved growth on heme
The main known target of ChrA is the hrtBA operon encoding a heme export system. To unravel the impact of the frame-shift mutation in chrS on its activation kinetics, hrtBA reporter assays were performed, using the reporter plasmid pJC1-PhrtBA-eyfp (
FIGURE 4

Mutational analysis provides insights in the hrtBA activating cascade of the evolved C. glutamicum 1.fs clone. Data represent the average of three biological replicates including standard deviations depicted as error bars. Cells were inoculated at a starting-OD600 of 1 in CGXII medium containing 2% glucose and the indicated amount of heme or iron. (A) Reporter assays visualizing hrtBA-expression using the plasmid pJC1-PhrtB-eyfp for transformation of the wild type (WT) (dark blue), the evolved clone 1.fs (dark orange), ΔchrS (light orange), ΔchrSΔhrrS (mid-blue), and chrS-Q191A (light blue). (B) Growth of the evolved clone (shades of orange) and the WT (shades of blue) natively or with each a deletion of hrrS. (C) Reporter assays for hrtB-expression using pJC1-PhrtB-eyfp upon deletion of hrrS in 1.fs (shades of orange) and the WT (shades of blue). (D) Growth of the wild type strain (shades of blue) and the evolved 1.fs strain (shades of orange) natively or possessing an amino acid exchange of the autophosphorylation histidine (H186) of ChrS. (E) Impact of ChrS truncation on growth; ChrS variants lacking only the CA-domain (light orange) or CA-domain and DHp-domain (dark orange) compared to the WT and 1.fs strain (shades of blue). TM, transmembrane domain; DHp, dimerization- and histidine-phosphotransfer-domain; CA, catalytic- and ATPase-domain.
In the following, we further analyzed the growth and hrtBA expression of the evolved clone in comparison to different mutant strains, including the deletion mutants ΔchrS and ΔchrSΔhrrS as well as the phosphatase mutant chrS-Q191A (Figures 3E, 4A;
Considering that the frame-shift mutation in chrS likely abolished the catalytic activity of the ChrS kinase, we wondered whether activation of ChrA is solely dependent on HrrS. Cross-talk between the kinases was previously described (
To investigate if the truncated version of ChrS still plays a role in the phosphotransfer to ChrA, we exchanged the conserved histidine residue at position 186 to an alanine (Figure 4D). This amino acid exchange did, however, not significantly influence the heme tolerance of strain 1.fs. Therefore, it can be assumed that the truncated version of ChrS does not participate in the phosphotransfer via autophosphorylation at the histidine 186.
Based on these results, we postulated that a lack of the catalytic activity of ChrS is beneficial for C. glutamicum heme tolerance. In line with this hypothesis, an in-frame deletion of the CA-domain of chrS led to improved growth on heme compared to the wild type. However, this strain did not reach comparably high heme tolerance like the 1.fs strain (Figure 4E). Upon additional deletion of the DHp-domain, this growth advantage was abolished suggesting the necessity of homo- or heterodimerization and/or interaction with ChrA.
ChrS and HrrS form heterodimers in vivo
In general, HKs act as homodimers. However, based on the fact that the autophosphorylation histidine residue of ChrS in the evolved mutant 1.fs was not relevant for its growth benefit on heme (Figure 4D), we aimed to investigate the homo- and heterodimerization properties of the C. glutamicum HKs ChrS and HrrS and the pseudokinase ChrS (ChrS_CA-fs).
To assess these protein-protein interactions between the respective monomers ChrS, HrrS and the truncated HK variant ChrS_CA-fs, we performed bacterial two-hybrid (BACTH) assays (Euromedex, Souffelweyersheim, France). The plate assays in Figure 5A as well as the quantitative β-galactosidase assay in Figure 5B show that homodimerization for both the native and the evolved HKs was observed when these fusion proteins were produced as C-terminal fusions in E. coli. The assays also revealed heterodimerization of the native HKs, while there was no significant evidence for heterodimerization of the truncated ChrS_CA-fs version. Similar results were observed when the proteins were produced as N-terminal fusions (Supplementary Figure 7).
FIGURE 5

Bacterial two-hybrid assays of interactions between ChrSA and HrrSA. (A) BACTH interactions between the histidine kinases ChrS, HrrS and the evolved ChrS variant (here ChrS_CA-fs) were analyzed as C-terminal fusions; results for N-terminal variants are shown in Supplementary Figure 7. Blueish color of the colonies indicate interaction, while white colonies indicate no interaction (Euromedex, Souffelweyersheim, France). First histidine kinase represents the T25-fusion, the second the T18-fusion. + = pKTN25-zip with pUT18-zip (leucine zipper, positive control), – = pKTN25 with pUT18 (negative control). (B) Quantitative analysis using a β-galactosidase assay. Triplicates were cultivated and treated according to
In a next set of experiments, we investigated the interaction between the sensor kinases and the RRs. Here, BACTH assays confirmed the interaction of ChrS with ChrA and HrrS with HrrA, as well as the cross-talk between ChrS and HrrA. Interaction was not observed for ChrS_CA-fs and neither ChrA nor HrrA within the β-galactosidase assay, although a slight signal appeared to be visible on plates (Figure 5). However, the assays did also not reveal the already reported cross-talk between HrrS and ChrA (
Heme-binding proteins contribute to heme tolerance
Within this study, we showed the crucial role of the HrtBA export system for C. glutamicum heme tolerance. To identify further potentially relevant factors, we performed a comparative transcriptome analysis of the ChrS-Ala245fs strain and C. glutamicum wild type (Table 2). As expected, the hrtBA operon showed significantly increased mRNA levels in the evolved clone (∼150-fold increase). Besides hrtBA, several other heme-related targets also showed increased expression levels, including the TCS chrSA itself, hmuO encoding heme oxygenase and the heme transport system hmuU. Quantitative PCR confirmed unaltered expression levels of hrrS (Supplementary Figure 4). Strikingly, all genes encoding known heme-binding proteins were significantly upregulated. Furthermore, many targets of the DtxR regulon were upregulated (
TABLE 2
| Category and cg gene number | Gene designation and description of product | mRNA ratioa | P-value |
| Heme-related genes | |||
| cg0468 | hmuU, hemin transport system, permease protein | 16.90 | 0.05 |
| cg2200 | chrA, two-component system, response regulator | 8.85 | 0.01 |
| cg2201 | chrS, two-component system, signal transduction histidine kinase | 43.60 | 0.04 |
| cg2202 | hrtB, ABC-type transport system, permease component | 170.13 | 0.00 |
| cg2204 | hrtA, ABC-type transport system, ATPase component | 151.54 | 0.00 |
| cg2445 | hmuO, heme oxygenase | 8.30 | 0.01 |
| Heme-binding proteins | |||
| cg0466 | htaA, secreted heme-transport associated protein | 5.80 | 0.01 |
| cg0467 | hmuT, hemin-binding periplasmic protein precursor | 21.23 | 0.06 |
| cg0470 | htaB, secreted heme transport-associated protein | 66.64 | 0.02 |
| cg0471 | htaC, secreted heme transport-associated protein | 16.87 | 0.04 |
| cg3156 | htaD, secreted heme transport-associated protein | 18.64 | 0.05 |
| DtxR regulon | |||
| cg0160 | Hypothetical protein cg0160 | 2.98 | 0.00 |
| cg1120 | ripA, transcriptional regulator of iron proteins, AraC family | 5.76 | 0.10 |
| cg1419 | Putative Na+-dependent transporter | 4.85 | 0.01 |
| cg1476 | thiC, thiamine biosynthesis protein ThiC | 2.48 | 0.04 |
| cg1695 | Putative plasmid maintenance system antidote protein | 0.34 | 0.05 |
| cg1930 | Putative secreted hydrolase | 5.68 | 0.01 |
| cg1930 | Putative secreted hydrolase | 5.68 | 0.01 |
| cg1931 | Putative secreted protein | 9.98 | 0.05 |
| cg1931 | Putative secreted protein | 9.98 | 0.05 |
| cg2311 | SAM-dependent methyltransferase | 3.47 | 0.00 |
| cg2444 | Hypothetical protein cg2444 | 4.67 | 0.01 |
| cg2782 | ftn, ferritin-like protein | 0.32 | 0.06 |
| cg2796 | MMGE/PRPD family protein | 11.01 | 0.00 |
| cg2962 | Uncharacterized enzyme involved in biosynthesis of extracellular polysaccharides | 6.69 | 0.02 |
| RipA regulon | |||
| cg0310 | katA, catalase | 0.14 | 0.00 |
| cg0445 | sdhC, succinate dehydrogenase | 0.36 | 0.02 |
| cg0446 | sdhA, succinate dehydrogenase | 0.40 | 0.00 |
| cg0447 | sdhB, succinate dehydrogenase | 0.45 | 0.00 |
| cg1343 | narH, probable respiratory nitrate reductase oxidoreduct | 0.49 | 0.03 |
| cg1344 | narG, nitrate reductase 2, alpha subunit | 0.30 | 0.00 |
| cg1487 | leuC, isopropylmalate isomerase large subunit | 0.29 | 0.01 |
| cg1737 | acn, aconitate hydratase | 0.29 | 0.01 |
| cg2636 | catA1, catechol 1,2-dioxygenase | 0.03 | 0.00 |
| cg3048 | pta, phosphate acetyltransferase | 0.24 | 0.00 |
Comparative transcriptome analysis of C. glutamicum wild type and C. glutamicum chrS-Ala245fs growing on 4 μM heme.
aExpression of selected genes given as the mRNA ratio of the evolved strain compared to the WT (>2-fold or < 0.5-fold, p-value < 0.05). Data represent the average of three biological replicates (for a complete list of up- and downregulated genes, see Supplementary Table 4).
To test whether the upregulation of heme-binding proteins could also contribute to heme tolerance mediated by heme sequestration, we further analyzed the impact of heme binding proteins by the construction of serial deletions (Figure 6). A strain lacking the heme binding proteins hmuT, htaA, htaB, htaC, and htaD showed wild typic growth at low (4 μM) and high (20 μM) heme levels. However, at moderate (10 μM) heme concentrations, the mutant showed a significant growth defect. These results suggested that heme sequestration via heme-binding proteins could promote heme tolerance at moderate levels.
FIGURE 6

Heme binding proteins contribute to heme tolerance. The C. glutamicum wild type (WT) (blue) and the heme binding deficient mutant ΔhmuTΔhtaAΔhtaBCΔhtaD strain (orange) were inoculated at a starting-OD600 of 1 in CGXII medium containing 2% glucose and either 4, 10, or 20 μM heme. Data represent the average of three biological replicates including standard deviations depicted as error bars.
Discussion
In this study, we pursued a laboratory evolution approach to adapt C. glutamicum ATCC 13032 to high heme levels. This ALE approach resulted in the isolation of strains harboring a frameshift mutation in the chrS HK gene yielding a catalytically inactive pseudokinase, which was shown to promote high heme tolerance of up to 100 μM. This effect could mainly be attributed to the strong upregulation of the heme exporter HrtBA and was strictly dependent on phosphotransfer via the non-cognate HK HrrS to the RR ChrA. Further mutational analysis confirmed that the conserved histidine residue of the ChrS pseudokinase (ChrS_CA-fs) was not involved in this phosphotransfer reaction.
Continuously high hrtBA expression levels observed in our evolved 1.fs strain are in agreement with a defect in ChrS phosphatase function of this strain. Remarkably, the evolved strains encoding the truncated pseudokinase variant ChrS showed significantly higher heme tolerance and higher hrtBA expression levels compared to a phosphatase deficient strain chrS-Q191A or a strain lacking chrS completely (ΔchrS) (
Our results indicate that the catalytically inactive version of ChrS promotes—directly or indirectly—the efficient phosphotransfer reaction from the paralogous HrrS to ChrA leading to the constitutive activation of hrtBA (due to the absence of ChrS phosphatase activity). Gene duplication is a powerful evolutionary driving force and the presence of paralogs has previously been shown to be beneficial for adaptations to new environmental conditions (
Pseudokinases are described as kinases lacking catalytic functions, but can contribute to signaling via functioning as allosteric modulators, dynamic scaffolds, or competitors of protein-protein interactions (
Although only a few reports on heterodimerization of HKs exist up to date (
Apart from the conserved heme export system HrtBA, transcriptome analysis gave further hints for additional players contributing to heme tolerance in C. glutamicum. Here, we observed a high upregulation of all genes encoding heme-binding proteins and a deletion mutant showed reduced tolerance to intermediate heme levels. In fact, heme-binding proteins could serve a detoxifying role via heme sequestration, as described throughout the literature for other organisms, including HbpC of Bartonella henselae (
Moreover, several transport systems are differently expressed in the evolved strains. Export of further toxic heme-related products, or even import of neutralizing compounds coping with H2O2, like e.g., described for the ribulose-5-phosphate 3-epimerase in Escherichia coli or the Mn(II) uptake system of Neisseria gonorrhoeae importing manganese (
The appearance of a ChrS pseudokinase was not yet described to occur naturally in corynebacterial strains (
Strains featuring an elevated heme tolerance are also highly interesting for the biotechnological production of heme, which is commercially produced for medical uses or the food sector for artificial meat products. Recent metabolic engineering efforts resulted in E. coli (
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
JF and AK conceived and designed the analysis and wrote the manuscript. AK performed the experiments and collected the data, contributed to data and analysis tools, and performed the analysis. Both authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Helmholtz Association (grant no. W2/W3-096).
Acknowledgments
We thank Ulrike Viets for whole genome sequencing, Tino Polen for help with the microarrays and upload to databases and Elvira Mukinovic for assisting with the qPCR experiments.
Conflict of interest
AK and JF were employed by the Forschungszentrum Jülich GmbH.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.997448/full#supplementary-material
References
1
AjiokaR. S.PhillipsJ. D.KushnerJ. P. (2006). Biosynthesis of heme in mammals.Biochim. Biophys. Acta1763723–736. 10.1016/j.bbamcr.2006.05.005
2
AndrewsS. C.RobinsonA. K.Rodríguez-QuiñonesF. (2003). Bacterial iron homeostasis.FEMS Microbiol. Rev.27215–237. 10.1016/S0168-6445(03)00055-X
3
AnzaldiL. L.SkaarE. P. (2010). Overcoming the heme paradox: Heme toxicity and tolerance in bacterial pathogens.Infect. Immun.784977–4989. 10.1128/IAI.00613-10
4
BaumgartM.LuderK.GroverS.GätgensC.BesraG. S.FrunzkeJ. (2013). IpsA, a novel LacI-type regulator, is required for inositol-derived lipid formation in Corynebacteria and Mycobacteria.BMC Biol.11:122. 10.1186/1741-7007-11-122
5
BibbL. A.KunkleC. A.SchmittM. P. (2007). The ChrA-ChrS and HrrA-HrrS signal transduction systems are required for activation of the hmuO promoter and repression of the hemA promoter in Corynebacterium diphtheriae.Infect. Immun.752421–2431. 10.1128/IAI.01821-06
6
BottM.BrockerM. (2012). Two-component signal transduction in Corynebacterium glutamicum and other corynebacteria: On the way towards stimuli and targets.Appl. Microbiol. Biotechnol.941131–1150. 10.1007/s00253-012-4060-x
7
BratlieM. S.JohansenJ.ShermanB. T.HuangD. W.LempickiR. A.DrabløsF. (2010). Gene duplications in prokaryotes can be associated with environmental adaptation.BMC Genomics11:588. 10.1186/1471-2164-11-588
8
BuelowD. R.RaivioT. L. (2010). Three (and more) component regulatory systems – auxiliary regulators of bacterial histidine kinases.Mol. Microbiol.75547–566. 10.1111/j.1365-2958.2009.06982.x
9
CapraE. J.LaubM. T. (2012). Evolution of two-component signal transduction systems.Annu. Rev. Microbiol.66325–347. 10.1146/annurev-micro-092611-150039
10
CasinoP.RubioV.MarinaA. (2009). Structural insight into partner specificity and phosphoryl transfer in two-component signal transduction.Cell139325–336. 10.1016/j.cell.2009.08.032
11
DavisJ. J.WattamA. R.AzizR. K.BrettinT.ButlerR.ButlerR. M.et al (2020). The PATRIC bioinformatics resource center: Expanding data and analysis capabilities.Nucleic Acids Res.48D606–D612. 10.1093/nar/gkz943
12
DuttaR.QinL.InouyeM. (1999). Histidine kinases: Diversity of domain organization.Mol. Microbiol.34633–640. 10.1046/j.1365-2958.1999.01646.x
13
EikmannsB. J.Thum-SchmitzN.EggelingL.LüdtkeK.-U.SahmH. (1994). Nucleotide sequence, expression and transcriptional analysis of the Corynebacterium glutamicum gltA gene encoding citrate synthase.Microbiology1401817–1828. 10.1099/13500872-140-8-1817
14
FrancisV. I.PorterS. L. (2019). Multikinase networks: Two-component signaling networks integrating multiple stimuli.Annu. Rev. Microbiol.73199–223. 10.1146/annurev-micro-020518-115846
15
FrancisV. I.WatersE. M.Finton-JamesS. E.GoriA.KadiogluA.BrownA. R.et al (2018). Multiple communication mechanisms between sensor kinases are crucial for virulence in Pseudomonas aeruginosa.Nat. Commun.9:2219. 10.1038/s41467-018-04640-8
16
GaoR.StockA. M. (2009). Biological insights from structures of two-component proteins.Annu. Rev. Microbiol.63133–154. 10.1146/annurev.micro.091208.073214
17
GeversD.VandepoeleK.SimillionC.Van de PeerY. (2004). Gene duplication and biased functional retention of paralogs in bacterial genomes.Trends Microbiol.12148–154. 10.1016/j.tim.2004.02.007
18
GibsonD. G.YoungL.ChuangR. Y.VenterJ. C.HutchisonC. A.IIISmithH. O. (2009). Enzymatic assembly of DNA molecules up to several hundred kilobases.Nat. Methods6343–345. 10.1038/nmeth.1318
19
GoodmanA. L.MerighiM.HyodoM.VentreI.FillouxA.LoryS. (2009). Direct interaction between sensor kinase proteins mediates acute and chronic disease phenotypes in a bacterial pathogen.Genes Dev.23249–259. 10.1101/gad.1739009
20
GriffithK. L.WolfR. E.Jr. (2002). Measuring beta-galactosidase activity in bacteria: Cell growth, permeabilization, and enzyme assays in 96-well arrays.Biochem. Biophys. Res. Commun.290397–402. 10.1006/bbrc.2001.6152
21
HentschelE.MackC.GatgensC.BottM.BrockerM.FrunzkeJ. (2014). Phosphatase activity of the histidine kinases ensures pathway specificity of the ChrSA and HrrSA two-component systems in Corynebacterium glutamicum.Mol. Microbiol.921326–1342. 10.1111/mmi.12633
22
HeyerA.GatgensC.HentschelE.KalinowskiJ.BottM.FrunzkeJ. (2012). The two-component system ChrSA is crucial for haem tolerance and interferes with HrrSA in haem-dependent gene regulation in Corynebacterium glutamicum.Microbiology1583020–3031. 10.1099/mic.0.062638-0
23
HorsburghM. J.WhartonS. J.KaravolosM.FosterS. J. (2002). Manganese: Elemental defence for a life with oxygen.Trends Microbiol.10496–501. 10.1016/S0966-842X(02)02462-9
24
IniestaA. A.HillsonN. J.ShapiroL. (2010). Cell pole-specific activation of a critical bacterial cell cycle kinase.Proc. Natl. Acad. Sci. U.S.A.1077012–7017. 10.1073/pnas.1001767107
25
ItoY.NakagawaS.KomagataA.Ikeda-SaitoM.ShiroY.NakamuraH. (2009). Heme-dependent autophosphorylation of a heme sensor kinase, ChrS, from Corynebacterium diphtheriae reconstituted in proteoliposomes.FEBS Lett.5832244–2248. 10.1016/j.febslet.2009.06.001
26
JingX.JawJ.RobinsonH. H.SchubotF. D. (2010). Crystal structure and oligomeric state of the RetS signaling kinase sensory domain.Proteins781631–1640. 10.1002/prot.22679
27
JumperJ.EvansR.PritzelA.GreenT.FigurnovM.RonnebergerO.et al (2021). Highly accurate protein structure prediction with AlphaFold.Nature596583–589. 10.1038/s41586-021-03819-2
28
KeilhauerC.EggelingL.SahmH. (1993). Isoleucine synthesis in Corynebacterium glutamicum: Molecular analysis of the ilvB-ilvN-ilvC operon.J. Bacteriol.1755595–5603. 10.1128/jb.175.17.5595-5603.1993
29
KensyF.ZangE.FaulhammerC.TanR.-K.BüchsJ. (2009). Validation of a high-throughput fermentation system based on online monitoring of biomass and fluorescence in continuously shaken microtiter plates.Microb. Cell Fact.8:31. 10.1186/1475-2859-8-31
30
KeppelM.DavoudiE.GätgensC.FrunzkeJ. (2018). Membrane topology and heme binding of the histidine kinases HrrS and ChrS in Corynebacterium glutamicum.Front. Microbiol.9:183. 10.3389/fmicb.2018.00183
31
KeppelM.HünnefeldM.FilipchykA.VietsU.DavoudiC. F.KrügerA.et al (2020). HrrSA orchestrates a systemic response to heme and determines prioritization of terminal cytochrome oxidase expression.Nucleic Acids Res.486547–6562. 10.1093/nar/gkaa415
32
KeppelM.PiepenbreierH.GatgensC.FritzG.FrunzkeJ. (2019). Toxic but tasty - temporal dynamics and network architecture of heme-responsive two-component signaling in Corynebacterium glutamicum.Mol. Microbiol.1111367–1381. 10.1111/mmi.14226
33
KimD. J.ForstS. (2001). Genomic analysis of the histidine kinase family in bacteria and archaea.Microbiology (Reading)1471197–1212. 10.1099/00221287-147-5-1197
34
KoY. J.KimM.YouS. K.ShinS. K.ChangJ.ChoiH. J.et al (2021). Animal-free heme production for artificial meat in Corynebacterium glutamicum via systems metabolic and membrane engineering.Metabolic Eng.66217–228. 10.1016/j.ymben.2021.04.013
35
KrügerA.KeppelM.SharmaV.FrunzkeJ. (2022). The diversity of heme sensor systems - heme-responsive transcriptional regulation mediated by transient heme protein interactions.FEMS Microbiol. Rev.46:fuac002. 10.1093/femsre/fuac002
36
KumarS.BandyopadhyayU. (2005). Free heme toxicity and its detoxification systems in human.Toxicol. Lett.157175–188. 10.1016/j.toxlet.2005.03.004
37
KungJ. E.JuraN. (2019). Prospects for pharmacological targeting of pseudokinases.Nat. Rev. Drug Discov.18501–526. 10.1038/s41573-019-0018-3
38
KwonA.ScottS.TaujaleR.YeungW.KochutK. J.EyersP. A.et al (2019). Tracing the origin and evolution of pseudokinases across the tree of life.Sci. Signal.12:eaav3810. 10.1126/scisignal.aav3810
39
KwonS. J.BoerA. L. D.PetriR.Schmidt-DannertC. (2003). High-Level production of porphyrins in metabolically engineered Escherichia coli: Systematic extension of a pathway assembled from overexpressed genes involved in heme biosynthesis.Appl. Environ. Microbiol.694875–4883. 10.1128/AEM.69.8.4875-4883.2003
40
LaubM. T.GoulianM. (2007). Specificity in two-component signal transduction pathways.Annu. Rev. Genet.41121–145. 10.1146/annurev.genet.41.042007.170548
41
LayerG.ReicheltJ.JahnD.HeinzD. W. (2010). Structure and function of enzymes in heme biosynthesis.Protein Sci.191137–1161. 10.1002/pro.405
42
LylesK. V.ThomasL. S.OuelletteC.CookL. C. C.EichenbaumZ. (2022). HupZ, a unique heme-binding protein, enhances group A Streptococcus fitness during mucosal colonization.Front. Cell. Infect. Microbiol.12:867963. 10.3389/fcimb.2022.867963
43
MaceP. D.MurphyJ. M. (2021). There’s more to death than life: Noncatalytic functions in kinase and pseudokinase signaling.J. Biol. Chem.296:100705. 10.1016/j.jbc.2021.100705
44
MascherT.HelmannJ. D.UndenG. (2006). Stimulus perception in bacterial signal-transducing histidine kinases.Microbiol. Mol. Biol. Rev.70910–938. 10.1128/MMBR.00020-06
45
NakamuraH.HisanoT.RahmanM. M.ToshaT.ShirouzuM.ShiroY. (2022). Structural basis for heme detoxification by an ATP-binding cassette-type efflux pump in gram-positive pathogenic bacteria.Proc. Natl. Acad. Sci. U.S.A.119:e2123385119. 10.1073/pnas.2123385119
46
NiebischA.BottM. (2001). Molecular analysis of the cytochrome bc1-aa3 branch of the Corynebacterium glutamicum respiratory chain containing an unusual diheme cytochrome c1.Arch. Microbiol.175282–294. 10.1007/s002030100262
47
PadmanabanG.VenkateswarV.RangarajanP. N. (1989). Haem as a multifunctional regulator.Trends Biochem. Sci.14492–496. 10.1016/0968-0004(89)90182-5
48
PeregoM.HochJ. A. (1996). Protein aspartate phosphatases control the output of two-component signal transduction systems.Trends Genet.1297–101. 10.1016/0168-9525(96)81420-X
49
PonkaP. (1999). Cell biology of heme.Am. J. Med. Sci.318241–256. 10.1016/S0002-9629(15)40628-7
50
RajuS.ShawA. S. (2015). What is the point of pseudokinases?Elife4:e07771. 10.7554/eLife.07771
51
ReisingerS. J.HuntworkS.ViollierP. H.RyanK. R. (2007). DivL performs critical cell cycle functions in Caulobacter crescentus independent of kinase activity.J. Bacteriol.1898308–8320. 10.1128/JB.00868-07
52
ReitererV.EyersP. A.FarhanH. (2014). Day of the dead: Pseudokinases and pseudophosphatases in physiology and disease.Trends Cell Biol.24489–505. 10.1016/j.tcb.2014.03.008
53
RodenJ. A.WellsD. H.ChomelB. B.KastenR. W.KoehlerJ. E. (2012). Hemin binding protein C is found in outer membrane vesicles and protects Bartonella henselae against toxic concentrations of hemin.Infect. Immun.80929–942. 10.1128/IAI.05769-11
54
SalvadoB.VilaprinyoE.KarathiaH.SorribasA.AlvesR. (2012). Two component systems: Physiological effect of a third component.PLoS One7:e31095. 10.1371/journal.pone.0031095
55
SambrookJ. F.RussellD. (2001). Molecular cloning: A laboratory manual (3-Volume Set), Vo. 1. New York, NY: Cold Spring Harbor Laboratory Press.
56
SchäferA.TauchA.JägerW.KalinowskiJ.ThierbachG.PühlerA. (1994). Small mobilizable multi-purpose cloning vectors derived from the Escherichia coli plasmids pK18 and pK19: Selection of defined deletions in the chromosome of Corynebacterium glutamicum.Gene14569–73. 10.1016/0378-1119(94)90324-7
57
SeibK. L.TsengH. J.McEwanA. G.ApicellaM. A.JenningsM. P. (2004). Defenses against oxidative stress in Neisseria gonorrhoeae and Neisseria meningitidis: Distinctive systems for different lifestyles.J. Infect. Dis.190136–147. 10.1086/421299
58
SkaarE. P.GasparA. H.SchneewindO. (2006). Bacillus anthracis IsdG, a heme-degrading monooxygenase.J. Bacteriol.1881071–1080. 10.1128/JB.188.3.1071-1080.2006
59
SobotaJ. M.ImlayJ. A. (2011). Iron enzyme ribulose-5-phosphate 3-epimerase in Escherichia coli is rapidly damaged by hydrogen peroxide but can be protected by manganese.Proc. Natl. Acad. Sci. U.S.A.1085402–5407. 10.1073/pnas.1100410108
60
StauffD. L.SkaarE. P. (2009b). The heme sensor system of Staphylococcus aureus.Contrib. Microbiol.16120–135. 10.1159/000219376
61
StauffD. L.SkaarE. P. (2009a). Bacillus anthracis HssRS signaling to HrtAB regulates heme resistance during infection.Mol. Microbiol.72763–778. 10.1111/j.1365-2958.2009.06684.x
62
StockA. M.RobinsonV. L.GoudreauP. N. (2000). Two-component signal transduction.Annu. Rev. Biochem.69183–215. 10.1146/annurev.biochem.69.1.183
63
StojiljkovicI.HantkeK. (1994). Transport of haemin across the cytoplasmic membrane through a haemin-specific periplasmic binding-protein-dependent transport system in Yersinia enterocolitica.Mol. Microbiol.13719–732. 10.1111/j.1365-2958.1994.tb00465.x
64
StojiljkovicI.KumarV.SrinivasanN. (1999). Non-iron metalloporphyrins: Potent antibacterial compounds that exploit haem/Hb uptake systems of pathogenic bacteria.Mol. Microbiol.31429–442. 10.1046/j.1365-2958.1999.01175.x
65
TomoniA.LeesJ.SantanaA. G.Bolanos-GarciaV. M.BastidaA. (2019). Pseudokinases: From allosteric regulation of catalytic domains and the formation of macromolecular assemblies to emerging drug targets.Catalysts9:778. 10.3390/catal9090778
66
TsokosC. G.PerchukB. S.LaubM. T. (2011). A dynamic complex of signaling proteins uses polar localization to regulate cell-fate asymmetry in Caulobacter crescentus.Dev. Cell20329–341. 10.1016/j.devcel.2011.01.007
67
van der RestM. E.LangeC.MolenaarD. (1999). A heat shock following electroporation induces highly efficient transformation of Corynebacterium glutamicum with xenogeneic plasmid DNA.Appl. Microbiol. Biotechnol.52541–545. 10.1007/s002530051557
68
VaradiM.AnyangoS.DeshpandeM.NairS.NatassiaC.YordanovaG.et al (2021). AlphaFold Protein Structure Database: Massively expanding the structural coverage of protein-sequence space with high-accuracy models.Nucleic Acids Res.50D439–D444. 10.1093/nar/gkab1061
69
WennerholdJ.BottM. (2006). The DtxR regulon of Corynebacterium glutamicum.J. Bacteriol.1882907–2918. 10.1128/JB.188.8.2907-2918.2006
70
WennerholdJ.KrugA.BottM. (2005). The AraC-type regulator RipA represses aconitase and other iron proteins from Corynebacterium under iron limitation and is itself repressed by DtxR*.J. Biol. Chem.28040500–40508. 10.1074/jbc.M508693200
71
WillettJ. W.CrossonS. (2017). Atypical modes of bacterial histidine kinase signaling.Mol. Microbiol.103197–202. 10.1111/mmi.13525
72
WolaninP. M.ThomasonP. A.StockJ. B. (2002). Histidine protein kinases: Key signal transducers outside the animal kingdom.Genome Biol.3:REVIEWS3013. 10.1186/gb-2002-3-10-reviews3013
73
ZhaoX. R.ChoiK. R.LeeS. Y. (2018). Metabolic engineering of Escherichia coli for secretory production of free haem.Nat. Catalysis1720–728. 10.1038/s41929-018-0126-1
Summary
Keywords
adaptive laboratory evolution (ALE), heme, pseudokinase, two-component system (TCS), histidine kinase, CA-domain
Citation
Krüger A and Frunzke J (2022) A pseudokinase version of the histidine kinase ChrS promotes high heme tolerance of Corynebacterium glutamicum. Front. Microbiol. 13:997448. doi: 10.3389/fmicb.2022.997448
Received
18 July 2022
Accepted
18 August 2022
Published
07 September 2022
Volume
13 - 2022
Edited by
Ulrike Kappler, The University of Queensland, Australia
Reviewed by
Miaomiao Shi, The University of Chicago, United States; Haibi Wang, Cornell University, United States
Updates

Check for updates
Copyright
© 2022 Krüger and Frunzke.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Julia Frunzke, j.frunzke@fz-juelich.de
This article was submitted to Microbial Physiology and Metabolism, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.