A decay effect of the growth rate associated with genome reduction in Escherichia coli

Bacterial genome reduction for optimal chassis of synthetic biology: a review

Bacteria with streamlined genomes, that harbor full functional genes for essential metabolic networks, are able to synthesize the desired products more effectively and thus have advantages as production platforms in industrial applications. To obtain streamlined chassis genomes, a large amount of effort has been made to reduce existing bacterial genomes. This work falls into two categories: rational and random reduction. The identification of essential gene sets and the emergence of various genome-deletion techniques have greatly promoted genome reduction in many bacteria over the past few decades. Some of the constructed genomes possessed desirable properties for industrial applications, such as: increased genome stability, transformation capacity, cell growth, and biomaterial productivity. The decreased growth and perturbations in physiological phenotype of some genome-reduced strains may limit their applications as optimized cell factories. This review presents an assessment of the advancements made to date in bacterial genome reduction to construct optimal chassis for synthetic biology, including: the identification of essential gene sets, the genome-deletion techniques, the properties and industrial applications of artificially streamlined genomes, the obstacles encountered in constructing reduced genomes, and the future perspectives.

Repurposing fusidic acid as an antimicrobial against enterococci with a low probability of resistance development

Drug repurposing constitutes a strategy to combat antimicrobial resistance, by using agents with known safety, pharmacokinetics, and pharmacodynamics. Previous studies have implemented new fusidic acid (FA) front-loading-dose regimens, allowing higher serum levels than those achievable with ordinary doses. As susceptibility breakpoints are affected by serum level, we evaluated the repurposing of FA as an antimicrobial product against enterococci. FA minimum inhibitory concentrations (MICs) against standard enterococci strains; Enterococcus faecalis ATCC 29212 and Enterococcus faecium ATCC 27270 were 2 and 4 µg/mL, respectively. The MIC against 98 enterococcal clinical isolates was ≤ 8 µg/mL; all would be susceptible if categorized according to recalculated breakpoints (≥ 16 µg/mL), based on the serum level achieved using the front-loading regimen. FA administration in vivo, using the BALB/c mouse infection model, significantly reduced bacterial burden by two to three log10 units in the liver and spleen of mice infected with vancomycin-susceptible and -resistant strains. Exposure of the standard enterococcal strains to increasing, but not fixed, FA concentrations resulted in resistant strains (MIC = 128 µg/mL), with thicker cell walls and slower growth rates. Only one mutation (M651I) was detected in the fusA gene of the resistant strain derived from serial passage of E. faecium ATCC 27270, which was retained in the revertant strain after passage in the FA-free medium. In conclusion, FA can be repurposed as an antimicrobial drug against enterococci with a low probability of mutational resistance development, and can be employed for treatment of infections attributable to vancomycin-resistant enterococci.

Multiplex gene knockout raises Ala-Gln production by Escherichia coli expressing amino acid ester acyltransferase

L-Alanyl-L-Glutamine (Ala-Gln) is a common parenteral nutritional supplement. In our previous study, the recombinant whole-cell catalyst Escherichia coli BL21(DE3) overexpressing α-amino acid ester acyltransferase (BPA) to produce Ala-Gln has high activity and has been applied to large-scale production experiments. However, the degradation of Ala-Gln is detected under prolonged incubation, and endogenous broad-spectrum dipeptidase may be the primary cause. In this study, a CRISPR-Cas9 method was used to target pepA, pepB, pepD, pepN, dpp, and dtp to knock out one or more target genes. The deletion combination was optimized, and a triple knockout strain BL21(DE3)-ΔpepADN was constructed. The degradation performance of the knockout chassis was measured, and the results showed that the degradation rate of Ala-Gln was alleviated by 48% compared with the control. On this basis, B_pADNPA (BPA-ΔpepADN) was built, and the production of Ala-Gln was 129% of the BPA's accumulation, proving that the ΔpepADN knockout conducive to the accumulation of dipeptide. This study will push forward the industrialization process of Ala-Gln production by whole-cell catalyst Escherichia coli expressing α-amino acid ester acyltransferase. KEY POINTS: • Endogenous dipeptidase knockout alleviates the degradation of Ala-Gln by the chassis • The balanced gene knockout combination is pepA, pepD, and pepN • The accumulation of Ala-Gln with B_pADNPA was 129% of the control.

