- Open Access
Thermotolerant genes essential for survival at a critical high temperature in thermotolerant ethanologenic Zymomonas mobilis TISTR 548
Biotechnology for Biofuelsvolume 10, Article number: 204 (2017)
High-temperature fermentation (HTF) technology is expected to reduce the cost of bioconversion of biomass to fuels or chemicals. For stable HTF, the development of a thermotolerant microbe is indispensable. Elucidation of the molecular mechanism of thermotolerance would enable the thermal stability of microbes to be improved.
Thermotolerant genes that are essential for survival at a critical high temperature (CHT) were identified via transposon mutagenesis in ethanologenic, thermotolerant Zymomonas mobilis TISTR 548. Surprisingly, no genes for general heat shock proteins except for degP were included. Cells with transposon insertion in these genes showed a defect in growth at around 39 °C but grew normally at 30 °C. Of those, more than 60% were found to be sensitive to ethanol at 30 °C, indicating that the mechanism of thermotolerance partially overlaps with that of ethanol tolerance in the organism. Products of these genes were classified into nine categories of metabolism, membrane stabilization, transporter, DNA repair, tRNA modification, protein quality control, translation control, cell division, and transcriptional regulation.
The thermotolerant genes of Escherichia coli and Acetobacter tropicalis that had been identified can be functionally classified into 9 categories according to the classification of those of Z. mobilis, and the ratio of thermotolerant genes to total genomic genes in Z. mobilis is nearly the same as that in E. coli, though the ratio in A. tropicalis is relatively low. There are 7 conserved thermotolerant genes that are shared by these three or two microbes. These findings suggest that Z. mobilis possesses molecular mechanisms for its survival at a CHT that are similar to those in E. coli and A. tropicalis. The mechanisms may mainly contribute to membrane stabilization, protection and repair of damage of macromolecules and maintenance of cellular metabolism at a CHT. Notably, the contribution of heat shock proteins to such survival seems to be very low.
Zymomonas mobilis is an efficient ethanologenic microbe that has been isolated from sugarcane or alcoholic beverages such as African palm wine, and it causes cider sickness and spoiling of beer . The organism bears an anaerobic catabolism via the Entner–Doudoroff pathway , which utilizes 1 mol of glucose to yield 2 mol of pyruvate, which is then decarboxylated to acetaldehyde and reduced to ethanol. Due to its strong metabolic activity and low ATP productivity compared to those of the Emden–Meyerhof pathway in the conventional ethanol producer yeast and high-yield ethanol production as a result of the Entner–Doudoroff pathway [1, 3] as well as the fact that the organism is generally regarded as being safe (GRAS) , Z. mobilis has been focused for its applications to production of useful materials including ethanol as a biofuel, oligosaccharides as food additives, and levan as a medicine [5, 6].
Since the ethanol fermentation process is exothermic [7, 8], ethanologenic microorganisms are exposed to heat stress in addition to other stresses including ethanol [9, 10]. Heat stress has an impact on their growth or viability [11, 12] to prevent fermentation, and the impact is enhanced in the presence of other inhibiting factors, i.e., low pH, high ethanol concentration, and high osmolarity [13,14,15,16,17,18]. Thus, thermotolerant Z. mobilis is thought to be beneficial for the production of useful materials. Z. mobilis TISTR 548 is a thermotolerant strain that can grow even at 39 °C [19,20,21], which is 5–10 °C higher than the optimum temperature for the same genus  and the same species [1, 23], and it can efficiently produce ethanol to an extent similar to that of ZM4 . However, information on the molecular mechanism of the thermotolerance of thermotolerant Z. mobilis is limited, though some heat shock proteins have been analyzed [24, 25].
Elucidation of the molecular mechanism of microbial survival at a critical high temperature (CHT) may be useful for the development of high-temperature fermentation systems, which have several advantages including reduction in cooling cost, saving of enzyme cost in simultaneous saccharification and fermentation or prevention of contamination of unfavorable microbes [26, 27]. We thus performed transposon mutagenesis of the thermotolerant Z. mobilis TISTR 548 to isolate thermosensitive mutants, each of which is defective of one of the so-called thermotolerant genes. The physiological functions of these genes allow us to decipher the molecular mechanism of its survival at a CHT. Moreover, we may be able to understand the general strategy of Gram-negative bacteria to cope with thermal stresses at their individual CHTs by comparison of the mechanism in Z. mobilis as α-proteobacteria with those of other bacteria, Escherichia coli as γ-proteobacteria and Acetobacter tropicalis as α-proteobacteria, that have been investigated [28, 29]. E. coli is intrinsically thermotolerant compared to general mesophilic microbes and used for production of useful materials like amino acids, hormones, or vaccines. Z. mobilis TISTR548 and A. tropicalis are thermotolerant and efficiently produces ethanol and acetic acid, respectively, at relatively high temperatures [19, 29]. Thus, the knowledge of the general strategy might be applicable for relatively thermosensitive mesophilic microbes that have been utilized for production of useful materials in fermentation companies.
Isolation of thermosensitive mutants by transposon mutagenesis in thermotolerant Z. mobilis
Thermotolerant Z. mobilis strain TISTR 548 was subjected to transposon mutagenesis via E. coli S17-1 harboring pSUP2021Tn10 as a donor strain for conjugal mating . The growth levels of about 8000 transconjugants obtained were compared on YPD plates at 30 and 39.5 °C, and thermosensitive ones that exhibited no or almost no growth at the high temperature were selected. They were subjected to repeated examination on YPD plates as a second screening and resultantly obtained 123 thermosensitive isolates were further subjected to the final screening in a YPD liquid medium under a static condition at 30 and 39.5 °C. Eventually, 38 isolates that exhibited defective or very weak growth in the liquid culture at the high temperatures were selected as thermosensitive mutants and were used for the following experiments.
The insertion site of Tn10 in the genome of each mutant was determined by thermal asymmetric interlaced (TAIL)-PCR followed by nucleotide sequencing. The genomic sequences flanking Tn10 were analyzed by using public databases to identify a disrupted gene. As a result, out of the 38 thermosensitive mutants, only 26 were found to have a Tn10 insertion in independent genes and 12 were overlapped (Additional file 1: Table S1). This overlapping suggests that the isolation of thermosensitive mutants was nearly saturated. The 26 thermosensitive mutants including 14 representatives showed impaired growth at 39 or 39.5 °C but a similar level of growth to that of the parental strain at 30 °C (Additional file 1: Figure S1).
The gene organization around each Tn10-inserted gene might cause a polar effect of the insertion on the transcription of a downstream gene(s) that is intrinsically transcribed by read-through from an upstream promoter(s). Such an organization was found in 12 of the 26 mutants (Additional file 1: Figure S2). The possibility of such polar effects was thus examined by RT-PCR with total RNA that had been prepared from cells grown at 30 and 39.5 °C (Additional file 1: Figure S3). The data suggest that all genes located downstream of the transposon-inserted genes are expressed at the same levels of expression as those in the parental strain. Therefore, it is thought that the thermosensitive phenotype of the 26 thermosensitive mutants is due to the disruption of each gene inserted by Tn10, not due to a polar effect on its downstream gene(s). Taken together, 26 independent thermosensitive mutants were obtained and thus 26 thermotolerant genes were identified in thermotolerant Z. mobilis TISTR 548.
Function and classification of thermotolerant genes in thermotolerant Z. mobilis
In order to know the physiological functions of thermotolerant genes, database searching was performed. As a result, out of the 26 thermotolerant genes, 24 genes were functionally annotated and classified into 9 categories of general metabolism, membrane stabilization, transporter, DNA repair, tRNA/rRNA modification, protein quality control, translation control, cell division, and transcriptional regulation (Table 1). The remaining 2 genes encode unknown proteins.
Group A consists of two genes related to general metabolism, ZZ6_0707 and ZZ6_1376, that encode glucose sorbosone dehydrogenase and 5, 10-methylenetetrahydrofolate reductase, respectively. The former oxidizes glucose or sorbosone and belongs to a family that possesses a beta-propeller fold. The best characterized in the family is soluble glucose dehydrogenase from Acinetobacter calcoaceticus, which oxidizes glucose to glucono-δ-lactone . The latter catalyzes the conversion of 5,10-methylenetetrahydrofolate, which is used for de novo thymidylate biosynthesis, to 5-methyltetrahydrofolate , which is used for methionine biosynthesis .
Group B is the largest group that consists of 12 genes related to membrane stabilization or membrane formation. Of these, ZZ6_1146 encodes glucosamine/fructose 6-phosphate aminotransferase, which is the first and rate-limiting enzyme in the hexosamine biosynthetic pathway and catalyzes the formation of glucosamine-6-phosphate using glutamine as an ammonia donor. This amino sugar is essential for the formation of a plethora of glycoconjugates for the peptidoglycan macromolecule in prokaryotes . ZZ6_0929 encodes glycosyltransferase group 1, which is involved in biosynthesis of the lipopolysaccharide (LPS) core . This enzyme has two putative conserved domains: one domain covering 94% of the protein is named GT1_mtfB_like. MtfB (mannosyltransferase B) in E. coli has been shown to direct growth of the O9-specific polysaccharide chain . The other covering 53% of the protein is named RfaB and is involved in assembly of the lipopolysaccharide core in E. coli . ZZ6_0923 encodes phospholipase D/transphosphatidylase possessing the domain of cardiolipin synthase, which catalyzes phosphatidyl group transfer from one phosphatidylglycerol molecule to another to form cardiolipin and glycerol . The cls − for a defective cardiolipin synthase that shows a low level of cardiolipin in phospholipid composition has been reported , and the cls gene may be related to membrane stabilization. ZZ6_1551 encodes squalene hopene cyclase, which is a key enzyme for hopanoid biosynthesis and cyclizes squalene to hopene . Hopanoids belong to a triterpene series widespread among prokaryotes and play roles in membrane stabilization. Several different hopanoid derivatives are present in Z. mobilis . ZZ6_1046 and ZZ6_1043 encode TolQ and TolB, respectively. Both proteins are components of the Tol–Pal (peptidoglycan-associated lipoprotein) system, which is involved in the maintenance of outer membrane stability . Tol proteins are located in the cell envelope and are thought to be involved in the integration of some outer membrane components such as porins and lipopolysaccharides . ZZ6_1254 encodes a protein-export membrane protein, SecD, in the Sec system, and mutations of the gene exhibit pleiotropic defects in protein export in E. coli . ZZ6_1477 encodes a preprotein import (inner membrane) translocase subunit, Tim44. In mitochondria, Tim44 is a component to anchor mHsp70 to the TIM23 channel and associates transiently with the TIM23 complex for import of matrix-localized proteins in mitochondria . ZZ6_0158 encodes an autotransporter secretion inner membrane protein, TamB, that forms a complex of the translocation and assembly module with the outer membrane protein, TamA. The complex functions in translocation of autotransporters across the outer membrane . ZZ6_1210 encodes a competence protein, ComEC, that is a DNA transformation transporter (DNA-T) core component (KEGG). Competent cells generally possess a DNA transport complex that is most likely composed of surface-exposed DNA receptors, which facilitate DNA translocation through the cell wall, membrane pores, and motor molecules that power DNA transport . ZZ6_0840 encodes a hypothetical transmembrane protein that possesses a zinc finger domain at its N-terminal portion and a Hid1 superfamily domain at its middle portion as putative conserved protein domains. Hid1 is a high-temperature-induced dauer-forming protein 1 with many putative transmembrane segments in Caenorhabditis elegans . ZZ6_0541 encodes a protein bearing an SH3-like domain (COG3807). There are many SH3-like domain-containing proteins , but the function of the domain has not been clarified yet except for SH3-like domain-dependent interaction between CheA and CheW .
Group C as transporter includes a single gene, ZZ6_1289, that encodes a putative Fe2+/Mn2+ transporter, which shares 58% identity to Fe2+/Mn2+ transporter pcl1 in Acetobacter pasteurianus.
Group D consists of genes for DNA repair. ZZ6_0616 encodes the DNA repair protein RadC. RadC functions specifically in recombination repair that is associated with a replication fork and is required for growth-medium-dependent repair of DNA double strand breaks in E. coli . ZZ6_0934 encodes XseA, a large subunit of exonuclease VII that is implicated in the resection of a nicked mismatched strand in a methyl-directed mismatch repair pathway . ZZ6_0681 encodes the DNA repair protein RadA. In E. coli, RadA is involved in recombination and recombination repair and is likely involved in the stabilization or processing of branched DNA molecules or blocked replication forks . radA mutants show a modest decrease in survival after UV or X-irradiation exposure .
Group E consists of one gene for tRNA/rRNA modification. ZZ6_0023 encodes SpoU, which is a tRNA/rRNA methyltransferase. This enzyme may contribute to stabilization of the structure of tRNA or ribosome . Analysis of the nucleoside modification pattern of tRNA, 16S rRNA, and 23S rRNA in E. coli has shown that the modified nucleoside 2′-O-methylguanosine, present in a subset of tRNAs at residue 18, is completely absent in the spoU mutant .
Group F genes are related to protein quality control. ZZ6_1659 encodes a Zn-dependent peptidase (peptidase with a M16 domain) (KEGG). The M16 family of zinc peptidases comprises a pair of homologous domains that form two halves of a ‘‘clam-shell’’ surrounding the active site, and closure of the clam-shell is required for proteolytic activity . ZZ6_0980 encodes the serine protease DegP, and the orthologue gene has been identified as a thermotolerant gene in E. coli and A. tropicalis [28, 29]. DegP is a chaperone/serine protease located in the periplasm and acts to remove damaged proteins [57, 58].
Group G consists of one gene for translation control. ZZ6_0702 encodes the ATP-dependent helicase HrpB, that acts as an RNA helicase. Some in this helicase group unwind RNA molecules with a 3′ to 5′ polarity . HrpA is an orthologue of HrpB involved in mRNA processing in E. coli. hrpA mutations in regions for predicted binding and hydrolysis of nucleotide triphosphate abolish the ability for mRNA processing .
Group H as cell division includes ZZ6_0979 for ParA/MinD-like ATPase. In E. coli, MinD activates a MinC-dependent mechanism responsible for the inactivation of potential division sites and renders the division inhibition system sensitive to MinE, which are required for correct placement of a division site . MinD binds ATP and bears ATPase activity. On the other hand, ParA is required for the equipartition of P1 plasmids during cell division .
Group I consists of one gene related to transcriptional regulation. ZZ6_0019 encodes the flavoprotein WrbA, that binds to the tryptophan repressor TrpR and functions as an accessory element in blocking the TrpR-specific transcriptional process . WrbA enhances the formation and/or stabilization of noncovalent complexes between TrpR holorepressor and its primary operator targets . WrbA also functions as an NAD(P)H/quinone oxidoreductase  and belongs to the family of multimeric flavodoxin-like proteins  as a new type (type IV) of NAD(P)H:quinone oxidoreductase, which protects cells against oxidative stress  and may prepare cells for long-term maintenance under stress conditions .
There are two genes that deviate from the 9 categories. ZZ6_0962 is named as a pseudogene but should have a crucial function at a high temperature as observed in this study. The pseudogene has an inserted transposon in the gene, but the contribution of the transposon to thermotolerance is unknown. ZZ6_0861 encodes a hypothetical small protein consisting of 82 amino acid residues.
Effect of supplemented MgCl2 on growth of thermosensitive mutants
Mg2+ is known to stabilize the outer membrane structure in cells by binding extracellularly  and the thermosensitive phenotype of mutants due to the disruption of genes for membrane stabilization is suppressed by the addition of MgCl2 at a CHT in E. coli . Thus, the effect of MgCl2 on growth of thermosensitive mutants in Z. mobilis was tested at its CHT.
Thermosensitive mutants and the parental strain were grown in YPD medium with or without 20 mM MgCl2 at 39.5 °C for 24 h under a static condition (Additional file 1: Figure S4; Table 1). The growth of 13 thermosensitive mutants was significantly improved by the supplementation of MgCl2, 120–260% of that of the parental strain. Eight of them were in Group B and have disrupted genes for membrane stabilization or membrane formation. These results suggest that Mg2+ stabilizes the membrane structure at a CHT and protects cells from heat, as has been proposed in E. coli.
Effect of ethanol stress on growth of thermosensitive mutants
Zymomonas mobilis as an efficient ethanol producer is often exposed to ethanol stress under fermentation conditions. The effect of exogenous ethanol on thermosensitive mutants was thus examined on YPD plates containing 2.0 or 2.5% ethanol at 30 °C. In consequence, about half of the thermosensitive mutants exhibited repressed growth in the presence of ethanol, less than 50% growth compared to that in the absence of ethanol (Table 1). Interestingly, most of the thermosensitive mutants that were classified into the membrane stabilization group exhibited sensitivity to ethanol stress, and most of the ethanol-sensitive mutants were classified into the group in which the thermosensitive growth phenotype was suppressed by the addition of MgCl2. Therefore, these results suggest that the mechanism of thermotolerance at a CHT partially overlaps with that of ethanol stress resistance and allows us to speculate that stabilization of the membrane structure is one of crucial points for ethanol tolerance.
In this study, we isolated 38 thermosensitive mutants by transposon mutagenesis and finally identified 26 thermotolerant genes that are required for survival at a CHT in thermotolerant Z. mobilis TISTR 548. Physiological functions and classification of these gene products may allow us to obtain a clue regarding the thermotolerance mechanism of this organism. The gene products were classified into 9 categories (Table 1). About half of them are related to membrane stabilization or membrane formation including enzymes for peptidoglycan or lipid biosynthesis and proteins for protein secretion systems. Most of these, genes for glucosamine/fructose 6-phosphate aminotransferase (ZZ6_1146), glycosyltransferase (ZZ6_0929), squalene hopene cyclase (ZZ6_1551), protein export membrane protein SecD (ZZ6_1254), autotransporter secretion inner membrane protein TamB (ZZ6_0158), competence protein ComEC (ZZ6_1210), hypothetical transmembrane protein (ZZ6_0840), and hypothetical transmembrane protein (ZZ6_0541) were found to be required for ethanol tolerance. Therefore, it is thought that membrane stabilization and maintenance are essential for survival at a CHT. Surprisingly, as found in E. coli , there was no heat shock protein in these thermotolerant gene products except for DegP, suggesting that not all heat shock proteins may be essential for survival under high temperatures. DegP, which functions in the periplasm as a chaperone at low temperatures and as a protease at high temperatures , is thought to play a role in the maintenance of homeostasis of the periplasm or membranes. In E. coli, groEL as an essential gene was induced at a CHT  and thus some heat shock proteins may be required under such an extreme condition.
Thermotolerant genes have also been identified in E. coli BW25113 and A. tropicalis SKU1100: 72 and 24 genes, respectively [28, 29; unpublished data]. The thermotolerant genes of the two microbes can be classified into 9 categories according to the classification of those of Z. mobilis, and the number and distribution of these genes are shown in Table 2. The ratios of thermotolerant genes to total genomic genes in Z. mobilis, E. coli, and A. tropicalis are 1.47, 1.68, and 0.70%, respectively. We do not know the reason why the ratio in A. tropicalis is relatively low. In the case of E. coli, a single-gene knockout library was used for screening thermosensitive mutants and thus almost all of the genes except for essential genes were examined. On the other hand, in the case of Z. mobilis and A. tropicalis, transposon mutagenesis was applied for screening thermosensitive mutants, and the ratios of the number of thermotolerant genes, for each of which two or more transposon-inserted mutants were isolated, to the total number of thermotolerant genes (Additional file 1: Table S1)  were 35 and 21%, respectively. Therefore, the low ratio of multiple mutants for the same gene in A. tropicalis suggests the possibility that there are still unidentified thermotolerant genes in A. tropicalis SKU1100. In all categories except for general metabolism, ratios of thermotolerant genes in Z. mobilis are closer to those in E. coli than those in A. tropicalis. Notably, Z. mobilis has a higher ratio of thermotolerant genes for membrane stabilization than the ratios in other two microbes: 46, 25, and 20% in Z. mobilis, E. coli, and A. tropicalis, respectively.
On the other hand, E. coli possesses several discriminating sets of thermotolerant genes, which are absent in the other two microbes: 4 genes (aceE, aceF, lpd, and lipA) for pyruvate metabolism, 3 genes (atpA, atpD, and atpG) for ATPase, 3 genes (cydB, yhcB, and cydD) for ubiquinol oxidase or its formation, and 3 genes (ubiE, ubiH, and ubiX) for ubiquinone biosynthesis in the category of general metabolism, 8 genes (gmhB, lpcA, rfaC, rfaD, afaE, rfaF, rfaG, and lpxL) for lipopolysaccharide biosynthesis and 5 genes (ydcL, yfdL, ynbE, nlpI, and ycdO) for peptidoglycan-associated lipoproteins or predicted lipoproteins in the category of membrane stability, 5 genes (dnaQ, holC, priA, ruvA, and ruvC) for DNA double-strand break repair in the category of DNA repair, and 6 genes (iscS, yheL, yheM, yheN, yhhP, and yccM) for a sulfur relay system in the category of tRNA modification [28; unpublished data]. Of these sets, genes for the lipopolysaccharide biosynthesis and the sulfur relay system are postulated to have been acquired by horizontal gene transfer . The genes in the 4 categories described above seem to contribute to specific strategies for thermotolerance in E. coli [28; some thermotolerant genes will be described elsewhere].
There are common thermotolerant genes or thermotolerant genes related to the same physiological function or pathway among the three microbes. In the category of protein quality control, the three microbes share degP and both Z. mobilis and A. tropicalis have a gene for Zn-dependent protease (ZZ6_1659 and ATPR_0429, respectively). In membrane stabilization, one gene related to hopanoid biosynthesis is present in Z. mobilis and A. tropicalis (shc and ATPR_1188, respectively) and two to three genes for the Tol-Pal system are present in Z. mobilis (tolQ and tolB) and E. coli (pal, tolQ and tolR). One gene related to MinC-dependent cell division inhibition in cell division is present in Z. mobilis and A. tropicalis (minD and minC, respectively), and wrbA in transcriptional regulation and nhaA for the Na+/H+ antiporter in transporters are shared by Z. mobilis and A. tropicalis. On the basis of the functions of these genes and combinations of other thermotolerant genes in each category, some common strategies for thermotolerance have emerged: in the category of membrane stabilization, synthesis or modification of peptidoglycan and maintenance of integrity for all three microbes, and hopanoid or lipid synthesis for Z. mobilis and A. tropicalis; in DNA repair, double-strand DNA repair, which may be accumulated at a CHT, for Z. mobilis and E. coli; tRNA modification, probably for a stable structure at such a high temperature, for Z. mobilis and E. coli; in chaperone and protease, removal of damaged proteins, especially by periplasmic serine protease DegP, for all three microbes; control of chromosome segregation for E. coli and A. tropicalis, and control of cell division for all three microbes; and in transcriptional regulation, Trp repressor-binding protein WrbA (still unclear why necessary) for Z. mobilis and A. tropicalis. In addition, import or export of some metal ions may be important probably for keeping homeostasis of some ions, export of toxic ions or maintenance of membrane potential.
At a CHT, several problems including protein unfolding or increase in membrane fluidity occur. Reactive oxygen species increase as the temperature increases , causing the damage of macromolecules including DNA [70, 71]. The requirement of genes for the 9 categories allows us to make speculations about various types of damage of membrane and proteins or about the abnormal structures of macromolecules including proteins, DNAs and RNAs at a CHT. Microbes would have thus acquired thermotolerant genes to overcome these problems. Moreover, it is assumed that these genes are involved in the response of cells to other stresses including osmotic stress or oxygen stress. In fact, Z. mobilis increases thermotolerance by the addition of sorbitol  and exhibits faster growth and higher ethanol production under a static condition than that under a shaking condition [19, unpublished]. Further experiments are required for clarifying this assumption.
The thermotolerant genes of thermotolerant ethanologenic Z. mobilis TISTR 548 have been identified. Comparison with thermotolerant genes in E. coli and A. tropicalis reveal that these genes of the three microbes can be classified into 9 categories and that there are common thermotolerant genes or thermotolerant genes related to the same physiological function or pathway among the three microbes, which suggest several common strategies, including membrane stabilization, protection and repair of macromolecules of proteins, DNAs and RNAs, and maintenance of cellular metabolism-like cell division, transcription or translation, for the three microbes to survive at CHT. Considering the genetic conversion of non-thermotolerant to thermotolerant bacteria, such strategies might be applicable.
A DNA sequencing Kit (ABI PRISM ® Terminator v 3.1 Cycle sequencing Kit) was obtained from Applied Biosystem Japan. Oligonucleotide primers were synthesized by Proligo Japan K.K. (Tokyo, Japan). Other chemicals were all of analytical grade and obtained from commercial sources.
Microorganisms and media
Zymomonas mobilis TISTR 548 [19, 20] and its derivatives were grown in YPD (3% glucose, 0.5% peptone, and 0.3% yeast extract) medium. E. coli S17-1 harboring pSUP2021 Tn10  was grown in LB (0.5% yeast extract, 1% NaCl, and 1% Bactotryptone) medium supplemented with 12.5 µg/ml of tetracycline.
Conjugation and transposon mutagenesis
Escherichia coli S17-1 harboring pSUP2021 Tn10 as a donor for conjugal mating was grown in LB medium containing 12.5 µg/ml of tetracycline under a shaking condition at 100 rpm at 37 °C. The recipient Z. mobilis TISTR 548 was grown in YPD medium under a static condition at 30 °C. Cells of both strains were grown to the mid-log phase, washed three times with LB medium, recovered by centrifugation at 5000 rpm for 1 min, and suspended in a small volume of LB medium. Both cell suspensions were then mixed at a ratio of donor and recipient of 3:2 and stood for 3 h at 30 °C. The suspensions were spotted on the surfaces of LB agar plates and incubated at 30 °C for 5 h. After the mating steps, cells were recovered, resuspended in a small volume of YPD medium, and spread on YPD agar plates containing 0.15% acetic acid and 12.5 µg/ml of tetracycline. Transconjugants (transposon-inserted mutants) that appeared on the plates after 3-day incubation at 30 °C were subjected to the following screening.
Screening of thermosensitive mutants
About 8000 transconjugants were subjected to the first screening in which they were grown at 30 and 39.5 °C on YPD agar plates. Transposon-inserted mutants that showed no or almost no growth on the plates at 39.5 °C were selected for the next screening. The second screening was performed under the same condition as that in the first screening. Selected mutants were then subjected to the last screening in which their thermosensitivity was examined in 2-ml liquid culture of YPD medium at 30 and 39.5 °C for 24 h under a static condition. Cell growth was determined by measuring cell turbidity at OD550. Mutants that showed a value at OD550 significantly less than that of the parent strain were selected and defined as thermosensitive mutants.
Examination of the effects of heat and ethanol stresses on growth of thermosensitive mutants
Thermosensitive mutants and the parental strain were pre-cultured in YPD medium under a static condition at 30 °C until a mid-log phase. For the heat stress experiment, the pre-cultured cells were serially diluted with YPD medium, spotted on YPD agar plates, and incubated at 30, 38, 39, and 39.5 °C for 40 h. For the ethanol stress experiment, the pre-cultured cells were serially diluted, spotted on YPD plates supplemented with 2.0 or 2.5% ethanol, and incubated at 30 °C for 40 h. Growth ability was examined in triplicate.
Effect of Mg2+ on growth of thermosensitive mutants
Thermosensitive mutants and the parental strain were pre-cultured in YPD medium under a static condition at 30 °C until a mid-log phase. The pre-cultured cells were inoculated in YPD medium with or without 20 mM MgCl2 and incubated at 39.5 °C for 24 h under a static condition. The experiments were performed more than 3 times. The significance of the effect of MgCl2 on cell growth was evaluated by a t test.
Identification of the transposon (Tn10)-inserted site in a thermosensitive mutant genome by TAIL-PCR followed by nucleotide sequencing
The Tn10-inserted site in the genome of each thermosensitive mutant was determined by TAIL-PCR  followed by nucleotide sequencing. The genomic DNA from thermosensitive mutants was isolated as described previously . The concentration of isolated genomic DNA was measured by using Nanodrop (Nanodrop Technologies, Wilmington, DE). TAIL-PCR was performed by using TaKaRa PCR Thermal Cycler Dice® mini (TaKaRa). Three specific primers for TAIL-PCR were TnISR-1 (GATCCTCTCGTTTGTTGCGGTCAGGCC) , TnISR-1.5 (AGGGCTGCTAAAGGAAGCGG) (this work) and TnISR-2 (ACGAAGCGCAAAGAGGAAGCAGG) , and an arbitrary degenerated primer was AD2 (GTNCGASWCANAWGTT) . The first PCR was carried out in a 50-µl mixture containing 10 ng of chromosome DNA, 5.0 µM TnISR-1, 25 µM AD2 primer, 500 µM each of dNTPs, 0.5 U PrimeSTAR (TaKaRa) and 1× buffer supplied for the enzyme. Two percent of the first PCR product was used as a template for the second PCR, which was performed using the same reaction mixture as that used for the first PCR except that TnISR-1.5 was used as a specific primer. The third PCR was also performed using the same reaction mixture as that used for the first PCR except that TnISR2 was used as a specific primer and the concentration of AD2 was reduced to 12.5 µM . The second or third PCR product was purified by using a PCR product purification kit (Qiagen) and subjected to nucleotide sequencing on an ABI PRISM 310 Genetic Analyzer (Applied Biosystems) or DNA Sequencer GenomeLab GeXP (Beckman Coulter). The sequencing reaction was performed with a BigDye® Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems) or a GenomeLab Dye Terminator Cycle Sequencing with Quick Start Kit (Beckman Coulter).
Zymomonas mobilis cells were grown in 50 ml of YPD medium under a static condition at 30 °C until exponential phase, and then the temperature was increased to 39.5 °C and the cultivation was continued for 8 min. As a control, the cultivation was continued for 8 min at 30 °C. Total RNA was prepared from these heat-stressed or not heat-stressed cells by the hot phenol method . RT-PCR analysis was performed using an mRNA-selective RT-PCR kit (TaKaRa) and primers (Additional file 1: Table S2) to examine the expression of immediate downstream genes of Tn10-inserted genes as described previously . The reverse transcription reaction was carried out at 42 °C for 15 min, followed by PCR at 85 °C for 1 min, 45 °C for 1 min, and extension at 72 °C for 1 min, using the two specific primers for each gene. After the completion of 15, 20, 25, and 30 cycles, the PCR products were analyzed by 0.9% agarose gel electrophoresis and stained with ethidium bromide . The relative amounts of RT-PCR products on the gel were compared by measuring the density of bands on the gel by using image J (https://imagej.nih.gov/ij/). Under our conditions, the RNA-selective RT-PCR was able to specifically detect mRNA because no band was observed when reverse transcriptase was omitted.
The intrinsic gene that was inserted by Tn10 in each thermotolerant mutant was confirmed to be a thermotolerant gene after analyses of the gene organization and/or expression of its downstream gene. Thermotolerant genes were then subjected to functional classification by bioinformatics analysis mainly according to the instructions of KEGG (http://www.genome.jp/kegg/), NCBI (http://www.ncbi.nlm.nih.gov/), Inter Pro (http://www.ebi.ac.uk/interpro/), and Uniprot (http://www.uniprot.org/). Protein type was analyzed by TMHMM (http://www.cbs.dtu.dk/services/TMHMM/). Homology searching and alignment were performed using BLAST . The Z. mobilis TISTR 548 thermotolerant genes were designed as ZZ6_XXXX according to Z. mobilis subsp. mobilis ATCC29191 because the genome sequence of TISTR 548 was found to be almost identical to that of ATCC29191 after draft sequencing (unpublished).
Thailand Institute of Scientific and Technological Research
generally regarded as being safe
critical high temperature
thermal asymmetric interlaced PCR
DNA transformation transporter
reduced form of nicotinamide adenine dinucleotide
reduced form of nicotinamide adenine dinucleotide phosphate
Swings J, De Ley J. The biology of Zymomonas. Bacteriol Rev. 1977;41:1–46.
Gibbs M, De Moss RD. Anaerobic dissimilation of C14. Labelled glucose fructose by Pseudomonas lindneri. J Biol Chem. 1954;207:689–94.
Seo J, Chong H, Park HS, Yoon K, Jung C, Kim JJ, et al. The genome sequence of the ethanologenic Bacterium Zymomonas mobilis ZM4. Nat Biotechnol. 2005;23:63–8.
Yang S, Pelletier DA, Lu TYS, Brown SD. The Zymomonas mobilis regulator hfq contributes to tolerance against multiple lignocellulosic pretreatment inhibitors. BMC Microbiol. 2010;135:1–11.
Calazans GMT, Lopes CE, Lima RMOC, Defranca FP. Antitumor activities of levans produced by Zymomonas mobilis strains. Biotechnol Lett. 1997;19:19–21.
Yoo SH, Yoon EJ, Cha J, Lee HG. Antitumor activity of levan polysaccharides from selected microorganisms. Int J Biol Macromol. 2004;34:37–41.
Uden VN, Duarte HDC. Effects of ethanol on the temperature profile of Saccharomyces cerevisiae. Z Allg Mikrobiol. 1981;21:743–50.
Ghose TK, Bandyopadhyay KK. Studies on immobilized Saccharomyces cerevisiae. II. Effect of temperature distribution on continuous rapid ethanol formation in molasses fermentation. Biotechnol Bioeng. 1982;24:797–804.
Attfield PV. Stress tolerance: the key to effective strains of industrial baker’s yeast. Nat Biotechnol. 1997;15:1351–7.
Wang Y, Gong L, Liang J, Zhang Y. Effects of alcohol on expressions of apoE in mice livers and brains. Wei Sheng Yan Jiu. 2007;36:737–40.
Basso LC, de Amorim HV, de Oliveira AJ, Lopes ML. Yeast selection for fuel ethanol production in Brazil. FEMS Yeast Res. 2008;8:1155–63.
Babiker MA, Banat A, Hoshida H, Ano A, Nonklang S, Akada R. High-temperature fermentation: how can processes for ethanol production at high temperatures become superior to the traditional process using mesophilic yeast? Appl Microbiol Biotechnol. 2010;85:861–7.
Piper PW. The heat shock and ethanol stress responses of yeast exhibit extensive similarity and functional overlap. FEMS Microbiol Lett. 1995;134:121–7.
Carmelo V, Santos R, Viegas CA, Sa´-Correia I. Modification of Saccharomyces cerevisiae thermotolerance following rapid exposure to acid stress. Int J Food Microbiol. 1998;42:225–30.
Ciani M, Beco L, Comitini F. Fermentation behaviour and metabolic interactions of multistarter wine yeast fermentations. Int J Food Microbiol. 2006;108:239–45.
Pizarro F, Varela C, Martabit C, Bruno C, Pe´rez-Correa JR, Agosin E. Coupling kinetic expressions and metabolic networks for predicting wine fermentations. Biotechnol Bioeng. 2007;98:986–98.
Coleman MC, Fish R, Block DE. Temperature-dependent kinetic model for nitrogen-limited wine fermentations. Appl Environ Microb. 2007;73:5875–84.
Gibson BR, Lawrence SJ, Leclaire JP, Powell CD, Smart KA. Yeast responses to stresses associated with industrial brewery handling. FEMS Microbiol Rev. 2007;31:535–69.
Sootsuwan K, Irie A, Murata M, Lertwattanasakul N, Thanonkeo P, Yamada M. Thermotolerant Zymomonas mobilis: comparison of ethanol fermentation capability with that of an efficient type strain. Open Biotechnol J. 2007;1:59–65.
Charoensuk K, Irie A, Lertwattanasakul N, Sootsuwan K, Thanonkeo P, Yamada M. Physiological importance of cytochrome c peroxidase in ethanologenic Thermotolerant Zymomonas mobilis. J Mol Microbiol Biotechnol. 2011;20:70–82.
Thanonkeo P, Thanonkeo S, Charoensuk K, Yamada M. Ethanol production from Jerusalem artichoke (Helianthus tuberosus L.) by Zymomonas mobilis TISTR 548. Afr J Biotechnol. 2011;10:10691–7.
Manaia CM, Moore ERB. Pseudomonas thermotolerans sp. nov., a thermotolerant species of the genus Pseudomonas sensustricto. Int J Syst Evol Microbiol. 2002;52:2203–9.
Saeki A, Theeragool G, Matsushita K, Toyama H, Lotong N, Adachi O. Development of thermotolerant acetic acid bacteria useful for vinegar fermentation at higher temperatures. Biosci Biotechnol Biochem. 1997;61:138–45.
Michel GPF, Starka J. Effect of ethanol and heat stresses on the protein pattern of Zymomonas mobilis. J Bacteriol. 1986;165:1040–2.
Thanonkeo P, Sootsuwan K, Leelavacharamas V, Yamada M. Cloning and transcriptional analysis of groES and groEL in ethanol-producing bacterium Zymomonas mobilis TISTR 548. Pak J Biol Sci. 2007;10:13–22.
Rodrussamee N, Lertwattanasakul N, Hirata K, Suprayogi, Limtong S, Kosaka T, Yamada M. Growth and ethanol fermentation ability on hexose and pentose sugars and glucose effect under various conditions in thermotolerant yeast Kluyveromyces marxianus. Appl Microbiol Biotechnol. 2011;90:1573–86.
Murata M, Nitiyon S, Lertwattanasakul N, Sootsuwan K, Kosaka T, Thanonkeo P, Limtong S, Yamada M. High-temperature fermentation technology for low-cost bioethanol. J Jpn Inst Energy. 2015;94:1154–212.
Murata M, Fujimoto H, Nishimura K, Charoensuk K, Nagamitsu H, Raina S, et al. Molecular strategy for survival at a critical high temperature in Eschierichia coli. PLoS ONE. 2011;6:e20063.
Soemphol W, Deeraksa A, Matsutani M, Yakushi T, Toyama H, Adachi O, et al. Global analysis of the genes involved in the thermotolerance mechanism of thermotolerant Acetobacter tropicalis SKU1100. Biosci Biotechnol Biochem. 2011;75:1921–8.
Deeraksa A, Moonmangmee S, Toyama HMY, Adachi O, Matsushita K. Characterization and spontaneous mutation of a novel gene, polE, involved in pellicle formation in Acetobacter tropicalis SKU1100. Microbiology. 2005;151:4111–20.
Oubrie A, Rozeboom HJ, Dijkstra BW. Active-site structure of the soluble quinoprotein glucose dehydrogenase complexed with methylhydrazine: a covalent cofactor-inhibitor complex. Proc Natl Acad Sci USA. 1999;96:11787–91.
Kim SH, Lee BR, Kim JN, Kim BG. NdgR, a common transcriptional activator for methionine and leucine biosynthesis in Streptomyces coelicolor. J Bacteriol. 2012;94:6837–46.
Badet-Denisot MA, Fernandez-Herrero LA, Berenguer J, Ooi T, Badet B. Characterization of l-glutamine:d-fructose-6-phosphate amidotransferase from an extreme thermophile Thermus thermophilus HB8. Arch Biochem Biophys. 1997;337:129–36.
Roncero C, Casadaban MJ. Genetic analysis of the genes involved in synthesis of the lipopolysaccharide core in Escherichia coli K-12: three operons in the rfa locus. J Bacteriol. 1992;174:3250–60.
Kido N, Torgov VI, Sugiyama T, Uchiya K, Sugihara H, Komatsu T, et al. Expression of the O9 polysaccharide of Escherichia coli: sequencing of the E. coli O9 rfb gene cluster, characterization of mannosyl transferases, and evidence for an ATP-binding cassette transport system. J Bacteriol. 1995;177:2178–87.
Pradel E, Parker CT, Schnaitman CA. Structures of the rfaB, rfaI, rfaJ, and rfaS genes of Escherichia coli K-12 and their roles in assembly of the lipopolysaccharide core. J Bacteriol. 1992;174:4736–45.
Tropp BE. Cardiolipin synthase from Escherichia coli. Biochim Biophys Acta. 1997;1348:192–200.
Shibuya I, Miyazaki C, Ohta A. Alteration of phospholipid composition by combined defects in phosphatidylserine and cardiolipin synthases and physiological consequences in Escherichia coli. J Bacteriol. 1985;161:1086–92.
Siedenburg G, Jendrossek D. Squalene-hopene cyclases. Appl Environ Microbiol. 2011;77:3905–15.
Hermans MA, Neuss B, Sahm H. Content and composition of hopanoids in Zymomonas mobilis under various growth conditions. J Bacteriol. 1991;173:5592–5.
Kampfenkel K, Braun V. Membrane topologies of the TolQ and TolR proteins of Escherichia coli: inactivation of TolQ by a missense mutation in the proposed first transmembrane segment. J Bacteriol. 1993;75:4485–91.
Ray MC, Germon P, Vianney A, Portalier R, Lazzaroni JC. Identification by genetic suppression of Escherichia coli TolB residues important for TolB–Pal interaction. J Bacteriol. 2000;182:821–4.
Gardel C, Benson S, Hunt J, Michaelis S, Beckwith J. secD, a new gene involved in protein export in Escherichia coli. J Bacteriol. 1987;169:1286–90.
Slutsky-Leiderman O, Marom M, Iosefson O, Levy R, Maoz S, Azem A. The interplay between components of the mitochondrial protein translocation motor studied using purified components. J Biol Chem. 2007;282:33935–44942.
Selkrig J, Mosbahi K, Webb CT, Belousoff MJ, Perry AJ, Wells TJ, et al. Discovery of an archetypal protein transport system in bacterial outer membranes. Nat Struct Mol Biol. 2012;19:506–10.
Chen I, Dubnau D. DNA uptake during bacterial transformation. Nat Rev Microbiol. 2004;2:241–9.
Ailion M, Thomas JH. Isolation and characterization of high-temperature-induced dauer formation mutants in Caenorhabiditis elegans. Genetics. 2003;165:127–44.
Whisstock JC, Lesk AM. SH3 domains in prokaryotes. Trends Biochem Sci. 1999;24:132–3.
Bilwes AM, Alex LA, Crane BR, Simon MI. Structure of CheA, a signal-transducing histidine kinase. Cell. 1999;96:131–41.
Saveson CJ, Lovett ST. Tandem repeat recombination induced by replication fork defects in Escherichia coli requires a novel factor, RadC. Genetics. 1999;152:5–13.
Harris RS, Ross KJ, Lombardo MJ, Rosenberg SM. Mismatch repair in Escherichia coli cells lacking single-strand exonucleases ExoI, ExoVII, and RecJ. J Bacteriol. 1998;180:989–93.
Beam CE, Saveson CJ, Lovett ST. Role for radA/sms in recombination intermediate processing in Escherichia coli. J Bacteriol. 2002;184:6836–44.
Sargentini NJ, Smith KC. Quantitation of the involvement of the recA, recB, recC, recF, recJ, recN, lexA, radA, radB, uvrD, and umuC genes in the repair of X-ray-induced DNA double-strand breaks in Escherichia coli. Radiat Res. 1986;107:58–72.
Decatur WA, Fournier MJ. rRNA modifications and ribosome function. Trends Biochem Sci. 2002;27:344–51.
Persson BC, Jäger G, Gustafsson C. The spoU gene of Escherichia coli, the fourth gene of the spoT operon, is essential for tRNA (Gm18) 2′-O-methyltransferase activity. Nucleic Acids Res. 1997;25:4093–147.
Aleshin AE, Gramatikova S, Hura GL, Bobkov A, Strongin AY, Stec B, et al. Crystal and solution structures of a prokaryotic M16B peptidase: an open and shut case. Structure. 2009;17:1465–75.
Lipinska B, Zylicz M, Georgopoulos C. The HtrA (DegP) protein, essential for Escherichia coli survival at high temperatures, is an endopeptidase. J Bacteriol. 1990;172:1791–7.
Meltzer M, Hasenbein S, Mamant N, Merdanovic M, Poepsel S, Hauske P, et al. Structure, function and regulation of the conserved serine proteases DegP and DegS of Escherichia coli. Res Microbiol. 2009;160:660–6.
Lee CG, Hurwitz J. A new RNA helicase isolated from HeLa cells that catalytically translocates in the 3′ to 5′ direction. J Biol Chem. 1992;267:4398–407.
Koo JT, Choe J, Moseley SL. HrpA, a DEAH-box RNA helicase, is involved in mRNA processing of a fimbrial operon in Escherichia coli. Mol Microbiol. 2004;52:1813–26.
de Boer PA, Crossley RE, Hand AR, Rothfield LI. The MinD protein is a membrane ATPase required for the correct placement of the Escherichia coli division site. EMBO J. 1991;10:4371–80.
Nordström K, Austin SJ. Mechanisms that contribute to the stable segregation of plasmids. Annu Rev Genet. 1989;23:37–69.
Yang W, Ni L, Somerville RL. A stationary-phase protein of Escherichia coli that affects the mode of association between the trp repressor protein and operator-bearing DNA. Proc Natl Acad Sci USA. 1993;90:5796–800.
Patridge EV, Ferry JG. WrbA from Escherichia coli and Archaeoglobus fulgidus is an NAD(P)H: quinone oxidoreductase. J Bacteriol. 2006;188:3498–506.
Grandori R, Khalifah P, Boice JA, Fairman R, Giovanielli K, Carey J. Biochemical characterization of WrbA, founding member of a new family of multimeric flavodoxin-like proteins. J Biol Chem. 1998;273:20960–6.
Chang DE, Smalley DJ, Conway T. Gene expression profiling of Escherichia coli growth transitions: an expanded stringent response model. Mol Microbiol. 2002;45:289–306.
Nikaido H. Molecular basis of bacterial outer membrane permeability revisited. Microbiol Mol Biol Rev. 2003;67:593–656.
Spiess C, Beil A, Ehrmann M. A temperature-dependent switch from chaperone to protease in a widely conserved heat shock protein. Cell. 1999;97:339–47.
Noor R, Murata M, Yamada M. Oxidative stress as a trigger for growth phase-specific sigmaE-dependent cell lysis in Escherichia coli. J Mol Microbiol Biotechnol. 2009;17:177–87.
Condon S. Responses of lactic acid bacteria to oxygen. FEMS Microbiol Rev. 1987;46:269–80.
Zagorski N, Imlay J. The chemistry behind oxidative damage. ASBMB Today. 2009;4:33–6.
Sootsuwan K, Thanonkeo P, Keeratirakha N, Thanonkeo S, Jaisil P, Yamada M. Sorbitol required for cell growth and ethanol production by Zymomonas mobilis under heat, ethanol, and osmotic stresses. Biotechnol Biofuels. 2013;6(1):180
Lui YG, Mitsukawa N, Oosumi T, Whittier R. Efficient isolation and mapping of Arabidopsis Thaliana T-DNA insert junctions by thermal asymmetric interlaced PCR. Plant J. 1995;8:457–63.
Sambrook J, Russell DW. Molecular cloning: a laboratory manual. 3rd ed. Cold Spring Harbour: Cold Spring Harbour Laboratory Press; 2001.
Aiba H, Adhya S, de Cromburgghe B. Evidence for two functional gal promoters in intact Escherichia coli cells. J Biol Chem. 1981;256:11905–10.
Tsunedomi R, Izu H, Kawai T, Matsushita K, Ferenci T, Yamada M. The activator of GntII genes for gluconate metabolism, GntH, exerts negative control of GntR-regulated GntI genes in Escherichia coli. J Bacteriol. 2003;185:1783–95.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Desiniotis A, Kouvelis VN, Davenport K, Bruce D, Detter C, Tapia R, et al. Complete genome sequence of the ethanol-producing Zymomonas mobilis subsp. mobilis centrotype ATCC 29191. J Bacteriol. 2012;194:5966–7.
Conceived and designed the experiments: PT, MM, MY. Performed the experiments: KC, TS, AT, MM. Analyzed the data: KC, TS, AT, MM, TK, PT, MY. Wrote the paper: KC, MM, MY. All authors read and approved the final manuscript.
We thank K. Matsushita, T. Yakushi, W. Soemphol, and N. Lertwattanasakul for their helpful discussion. This work was supported by The Core to Core Program A. Advanced Research Networks, which was granted by the Japan Society for the Promotion of Science, the National Research Council of Thailand, Ministry of Science and Technology in Vietnam, National Univ. of Laos, Univ. of Brawijaya and Beuth Univ. of Applied Science Berlin, and supported by the Program for Promotion of Basic Research Activities for Innovative Biosciences, which was granted by Japan Science and Technology Agency.
The authors declare that they have no competing interests.
Current Funding Sources is 16H02485. Recipient person is Mamoru Yamada, Ph. D.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.