Mutation accumulation under UV radiation in Escherichia coli

Expanding the focus of the One Health concept: links between the Earth-system processes of the planetary boundaries framework and antibiotic resistance

The scientific community warns that our impact on planet Earth is so acute that we are crossing several of the planetary boundaries that demarcate the safe operating space for humankind. Besides, there is mounting evidence of serious effects on people's health derived from the ongoing environmental degradation. Regarding human health, the spread of antibiotic resistant bacteria is one of the most critical public health issues worldwide. Relevantly, antibiotic resistance has been claimed to be the quintessential One Health issue. The One Health concept links human, animal, and environmental health, but it is frequently only focused on the risk of zoonotic pathogens to public health or, to a lesser extent, the impact of contaminants on human health, i.e., adverse effects on human health coming from the other two One Health "compartments". It is recurrently claimed that antibiotic resistance must be approached from a One Health perspective, but such statement often only refers to the connection between the use of antibiotics in veterinary practice and the antibiotic resistance crisis, or the impact of contaminants (antibiotics, heavy metals, disinfectants, etc.) on antibiotic resistance. Nonetheless, the nine Earth-system processes considered in the planetary boundaries framework can be directly or indirectly linked to antibiotic resistance. Here, some of the main links between those processes and the dissemination of antibiotic resistance are described. The ultimate goal is to expand the focus of the One Health concept by pointing out the links between critical Earth-system processes and the One Health quintessential issue, i.e., antibiotic resistance.

Exploring the links between SOS response, mutagenesis, and resistance during the recovery period

Although the mechanistic connections between SOS-induced mutagenesis and antibiotic resistance are well established, our current understanding of the impact of SOS response levels, recovery durations, and transcription/translation activities on mutagenesis remains relatively limited. In this study, when bacterial cells were exposed to mutagens like ultraviolet light for defined time intervals, a compelling connection between the rate of mutagenesis and the RecA-mediated SOS response levels became evident. Our observations also indicate that mutagenesis primarily occurs during the subsequent recovery phase following the removal of the mutagenic agent. When transcription/translation was inhibited or energy molecules were depleted at the onset of treatment or during the early recovery phase, there was a noticeable decrease in SOS response activation and mutagenesis. However, targeting these processes later in the recovery phase does not have the same effect in reducing mutagenesis, suggesting that the timing of inhibiting transcription/translation or depleting energy molecules is crucial for their efficacy in reducing mutagenesis. Active transcription, translation, and energy availability within the framework of SOS response and DNA repair mechanisms appear to be conserved attributes, supported by their consistent manifestation across diverse conditions, including the use of distinct mutagens such as fluoroquinolones and various bacterial strains.

Competition and evolutionary selection among core regulatory motifs in gene expression control

Gene products that are beneficial in one environment may become burdensome in another, prompting the emergence of diverse regulatory schemes that carry their own bioenergetic cost. By ensuring that regulators are only expressed when needed, we demonstrate that autoregulation generally offers an advantage in an environment combining mutation and time-varying selection. Whether positive or negative feedback emerges as dominant depends primarily on the demand for the target gene product, typically to ensure that the detrimental impact of inevitable mutations is minimized. While self-repression of the regulator curbs the spread of these loss-of-function mutations, self-activation instead facilitates their propagation. By analyzing the transcription network of multiple model organisms, we reveal that reduced bioenergetic cost may contribute to the preferential selection of autoregulation among transcription factors. Our results not only uncover how seemingly equivalent regulatory motifs have fundamentally different impact on population structure, growth dynamics, and evolutionary outcomes, but they can also be leveraged to promote the design of evolutionarily robust synthetic gene circuits.

Construction of a Compact Array of Microplasma Jet Devices and Its Application for Random Mutagenesis of Rhodosporidium toruloides

A small and efficient DNA mutation-inducing machine was constructed with an array of microplasma jet devices (7 × 1) that can be operated at atmospheric pressure for microbial mutagenesis. Using this machine, we report disruption of a plasmid DNA and generation of mutants of an oleaginous yeast Rhodosporidium toruloides. Specifically, a compact-sized microplasma channel (25 × 20 × 2 mm3) capable of generating an electron density of greater than 1013 cm-3 was constructed to produce reactive species (N2*, N2+, O, OH, and Hα) under helium atmospheric conditions to induce DNA mutagenesis. The length of microplasma channels in the device played a critical role in augmenting both the volume of plasma and the concentration of reactive species. First, we confirmed that microplasma treatment can linearize a plasmid by creating nicks in vitro. Second, we treated R. toruloides cells with a jet device containing 7 microchannels for 5 min; 94.8% of the treated cells were killed, and 0.44% of surviving cells showed different colony colors as compared to their parental colony. Microplasma-based DNA mutation is energy-efficient and can be a safe alternative for inducing mutations compared to conventional methods using toxic mutagens. This compact and scalable device is amenable for industrial strain improvement involving large-scale mutagenesis.

Light tolerance of extended spectrum β-lactamase producing Escherichia coli strains after repetitive exposure to far-UVC and blue LED light

AIMS

The aim of this study was to investigate dual far-UVC (Ultraviolet-C) (222 nm) and blue LED (Light Emitting Diode) (405 nm) light on the inactivation of extended spectrum β-lactamase-producing Escherichia coli (ESBL-Ec) and to determine if repetitive exposure to long pulses of light resulted in changes to light tolerance, and antibiotic susceptibility.

METHODS AND RESULTS