A unified and simple medium for growing model methanogens

Apart from their archetypic use in anaerobic digestion (AD) methanogenic archaea are targeted for a wide range of applications. Using different methanogenic archaea for one specific application requires the optimization of culture media to enable the growth of different strains under identical environmental conditions, e.g., in microbial electrochemical technologies (MET) for (bio)electromethanation. Here we present a new culture medium (BFS01) adapted from the DSM-120 medium by omitting resazurin, yeast extract, casitone, and using a low salt concentration, that was optimized for Methanosarcina barkeri, Methanobacterium formicicum, and Methanothrix soehngenii. The aim was to provide a medium for follow-up co-culture studies using specific methanogens and Geobacter spp. dominated biofilm anodes. All three methanogens showed growth and activity in the BFS01 medium. This was demonstrated by estimating the specific growth rates ( μ ) and doubling times ( t d ) of each methanogen. The μ and t d based on methane accumulation in the headspace showed values consistent with literature values for M. barkeri and M. soehngenii. However, μ and t d based on methane accumulation in the headspace differed from literature data for M. formicicum but still allowed sufficient growth. The lowered salt concentration and the omission of chemically complex organic components in the medium may have led to the observed deviation from μ and t d for M. formicicum as well as the changed morphology. 16S rRNA gene-based amplicon sequencing and whole genome nanopore sequencing further confirmed purity and species identity.

Machine learning-assisted medium optimization revealed the discriminated strategies for improved production of the foreign and native metabolites

The composition of medium components is crucial for achieving the best performance of synthetic construction in genetically engineered cells. Which and how medium components determine the performance, e.g., productivity, remain poorly investigated. To address the questions, a comparative survey with two genetically engineered Escherichia coli strains was performed. As a case study, the strains carried the synthetic pathways for producing the aromatic compounds of 4-aminophenylalanine (4APhe) or tyrosine (Tyr), common in the upstream but differentiated in the downstream metabolism. Bacterial growth and compound production were examined in hundreds of medium combinations that comprised 48 pure chemicals. The resultant data sets linking the medium composition to bacterial growth and production were subjected to machine learning for improved production. Intriguingly, the primary medium components determining the production of 4PheA and Tyr were differentiated, which were the initial resource (glucose) of the synthetic pathway and the inducer (IPTG) of the synthetic construction, respectively. Fine-tuning of the primary component significantly increased the yields of 4APhe and Tyr, indicating that a single component could be crucial for the performance of synthetic construction. Transcriptome analysis observed the local and global changes in gene expression for improved production of 4APhe and Tyr, respectively, revealing divergent metabolic strategies for producing the foreign and native metabolites. The study demonstrated that ML-assisted medium optimization could provide a novel point of view on how to make the synthetic construction meet the designed working principle and achieve the expected biological function.

Random genome reduction coupled with polyhydroxybutyrate biosynthesis to facilitate its accumulation in Escherichia coli

