Skip to main content

Metabolic adaptations of Escherichia coli to extended zinc exposure: insights into tricarboxylic acid cycle and trehalose synthesis

Abstract

Balanced bacterial metabolism is essential for cell homeostasis and growth and can be impacted by various stress factors. In particular, bacteria exposed to metals, including the nanoparticle form, can significantly alter their metabolic processes. It is known that the extensive and intensive use of food and feed supplements, including zinc, in human and animal nutrition alters the intestinal microbiota and this may negatively impact the health of the host. This study examines the effects of zinc (zinc oxide and zinc oxide nanoparticles) on key metabolic pathways of Escherichia coli. Transcriptomic and proteomic analyses along with quantification of intermediates of tricarboxylic acid (TCA) were employed to monitor and study the bacterial responses. Multi-omics analysis revealed that extended zinc exposure induced mainly oxidative stress and elevated expression/production of enzymes of carbohydrate metabolism, especially enzymes for synthesis of trehalose. After the zinc withdrawal, E. coli metabolism returned to a baseline state. These findings shed light on the alteration of TCA and on importance of trehalose synthesis in metal-induced stress and its broader implications for bacterial metabolism and defense and consequently for the balance and health of the human and animal microbiome.

Peer Review reports

Introduction

Zinc oxide (ZnO) and zinc oxide nanoparticles (ZnONPs) are inorganic compounds known for their stable chemical and physical properties, easy synthesis, and applications in various fields [1]. They are commonly used in agriculture as pesticides and fertilizers as well as food supplements for livestock animals. Due to their antimicrobial properties, they also show promise in medical applications, such as wound dressings, bandages, and disinfectants [2,3,4]. Recently, ZnONPs have received increased attention due to their effectiveness against antibiotic resistant bacteria and their low ecological footprint [5]. Their antimicrobial activity is attributed to the release of Zn2+ ions and the production of reactive oxygen species (ROS) [6, 7] leading to damage of bacterial cell membranes, leakage of cellular contents, and cell death [8, 9]. However, a long-term exposure of bacteria to ZnO and ZnONPs has been also linked to the development of antimicrobial resistance and virulence [10,11,12,13].

In our previous study, we demonstrated elevated minimum inhibitory concentrations (MICs) to cephalosporins and aminoglycosides in E. coli subjected to the ZnO treatment for 20 subcultures (ZnO20) and this resistance remained for at least 40 subcultures (ZnO40) (15). The highest resistance increase was detected to cefazolin, gentamicin and amikacin with the MICs at the clinical breakpoint. Interestingly, an extended treatment with ZnO nanoparticles (ZnONPs20 and ZnONPs40) only resulted in the higher MIC to chloramphenicol. Furthermore, the elevated biofilm formation was observed under the ZnO40 treatment in comparison to that of control (C40). In contrast, ZnONPs40 treatment resulted in decreased biofilm formation when compared to that of C40. When the zinc treatment was ceased after the 20th subculture and sub-culturing continued until the 40th subculture without zinc (ZnO20 + 20, ZnONPs20 + 20), the antibiotic resistance and virulence phenotypes returned to their original state. Zinc treatments also caused alterations in the bacterial growth and cell morphology [14]. All of these effects are intricately tied to alterations of bacterial metabolism, that can trigger changes in the expression/production of virulence and resistance genes/proteins and vice versa [15,16,17]. Previous reports have indicated that bacteria under various stresses increase production of protection molecules, including carbohydrates (trehalose, glycogen, cellulose) and proteins/enzymes (efflux pumps, DNA-repairing enzymes, heat-shock or cold-shock proteins and more) [18, 19].

It is becoming evident that ZnO and ZnONPs can alter bacterial basal metabolism resulting in changes in bacterial resistance and virulence, yet research on this topic remains limited. In this study, we analyzed major metabolic pathways of E. coli exposed to zinc through multi-omic approaches with the aim to understand how the extended zinc exposure influences gene expression and protein abundance. We also analyzed intermediate compounds in the tricarboxylic acid (TCA) cycle to gain deeper insights into metabolic changes. This study investigates the interplay between zinc and bacterial metabolism and provides new insights for understanding of bacterial responses to environmental stress and the development of antibiotic resistance and virulence mechanisms.

Results

Transcriptomic and proteomic analysis

Differentially expressed genes (DEGs) and differentially abundant proteins (DAPs) in ZnO40, ZnONPs40, ZnO20 + 20, and ZnONPs20 + 20 treatments (Figs. 1 and 2) were studied in comparison to that of untreated control (C40). The study focused on DEGs and DAPs with statistical significance (p ˂ 0.05) and notable changes in expression (log2 fold change > 1.5 or ˂ -1.5). Subsequently, these DEGs and DAPs were further selected based on their gene ontology (GO) terms for carbohydrate, amino acid and fatty acid metabolism. Initially, the number of DEGs (Fig. 1A) and DAPs (Fig. 2A) across the three metabolism categories were analyzed in all treatments. Subsequently, the top 10 GO pathways with the highest rich factor of DEGs (Fig. 1B) and DAPs (Fig. 2B) were identified across all treatments in each metabolism category.

Fig. 1
figure 1

Transcriptome results. (A) Number of significant DEGs in ZnO40, ZnONPs40 ZnO20 + 20, and ZnONPs20 + 20 treatments connected to carbohydrate, amino acid, and fatty acid metabolism according to GO terms. (B) Top 10 GO pathways in ZnO40, ZnONPs40, ZnO20 + 20, and ZnONPs20 + 20 treatments with the highest rich factors connected to (a) carbohydrate metabolism (b) amino acid metabolism (c) fatty acid metabolism

Fig. 2
figure 2

Proteome results. (A) Number of significant DAPs in ZnO40, ZnONPs40, ZnO20 + 20, and ZnONPs20 + 20 treatments connected to carbohydrate, amino acid, and fatty acid metabolism according to GO terms (B) Top 10 GO pathways in ZnO40, ZnONPs40, ZnO20 + 20, and ZnONPs20 + 20 treatments with the highest rich factors connected to (a) carbohydrate metabolism (b) amino acid metabolism (c) fatty acid metabolism