Antimicrobial efficiency of dual and individual light wavelengths and development of light tolerance in E. coli was evaluated through a spread plate method after exposure to light at 25 cm. Dual light exposure for 30 min resulted in a 5-6 log10 CFU mL-1 reduction in two ESBL-Ec and two antibiotic-sensitive control E. coli strains. The overall inhibition achieved by dual light treatment was always greater than the combined reductions (log10 CFU) observed from exposure to individual light wavelengths (combined 222-405 nm), indicating a synergistic relationship between blue LED and far-UVC light when used together. Repetitive long pulses of dual and individual far-UVC light exposure resulted in light tolerance in two ESBL-Ec strains but not the antibiotic-sensitive E. coli strains. Subsequent passages of repetitive light-treated ESBL-Ec strains continued to exhibit light tolerance. Antibiotic susceptibility was determined through a standard disk diffusion method. No changes were observed in the antibiotic susceptibility profiles for any of the four strains after exposure to either dual or individual wavelengths.

CONCLUSIONS

Dual light exposure was effective in the disinfection of ESBL-Ec in solution; however, antibiotic-resistant E. coli were able to develop light tolerance after repetitive exposure to light.

Model selection and robust inference of mutational signatures using Negative Binomial non-negative matrix factorization

BACKGROUND

The spectrum of mutations in a collection of cancer genomes can be described by a mixture of a few mutational signatures. The mutational signatures can be found using non-negative matrix factorization (NMF). To extract the mutational signatures we have to assume a distribution for the observed mutational counts and a number of mutational signatures. In most applications, the mutational counts are assumed to be Poisson distributed, and the rank is chosen by comparing the fit of several models with the same underlying distribution and different values for the rank using classical model selection procedures. However, the counts are often overdispersed, and thus the Negative Binomial distribution is more appropriate.

RESULTS

We propose a Negative Binomial NMF with a patient specific dispersion parameter to capture the variation across patients and derive the corresponding update rules for parameter estimation. We also introduce a novel model selection procedure inspired by cross-validation to determine the number of signatures. Using simulations, we study the influence of the distributional assumption on our method together with other classical model selection procedures. We also present a simulation study with a method comparison where we show that state-of-the-art methods are highly overestimating the number of signatures when overdispersion is present. We apply our proposed analysis on a wide range of simulated data and on two real data sets from breast and prostate cancer patients. On the real data we describe a residual analysis to investigate and validate the model choice.

CONCLUSIONS

With our results on simulated and real data we show that our model selection procedure is more robust at determining the correct number of signatures under model misspecification. We also show that our model selection procedure is more accurate than the available methods in the literature for finding the true number of signatures. Lastly, the residual analysis clearly emphasizes the overdispersion in the mutational count data. The code for our model selection procedure and Negative Binomial NMF is available in the R package SigMoS and can be found at https://github.com/MartaPelizzola/SigMoS .

Abiotic selection of microbial genome size in the global ocean

Strong purifying selection is considered a major evolutionary force behind small microbial genomes in the resource-poor photic ocean. However, very little is currently known about how the size of prokaryotic genomes evolves in the global ocean and whether patterns reflect shifts in resource availability in the epipelagic and relatively stable deep-sea environmental conditions. Using 364 marine microbial metagenomes, we investigate how the average genome size of uncultured planktonic prokaryotes varies across the tropical and polar oceans to the hadal realm. We find that genome size is highest in the perennially cold polar ocean, reflecting elongation of coding genes and gene dosage effects due to duplications in the interior ocean microbiome. Moreover, the rate of change in genome size due to temperature is 16-fold higher than with depth up to 200 m. Our results demonstrate how environmental factors can influence marine microbial genome size selection and ecological strategies of the microbiome.

High-throughput approaches to understand and engineer bacteriophages

Bacteriophage research has been vital to fundamental aspects of modern biology. Advances in metagenomics have revealed treasure troves of new and uncharacterized bacteriophages ('phages') that are not yet understood. However, our ability to find new phages has outpaced our understanding of how sequence encodes function in phages. Traditional approaches for characterizing phages are limited in scale and face hurdles in determining how changes in sequence drive function. We describe powerful emerging technologies that can be used to clarify sequence-function relationships in phages through high-throughput genome engineering. Using these approaches, up to 105 variants can be characterized through pooled selection experiments and deep sequencing. We describe caveats when using these tools and provide examples of basic science and engineering goals that are pursuable using these approaches.

Genome-wide association study of Klebsiella pneumoniae identifies variations linked to carbapenems resistance

Klebsiella pneumoniae (KP) is one of the microorganisms that can acquire carbapenem-resistance (CR), and few antimicrobial therapy options exist for infections caused by Carbapenem-Resistant KP (CRKP). In recent years, with the increase of carbapenem resistance rates, treating CRKP has become a serious public health threat in clinical practice. We have collected 2,035 clinical KP isolates from a tertiary hospital in China. Whole genome sequencing data coupled with their binary antimicrobial susceptibility testing data were obtained to conduct the genome-wide association study using a bayesian-based method, including single nucleotide polymorphisms (SNPs) and genes. We identified 28 and 37 potential maker genes associated with imipenem and meropenem resistance, respectively. Among which 19 of them were selected in both drugs by genome-wide association study (GWAS), 11 genes among them were simultaneously validated in independent datasets. These genes were likely related to biofilm formation, efflux pump, and DNA repairing. Moreover, we identified 13 significant CR related SNPs in imipenem or meropenem, with one SNP located in the non-coding region and validated in the independent datasets. Our study indicates complex mechanisms of carbapenems resistance and further investigation of CRKP-related factors are warranted to better understand their contributions to carbapenems resistance. These identified biomarkers may provide targets for future drug interventions or treatments.

Global coordination of the mutation and growth rates across the genetic and nutritional variety in Escherichia coli