Genome reduction has been emerged as a powerful tool to construct ideal chassis for synthetic biology. Random genome reduction couple genomic deletion with growth and has the potential to construct optimum genome for a given environment. Recently, we developed a transposon-mediated random deletion (TMRD) method that allows the random and continuous reduction of Escherichia coli genome. Here, to prove its ability in constructing optimal cell factories, we coupled polyhydroxybutyrate (PHB) accumulation with random genome reduction and proceeded to reduce the E. coli genome. Five mutants showed high biomass and PHB yields were selected from 18 candidates after ten rounds of genome reduction. And eight or nine genomic fragments (totally 230.1–270.0 Kb) were deleted in their genomes, encompassing 4.95%–5.82% of the parental MG1655 genome. Most mutants displayed better growth, glucose utilization, protein expression, and significant increase of electroporation efficiency compared with MG1655. The PHB content and concentration enhanced up to 13.3%–37.2% and 60.2%–102.9% when batch fermentation was performed in M9-glucose medium using the five mutants. Particularly, in mutant H16, lacking 5.28% of its genome, the increase of biomass and PHB concentration were more than 50% and 100% compared with MG1655, respectively. This work expands the strategy for creating streamlined chassis to improve the production of high value-added products.

Contribution of the genomic and nutritional differentiation to the spatial distribution of bacterial colonies

Colony growth is a common phenomenon of structured populations dispersed in nature; nevertheless, studies on the spatial distribution of colonies are largely insufficient. Here, we performed a systematic survey to address the questions of whether and how the spatial distribution of colonies was influenced by the genome and environment. Six Escherichia coli strains carrying either the wild-type or reduced genomes and eight media of varied nutritional richness were used to evaluate the genomic and environmental impacts, respectively. The genome size and nutritional variation contributed to the mean size and total area but not the variation and shape of size distribution of the colonies formed within the identical space and of equivalent spatial density. The spatial analysis by means of the Voronoi diagram found that the Voronoi correlation remained nearly constant in common, in comparison to the Voronoi response decreasing in correlation to genome reduction and nutritional enrichment. Growth analysis at the single colony level revealed positive correlations of the relative growth rate to both the maximal colony size and the Voronoi area, regardless of the genomic and nutritional variety. This result indicated fast growth for the large space assigned and supported homeostasis in the Voronoi correlation. Taken together, the spatial distribution of colonies might benefit efficient clonal growth. Although the mechanisms remain unclear, the findings provide quantitative insights into the genomic and environmental contributions to the growth and distribution of spatially or geographically isolated populations.

Investigating the effect of nanoparticle on phenanthrene biodegradation by Labedella gwakjiensis strain KDI

Polycyclic aromatic hydrocarbons (PAHs), as persistent organic contaminants, are a major source of concern due to their toxic effect on ecosystems and human health. This study attempted to isolate halotolerant PAHs degrading bacteria from saline oil-contaminated soils. Among the isolates, strain KDI with the highest 16S rRNA gene sequence similarity to Labedella gwakjiensis was able to reduce surface tension (ST) from 65.42 to 26.60 mN m^-1 and increase the emulsification index to 81.04%, as a result of significant biosurfactant production. Response Surface Methodology (RSM) analysis was applied to optimize the factors, i.e. PAHs concentration and NaCl concentration as well as to determine the effect of these important variables on PAHs biodegradation. The Carbon Quantum Dots. Iron Oxide (CQDs.Fe₃O₄) nanoparticles were characterized by several popular analytical techniques, after which the effect of CQD.Fe₃O₄ nanoparticles on biodegradation was examined. PAHs biodegradation rate and efficiency of strain KDI to degrade PHE in the presence of CQD.Fe₃O₄ nanoparticles was analyzed by GC. According to the results during biodegradation both the concentration of PAHs and the amount of NaCl were effective. The biodegradation rate significantly increased in the presence of CQD.Fe₃O₄. The highest biodegradation of PHE occurred in the presence of 0.5 g/L of CQD.Fe₃O₄ which was 63.63% and 81.77% after 48 and 72 h of incubation. To the best of our knowledge, this is the first report on optimization of PAHs concentration and salinity by RSM and nanobioremediation of PHE using a bacterial strain in the presence of CQD.Fe₃O₄ nanoparticles.

Reduction of the Bacterial Genome by Transposon-Mediated Random Deletion