Extended ZnO/ZnONPs treatments mainly altered transcripts in carbohydrate metabolism

Bacterial strains subjected to extended zinc treatments exhibited the most significant changes in the number of DEGs in carbohydrate metabolism (35 DEGs in ZnO40 and 31 DEGs in ZnONPs40). This was followed by changes in amino acid metabolism (20 DEGs in ZnO40 and 13 DEGs in ZnONPs40). In both metabolism categories, these treatments resulted primarily in the up-regulation of transcripts. In contrast, the ZnONPs20 + 20 treatment led mainly to down-regulation of DEGs, especially in metabolism of carbohydrates (21 DEGs), followed by amino acid metabolism (18 DEGs), and fatty acid metabolism (9 DEGs). In all treatments, the number of DEGs involved in fatty acid metabolism was the least affected. Notably, the number of DEGs of ZnO20 + 20 treatment was similar to that of C40 (Fig. 1A).

We then focused on identification of the top 10 enriched GO pathways associated with the altered DEGs across all treatments within each metabolism category (Fig. 1B). The complete list of GO pathway terms, along with its rich factor values for the DEGs across all treatments, is presented in Figure S1. In the carbohydrate metabolism (Fig. 1Ba), the most significantly influenced DEGs pertained to trehalose metabolism; including trehalose-phosphatase activity [GO:0004805], trehalose metabolism in response to cold stress [GO:0070415], and the trehalose biosynthetic process [GO:0005992]. This was particularly notable in strains subjected to ZnO40 and ZnONPs40 treatments. The phosphatase activity of trehalose was also affected in E. coli exposed to the ZnO20 + 20 treatment. Additionally, transcripts involved in glycogen catabolism were impacted under ZnONPs40 and ZnONPs20 + 20 treatments. Furthermore, the biosynthetic processes of mannosylglycerate were altered by ZnONPs40, while mannosyl-3-phosphoglycerate phosphatase activity was affected by ZnO40.

In the amino acid metabolism (Fig. 1Bb), the most affected DEGs were involved in catabolism of arginine [GO:0019545, GO:0019544], L-lysine [GO:0019477] in ZnO40 and ZnONPs40, and catabolism of ornithine [GO:0006593] in ZnO40. Under the ZnONPs20 + 20 treatment, the most impacted DEG were involved in the transport of S-methylmethionine [GO:001586], beta-alanine [GO:0001762], and modification of peptidyl-methionine [GO:0018206]. In contrast, the DEGs involved in fatty acid metabolism (Fig. 1Bc) were the least affected by zinc exposure. In the ZnONPs20 + 20 treatment, 40% of DEGs (rich factor 0.4) (Fig. 1Bc) were altered in the monounsaturated fatty acid biosynthetic process [GO:1903966]. Similarly, in the ZnO40 treatment, 40% of DEGs (rich factor 0.4) were affected in the intermembrane phospholipid transfer [GO:0120010] (Fig. 1Bc).

Overall, among all metabolism categories, the ZnO20 + 20 treatment had the least impact on the gene expression, indicating the bacterial ability to return to its basal metabolic state prior to zinc treatments.

Extended ZnO/ZnONPs treatments mainly altered proteins in carbohydrate and amino acid metabolism

The bacterial strains subjected to extended zinc treatments exhibited the most significant changes in the number of DAPs involved in carbohydrate metabolism, with 24 DAPs in ZnO40 and 25 DAPs in ZnONPs40. The impact on amino acid metabolism was also notable with 25 DAPs in ZnO40 and 22 DAPs in ZnONPs40. In both metabolism categories, the majority of proteins was up-regulated. Across all treatments, the number of DAPs involved in fatty acid metabolism was the least affected. Overall, a higher number of altered DAPs in fatty acid metabolism was observed in strains with the withdrawal of zinc treatments (ZnO20 + 20 and ZnONPs20 + 20) (Fig. 2A).

The study then focused on identification of the top 10 enriched GO pathways associated with affected DAPs across all treatments within each metabolism category (Fig. 2B). A comprehensive list of GO pathway terms, along with rich factor values for the DAPs across all treatments, is in Figure S2. In the carbohydrate metabolism (Fig. 2Ba), the catabolic process of sorbitol [GO:0006062] was the most affected under the ZnONPs40 and ZnONPs20 + 20 treatments. In bacterial strains with extended zinc treatment, there was influence of 50% of DAPS (rich factor 0.5) (Fig. 2Bc) involved in trehalose metabolism; this included trehalose metabolism in response to cold stress [GO:0070415], trehalose biosynthetic process [GO:0005992] and trehalose transport [GO:0015771]. In all treatments, there was influence of plasma membrane pyruvate transport [GO:0006849], fructose metabolic process [GO:0006000] and the catabolic process of 6-sulfoquinovose(1-) [GO:1902777] and 6-sulfoquinovose(1-) to glycerole phosphate and 3-sulfolactaldehyde [GO:0061720].

In amino acid metabolism (Fig. 2Bb), the ZnONPs40 treatment primarily affected DAPs, especially in pathways, such as peptidyl-cysteine S-nitrosylation [GO: 0018119] and peptidyl-cysteine S-trans-nitrosylation [GO:0035606]. Additionally, along with ZnO40 treatment, there was an influence of DAPs in the transport of tyrosine [GO:0015828], p-aminobenzoyl-glutamate [GO:0015814], p-aminobenzoyl-glutamate transmembrane [GO:1902604] and tryptophan [GO:0015827]. This also included peptidyl-S-carbamoyl-L-cysteine dehydration [GO:0046892] and the catabolic processes of tryptophan [GO:0006569] and glutamate [GO:0006538].

Proteins involved in fatty acid metabolism showed the least change. The greatest changes were observed in the fatty acid beta-oxidation using acyl-CoA dehydrogenase [GO:0033539] under ZnO40 treatment across all treatments and in the anaerobic glycerol catabolic process [GO:0019588] under the ZnONPs40 and ZnO20 + 20 treatments.

Integration of transcriptomic and proteomic data