Fitness and mutability are the primary traits of living organisms for adaptation and evolution. However, their quantitative linkage remained largely deficient. Whether there is any general relationship between the two features and how genetic and environmental variables influence them remained unclear and were addressed here. The mutation and growth rates of an assortment of Escherichia coli strain collections, including the wild-type strains and the genetically disturbed strains of either reduced genomes or deletion of the genes involved in the DNA replication fidelity, were evaluated in various media. The contribution of media to the mutation and growth rates was differentiated depending on the types of genetic disturbance. Nevertheless, the negative correlation between the mutation and growth rates was observed across the genotypes and was common in all media. It indicated the comprehensive association of the correlated mutation and growth rates with the genetic and medium variation. Multiple linear regression and support vector machine successfully predicted the mutation and growth rates and the categories of genotypes and media, respectively. Taken together, the study provided a quantitative dataset linking the mutation and growth rates, genotype, and medium and presented a simple and successful example of predicting bacterial growth and mutability by data-driven approaches.

Purifying selection enduringly acts on the sequence evolution of highly expressed proteins in Escherichia coli

The evolutionary speed of a protein sequence is constrained by its expression level, with highly expressed proteins evolving relatively slowly. This negative correlation between expression levels and evolutionary rates (known as the E–R anticorrelation) has already been widely observed in past macroevolution between species from bacteria to animals. However, it remains unclear whether this seemingly general law also governs recent evolution, including past and de novo, within a species. However, the advent of genomic sequencing and high-throughput phenotyping, particularly for bacteria, has revealed fundamental gaps between the 2 evolutionary processes and has provided empirical data opposing the possible underlying mechanisms which are widely believed. These conflicts raise questions about the generalization of the E–R anticorrelation and the relevance of plausible mechanisms. To explore the ubiquitous impact of expression levels on molecular evolution and test the relevance of the possible underlying mechanisms, we analyzed the genome sequences of 99 strains of Escherichia coli for evolution within species in nature. We also analyzed genomic mutations accumulated under laboratory conditions as a model of de novo evolution within species. Here, we show that E–R anticorrelation is significant in both past and de novo evolution within species in E. coli. Our data also confirmed ongoing purifying selection on highly expressed genes. Ongoing selection included codon-level purifying selection, supporting the relevance of the underlying mechanisms. However, the impact of codon-level purifying selection on the constraints in evolution within species might be smaller than previously expected from evolution between species.

Genetic variation along an altitudinal gradient in the Phytophthora infestans effector gene Pi02860

Effector genes, together with climatic and other environmental factors, play multifaceted roles in the development of plant diseases. Understanding the role of environmental factors, particularly climate conditions affecting the evolution of effector genes, is important for predicting the long-term value of the genes in controlling agricultural diseases. Here, we collected Phytophthora infestans populations from five locations along a mountainous hill in China and sequenced the effector gene Pi02860 from >300 isolates. To minimize the influence of other ecological factors, isolates were sampled from the same potato cultivar on the same day. We also expressed the gene to visualise its cellular location, assayed its pathogenicity and evaluated its response to experimental temperatures. We found that Pi02860 exhibited moderate genetic variation at the nucleotide level which was mainly generated by point mutation. The mutations did not change the cellular location of the effector gene but significantly modified the fitness of P. infestans. Genetic variation and pathogenicity of the effector gene were positively associated with the altitude of sample sites, possibly due to increased mutation rate induced by the vertical distribution of environmental factors such as UV radiation and temperature. We further found that Pi02860 expression was regulated by experimental temperature with reduced expression as experimental temperature increased. Together, these results indicate that UV radiation and temperature are important environmental factors regulating the evolution of effector genes and provide us with considerable insight as to their future sustainable action under climate and other environmental change.

Role of environmental stresses in elevating resistance mutations in bacteria: Phenomena and mechanisms

Mutations are an important origin of antibiotic resistance in bacteria. While there is increasing evidence showing promoted resistance mutations by environmental stresses, no retrospective research has yet been conducted on this phenomenon and its mechanisms. Herein, we summarized the phenomena of stress-elevated resistance mutations in bacteria, generalized the regulatory mechanisms and discussed the environmental and human health implications. It is shown that both chemical pollutants, such as antibiotics and other pharmaceuticals, biocides, metals, nanoparticles and disinfection byproducts, and non-chemical stressors, such as ultraviolet radiation, electrical stimulation and starvation, are capable of elevating resistance mutations in bacteria. Notably, resistance mutations are more likely to occur under sublethal or subinhibitory levels of these stresses, suggesting a considerable environmental concern. Further, mechanisms for stress-induced mutations are summarized in several points, namely oxidative stress, SOS response, DNA replication and repair systems, RpoS regulon and biofilm formation, all of which are readily provoked by common environmental stresses. Given bacteria in the environment are confronted with a variety of unfavorable conditions, we propose that the stress-elevated resistance mutations are a universal phenomenon in the environment and represent a nonnegligible risk factor for ecosystems and human health. The present review identifies a need for taking into account the pollutants' ability to elevate resistance mutations when assessing their environmental and human health risks and highlights the necessity of including resistance mutations as a target to prevent antibiotic resistance evolution.

Whole genome analysis of UV-induced mutagenesis in Caulobacter crescentus

UV-induced mutagenesis is, to greater extent, a phenomenon dependent on translesion synthesis (TLS) and regulated by the SOS response in bacteria. Caulobacter crescentus, like many bacterial species, employs the ImuABC (ImuAB DnaE2) pathway in TLS. To have a better understanding of the characteristics of UV-induced mutagenesis in this organism, we performed a whole genome analysis of mutations present in survivors after an acute UVC exposure (300 J/m²). We found an average of 3.2 mutations/genome in irradiated samples, distributed in a mutational spectrum consisting exclusively of base substitutions, including tandem mutations. Although limited in conclusions by the small number of mutations identified, our study points to the feasibility of using whole-genome sequencing to study mutagenesis occurring in experiments involving a single acute exposure to genotoxic agents.