Genome reduction is an important strategy in synthetic biology for constructing functional chassis cells or minimal genomes. However, the limited knowledge of complex gene functions and interactions makes genome reduction by rational design encounter a bottleneck. Here, we present an iterative and random genome reduction method for Escherichia coli, named "transposon-mediated random deletion (TMRD)". TMRD generates random double-strand breaks (DSBs) in the genome by combining Tn5 transposition with the CRISPR/Cas9 system and allows genomic deletions of various sizes at random positions during DSB repair through the intracellular alternative end-joining mechanism. Using E. coli MG1655 as the original strain, a pool of cells with multiple random genomic deletions were obtained after five reduction cycles. The growth rates of the obtained cells were comparable to that of MG1655, while the electroporation efficiency increased by at least 2 magnitudes. TMRD can generate a small E. coli library carrying multiple and random genomic deletions while enriching the cells with environmental fitness in the population. TMRD has the potential to be widely applied in the construction of minimal genomes or chassis cells for metabolic engineering.

The rise and the fall of a Pseudomonas aeruginosa endemic lineage in a hospital

The biological features that allow a pathogen to survive in the hospital environment are mostly unknown. The extinction of bacterial epidemics in hospitals is mostly attributed to changes in medical practice, including infection control, but the role of bacterial adaptation has never been documented. We analysed a collection of Pseudomonas aeruginosa isolates belonging to the Besançon Epidemic Strain (BES), responsible for a 12year nosocomial outbreak, using a genotype-to-phenotype approach. Bayesian analysis estimated the emergence of the clone in the hospital 5 years before its opening, during the creation of its water distribution network made of copper. BES survived better than the reference strains PAO1 and PA14 in a copper solution due to a genomic island containing 13 metal-resistance genes and was specifically able to proliferate in the ubiquitous amoeba Vermamoeba vermiformis. Mutations affecting amino-acid metabolism, antibiotic resistance, lipopolysaccharide biosynthesis, and regulation were enriched during the spread of BES. Seven distinct regulatory mutations attenuated the overexpression of the genes encoding the efflux pump MexAB-OprM over time. The fitness of BES decreased over time in correlation with its genome size. Overall, the resistance to inhibitors and predators presumably aided the proliferation and propagation of BES in the plumbing system of the hospital. The pathogen further spread among patients via multiple routes of contamination. The decreased prevalence of patients infected by BES mirrored the parallel and convergent genomic evolution and reduction that affected bacterial fitness. Along with infection control measures, this may have participated in the extinction of BES in the hospital setting.

Chlorhexidine leads to the evolution of antibiotic-resistant Pseudomonas aeruginosa

Antimicrobial resistance is a major public-health concern. We evaluate chlorhexidine role in selection of resistant Pseudomonas aeruginosa mutants and their antibiotic cross-resistance. Mutation frequency and mutation rate after short-term exposure to sub-inhibitory concentrations of chlorhexidine were compared to those after spontaneous chlorhexidine-exposure, in P. aeruginosa PAO1 strain. Chlorhexidine-resistant mutants were generated, either by serial passage in increasing chlorhexidine concentrations or by single exposure to lethal chlorhexidine concentration. The generated mutants were tested for cross-resistance to different antibiotics, by determination of minimum inhibitory concentrations (MIC). The accompanied phenotypic changes in membrane permeability, outer membrane proteins (OMP), and efflux function were evaluated. The effect of exposure to chlorhexidine on MexAB-OprM, MexEF-oprN, and MexXY efflux pumps expression was investigated. No significant change was recorded between the mutation frequencies and mutation rates after short-term exposure to sub-inhibitory concentrations of chlorhexidine and after spontaneous chlorhexidine-exposure, in P. aeruginosa PAO1 strain. Twelve stable mutants, with ≥ eight-fold increase in chlorhexidine MIC, were generated. Several mutants showed increase in the MIC of colistin, cefepime, ceftazidime, meropenem, ciprofloxacin, and amikacin; seven mutants expressed meropenem cross-resistance. This was accompanied by decreased outer membrane permeability and changes in OMP. Using efflux pump inhibitor, chlorhexidine resistance was reverted in most isolates. Exposure to sub-inhibitory concentration of chlorhexidine induced the expression of MexXY efflux pump. Some resistant mutants had overexpressed MexXY efflux pump. Chlorhexidine can select P. aeruginosa strains with antibiotic cross-resistance. This necessitates implementing special protocols for chlorhexidine use and re-evaluation of its benefit versus risk in personal-care products.