The 3D graphics in Fig. 3A illustrate the common GO pathway terms for DEGs (Fig. 3Aa) and DAPs (Fig. 3Ab), with the corresponding rich factor values for both, transcriptomic (Fig. 3Aa) and proteomic data (Fig. 3Ab) shown on the z-axis (the y-axis numbers for each treatment are in Fig. 3Ac). Among the treatments, the most significant influence of DEGs and DAPs was observed under the ZnO40 and ZnONPs40 treatments, especially within carbohydrate metabolism, focusing on trehalose metabolism pathways: trehalose biosynthetic process [GO:0005992] and trehalose metabolism in response to cold stress [GO:0070415]. In these pathways, the genes encoding trehalose-6-phosphate synthase (otsA) and trehalose-6-phosphate phosphatase (otsB), as well as the OtsB protein, were significantly affected, as evidenced by smaller rich factor values for the proteomic data (Fig. 3B). Up-regulation was observed for the genes otsA and otsB and protein OtsB (Fig. 3B).

Fig. 3
figure 3

Integration of transcriptome and proteome results for amino acid, carbohydrate, and fatty acid metabolism Aa. Transcriptome (TR) and Ab. proteome (PR) results presenting rich factor values for Ac. integrated GO pathways of treated strains ZnO40, ZnONPs40, ZnO20 + 20, and ZnONPs20 + 20 B. DEGs and DAPs for common GO terms related to amino acid (green), carbohydrate (orange), and fatty acid (blue) metabolism in ZnO40, ZnONPs40, ZnO20 + 20, and ZnONPs20 + 20 strains from integrated GO pathways

In amino acid metabolism, the main influence of ZnO40 and ZnONPs40 treatments was detected in the L-lysine catabolic process [GO:0019477], affecting the expression of genes encoding putrescine aminotransferase (ygjG), L-2-hydroxyglutarate dehydrogenase (lhgO), and Glutarate L-2-Hydroxylase (csiD), along with the abundance of the CsiD protein, supporting the smaller rich factor values for the proteomic data. An up-regulation was detected for genes ygjG, lhgO, csiD and the protein CsiD (Fig. 3B). Rich factor values in the other treatments did not reveal significant changes.

The log2 fold change values of all DEGs and DAPs in common GO pathway terms of carbohydrate, amino acid and fatty acid metabolisms for all treatments are presented in the heatmap (Fig. 3B). A summary of GO pathway terms with linked expressed genes, abundant proteins, their accession numbers and log2 fold change values for each treatment are in Table S1. Upon comparing the effects of extended ZnO treatment to its nanoparticle form (ZnONPs), noticeable differences in log2 fold change values emerged for several DEGs (fabB, fadE, lpxC and lpxP) and DAPs (AgaS, BglF, MalT, RpiB, RstB, FadE, and LpxP), all of which play roles in carbohydrate, amino acid, and fatty acid metabolic pathways.

Integration of metabolic, transcriptomic, and proteomic data of TCA cycle intermediates

Figure 4 presents a visual representation of the detected changes in the production of TCA cycle metabolites, gene expression, and protein abundance under various Zn treatments. Color bars indicate the production level of metabolites, with stars denoting significance based on the p-values (* p-value ˂ 0.05, ** p-value ˂ 0.01, *** p-value ˂ 0.001, **** p-value ˂ 0.0001). Additionally, the expression and abundance of enzymes involved in this cycle are depicted in the accompanying tables.

Fig. 4
figure 4

Influence of TCA cycle of E. coli by ZnO/ZnONPs: Production of metabolites related to TCA cycle [ng] for all tested strains are visualized in bar graphs: Control C40 (gray), ZnO40 (violet), ZnONPs40 (pink), ZnO20 + 20 (green), and ZnONPs20 + 20 (blue) with error bars and asterisks indicating the significance according to p-value (* p-value ˂ 0.05, ** p-value ˂ 0.01, *** p-value ˂ 0.001, **** p-value ˂ 0.0001). To compare the treatments (ZnO40, ZnONPs40, ZnO20 + 20, and ZnONPs20 + 20) with the control (C40), an unpaired t-test was performed, and the significance is represented by asterisk symbols. Furthermore, this figure includes log2 fold change values of enzymes involved in this cycle, presented in tables of DEGs on the left and DAPs on the right. Only DEGs or DAPs whose significance was p ˂ 0.05 and log2 fold change > 1.5 and ˂ -1.5, were chosen for depiction. Up-regulated genes/proteins are marked as up arrow, down-regulated genes/proteins are marked as down arrow, hyphens present not changed expression/abundance compared to that of C40, and X marks no identification of DEG/DAP for treatment

Our analysis of TCA cycle intermediates revealed the increased production of pyruvate, citric acid, glutamic acid, ketoglutarate, and succinate in strains subjected to ZnO40 and ZnONPs40 compared to that of control (C40). Notably, E. coli exposed to ZnONPs demonstrated the higher production of pyruvate, citric acid, glutamic acid, and α-ketoglutarate than that under ZnO40. Both extended treatments resulted in a significant decrease in isocitrate, fumarate, and malate. In strains with the zinc withdrawal after 20 subcultures (ZnO20 + 20 and ZnONPs20 + 20), TCA cycle intermediate production remained unchanged relative to C40, with the exception of decrease in isocitrate, citric acid, and lactate under ZnONPs20 + 20.

Figure 4 also illustrates the log2 fold change values of DEGs and DAPs of enzymes involved in the TCA cycle, presented as a heatmap. This map provides a comprehensive view of the gene expression and protein abundance level changes within the TCA cycle, highlighting DEGs and DAPs with significant values (p ˂ 0.05, log2FC > 1.5 and ˂ -1.5). The most notable changes were observed in the expression of genes encoding pyruvate dehydrogenase (poxB) and pyruvate kinase I (pykF) and in the abundance of PoxB in extended zinc treatments. Additionally, the lower abundance of isocitrate dehydrogenase kinase/phosphatase (AceK) was detected in bacteria under the ZnO40 treatment, as well as that in ZnO20 + 20 and ZnONPs20 + 20. A heatmap representing all DEGs and DAPs of the TCA enzymes is provided in Figure S3. When comparing the impact of ZnO to that of ZnONPs on E. coli through changes in TCA cycle intermediate production, only minor differences were detected. Specifically, ZnONPs resulted in slightly higher production of α-ketoglutarate, pyruvate, citric acid and glutamic acid compared to ZnO (Figure S3).