The effect of weak clonal interference on average fitness trajectories in the presence of macroscopic epistasis

Adaptation dynamics on fitness landscapes is often studied theoretically in the strong-selection, weak-mutation regime. However, in a large population, multiple beneficial mutants can emerge before any of them fixes in the population. Competition between mutants is known as clonal interference, and while it is known to slow down the rate of adaptation (when compared to the strong-selection, weak-mutation model with the same parameters), how it affects the shape of long-term fitness trajectories in the presence of epistasis is an open question. Here, by considering how changes in fixation probabilities arising from weak clonal interference affect the dynamics of adaptation on fitness-parameterized landscapes, we find that the change in the shape of fitness trajectory arises only through changes in the supply of beneficial mutations (or equivalently, the beneficial mutation rate). Furthermore, a depletion of beneficial mutations as a population climbs up the fitness landscape can speed up the rescaled fitness trajectory (where adaptation speed is measured relative to its value at the start of the experiment), while an enhancement of the beneficial mutation rate does the opposite of slowing it down. Our findings suggest that by carrying out evolution experiments in both regimes (with and without clonal interference), one could potentially distinguish the different sources of macroscopic epistasis (fitness effect of mutations vs change in fraction of beneficial mutations).

DNA Damage Protection for Enhanced Bacterial Survival Under Simulated Low Earth Orbit Environmental Conditions in Escherichia coli

Some organisms have shown the ability to naturally survive in extreme environments, even outer space. Some of these have natural mechanisms to resist severe DNA damage from conditions such as ionizing and non-ionizing radiation, extreme temperatures, and low pressures or vacuum. A good example can be found in Deinococcus radiodurans, which was exposed to severe conditions such as those listed in the Exposure Facility of the International Space Station (ISS) for up to three years. Another example are tardigrades (Ramazzottius varieornatus) which are some of the most resilient animals known. In this study, the survival under simulated Low earth Orbit (LEO) environmental conditions was tested in Escherichia coli. The radiation resistance of this bacteria was enhanced using the Dsup gene from R. varieornatus, and two more genes from D. radiodurans involved in DNA damage repair, RecA and uvrD. The enhanced survival to wide ranges of temperatures and low pressures was then tested in the new strains. This research constitutes a first step in the creation of new bacterial strains engineered to survive severe conditions and adapting existing species for their survival in remote environments, including extra-terrestrial habitats. These strains could be key for the development of environments hospitable to life and could be of use for ecological restoration and space exploration. In addition, studying the efficacy and the functioning of the DNA repair mechanisms used in this study could be beneficial for medical and life sciences engineering.

Honey-based trap for Pseudomonas: a sustainable prototype for water disinfection

This paper introduces a novel prototype for the removal of Pseudomonas from water samples. Bacterial cells have the tendency to get attracted towards specific chemicals (chemotaxis); a 'honey-based trap' (henceforth, addressed as 'honey-trap') strip was conceptualized by integrating a combination of serine, pseudomonas-specific chemoattractant and honey to attract and inhibit the bacteria in situ. Honey, a natural antimicrobial agent, has garnered the attention as an effective inhibitor for Pseudomonal biofilms and wound infections. Dipping serine side of the strip attracted bacteria towards honey-trap, whereby the porous nature of the strip facilitated the 'trapping' and subsequent diffusion of the bacterial cells towards honey-adsorbed end of the strip. This 'honey-trap' reportedly leads to the targeted elimination of Pseudomonas, hence facilitating its removal. The percentage efficacy of this 'honey-trap' device is 96% with a log reduction equivalent to 1.6 within a time frame of 2 h. Pseudomonas aeruginosa, although, not a natural contaminant of potable water, enters circulation due to improperly maintained plumbing fixtures and storage facilities. Honey-trap strip is an easy to use, biodegradable and cost-effective sustainable solution, and thus a scaled-up version of this device may enable substantial improvement in quality of potable water. Schematics showing the preparation and working of the Pseudomonas Honey-trap. Serine as an attractant and honey as an inhibitor was absorbed on filter strips (HT) for use. The strip was dipped in culture from serine end. After different time period of incubation, difference in bacterial load was confirmed by measuring the electrical conductivity and OD600nm of the culture. Additionally, inhibitory effect of HS was confirmed by placing the strip incubated with culture on agar plates and differences in bacterial lawn were monitored. Removal of bacterial cells from the suspension was also confirmed using absorption spectroscopy.

Genotoxic Agents Produce Stressor-Specific Spectra of Spectinomycin Resistance Mutations Based on Mechanism of Action and Selection in Bacillus subtilis

Mutagenesis is integral for bacterial evolution and the development of antibiotic resistance. Environmental toxins and stressors are known to elevate the rate of mutagenesis through direct DNA toxicity, known as stress-associated mutagenesis, or via a more general stress-induced process that relies on intrinsic bacterial pathways. Here, we characterize the spectra of mutations induced by an array of different stressors using high-throughput sequencing to profile thousands of spectinomycin-resistant colonies of Bacillus subtilis. We found 69 unique mutations in the rpsE and rpsB genes, and that each stressor leads to a unique and specific spectrum of antibiotic-resistance mutations. While some mutations clearly reflected the DNA damage mechanism of the stress, others were likely the result of a more general stress-induced mechanism. To determine the relative fitness of these mutants under a range of antibiotic selection pressures, we used multistrain competitive fitness experiments and found an additional landscape of fitness and resistance. The data presented here support the idea that the environment in which the selection is applied (mutagenic stressors that are present), as well as changes in local drug concentration, can significantly alter the path to spectinomycin resistance in B. subtilis.

Emerging trends in the photodynamic inactivation (PDI) applied to the food decontamination