Biochar Mitigates N2O Emission of Microbial Denitrification through Modulating Carbon Metabolism and Allocation of Reducing Power

To elucidate the direct effects of biochar on denitrification metabolism at the cellular level, the global response of model denitrifying soil bacterium (Paracoccus denitrificans) to biochar addition was investigated by physiological, proteomic, and metabolomics analyses. The enhancement effect on denitrification was positively correlated with its pyrolysis temperatures (300-500 °C) and dosages (0.1-1%), regardless of precursors [corn straw (CS) or wheat straw). Moreover, the stimulating effect of CS biochar made at 500 °C (CS-500) was mainly attributed to the bulk particles rather than the released soluble compounds. Without direct contact with cells, bulk CS-500 particles might directly modulate the carbon metabolism by the adsorption of extracellular metabolites. Since carbon flux to storage was shifted to oxidative catabolism and growth assimilation, more share of the produced reducing power was used for nitrogen reduction. Meanwhile, except for nitrate reductase, both protein expressions and enzyme activities of nitrite reductase, nitric oxide reductase, and nitrous oxide reductase were up-regulated. Accordingly, the accumulation of N₂O was reduced by 98% due to the optimized electron distribution among denitrifying enzymes. Eventually, the growth rate of Pc. denitrificans enhanced because of the improved energy utilization efficiency. These results updated the regulation mechanism of biochar on denitrification metabolism and N₂O mitigation.

Correlation between the spatial distribution and colony size was common for monogenetic bacteria in laboratory conditions

Background

Geographically separated population growth of microbes is a common phenomenon in microbial ecology. Colonies are representative of the morphological characteristics of this structured population growth. Pattern formation by single colonies has been intensively studied, whereas the spatial distribution of colonies is poorly investigated.

Results

The present study describes a first trial to address the questions of whether and how the spatial distribution of colonies determines the final colony size using the model microorganism Escherichia coli, colonies of which can be grown under well-controlled laboratory conditions. A computational tool for image processing was developed to evaluate colony density, colony size and size variation, and the Voronoi diagram was applied for spatial analysis of colonies with identical space resources. A positive correlation between the final colony size and the Voronoi area was commonly identified, independent of genomic and nutritional differences, which disturbed the colony size and size variation.

Conclusions

This novel finding of a universal correlation between the spatial distribution and colony size not only indicated the fair distribution of spatial resources for monogenetic colonies growing with identical space resources but also indicated that the initial localization of the microbial colonies decided by chance determined the fate of the subsequent population growth. This study provides a valuable example for quantitative analysis of the complex microbial ecosystems by means of experimental ecology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02180-8.

Targeted Genome Reduction of Pseudomonas aeruginosa Strain PAO1 Led to the Development of Hypovirulent and Hypersusceptible rDNA Hosts

Pseudomonas aeruginosa is a human opportunistic pathogen responsible for nosocomial infections, which is largely used as a model organism to study antibiotic resistance and pathogenesis. As other species of the genus, its wide metabolic versatility appears to be attractive to study biotechnological applications. However, its natural resistance to antibiotics and its capacity to produce a wide range of virulence factors argue against its biotechnological potential. By reducing the genome of the reference strain PAO1, we explored the development of four hypovirulent and hypersusceptible recombinant DNA hosts (rDNA hosts). Despite deleting up to 0.8% of the core genome, any of the developed strains presented alterations of fitness when cultured under standard laboratory conditions. Other features such as antibiotic susceptibility, cytotoxicity, in vivo pathogenesis, and expression of heterologous peptides were also explored to highlight the potential applications of these models. This work stands as the first stage of the development of a safe-platform strain of Pseudomonas aeruginosa that will be further optimized for biotechnological applications.