The expression of genes encoding enzymes participating in the TCA cycle is regulated via different transcriptional factors. The major transcriptional factors regulating the TCA genes are listed in Table S2 showing log2 fold changes of expressed genes and abundant proteins with their adjusted p-values.

Discussion

Bacterial metabolism comprises a complex interconnected metabolic network that enables bacteria to process a variety of organic and inorganic compounds essential for the growth and the maintenance of cellular homeostasis [20]. A variety of environmental factors, such as temperature, pH, radiation, nutrient limitation, and oxidative or osmotic stress can lead to the disruption of various metabolic processes [21, 22]. In response to such stressors, especially from heavy metals, bacteria developed a range of adaptive mechanisms. Although these mechanisms are critical for bacterial survival under stress and result to changes in resistance and virulence, they are still not well understood [23, 24].

Understanding the metabolic networks and their regulation is crucial for gaining insights into bacterial physiology and adaptability to changing environmental conditions. Here, we focused on three major metabolic pathways of E. coli exposed to zinc that play a critical role in energy production, biosynthesis, and cellular homeostasis [25,26,27]. Extended zinc treatments (ZnO40, ZnONPs40) of E. coli resulted in alterations in gene expression and protein abundance in metabolism of carbohydrates, proteins, and lipids. Transcriptomic data revealed an increased number of DEGs, involved primarily in carbohydrate metabolism. A similar trend was observed in the proteomic data along with an increased number of DAPs in amino acid metabolism.

It is known that ZnO and ZnONPs induce the production of reactive oxygen species (ROS) in E. coli, such as hydrogen peroxide (H2O2), superoxide anions (O2−.), and hydroxyl radicals (.OH) leading to oxidative stress and potential damage to cellular components, including nucleic acids, amino acids, proteins, and lipids in cell membranes [28, 29]. Genes encoding the stress-protective enzymes (superoxide dismutase and catalase) are regulated by transcription factors through two major oxidative stress regulons, OxyR and SoxRS, as well as the general stress regulon RpoS [30, 31]. RpoS, known for its involvement in responses to multiple stresses [32], directly or indirectly regulates the expression of a large set of genes that aid bacterial adaptation and survival under adverse conditions [33]. This may explain the increased number of DEGs and DAPs linked to carbohydrate and amino acid metabolism in bacteria exposed to ZnO and ZnONPs resulting in stress alleviation. Furthermore, RpoS has been implicated in the regulation of otsA and otsB, that are involved in trehalose biosynthesis [34, 35]. Increased levels of OtsA and OtsB were observed at transcriptome and proteome of extended zinc treatments, corroborated by our previous study, where the up-regulation of otsB was confirmed in both treatments by using real-time qPCR [14]. These combined results indicate the elevated trehalose synthesis.

Trehalose is known to protect membranes and biomacromolecules against abiotic stresses, including osmotic stress, heat stress, desiccation, freezing, and other environmental conditions [36,37,38]. In our study, trehalose was likely synthesized as osmoprotectant to stabilize effects of zinc on proteins and lipids in the cell membranes and to eliminate ROS. Several studies have shown that plants and yeasts can produce trehalose in response to mild heat shock and utilize it to scavenge oxygen radicals and protect cellular proteins from damage [39, 40]. Furthermore, some plants have been found to deal with stress from heavy metals such as cadmium, zinc, and nickel through trehalose production [41, 42]. Exogenously added trehalose has also been shown to enhance stress tolerance and yield in various plant species, indicating its potential as a protective agent against metal-induced damage [43]. However, the role of trehalose in bacteria under stress from heavy metals has not been explored. One study demonstrated that Rhizobium sp. have the capacity to accumulate trehalose and this may play an important role in mediating stress responses in interactions with plants [44].

An extended treatment of E. coli with ZnO nanoparticles (ZnONPs40) led to increased catabolism of glycogen, D-allose, and sorbitol, while ZnO40 treatment enhanced the catabolism of cellobiose. Subsequently, glycolysis likely became a great source of pyruvate for the TCA cycle, generating a substantial amount of energy. The heightened production of pyruvate could be also attributed to the catabolism of 6-sulfoquinovose and the activity of deoxyribose phosphate that were notably elevated under extended ZnO/ZnONPs treatments. Furthermore, enzymes involved in the pentose phosphate shunt pathway were elevated and likely contributed to the production of pyruvate. Consequently, more pyruvate can enter the TCA cycle and this is corroborated by elevated intermediates in the TCA cycle in cells under ZnO40 and ZnONPs40 treatments. The TCA pathway plays a pivotal role in the catabolism of organic molecules in the presence of oxygen to harvest the energy for cell growth and division [45]. Increased production of the enzyme PoxB facilitates the conversion of pyruvate into acetyl-CoA. Extended zinc treatments likely augmented the catabolism of L-lysine, leading to a higher production of acetyl-CoA. This potentially contributed to an elevated abundance of citric acid in E. coli after the long-term exposure to ZnO and ZnONPs.

The most notable change in our analysis was in the production of α-ketoglutarate. Under ZnO40, the gene aceK coding for the protein responsible for activating isocitrate dehydrogenase [46] was up-regulated, potentially leading to an increased production of α-ketoglutarate. It has been demonstrated that high concentrations of zinc can alter the activity of isocitrate dehydrogenase by phosphorylation process resulting in its inactivity [47], which could lead to the synthesis of α-ketoglutarate from the L-glutamine amino acid pathway. Higher levels of NADH and GTP in strains with extended zinc treatments were likely generated by the conversion of α-ketoglutarate to succinyl-CoA and its subsequent conversion to succinate. Succinate levels were also elevated possibly due to its conversion from isocitrate, facilitated by the up-regulation of aceA and aceB in the glyoxylate pathway, as observed in ZnO40 treatment. Additionally, an increased production of succinate in the ZnONPs40 treatment could have been caused by the down-regulation of the sdhB and sdhD genes, coding for succinate: quinone oxidoreductase, responsible for the conversion of succinate to fumarate. Since the production of fumarate and malate decreased in these treatments, succinate was likely redirected to different pathways, such as the propionate synthesis pathway [48]. The oxidation of this short-chain fatty acid is then used as a source of carbons or energy [49]. High levels of intermediate compounds in the TCA cycle of E. coli can provide a greater amount of energy-rich molecules such as GTP and NADH, which can be converted to ATP to support anabolic processes, leading to cell protection, growth, and proliferation [50]. The elevated levels of TCA intermediates indicate a higher production of energy-rich molecules in cells under an extended zinc exposure.