The food and drink manufacturing industry is constantly seeking for alternative sanitation and disinfection systems that may achieve the same antimicrobial efficiency of conventional chemical sanitisers and at the same time be convenient in terms of energy and water savings. A candidate technology for this purpose is the use of light in combination with photosensitisers (PS) to generate a bioactive effect against microbial agents in a process defined as photodynamic inactivation (PDI). This technology can be applied to the food processing of different food matrices to reduce the microbial load of foodborne pathogens such as bacteria, fungi, viruses and protozoa. Also, the PDI can be exploited to increase the shelf-life period of food by inactivation of spoiling microbes. This review analyses new developments in the last five years for PDI systems applied to the food decontamination from foodborne pathogens. The photosensitisation mechanisms and methods are reported to introduce the applied technology against microbial targets in food matrices. Recent blue light emitting diodes (LED) lamp systems for the PDI mediated by endogenous PS are discussed as well PDI technologies with the use of exogenous PS from plant sources such as curcumin and porphyrin-based molecules. The updated overview of the most recent developments in the PDI technology both in wavelengths and employed PS will provide further points of analysis for the advancement of the research on new competitive and effective disinfection systems in the food industry.

Spatiotemporal Manipulation of the Mismatch Repair System of Pseudomonas putida Accelerates Phenotype Emergence

The development of complex phenotypes in industrially relevant bacteria is a major goal of metabolic engineering, which encompasses the implementation of both rational and random approaches. In the latter case, several tools have been developed toward increasing mutation frequencies, yet the precise control of mutagenesis processes in cell factories continues to represent a significant technical challenge. Pseudomonas species are endowed with one of the most efficient DNA mismatch repair (MMR) systems found in the bacterial domain. Here, we investigated if the endogenous MMR system could be manipulated as a general strategy to artificially alter mutation rates in Pseudomonas species. To bestow a conditional mutator phenotype in the platform bacterium Pseudomonas putida, we constructed inducible mutator devices to modulate the expression of the dominant-negative mutL^E36K allele. Regulatable overexpression of mutL^E36K in a broad-host-range, easy-to-cure plasmid format resulted in a transitory inhibition of the MMR machinery, leading to a significant increase (up to 438-fold) in DNA mutation frequencies and a heritable fixation of mutations in the genome. Following such an accelerated mutagenesis-followed by selection approach, three phenotypes were successfully evolved: resistance to antibiotics streptomycin and rifampicin (either individually or combined) and reversion of a synthetic uracil auxotrophy. Thus, these mutator devices could be applied to accelerate the evolution of metabolic pathways in long-term evolutionary experiments, alternating cycles of (inducible) mutagenesis coupled to selection schemes toward the desired phenotype(s).

Biodiversity-based development and evolution: the emerging research systems in model and non-model organisms

Evolutionary developmental biology, or Evo-Devo for short, has become an established field that, broadly speaking, seeks to understand how changes in development drive major transitions and innovation in organismal evolution. It does so via integrating the principles and methods of many subdisciplines of biology. Although we have gained unprecedented knowledge from the studies on model organisms in the past decades, many fundamental and crucially essential processes remain a mystery. Considering the tremendous biodiversity of our planet, the current model organisms seem insufficient for us to understand the evolutionary and physiological processes of life and its adaptation to exterior environments. The currently increasing genomic data and the recently available gene-editing tools make it possible to extend our studies to non-model organisms. In this review, we review the recent work on the regulatory signaling of developmental and regeneration processes, environmental adaptation, and evolutionary mechanisms using both the existing model animals such as zebrafish and Drosophila, and the emerging nonstandard model organisms including amphioxus, ascidian, ciliates, single-celled phytoplankton, and marine nematode. In addition, the challenging questions and new directions in these systems are outlined as well.

Can multidrug-resistant organisms become resistant to ultraviolet (UV) light following serial exposures? Characterization of post-UV genomic changes using whole-genome sequencing

OBJECTIVES

No-touch disinfection systems like xenon- or mercury-based ultraviolet (UV) are now commonly being used for hospital room disinfection. However, serial exposure to UV light can potentially lead to the development of bacterial resistance. We sought to determine whether UV resistance develops due to serial exposure to UV light using 3 epidemiologically important multidrug-resistant microbial strains.

METHODS

Methicillin-resistant Staphylococcus aureus (MRSA), carbapenemase-producing Klebsiella pneumoniae (KPC) and metallo-β-lactamase-producing Klebsiella pneumoniae (MBL) were serially exposed to 25 growth-irradiation cycles of UV produced by a xenon-based UV (Xe-UV) lamp for 5 minutes or a mercury-based UV (Hg-UV) lamp for 10 minutes. After each UV exposure cycle, the surviving colony-forming units (CFUs) were measured and compared with the initial inoculum of each cycle for each strain, respectively.

RESULTS

In each cycle, ˜1-10 million of MRSA, KPC, and MBL were used to test the effect of UV irradiation. Postexposure colony counts remained low (3-100 colonies) throughout the 25 serial exposures to both xenon- and mercury-based UV. The log-kill rate after each exposure showed no changes following UV disinfection by Xe-UV. The MRSA log-kill rate increased after repeated exposure to Hg-UV unlike KPC and MBL K. pneumoniae, which did not change. Whole-genome sequencing (WGS) analyses performed on these 3 strains demonstrated no significant genetic changes after multiple UV irradiation cycles.

CONCLUSIONS

Exposure of multidrug-resistant bacteria to UV produced from 2 different UV sources did not engender UV resistance after 25 serial exposures, as demonstrated by WGS analysis; thus, UV disinfection is unlikely to generate UV-resistant hospital flora.

High-throughput laboratory evolution reveals evolutionary constraints in Escherichia coli