Investigating the release of ZnO nanoparticles from cement mortars on microbiological models

Incorporating zinc oxide nanoparticles (ZnO NPs) into cement mortars may provide additional functions, e.g., self-cleaning and antibacterial or electroconductive ability. However, these NPs are also known for their potential toxicity. During the life cycle of a cement mortar, various abrasive forces cause the release of admixtures to the natural environment. The effect of the released NPs on model microorganisms has not been extensively studied. Previous studies have shown that nanomaterials may affect various microorganisms’ physiological responses, including changes in metabolic activity, biofilming, or growth rate. In this study, we have focused on evaluating the response of model microorganisms, Escherichia coli, Pseudomonas aeruginosa, Staphylococcus aureus, and Candida albicans, towards ZnO nanoparticles released from cement mortars in different deterioration scenarios. The addition of ZnO nanoparticles to cement mortars had a noticeable effect on impeding the strength development. We have also detected that depending on the deterioration scenario, the release of ZnO nanoparticles was varied. Our studies have also shown that even though the release of nanoform ZnO could be limited by poor dispersion or the used filtration technique, the eluates have caused slight but statistically significant changes in the physiological features of studied microorganisms showing relatively low toxicity.

Single Mathematical Parameter for Evaluation of the Microorganisms' Growth as the Objective Function in the Optimization by the DOE Techniques

The cultivation of bacteria sets a ground for studying biological processes in many scientific disciplines. The development of the bacterial population is commonly described with three factors that can be used to evaluate culture conditions. However, selecting only one of them for the optimization protocol is rather problematic and may lead to unintended errors. Therefore, we proposed a novel mathematical approach to obtain a single factor that could be used as the objective function to evaluate the whole growth dynamic and support the optimization of the biomass production process. The sigmoidal-shape curve, which is the commonly used function to plot the amount of biomass versus time, was the base for the mathematical analysis. The key process parameters, such as maximal specific growth rate and lag-phase duration were established with the use of mathematical coefficients of the model curve and combined to create the single growth parameter. Moreover, this parameter was used for the exemplary optimization of the cultivation conditions of Klebsiella pneumoniae that was cultured to be further used in the production of lytic bacteriophages. The proposed growth parameter was successfully validated and used to calculate the optimal process temperature of the selected bacterial strain. The obtained results indicated that the proposed mathematical approach could be effortlessly adapted for a precise evaluation of growth curves.

Clustering of Bacterial Growth Dynamics in Response to Growth Media by Dynamic Time Warping

Bacterial growth curves, representing population dynamics, are still poorly understood. The growth curves are commonly analyzed by model-based theoretical fitting, which is limited to typical S-shape fittings and does not elucidate the dynamics in their entirety. Thus, whether a certain growth condition results in any particular pattern of growth curve remains unclear. To address this question, up-to-date data mining techniques were applied to bacterial growth analysis for the first time. Dynamic time warping (DTW) and derivative DTW (DDTW) were used to compare the similarity among 1015 growth curves of 28 Escherichia coli strains growing in three different media. In the similarity evaluation, agglomerative hierarchical clustering, assessed with four statistic benchmarks, successfully categorized the growth curves into three clusters, roughly corresponding to the three media. Furthermore, a simple benchmark was newly proposed, providing a highly improved accuracy (~99%) in clustering the growth curves corresponding to the growth media. The biologically reasonable categorization of growth curves suggested that DTW and DDTW are applicable for bacterial growth analysis. The bottom-up clustering results indicate that the growth media determine some specific patterns of population dynamics, regardless of genomic variation, and thus have a higher priority of shaping the growth curves than the genomes do.