Overall, the transcriptomic and the proteomic data suggest that cells under ZnO40 and ZnONPs treatments divert energy from catabolic pathways to gluconeogenesis. This process likely leads to the synthesis of trehalose and other carbohydrates. The role of the TCA cycle in providing ATP for trehalose synthesis has been shown previously [50]. Consequently, the elevated trehalose might inactivate the oxygen radicals and acts as a stress protectant for lipids and proteins in the cell periplasm and the cytoplasmic membrane, which could be otherwise damaged [51]. A recent study revealed that in E. coli, trehalose functions as a chemical chaperone, stabilizing denatured proteins and facilitates their refolding. Trehalose may also protect proteins against aggregation by acting as a metabolite that indirectly counteracts harmful protein acetylation [52]. When the stress is relieved, E. coli can rapidly degrade trehalose and revert to basic metabolism. Degradation of trehalose involves its transfer from the periplasm to the cytoplasm, a process facilitated by the TreB protein [53], which was indeed elevated in cells under ZnO40 and ZnONPs40. After transport, trehalose is phosphorylated and degraded by trehalases TreC and TreF, into two glucose molecules, whose expression/abundance was not changed in comparison to that of control [53, 54].

It is very likely that under extended ZnO and ZnONPs treatments, E. coli experienced oxidative stress, which triggered the up-regulation of the otsAB system and the subsequent trehalose synthesis. This response provided protection against the damaging effects of zinc ions and allowed bacteria to adapt and survive. The higher energy content in cells exposed to ZnONPs40 was likely utilized for the increased trehalose production in comparison to that of ZnO40. Upon reversal ZnO and ZnONPs treatments, E. coli returned to a metabolic state similar to that of control C40. Future studies need to focus on detection and measurement of trehalose levels as well as to generate mutant strains with inactivated trehalose synthesis to confirm our findings from multi-omics approach.

In conclusion, extended zinc treatments, in the form of ZnO and ZnONPs, greatly impacted the carbohydrate metabolic network of E. coli. Transcriptomic and proteomic analyses revealed elevated expression/abundance of enzymes involved in glycolysis and the TCA cycle which was further confirmed by analysis of TCA intermediate compounds. These metabolic shifts are indicative of the bacterial response to stress as it generates the necessary energy for protection and adaptation. Importantly, the restoration of the metabolic state after withdrawal of zinc points out the bacterial capacity to adapt and recover. Our findings provide new and comprehensive insights into the metabolic adaptations of E. coli to extended zinc exposure and indicate the role of trehalose in protection from oxidative stress induced by a heavy metal and its nanoparticles. This is the first study indicating the role of TCA cycle and trehalose in protection from stress induced by a heavy metal. Considering the extent of zinc use in human and animal nutrition, our understanding the molecular mechanisms of bacterial responses is essential for elucidating the impact of metal exposure on bacterial metabolism and for developing strategies for maintenance of balanced intestinal microbiomes and combating bacterial resistance and virulence.

Materials and methods

Sub-culturing of E. Coli strain with zinc oxide and zinc oxide nanoparticles and its processing

In our previous study, we performed the sub-culturing procedure of E. coli ATCC 25,922 strain with sub-inhibitory doses of ZnO and ZnONPs in heptaplicate (n = 7) for each strain [14]. Briefly, bacteria E. coli ATCC 25,922 were cultivated on 5% Columbia blood agar (pH 7.3 ± 0.2; Lab Media Servis, Czech Republic) at 37 °C. Overnight cultures were diluted in double-concentrated Mueller-Hinton (MH) broth (pH 7.3 ± 0.1; Sigma Aldrich, USA) to an optical density at 600 nm (OD600) of 0.08–0.13 AU and then further diluted 100×. It was then incubated with ZnO/ZnONPs in 96 well microplates at 37 °C for 24 h with shaking at 120 rpm. Next day, the absorbance (620 nm) was measured using the MultiScan EX Microplate Photometer (Thermo Fisher, Germany) and the minimum inhibitory concentrations and sub-inhibitory concentrations were determined. A loopful (10 µl) of cells from the well containing the first concentration below the MIC (sub-inhibitory concentration) was transferred onto the surface of 5% Columbia blood agar. This entire process was repeated 19 times, resulting in 20 sub-cultures across seven independent replicates, with MIC measurements taken. After reaching the 20th sub-culture, the treated bacteria were subdivided into two lines and the experiment continued for additional 20 cycles with ZnO/ZnONPs (ZnO40, ZnONPs40) or mock (ZnO20 + 20, ZnONPs20 + 20). In parallel, control bacteria were incubated under the same conditions using PBS. Final sub-cultures reaching 40 passages (C40, ZnO40, ZnONPs40, ZnO20 + 20, ZnONPs20 + 20) were preserved through freezing after 24 h incubation on 5% Columbia blood agar. Specifically, bacteria were centrifuged, and the resulting pellet was washed three times with PBS. Each bacterial pellet was stored in cryo-medium at – 80 °C until further use for transcriptomic, proteomic approach and quantitative analysis of TCA cycle intermediates. This process followed the protocols detailed in our previous study [14]. Subsequently, bacterial pellets were utilized for isolation of RNA and proteins, as outlined in our earlier work. Furthermore, in the current study, metabolites were also extracted from the bacterial pellets.

Extraction of metabolites and analysis of TCA cycle intermediates