Understanding the constraints that shape the evolution of antibiotic resistance is critical for predicting and controlling drug resistance. Despite its importance, however, a systematic investigation of evolutionary constraints is lacking. Here, we perform a high-throughput laboratory evolution of Escherichia coli under the addition of 95 antibacterial chemicals and quantified the transcriptome, resistance, and genomic profiles for the evolved strains. Utilizing machine learning techniques, we analyze the phenotype–genotype data and identified low dimensional phenotypic states among the evolved strains. Further analysis reveals the underlying biological processes responsible for these distinct states, leading to the identification of trade-off relationships associated with drug resistance. We also report a decelerated evolution of β-lactam resistance, a phenomenon experienced by certain strains under various stresses resulting in higher acquired resistance to β-lactams compared to strains directly selected by β-lactams. These findings bridge the genotypic, gene expression, and drug resistance gap, while contributing to a better understanding of evolutionary constraints for antibiotic resistance.

Understanding evolutionary constraints in antibiotic resistance is crucial for prediction and control. Here, the authors use high-throughput laboratory evolution of Escherichia coli alongside machine learning to identify trade-off relationships associated with drug resistance.

In situ characterisation and manipulation of biological systems with Chi.Bio

The precision and repeatability of in vivo biological studies is predicated upon methods for isolating a targeted subsystem from external sources of noise and variability. However, in many experimental frameworks, this is made challenging by nonstatic environments during host cell growth, as well as variability introduced by manual sampling and measurement protocols. To address these challenges, we developed Chi.Bio, a parallelised open-source platform that represents a new experimental paradigm in which all measurement and control actions can be applied to a bulk culture in situ. In addition to continuous-culturing capabilities, it incorporates tunable light outputs, spectrometry, and advanced automation features. We demonstrate its application to studies of cell growth and biofilm formation, automated in silico control of optogenetic systems, and readout of multiple orthogonal fluorescent proteins in situ. By integrating precise measurement and actuation hardware into a single low-cost platform, Chi.Bio facilitates novel experimental methods for synthetic, systems, and evolutionary biology and broadens access to cutting-edge research capabilities.

Complete Genomic Analysis of Enterococcus faecium Heat-Resistant Strain Developed by Two-Step Adaptation Laboratory Evolution Method

Stress resistance is an important trait expected of lactic acid bacteria used in food manufacturing. Among the various sources of stress, high temperature is a key factor that interrupts bacterial growth. In this regards, constant efforts are made for the development of heat-resistant strains, but few studies were done accompanying genomic analysis to identify the causal factors of the resistance mechanisms. Furthermore, it is also thought that tolerance to multiple stresses are equally important. Herein, we isolated one Enterococcus faecium strain named BIOPOP-3 and completed a full-length genome sequence. Using this strain, a two-step adaptive laboratory evolution (ALE) method was applied to obtain a heat-resistant strain, BIOPOP-3 ALE. After sequencing the whole genome, we compared the two full-length sequences and identified one non-synonymous variant and four indel variants that could potentially confer heat resistance, which were technically validated by resequencing. We experimentally verified that the evolved strain was significantly enhanced in not only heat resistance but also acid and bile resistance. We demonstrated that the developed heat-resistant strain can be applied in animal feed manufacturing processes. The multi-stress-resistant BIOPOP-3 ALE strain developed in this study and the two-step ALE method are expected to be widely applied in industrial and academic fields. In addition, we expect that the identified variants which occurred specifically in heat-resistant strain will enhance molecular biological understanding and be broadly applied to the biological engineering field.

Changes in Cell Membrane Fatty Acid Composition of Streptococcus thermophilus in Response to Gradually Increasing Heat Temperature

In this study, a method of heat adaptation was implemented in an attempt to increase the upper thermal threshold of two Streptococcus thermophilus found in South Korea and identified the alterations in membrane fatty acid composition to adaptive response to heat. In order to develop heat tolerant lactic acid bacteria, heat treatment was continuously applied to bacteria by increasing temperature from 60°C until the point that no surviving cell was detected. Our results indicated significant increase in heat tolerance of heat-adapted strains compared to the wild type (WT) strains. In particular, the survival ratio of basically low heat-tolerant strain increased even more. In addition, the strains with improved heat tolerance acquired cross protection, which improved their survival ratio in acid, bile salts and osmotic conditions. A relation between heat tolerance and membrane fatty acid composition was identified. As a result of heat adaptation, the ratio of unsaturated to saturated fatty acids (UFA/SFA) and C18:1 relative concentration were decreased. C6:0 in only heatadapted strains and C22:0 in only the naturally high heat tolerant strain were detected. These results support the hypothesis, that the consequent increase of SFA ratio is a cellular response to environmental stresses such as high temperatures, and it is able to protect the cells from acid, bile salts and osmotic conditions via cross protection. This study demonstrated that the increase in heat tolerance can be utilized as a mean to improve bacterial tolerance against various environmental stresses.

Decoding the evolution and transmissions of the novel pneumonia coronavirus (SARS-CoV-2 / HCoV-19) using whole genomic data

The outbreak of COVID-19 started in mid-December 2019 in Wuhan, China. Up to 29 February 2020, SARS-CoV-2 (HCoV-19 / 2019-nCoV) had infected more than 85 000 people in the world. In this study, we used 93 complete genomes of SARS-CoV-2 from the GISAID EpiFlu ^TM database to investigate the evolution and human-to-human transmissions of SARS-CoV-2 in the first two months of the outbreak. We constructed haplotypes of the SARS-CoV-2 genomes, performed phylogenomic analyses and estimated the potential population size changes of the virus. The date of population expansion was calculated based on the expansion parameter tau ( τ) using the formula t= τ/2 u. A total of 120 substitution sites with 119 codons, including 79 non-synonymous and 40 synonymous substitutions, were found in eight coding-regions in the SARS-CoV-2 genomes. Forty non-synonymous substitutions are potentially associated with virus adaptation. No combinations were detected. The 58 haplotypes (31 found in samples from China and 31 from outside China) were identified in 93 viral genomes under study and could be classified into five groups. By applying the reported bat coronavirus genome (bat-RaTG13-CoV) as the outgroup, we found that haplotypes H13 and H38 might be considered as ancestral haplotypes, and later H1 was derived from the intermediate haplotype H3. The population size of the SARS-CoV-2 was estimated to have undergone a recent expansion on 06 January 2020, and an early expansion on 08 December 2019. Furthermore, phyloepidemiologic approaches have recovered specific directions of human-to-human transmissions and the potential sources for international infected cases.