Frozen bacterial pellets were resuspended in a cold solution containing 0.1% formic acid, 40% methanol, and 40% acetonitrile (v/v) in pentaplicate (n = 5) for each strain. Metabolites were isolated via the extraction protocol for TCA cycle products, as outlined by Carvalho et al. (2019) [55] using glass beads (Benchmark Scientific, Sayreville, NJ, USA). Prior to extraction, bacterial pellets were pre-homogenized (4500 rpm, 60s) using the MagnaLyzer and subsequently placed on ice for 30 s. The supernatant was then stored at − 80 °C until further analysis. This process was repeated four times, adding the cold solution to the remaining pellets, resulting in four frozen supernatants for each treatment. The samples were then centrifuged for 15 min at 14,000 rpm and 4 °C. Subsequently, 500 µl of the supernatant was vacuum dried and reconstituted in 200 µl of the mobile phase (1% methanol acidified with 0.2% formic acid). Then, 5.0 µl of the samples were analyzed by the LC-MS using reversed-phase chromatography. The analysis was conducted on the high-performance chromatographic system Agilent Technologies 1260 (Waldbronn, Germany) coupled to the LTQ VelosPro Orbitrap Elite hybrid mass spectrometer (Thermo Scientific, CA, USA). Chromatographic separation was achieved suing a C18 Kinetex column (150 × 2.1 mm with 5 μm particles) with a flow rate of 0.2 ml/min. The gradient of mobile phases A (0.2% formic acid in water) and B (methanol) was applied as follows: 1% of B at 1 min, 5% at 4 min, 20% at 6 min, 100% at 8 min. The spectra were evaluated using Xcalibur software v2.2 (Thermo Scientific, CA, USA). A constant amount of 13C-labeled L-lactic acid (Sigma Aldrich) was used as an internal standard, with 1 µl of the internal standard solution mixed with 99 µl of the sample/standard solution.

Transcriptome and proteome analysis

The extraction and processing of biological materials (RNA and proteins) from bacterial pellets to obtain transcriptomic and proteomic data, were described in detail in our previous study [14]. Differential expression analysis was always performed from zinc treatment (ZnO40, ZnONPs40, ZnO20 + 20 and ZnONPs20 + 20) compared to control (C40). In this study we worked only with those DEGs and DAPs whose p-value ˂ 0.05. For transcriptome, adjusted p-value threshold (˂ 0.05) was calculated through the false discovery rate (FDR). This was performed in Geneious Prime 2021.0.1. by DESeq2 package [56].

Gene ontology and rich factor calculation

The significantly regulated genes/proteins were subjected to gene enrichment analysis using the aGOtool within the UniProt keyword classification system for each treatment [57]. Generated GO terms for biological processes were then organized into appropriate metabolic categories, including carbohydrate, amino acid, and fatty acid metabolism. For each DEGs and DAPs, a rich factor was calculated. The rich factor represents the ratio of the number of genes/proteins in a specific pathway to all genes/proteins of that specific pathway. A higher rich factor indicates a greater degree of enrichment within the pathway. The total number of genes/proteins participating in each specific pathway was obtained from the EcoCyc database: https://biocyc.org/ECOLI/class-tree?object=Gene-Ontology-Terms.

Statistical analysis and data visualization

Metabolic data of TCA intermediates were processed, standard deviations were calculated, and significance in comparison to that of control strain was determined using an unpaired t-test in GraphPad Prism 8.0.1. (GraphPad Software, CA, USA). In case of transcriptome and proteome, the p-value was calculated by background-based t-test as mentioned in our previous study [14]. All graphics for the transcriptomic, proteomic and metabolomic data were generated using the ggplot2 in RStudio [58], while 3D graphics were visualized using the matplotlib in Python [59].

Data availability