Decoding the evolution and transmissions of the novel pneumonia coronavirus (SARS-CoV-2 / HCoV-19) using whole genomic data

The outbreak of COVID-19 started in mid-December 2019 in Wuhan, China. Up to 29 February 2020, SARS-CoV-2 (HCoV-19 / 2019-nCoV) had infected more than 85 000 people in the world. In this study, we used 93 complete genomes of SARS-CoV-2 from the GISAID EpiFlu TM database to investigate the evolution and human-to-human transmissions of SARS-CoV-2 in the first two months of the outbreak. We constructed haplotypes of the SARS-CoV-2 genomes, performed phylogenomic analyses and estimated the potential population size changes of the virus. The date of population expansion was calculated based on the expansion parameter tau ( τ) using the formula t= τ/2 u. A total of 120 substitution sites with 119 codons, including 79 non-synonymous and 40 synonymous substitutions, were found in eight coding-regions in the SARS-CoV-2 genomes. Forty non-synonymous substitutions are potentially associated with virus adaptation. No combinations were detected. The 58 haplotypes (31 found in samples from China and 31 from outside China) were identified in 93 viral genomes under study and could be classified into five groups. By applying the reported bat coronavirus genome (bat-RaTG13-CoV) as the outgroup, we found that haplotypes H13 and H38 might be considered as ancestral haplotypes, and later H1 was derived from the intermediate haplotype H3. The population size of the SARS-CoV-2 was estimated to have undergone a recent expansion on 06 January 2020, and an early expansion on 08 December 2019. Furthermore, phyloepidemiologic approaches have recovered specific directions of human-to-human transmissions and the potential sources for international infected cases.

Phenotypic delay in the evolution of bacterial antibiotic resistance: Mechanistic models and their implications

Phenotypic delay-the time delay between genetic mutation and expression of the corresponding phenotype-is generally neglected in evolutionary models, yet recent work suggests that it may be more common than previously assumed. Here, we use computer simulations and theory to investigate the significance of phenotypic delay for the evolution of bacterial resistance to antibiotics. We consider three mechanisms which could potentially cause phenotypic delay: effective polyploidy, dilution of antibiotic-sensitive molecules and accumulation of resistance-enhancing molecules. We find that the accumulation of resistant molecules is relevant only within a narrow parameter range, but both the dilution of sensitive molecules and effective polyploidy can cause phenotypic delay over a wide range of parameters. We further investigate whether these mechanisms could affect population survival under drug treatment and thereby explain observed discrepancies in mutation rates estimated by Luria-Delbrück fluctuation tests. While the effective polyploidy mechanism does not affect population survival, the dilution of sensitive molecules leads both to decreased probability of survival under drug treatment and underestimation of mutation rates in fluctuation tests. The dilution mechanism also changes the shape of the Luria-Delbrück distribution of mutant numbers, and we show that this modified distribution provides an improved fit to previously published experimental data.

The remarkable UV light invulnerability of thymine GNA dinucleotides

We herein demonstrate the UV resistance of glycol nucleic acid (GNA) dinucleotides. This resistance sustains the hypothesis of GNA as a nucleic acid prebiotic ancestor on early Earth, a time of intense solar UV light. Such photorobustness, due to the absence of intrastrand base stacking, could offer an opportunity for nanodevice development requiring challenging UV conditions.

Inactivation of a Mismatch-Repair System Diversifies Genotypic Landscape of Escherichia coli During Adaptive Laboratory Evolution

Adaptive laboratory evolution (ALE) is used to find causal mutations that underlie improved strain performance under the applied selection pressure. ALE studies have revealed that mutator populations tend to outcompete their non-mutator counterparts following the evolutionary trajectory. Among them, mutS-inactivated mutator cells, characterize d by a dysfunctional methyl-mismatch repair system, are frequently found in ALE experiments. Here, we examined mutS inactivation as an approach to facilitate ALE of Escherichia coli. The wild-type E. coli MG1655 and mutS knock-out derivative (ΔmutS) were evolved in parallel for 800 generations on lactate or glycerol minimal media in a serial-transfer experiment. Whole-genome re-sequencing of each lineage at 100-generation intervals revealed that (1) mutations emerge rapidly in the ΔmutS compared to in the wild-type strain; (2) mutations were more than fourfold higher in the ΔmutS strain at the end-point populations compared to the wild-type strain; and (3) a significant number of random mutations accumulated in the ΔmutS strains. We then measured the fitness of the end-point populations on an array of non-adaptive carbon sources. Interestingly, collateral fitness increases on non-adaptive carbon sources were more pronounced in the ΔmutS strains than the parental strain. Fitness measurement of single mutants revealed that the collateral fitness increase seen in the mutator lineages can be attributed to a pool of random mutations. Together, this study demonstrates that short-term mutator ALE extensively expands possible genotype space, resulting in versatile bacteria with elevated fitness levels across various carbon sources.

The emergence of adaptive laboratory evolution as an efficient tool for biological discovery and industrial biotechnology

Harnessing the process of natural selection to obtain and understand new microbial phenotypes has become increasingly possible due to advances in culturing techniques, DNA sequencing, bioinformatics, and genetic engineering. Accordingly, Adaptive Laboratory Evolution (ALE) experiments represent a powerful approach both to investigate the evolutionary forces influencing strain phenotypes, performance, and stability, and to acquire production strains that contain beneficial mutations. In this review, we summarize and categorize the applications of ALE to various aspects of microbial physiology pertinent to industrial bioproduction by collecting case studies that highlight the multitude of ways in which evolution can facilitate the strain construction process. Further, we discuss principles that inform experimental design, complementary approaches such as computational modeling that help maximize utility, and the future of ALE as an efficient strain design and build tool driven by growing adoption and improvements in automation.

Low genetic variation is associated with low mutation rate in the giant duckweed

Mutation rate and effective population size (N_e) jointly determine intraspecific genetic diversity, but the role of mutation rate is often ignored. Here we investigate genetic diversity, spontaneous mutation rate and N_e in the giant duckweed (Spirodela polyrhiza). Despite its large census population size, whole-genome sequencing of 68 globally sampled individuals reveals extremely low intraspecific genetic diversity. Assessed under natural conditions, the genome-wide spontaneous mutation rate is at least seven times lower than estimates made for other multicellular eukaryotes, whereas N_e is large. These results demonstrate that low genetic diversity can be associated with large-N_e species, where selection can reduce mutation rates to very low levels. This study also highlights that accurate estimates of mutation rate can help to explain seemingly unexpected patterns of genome-wide variation.

While the role of effective population size (N_e) in explaining variation in genetic diversity has received much attention, the role of spontaneous mutation rate is largely ignored. Here, Xu et al. show that giant duckweed has a high N_e yet low genetic diversity, likely due to its low mutation rate.

Strategies of Phytophthora infestans adaptation to local UV radiation conditions

Expected global changes in environmental conditions underline the need for a better understanding of genetic variation in ecological traits and their strategies of adaptation to the stresses. In this study, evolutionary mechanisms and processes of UV adaptation in plant pathogens were investigated by combining statistical genetics, physiological assays, and common garden experiment approaches in an assessment of the potato late blight pathogen, Phytophthora infestans, sampled from various geographic locations in China. We found spatial divergence caused by diversifying selection in UV tolerance among P. infestans populations. Local UV radiation was the driving force of selection as indicated by a positive correlation between UV tolerance in P. infestans populations and the altitude of collection sites. Plasticity accounted for 68% of population variation while heritability was negligible, suggesting temporary changes in gene expression and/or enzymatic activity play a more important role than permanent modification of gene structure in the evolution of UV adaptation. This adaptation strategy may explain the lack of fitness penalty observed in genotypes with higher UV tolerance.

Reframing gene essentiality in terms of adaptive flexibility

BACKGROUND

Essentiality assays are important tools commonly utilized for the discovery of gene functions. Growth/no growth screens of single gene knockout strain collections are also often utilized to test the predictive power of genome-scale models. False positive predictions occur when computational analysis predicts a gene to be non-essential, however experimental screens deem the gene to be essential. One explanation for this inconsistency is that the model contains the wrong information, possibly an incorrectly annotated alternative pathway or isozyme reaction. Inconsistencies could also be attributed to experimental limitations, such as growth tests with arbitrary time cut-offs. The focus of this study was to resolve such inconsistencies to better understand isozyme activities and gene essentiality.

RESULTS

In this study, we explored the definition of conditional essentiality from a phenotypic and genomic perspective. Gene-deletion strains associated with false positive predictions of gene essentiality on defined minimal medium for Escherichia coli were targeted for extended growth tests followed by population sequencing and transcriptome analysis. Of the twenty false positive strains available and confirmed from the Keio single gene knock-out collection, 11 strains were shown to grow with longer incubation periods making these actual true positives. These strains grew reproducibly with a diverse range of growth phenotypes. The lag phase observed for these strains ranged from less than one day to more than 7 days. It was found that 9 out of 11 of the false positive strains that grew acquired mutations in at least one replicate experiment and the types of mutations ranged from SNPs and small indels associated with regulatory or metabolic elements to large regions of genome duplication. Comparison of the detected adaptive mutations, modeling predictions of alternate pathways and isozymes, and transcriptome analysis of KO strains suggested agreement for the observed growth phenotype for 6 out of the 9 cases where mutations were observed.

CONCLUSIONS

Longer-term growth experiments followed by whole genome sequencing and transcriptome analysis can provide a better understanding of conditional gene essentiality and mechanisms of adaptation to such perturbations. Compensatory mutations are largely reproducible mechanisms and are in agreement with genome-scale modeling predictions to loss of function gene deletion events.

Development of an Automated UV Irradiation Device for Microbial Cell Culture

Ultraviolet (UV) mutagenesis is a widely used technique to increase bacterial mutation rates in laboratory experiments. UV mutagenesis requires fine regulation of UV dose, because the number of dead cells increases exponentially as the dose increases. Ignoring this hazard can cause extinction of UV-exposed populations. Therefore, an automated system that cooperatively conducts both growth measurement and UV irradiation is needed for efficient UV mutagenesis experiments. To address this task, we constructed an automated UV irradiation device for microbial cell culture. This device can measure cell density and irradiate the bacterial cells with UV light automatically according to the state of cell growth. We demonstrated that this growth feedback control avoided extinction and enabled accumulation of mutations in bacterial genomes at a rapid rate for a long period. Whole-genome sequencing revealed the high accumulation rate, neutrality, and spectrum of UV-induced mutations. These characteristics were all consistent with those obtained by manual UV irradiation. These results indicate that our automated device is useful in accelerating mutation accumulation over a long duration.