Data that support the findings of this study are presented in the main article and Supplementary Information files. Raw data generated from the RNA-Seq experiment of E. coli strains were deposited in the SRA archive of NCBI database under the BioProject: PRJNA915319 (https://www.ncbi.nlm.nih.gov/sra/PRJNA915319). The mass spectrometry proteomics data were deposited to the ProteomeXchange Consortium via the PRIDE (10.1093/nar/gky1106) partner repository with the data set identifier PXD027925.

Abbreviations

DAP:

differential abundance of protein

DEG:

differential expression of gene

ROS:

reactive oxygen species

TCA:

tricarboxylic acid

ZnO:

Zinc oxide

ZnONPs:

Zinc oxide nanoparticles

ZnO40:

Escherichia coli sub–cultured with sub–inhibitory doses of ZnO till 40th subculture

ZnONPs40:

Escherichia coli sub–cultured with sub–inhibitory doses of ZnONPs till 40th subculture

ZnO20 + 20:

Escherichia coli sub–cultured with sub–inhibitory doses of ZnO till 20th subculture and then subsequently sub–cultured without ZnO until 40th subculture

ZnONPs20 + 20:

Escherichia coli sub–cultured with sub–inhibitory doses of ZnONPs till 20th subculture and then subsequently sub–cultured without ZnONPs until 40th subculture

References

  1. Quadri TW, et al. Zinc oxide nanocomposites of selected polymers: synthesis, characterization, and corrosion inhibition studies on mild steel in HCl solution. ACS Omega. 2017;2(11):8421–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. García-Gómez C, et al. Study of Zn availability, uptake, and effects on earthworms of zinc oxide nanoparticle versus bulk applied to two agricultural soils: acidic and calcareous. Chemosphere. 2020;239:124814.

    Article  PubMed  Google Scholar 

  3. Mendes CR, et al. Antibacterial action and target mechanisms of zinc oxide nanoparticles against bacterial pathogens. Sci Rep. 2022;12(1):2658.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Shi T, et al. Inventories of heavy metal inputs and outputs to and from agricultural soils: a review. Ecotoxicol Environ Saf. 2018;164:118–24.

    Article  CAS  PubMed  Google Scholar 

  5. Additives EPo, Feed A. Scientific opinion on the potential reduction of the currently authorised maximum zinc content in complete feed. Efsa J. 2014;12(5):3668.

    Google Scholar 

  6. Espitia PJP, et al. Zinc oxide nanoparticles: synthesis, antimicrobial activity and food packaging applications. Food Bioprocess Technol. 2012;5:1447–64.

    Article  CAS  Google Scholar 

  7. Sawai J. Quantitative evaluation of antibacterial activities of metallic oxide powders (ZnO, MgO and CaO) by conductimetric assay. J Microbiol Methods. 2003;54(2):177–82.

    Article  CAS  PubMed  Google Scholar 

  8. Jones N, et al. Antibacterial activity of ZnO nanoparticle suspensions on a broad spectrum of microorganisms. FEMS Microbiol Lett. 2008;279(1):71–6.

    Article  CAS  PubMed  Google Scholar 

  9. Sinha R, et al. Interaction and nanotoxic effect of ZnO and ag nanoparticles on mesophilic and halophilic bacterial cells. Bioresour Technol. 2011;102(2):1516–20.

    Article  CAS  PubMed  Google Scholar 

  10. Brown LR et al. Increased zinc availability enhances initial aggregation and biofilm formation of Streptococcus pneumoniae. Front Cell Infect Microbiol, 2017: p. 233.

  11. Cui H, Smith AL. Impact of engineered nanoparticles on the fate of antibiotic resistance genes in wastewater and receiving environments: a comprehensive review. Environ Res. 2022;204:112373.

    Article  CAS  PubMed  Google Scholar 

  12. Kociova S, et al. Zinc phosphate-based nanoparticles as alternatives to zinc oxide in diet of weaned piglets. J Anim Sci Biotechnol. 2020;11(1):1–16.

    Article  Google Scholar 

  13. Wu T, et al. Zinc exposure promotes commensal-to-pathogen transition in Pseudomonas aeruginosa leading to mucosal inflammation and illness in mice. Int J Mol Sci. 2021;22(24):13321.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Rihacek M et al. Zinc effects on bacteria: insights from < i > Escherichia coli by multi-omics approach. mSystems. 0(0): p. e00733–23.

  15. Martínez JL. Bacterial pathogens: from natural ecosystems to human hosts. Environ Microbiol. 2013;15(2):325–33.

    Article  PubMed  Google Scholar 

  16. Martínez JL, Rojo F. Metabolic regulation of antibiotic resistance. FEMS Microbiol Rev. 2011;35(5):768–89.

    Article  PubMed  Google Scholar 

  17. Schroeder M, Brooks BD, Brooks AE. The complex relationship between virulence and antibiotic resistance. Genes. 2017;8(1):39.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Khan R, et al. Bacterial polysaccharides—A big source for prebiotics and therapeutics. Front Nutr. 2022;9:1031935.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Święciło A, Zych-Wężyk I. Bacterial stress response as an adaptation to life in a soil environment. Pol J Environ Stud, 2013. 22(6).

  20. Jurtshuk P Jr. Bacterial metabolism. 2011.

  21. Ron EZ. Bacterial stress response. The prokaryotes, 2006: pp. 1012–1027.

  22. Wood JM. Osmosensing by bacteria. Science’s STKE. 2006;2006(357):pe43–43.

    PubMed  Google Scholar 

  23. Belloch C, et al. Fermentative stress adaptation of hybrids within the Saccharomyces sensu stricto complex. Int J Food Microbiol. 2008;122(1–2):188–95.

    Article  CAS  PubMed  Google Scholar 

  24. Valls M, De Lorenzo V. Exploiting the genetic and biochemical capacities of bacteria for the remediation of heavy metal pollution. FEMS Microbiol Rev. 2002;26(4):327–38.

    Article  CAS  PubMed  Google Scholar 

  25. Bender DA. Amino acid metabolism. Wiley; 2012.

  26. Fraenkel D, Vinopal R. Carbohydrate metabolism in bacteria. Annual Reviews Microbiol. 1973;27(1):69–100.

    Article  CAS  Google Scholar 

  27. Fulco AJ. Fatty acid metabolism in bacteria. Prog Lipid Res. 1983;22(2):133–60.

    Article  CAS  PubMed  Google Scholar 

  28. Imlay JA. The molecular mechanisms and physiological consequences of oxidative stress: lessons from a model bacterium. Nat Rev Microbiol. 2013;11(7):443–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Matuła K, et al. Phenotypic plasticity of Escherichia coli upon exposure to physical stress induced by ZnO nanorods. Sci Rep. 2019;9(1):8575.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Seo S, Kim D, Szubin R. Palsson, and Bernhard, O.(2015). Genome-wide reconstruction of OxyR and SoxRS transcriptional regulatory networks under oxidative stress in Escherichia coli K-12 MG1655. Cell Rep. 12(8): pp. 1289–1299.

  31. Xiao X, Wu Z-C, Chou K-C. A multi-label classifier for predicting the subcellular localization of gram-negative bacterial proteins with both single and multiple sites. PLoS ONE. 2011;6(6):e20592.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Tramonti A, De Canio M, De Biase D. GadX/GadW-dependent regulation of the Escherichia coli acid fitness island: transcriptional control at the gady–gadw divergent promoters and identification of four novel 42 bp GadX/GadW‐specific binding sites. Mol Microbiol. 2008;70(4):965–82.

    Article  CAS  PubMed  Google Scholar 

  33. Weber H, et al. Genome-wide analysis of the general stress response network in Escherichia coli: σS-dependent genes, promoters, and sigma factor selectivity. J Bacteriol. 2005;187(5):1591–603.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Gibson RP, et al. Insights into trehalose synthesis provided by the structure of the retaining glucosyltransferase OtsA. Chem Biol. 2002;9(12):1337–46.

    Article  CAS  PubMed  Google Scholar 

  35. Hengge R. Proteolysis of σS (RpoS) and the general stress response in Escherichia coli. Res Microbiol. 2009;160(9):667–76.

    Article  CAS  PubMed  Google Scholar 

  36. De Virgilio C, et al. The role of trehalose synthesis for the acquisition of thermotolerance in yeast: I. Genetic evidence that trehalose is a thermoprotectant. Eur J Biochem. 1994;219(1–2):179–86.

    Article  PubMed  Google Scholar 

  37. Hengge-Aronis R, et al. Trehalose synthesis genes are controlled by the putative sigma factor encoded by rpoS and are involved in stationary-phase thermotolerance in Escherichia coli. J Bacteriol. 1991;173(24):7918–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Van Laere A. Trehalose, reserve and/or stress metabolite? FEMS Microbiol Lett. 1989;63(3):201–9.

    Article  Google Scholar 

  39. Benaroudj N, Goldberg AL. Trehalose accumulation during cellular stress protects cells and cellular proteins from damage by oxygen radicals. J Biol Chem. 2001;276(26):24261–7.

    Article  CAS  PubMed  Google Scholar 

  40. Luo Y, Li W-M, Wang W. Trehalose: protector of antioxidant enzymes or reactive oxygen species scavenger under heat stress? Environ Exp Bot. 2008;63(1–3):378–84.

    Article  CAS  Google Scholar 

  41. Garg N, Saroy K. Interactive effects of polyamines and arbuscular mycorrhiza in modulating plant biomass, N2 fixation, ureide, and trehalose metabolism in Cajanus cajan (L.) Millsp. Genotypes under nickel stress. Environ Sci Pollut Res. 2020;27(3):3043–64.

    Article  CAS  Google Scholar 

  42. Garg N, Singh S. Mycorrhizal inoculations and silicon fortifications improve rhizobial symbiosis, antioxidant defense, trehalose turnover in pigeon pea genotypes under cadmium and zinc stress. Plant Growth Regul. 2018;86(1):105–19.

    Article  CAS  Google Scholar 

  43. Rehman S, et al. Exogenously applied trehalose augments cadmium stress tolerance and yield of mung bean (Vigna radiata L.) grown in soil and hydroponic systems through reducing cd uptake and enhancing photosynthetic efficiency and antioxidant defense systems. Plants. 2022;11(6):822.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Streeter JG. Accumulation of alpha, alpha-trehalose by Rhizobium bacteria and bacteroids. J Bacteriol. 1985;164(1):78–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Meyer FM, et al. Physical interactions between tricarboxylic acid cycle enzymes in Bacillus subtilis: evidence for a metabolon. Metab Eng. 2011;13(1):18–27.

    Article  CAS  PubMed  Google Scholar 

  46. Cozzone AJ. Regulation of acetate metabolism by protein phosphorylation in enteric bacteria. Annu Rev Microbiol. 1998;52(1):127–64.

    Article  CAS  PubMed  Google Scholar 

  47. Buchachenko AL, et al. A specific role of magnetic isotopes in biological and ecological systems. Physics and biophysics beyond. Prog Biophys Mol Biol. 2020;155:1–19.

    Article  CAS  PubMed  Google Scholar 

  48. Louis P, Flint HJ. Formation of propionate and butyrate by the human colonic microbiota. Environ Microbiol. 2017;19(1):29–41.

    Article  CAS  PubMed  Google Scholar 

  49. Textor S, et al. Propionate oxidation in Escherichia coli: evidence for operation of a methylcitrate cycle in bacteria. Arch Microbiol. 1997;168:428–36.

    Article  CAS  PubMed  Google Scholar 

  50. Jurtshuk P. Bacterial metabolism. Medical microbiology, 1996. 4.

  51. Singer MA, Lindquist S. Multiple effects of trehalose on protein folding in vitro and in vivo. Mol Cell. 1998;1(5):639–48.

    Article  CAS  PubMed  Google Scholar 

  52. Moruno Algara M, et al. Trehalose protects Escherichia coli against carbon stress manifested by protein acetylation and aggregation. Mol Microbiol. 2019;112(3):866–80.

    Article  CAS  PubMed  Google Scholar 

  53. Steen JA, et al. The trehalose phosphotransferase system (PTS) in E. Coli W can transport low levels of sucrose that are sufficient to facilitate induction of the csc sucrose catabolism operon. PLoS ONE. 2014;9(2):e88688.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Rimmele M, Boos W. Trehalose-6-phosphate hydrolase of Escherichia coli. J Bacteriol. 1994;176(18):5654–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Carvalho SM, et al. Metabolomics of Escherichia coli treated with the antimicrobial carbon monoxide-releasing molecule CORM-3 reveals tricarboxylic acid cycle as major target. Antimicrob Agents Chemother. 2019;63(10). https://doi.org/10.1128/aac. 00643 – 19.

  56. Anders S, Huber W. Differential expression analysis for sequence count data. Nat Precedings, 2010: p. 1–1.

  57. Consortium U. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2018;46(5):2699.

    Article  Google Scholar 

  58. Team RC. R: A language and environment for statistical computing. 2022.

  59. Van Rossum G, Drake FL. Python reference manual. Volume 111. Centrum voor Wiskunde en Informatica Amsterdam; 1995.

Download references

Acknowledgements

The study was financially supported by ERDF “Multidisciplinary research to increase application potential of nanomaterials in agricultural practice” (No. CZ.03.1.01/0.0/0.0/16_025/0007314).

Funding

ERDF “Multidisciplinary research to increase application potential of nanomaterials in agricultural practice” (No. CZ.03.1.01/0.0/0.0/16_025/0007314).

Author information

Authors and Affiliations

Authors

Contributions

All authors made a significant contribution to the reported work. M.R. processed transcriptomic and metabolomic data, created graphics, and wrote the manuscript. L.K., T.R. and D.S.S. processed the samples and analyzed the metabolic data. T.F. analyzed proteomic data. V.A. assisted with the manuscript revisions. K.C. and L.Z. designed the project, supervised the experiments, and co-wrote the manuscript. The final text has been reviewed and approved by all authors.

Corresponding author

Correspondence to Kristyna Cihalova.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rihacek, M., Kosaristanova, L., Fialova, T. et al. Metabolic adaptations of Escherichia coli to extended zinc exposure: insights into tricarboxylic acid cycle and trehalose synthesis. BMC Microbiol 24, 384 (2024). https://doi.org/10.1186/s12866-024-03463-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12866-024-03463-6

Keywords