Rates of transposition in Escherichia coli

A metagenomics pipeline reveals insertion sequence-driven evolution of the microbiota

Insertion sequence (IS) elements are mobile genetic elements in bacterial genomes that support adaptation. We developed a database of IS elements coupled to a computational pipeline that identifies IS element insertions in the microbiota. We discovered that diverse IS elements insert into the genomes of intestinal bacteria regardless of human host lifestyle. These insertions target bacterial accessory genes that aid in their adaptation to unique environmental conditions. Using IS expansion in Bacteroides, we show that IS activity leads to the insertion of "hot spots" in accessory genes. We show that IS insertions are stable and can be transferred between humans. Extreme environmental perturbations force IS elements to fall out of the microbiota, and many fail to rebound following homeostasis. Our work shows that IS elements drive bacterial genome diversification within the microbiota and establishes a framework for understanding how strain-level variation within the microbiota impacts human health.

Insertion Sequences Determine Plasmid Adaptation to New Bacterial Hosts

Plasmids facilitate the vertical and horizontal spread of antimicrobial resistance genes between bacteria. The host range and adaptation of plasmids to new hosts determine their impact on the spread of resistance. In this work, we explore the mechanisms driving plasmid adaptation to novel hosts in experimental evolution. Using the small multicopy plasmid pB1000, usually found in Pasteurellaceae, we studied its adaptation to a host from a different bacterial family, Escherichia coli. We observed two different mechanisms of adaptation. One mechanism is single nucleotide polymorphisms (SNPs) in the origin of replication (oriV) of the plasmid, which increase the copy number in E. coli cells, elevating the stability, and resistance profile. The second mechanism consists of two insertion sequences (ISs), IS1 and IS10, which decrease the fitness cost of the plasmid by disrupting an uncharacterized gene on pB1000 that is harmful to E. coli. Both mechanisms increase the stability of pB1000 independently, but only their combination allows long-term maintenance. Crucially, we show that the mechanisms have a different impact on the host range of the plasmid. SNPs in oriV prevent the replication in the original host, resulting in a shift of the host range. In contrast, the introduction of ISs either shifts or expands the host range, depending on the IS. While IS1 leads to expansion, IS10 cannot be reintroduced into the original host. This study gives new insights into the relevance of ISs in plasmid-host adaptation to understand the success in spreading resistance. IMPORTANCE ColE1-like plasmids are small, mobilizable plasmids that can be found across at least four orders of Gammaproteobacteria and are strongly associated with antimicrobial resistance genes. Plasmid pB1000 carries the gene bla_ROB-1, conferring high-level resistance to penicillins and cefaclor. pB1000 has been described in various species of the family Pasteurellaceae, for example, in Haemophilus influenzae, which can cause diseases such as otitis media, meningitis, and pneumonia. To understand the resistance spread through horizontal transfer, it is essential to study the mechanisms of plasmid adaptation to novel hosts. In this work we identify that a gene from pB1000, which encodes a peptide that is toxic for E. coli, and the low plasmid copy number (PCN) of pB1000 in E. coli cells are essential targets in the described plasmid-host adaptation and therefore limit the spread of pB1000-encoded bla_ROB-1. Furthermore, we show how the interplay of two adaptation mechanisms leads to successful plasmid maintenance in a different bacterial family.

A population genetics theory for piRNA-regulated transposable elements

Transposable elements (TEs) are self-reproducing selfish DNA sequences that can invade the genome of virtually all living species. Population genetics models have shown that TE copy numbers generally reach a limit, either because the transposition rate decreases with the number of copies (transposition regulation) or because TE copies are deleterious, and thus purged by natural selection. Yet, recent empirical discoveries suggest that TE regulation may mostly rely on piRNAs, which require a specific mutational event (the insertion of a TE copy in a piRNA cluster) to be activated - the so-called TE regulation "trap model". We derived new population genetics models accounting for this trap mechanism, and showed that the resulting equilibria differ substantially from previous expectations based on a transposition-selection equilibrium. We proposed three sub-models, depending on whether or not genomic TE copies and piRNA cluster TE copies are selectively neutral or deleterious, and we provide analytical expressions for maximum and equilibrium copy numbers, as well as cluster frequencies for all of them. In the full neutral model, the equilibrium is achieved when transposition is completely silenced, and this equilibrium does not depend on the transposition rate. When genomic TE copies are deleterious but not cluster TE copies, no long-term equilibrium is possible, and active TEs are eventually eliminated after an active incomplete invasion stage. When all TE copies are deleterious, a transposition-selection equilibrium exists, but the invasion dynamics is not monotonic, and the copy number peaks before decreasing. Mathematical predictions were in good agreement with numerical simulations, except when genetic drift and/or linkage disequilibrium dominates. Overall, the trap-model dynamics appeared to be substantially more stochastic and less repeatable than traditional regulation models.

Targeted Screening for Spontaneous Insertion Mutations in a Lactic Acid Bacterium, Tetragenococcus halophilus

Studies on the microorganisms used in food production are of interest because microbial genotypes are reflected in food qualities such as taste, flavor, and yield. However, several microbes are nonmodel organisms, and their analysis is often limited by the lack of genetic tools. Tetragenococcus halophilus, a halophilic lactic acid bacterium used in soy sauce fermentation starter culture, is one such microorganism. The lack of DNA transformation techniques for T. halophilus makes gene complementation and disruption assays difficult. Here, we report that the endogenous insertion sequence ISTeha4, belonging to the IS4 family, is translocated at an extremely high frequency in T. halophilus and causes insertional mutations at various loci. We developed a method named targeting spontaneous insertional mutations in genomes (TIMING), which combines high-frequency insertional mutations and efficient PCR screening, enabling the isolation of gene mutants of interest from a library. The method provides a reverse genetics and strain improvement tool, does not require the introduction of exogenous DNA constructs, and enables the analysis of nonmodel microorganisms lacking DNA transformation techniques. Our results highlight the important role of insertion sequences as a source of spontaneous mutagenesis and genetic diversity in bacteria. IMPORTANCE Genetic and strain improvement tools to manipulate a gene of interest are required for the nontransformable lactic acid bacterium Tetragenococcus halophilus. Here, we demonstrate that an endogenous transposable element, ISTeha4, is transposed into the host genome at an extremely high frequency. A genotype-based and non-genetically engineered screening system was constructed to isolate knockout mutants using this transposable element. The method described enables a better understanding of the genotype-phenotype relationship and serves as a tool to develop food-grade-appropriate mutants of T. halophilus.

Improving the Editing Efficiency of CRISPR-Cas9 by Reducing the Generation of Escapers Based on the Surviving Mechanism

Due to the high specificity in targeting DNA and highly convenient programmability, CRISPR-Cas-based antimicrobials applied for eliminating specific strains such as antibiotic-resistant bacteria in the microbiome were gradually developed. However, the generation of escapers makes the elimination efficiency far lower than the acceptable rate (10^-8) recommended by the National Institutes of Health. Here, a systematic study was carried out providing insight into the escaping mechanisms in Escherichia coli, and strategies for reducing the escapers were devised accordingly. We first showed an escape rate of 10^-5-10^-3 in E. coli MG1655 under the editing of pEcCas/pEcgRNA established previously. Detailed analysis of the escapers obtained at ligA site in E. coli MG1655 uncovered that the disruption of cas9 was the main cause of the generation of survivors, notably the frequent insertion of IS5. Hence, the sgRNA was next designed to target the "perpetrator" IS5, and subsequently the killing efficiency was improved 4-fold. Additionally, the escape rate in IS-free E. coli MDS42 was also tested at the ligA site, ∼10-fold decrease compared with MG1655, but the disruption of cas9 was still observed in all survivors manifested in the form of frameshifts or point mutations. Thus, we optimized the tool itself by increasing the copy number of cas9 to retain some cas9 that still has the correct DNA sequence. Fortunately, the escape rates dropped below 10^-8 at 9 of the 16 tested genes. Furthermore, the λ-Red recombination system was added to generate the pEcCas-2.0, and a 100% gene deletion efficiency was achieved at genes cadA, maeB, and gntT in MG1655, whereas those genes were edited with low efficiency previously. Last, the application of pEcCas-2.0 was then expanded to the E. coli B strain BL21(DE3) and W strain ATCC9637. This study reveals the mechanism of E. coli surviving Cas9-mediated death, and a highly efficient editing tool is established based on the mechanism, which will accelerate the further application of CRISPR-Cas.

Ancient Darwinian replicators nested within eubacterial genomes

Integrative mobile genetic elements (MGEs), such as transposons and insertion sequences, propagate within bacterial genomes, but persistence times in individual lineages are short. For long-term survival, MGEs must continuously invade new hosts by horizontal transfer. Theoretically, MGEs that persist for millions of years in single lineages, and are thus subject to vertical inheritance, should not exist. Here we draw attention to an exception - a class of MGE termed REPIN. REPINs are non-autonomous MGEs whose duplication depends on non-jumping RAYT transposases. Comparisons of REPINs and typical MGEs show that replication rates of REPINs are orders of magnitude lower, REPIN population size fluctuations correlate with changes in available genome space, REPIN conservation depends on RAYT function, and REPIN diversity accumulates within host lineages. These data lead to the hypothesis that REPINs form enduring, beneficial associations with eubacterial chromosomes. Given replicative nesting, our hypothesis predicts conflicts arising from the diverging effects of selection acting simultaneously on REPINs and host genomes. Evidence in support comes from patterns of REPIN abundance and diversity in two distantly related bacterial species. Together this bolsters the conclusion that REPINs are the genetic counterpart of mutualistic endosymbiotic bacteria.

Molecular Factors and Mechanisms Driving Multidrug Resistance in Uropathogenic Escherichia coli-An Update

The rapid emergence of multidrug-resistant (MDR) bacteria indisputably constitutes a major global health problem. Pathogenic Escherichia coli are listed among the most critical group of bacteria that require fast development of new antibiotics and innovative treatment strategies. Among harmful extraintestinal Enterobacteriaceae strains, uropathogenic E. coli (UPEC) pose a significant health threat. UPEC are considered the major causative factor of urinary tract infection (UTI), the second-most commonly diagnosed infectious disease in humans worldwide. UTI treatment places a substantial financial burden on healthcare systems. Most importantly, the misuse of antibiotics during treatment has caused selection of strains with the ability to acquire MDR via miscellaneous mechanisms resulting in gaining resistance against many commonly prescribed antibiotics like ampicillin, gentamicin, cotrimoxazole and quinolones. Mobile genetic elements (MGEs) such as transposons, integrons and conjugative plasmids are the major drivers in spreading resistance genes in UPEC. The co-occurrence of various bacterial evasion strategies involving MGEs and the SOS stress response system requires further research and can potentially lead to the discovery of new, much-awaited therapeutic targets. Here, we analyzed and summarized recent discoveries regarding the role, mechanisms, and perspectives of MDR in the pathogenicity of UPEC.

Transposable Element Insertions into the Escherichia coli Polysialic Acid Gene Cluster Result in Resistance to the K1F Bacteriophage

Reviewing the genetics underlying the arms race between bacteria and bacteriophages can offer an interesting insight into the development of bacterial resistance and phage co-evolution. This study shows how the natural development of resistances to the K1F bacteriophage, a phage which targets the K1 capsule of pathogenic Escherichia coli, can come about through insertion sequences (IS). Of the K1F resistant mutants isolated, two were of particular interest. The first of these showed full resistance to K1F and was found to have disruptions to kpsE, the product of which is involved in polysialic acid translocation. The second, after showing an initial susceptibility to K1F which then developed to full resistance, had disruptions to neuC, a gene involved in one of the early steps of polysialic acid biosynthesis. Both of these mutations came with a fitness cost and produced considerable phenotypic differences in the completeness and location of the K1 capsule when compared with the wild type. Sequential treatment of these two K1F resistant mutants with T7 resulted in the production of a variety of isolates, many of which showed a renewed susceptibility to K1F, indicating that these insertion sequence mutations are reversible, as well as one isolate that developed resistance to both phages.

IMPORTANCE Bacteriophages have many potential uses in industry and the clinical environment as an antibacterial control measure. One of their uses, phage therapy, is an appealing alternative to antibiotics due to their high specificity. However, as with the rise in antimicrobial resistance (AMR), it is critical to improve our understanding of how resistance develops against these viral agents. In the same way as bacteria will evolve and mutate antibiotic receptors so they can no longer be recognized, resistance to bacteriophages can come about via mutations to phage receptors, preventing phage binding and infection. We have shown that Escherichia coli will become resistant to the K1F bacteriophage via insertion element reshufflings causing null mutations to elements of the polysialic acid biosynthetic cluster. Exposure to the T7 bacteriophage then resulted in further changes in the position of these IS elements, further altering their resistance and sensitivity profiles.

Frequency, composition and mobility of Escherichia coli-derived transposable elements in holdings of plasmid repositories

By providing the scientific community with uniform and standardized resources of consistent quality, plasmid repositories play an important role in enabling scientific reproducibility. Plasmids containing insertion sequence elements (IS elements) represent a challenge from this perspective, as they can change the plasmid structure and function. In this study, we conducted a systematic analysis of a subset of plasmid stocks distributed by plasmid repositories (The Arabidopsis Biological Resource Center and Addgene) which carry unintended integrations of bacterial mobile genetic elements. The integration of insertion sequences was most often found in, but not limited to, pBR322‐derived vectors, and did not affect the function of the specific plasmids. In certain cases, the entire stock was affected, but the majority of the stocks tested contained a mixture of the wild‐type and the mutated plasmids, suggesting that the acquisition of IS elements likely occurred after the plasmids were acquired by the repositories. However, comparison of the sequencing results of the original samples revealed that some plasmids already carried insertion mutations at the time of donation. While an extensive BLAST analysis of 47 877 plasmids sequenced from the Addgene repository uncovered IS elements in only 1.12%, suggesting that IS contamination is not widespread, further tests showed that plasmid integration of IS elements can propagate in conventional Escherichia coli hosts over a few tens of generations. Use of IS‐free E. coli hosts prevented the emergence of IS insertions as well as that of small indels, suggesting that the use of IS‐free hosts by donors and repositories could help limit unexpected and unwanted IS integrations into plasmids.

Bacterial resistance to CRISPR-Cas antimicrobials

In the age of antibiotic resistance and precise microbiome engineering, CRISPR-Cas antimicrobials promise to have a substantial impact on the way we treat diseases in the future. However, the efficacy of these antimicrobials and their mechanisms of resistance remain to be elucidated. We systematically investigated how a target E. coli strain can escape killing by episomally-encoded CRISPR-Cas9 antimicrobials. Using Cas9 from Streptococcus pyogenes (SpCas9) we studied the killing efficiency and resistance mutation rate towards CRISPR-Cas9 antimicrobials and elucidated the underlying genetic alterations. We find that killing efficiency is not correlated with the number of cutting sites or the type of target. While the number of targets did not significantly affect efficiency of killing, it did reduce the emergence of chromosomal mutations conferring resistance. The most frequent target of resistance mutations was the plasmid-encoded SpCas9 that was inactivated by bacterial genome rearrangements involving translocation of mobile genetic elements such as insertion elements. This resistance mechanism can be overcome by re-introduction of an intact copy of SpCas9. The work presented here provides a guide to design strategies that reduce resistance and improve the activity of CRISPR-Cas antimicrobials.

Genome Plasticity by Insertion Sequences Learned From a Case of Radiation-Resistant Bacterium Deinococcus geothermalis

The genome of the radiation-resistant bacterium Deinococcus geothermalis contains 19 types of insertion sequences (ISs), including 93 total transposases (Tpases) in 73 full-length ISs from the main chromosome and 2 mega plasmids. In this study, 68 ISs from the D. geothermalis genome were extracted to implicate the earlier genome before its mutation by transposition of ISs. The total size of eliminated ISs from genome was 78.85 kb. From these in silico corrections of mutation by the ISs, we have become aware of some bioinformatics factualness as follows: (1) can reassemble the disrupted genes if the exact IS region was eliminated, (2) can configure the schematic clustering of major DDE type Tpases, (3) can determine IS integration order across multiple hot spots, and (4) can compare genetic relativeness by the partial synteny analysis between D. geothermalis and Deinococcus strain S9. Recently, we found that several IS elements actively transferred to other genomic sites under hydrogen peroxide-induced oxidative stress conditions, resulting in the inactivation of functional genes. Therefore, the single species genome’s mobilome study provides significant support to define bacterial genome plasticity and molecular evolution from past and present progressive transposition events.

Evolutionary mechanisms that determine which bacterial genes are carried on plasmids

The evolutionary pressures that determine the location (chromosomal or plasmid‐borne) of bacterial genes are not fully understood. We investigate these pressures through mathematical modeling in the context of antibiotic resistance, which is often found on plasmids. Our central finding is that gene location is under positive frequency‐dependent selection: the higher the frequency of one form of resistance compared to the other, the higher its relative fitness. This can keep moderately beneficial genes on plasmids, despite occasional plasmid loss. For these genes, positive frequency dependence leads to a priority effect: whichever form is acquired first—through either mutation or horizontal gene transfer—has time to increase in frequency and thus becomes difficult to displace. Higher rates of horizontal transfer of plasmid‐borne than chromosomal genes therefore predict moderately beneficial genes will be found on plasmids. Gene flow between plasmid and chromosome allows chromosomal forms to arise, but positive frequency‐dependent selection prevents these from establishing. Further modeling shows that this effect is particularly pronounced when genes are shared across a large number of species, suggesting that antibiotic resistance genes are often found on plasmids because they are moderately beneficial across many species. We also revisit previous theoretical work—relating to the role of local adaptation in explaining gene location and to plasmid persistence—in light of our findings.

Targeted-antibacterial-plasmids (TAPs) combining conjugation and CRISPR/Cas systems achieve strain-specific antibacterial activity

The global emergence of drug-resistant bacteria leads to the loss of efficacy of our antibiotics arsenal and severely limits the success of currently available treatments. Here, we developed an innovative strategy based on targeted-antibacterial-plasmids (TAPs) that use bacterial conjugation to deliver CRISPR/Cas systems exerting a strain-specific antibacterial activity. TAPs are highly versatile as they can be directed against any specific genomic or plasmid DNA using the custom algorithm (CSTB) that identifies appropriate targeting spacer sequences. We demonstrate the ability of TAPs to induce strain-selective killing by introducing lethal double strand breaks (DSBs) into the targeted genomes. TAPs directed against a plasmid-born carbapenem resistance gene efficiently resensitise the strain to the drug. This work represents an essential step toward the development of an alternative to antibiotic treatments, which could be used for in situ microbiota modification to eradicate targeted resistant and/or pathogenic bacteria without affecting other non-targeted bacterial species.

Local genetic context shapes the function of a gene regulatory network

Gene expression levels are influenced by multiple coexisting molecular mechanisms. Some of these interactions such as those of transcription factors and promoters have been studied extensively. However, predicting phenotypes of gene regulatory networks (GRNs) remains a major challenge. Here, we use a well-defined synthetic GRN to study in Escherichia coli how network phenotypes depend on local genetic context, i.e. the genetic neighborhood of a transcription factor and its relative position. We show that one GRN with fixed topology can display not only quantitatively but also qualitatively different phenotypes, depending solely on the local genetic context of its components. Transcriptional read-through is the main molecular mechanism that places one transcriptional unit (TU) within two separate regulons without the need for complex regulatory sequences. We propose that relative order of individual TUs, with its potential for combinatorial complexity, plays an important role in shaping phenotypes of GRNs.

Evolution along the parasitism-mutualism continuum determines the genetic repertoire of prophages

Integrated into their bacterial hosts' genomes, prophage sequences exhibit a wide diversity of length and gene content, from highly degraded cryptic sequences to intact, functional prophages that retain a full complement of lytic-function genes. We apply three approaches-bioinformatics, analytical modelling and computational simulation-to understand the diverse gene content of prophages. In the bioinformatics work, we examine the distributions of over 50,000 annotated prophage genes identified in 1384 prophage sequences, comparing the gene repertoires of intact and incomplete prophages. These data indicate that genes involved in the replication, packaging, and release of phage particles have been preferentially lost in incomplete prophages, while tail fiber, transposase and integrase genes are significantly enriched. Consistent with these results, our mathematical and computational approaches predict that genes involved in phage lytic function are preferentially lost, resulting in shorter prophages that often retain genes that benefit the host. Informed by these models, we offer novel hypotheses for the enrichment of integrase and transposase genes in cryptic prophages. Overall, we demonstrate that functional and cryptic prophages represent a diversity of genetic sequences that evolve along a parasitism-mutualism continuum.

Short and long-read ultra-deep sequencing profiles emerging heterogeneity across five platform Escherichia coli strains

Reprogramming organisms for large-scale bioproduction counters their evolutionary objectives of fast growth and often leads to mutational collapse of the engineered production pathways during cultivation. Yet, the mutational susceptibility of academic and industrial Escherichia coli bioproduction host strains are poorly understood. In this study, we apply 2^nd and 3^rd generation deep sequencing to profile simultaneous modes of genetic heterogeneity that decimate engineered biosynthetic production in five popular E. coli hosts BL21(DE3), TOP10, MG1655, W, and W3110 producing 2,3-butanediol and mevalonic acid. Combining short-read and long-read sequencing, we detect strain and sequence-specific mutational modes including single nucleotide polymorphism, inversion, and mobile element transposition, as well as complex structural variations that disrupt the integrity of the engineered biosynthetic pathway. Our analysis suggests that organism engineers should avoid chassis strains hosting active insertion sequence (IS) subfamilies such as IS1 and IS10 present in popular E. coli TOP10. We also recommend monitoring for increased mutagenicity in the pathway transcription initiation regions and recombinogenic repeats. Together, short and long sequencing reads identified latent low-frequency mutation events such as a short detrimental inversion within a pathway gene, driven by 8-bp short inverted repeats. This demonstrates the power of combining ultra-deep DNA sequencing technologies to profile genetic heterogeneities of engineered constructs and explore the markedly different mutational landscapes of common E. coli host strains. The observed multitude of evolving variants underlines the usefulness of early mutational profiling for new synthetic pathways designed to sustain in organisms over long cultivation scales.

Active Transposition of Insertion Sequences by Oxidative Stress in Deinococcus geothermalis

Radiation-resistant bacterium Deinococcus geothermalis has a total of 73 insertion sequences (ISs) in genomes, and some of them are actively transposed to other loci with replicative mode due to oxidative stress of hydrogen peroxide treatment. Here, we detected two transposition events in wild-type (WT) strain and LysR family member gene disrupted strain (Δdgeo_2840). Similar to our previous report (Lee et al., 2019), phytoene desaturase (dgeo_0524), a key enzyme of carotenoid biosynthesis, was disrupted by the integration of IS element, thereby detected a single phenotypically non-pigmented colony in each WT and Δdgeo_2840 strain. Two separate types of IS element have been integrated into non-pigmented clones: ISDge11 for WT and ISDge6 for Δdgeo_2840 strain. Surprisingly, Δdgeo_2840 mutant strain revealed higher resistance to oxidative stress than WT strain at late exponential growth phase. From the qRT-PCR analysis, OxyR (dgeo_1888) was highly up-regulated to 30-fold by oxidative stress through hydrogen peroxide treatment in both WT and Δdgeo_2840 mutant strains. However, the oxidative stress response enzyme, catalase or superoxide dismutase, was not significantly induced by overexpressed OxyR. Thus, a putative LysR family regulator Dgeo_2840 controlled the expression of ISDge6 type transposase and the induction of OxyR under oxidative condition. There is LysR family DNA-binding protein dependent active transposition of specific type IS and the up-regulated OxyR has not positively controlled ROS scavenger enzymes in D. geothermalis.

Multicopy Chromosomal Integration Using CRISPR-Associated Transposases

Controlling the copy number of gene expression cassettes is an important strategy to engineer bacterial cells into high-efficiency biocatalysts. Current strategies mostly use plasmid vectors, but multicopy plasmids are often genetically unstable, and their copy numbers cannot be precisely controlled. The integration of expression cassettes into a bacterial chromosome has advantages, but iterative integration is laborious, and it is challenging to obtain a library with varied gene doses for phenotype characterization. Here, we demonstrated that multicopy chromosomal integration using CRISPR-associated transposases (MUCICAT) can be achieved by designing a crRNA to target multicopy loci or a crRNA array to target multiple loci in the Escherichia coli genome. Within 5 days without selection pressure, E. coli strains carrying cargos with successively increasing copy numbers (up to 10) were obtained. Recombinant MUCICAT E. coli containing genomic multicopy glucose dehydrogenase expression cassettes showed 2.6-fold increased expression of this important industrial enzyme compared to E. coli harboring the conventional protein-expressing plasmid pET24a. Successful extension of MUCICAT to Tatumella citrea further demonstrated that MUCICAT may be generally applied to many bacterial species.

Transposable element persistence via potential genome-level ecosystem engineering

BACKGROUND

The nuclear genomes of eukaryotes vary enormously in size, with much of this variability attributable to differential accumulation of transposable elements (TEs). To date, the precise evolutionary and ecological conditions influencing TE accumulation remain poorly understood. Most previous attempts to identify these conditions have focused on evolutionary processes occurring at the host organism level, whereas we explore a TE ecology explanation.

RESULTS

As an alternative (or additional) hypothesis, we propose that ecological mechanisms occurring within the host cell may contribute to patterns of TE accumulation. To test this idea, we conducted a series of experiments using a simulated asexual TE/host system. Each experiment tracked the accumulation rate for a given type of TE within a particular host genome. TEs in this system had a net deleterious effect on host fitness, which did not change over the course of experiments. As one might expect, in the majority of experiments TEs were either purged from the genome or drove the host population to extinction. However, in an intriguing handful of cases, TEs co-existed with hosts and accumulated to very large numbers. This tended to occur when TEs achieved a stable density relative to non-TE sequences in the genome (as opposed to reaching any particular absolute number). In our model, the only way to maintain a stable density was for TEs to generate new, inactive copies at a rate that balanced with the production of active (replicating) copies.

CONCLUSIONS

From a TE ecology perspective, we suggest this could be interpreted as a case of ecosystem engineering within the genome, where TEs persist by creating their own "habitat".

Transposition of Insertion Sequences was Triggered by Oxidative Stress in Radiation-Resistant Bacterium Deinococcus geothermalis

During an oxidative stress-response assay on a putative Dps-like gene-disrupted Δdgeo_0257 mutant strain of radiation-resistant bacterium Deinococcus geothermalis, a non-pigmented colony was observed among the normal reddish color colonies. This non-pigmented mutant cell subsequently displayed higher sensitivity to H₂O₂. While carotenoid has a role in protecting as scavenger of reactive oxygen species the reddish wild-type strain from radiation and oxidative stresses, it is hypothesized that the carotenoid biosynthesis pathway has been disrupted in the mutant D. geothermalis cell. Here, we show that, in the non-pigmented mutant cell of interest, phytoene desaturase (Dgeo_0524, crtI), a key enzyme in carotenoid biosynthesis, was interrupted by transposition of an ISDge7 family member insertion sequence (IS) element. RNA-Seq analysis between wild-type and Δdgeo_0257 mutant strains revealed that the expression level of ISDge5 family transposases, but not ISDge7 family members, were substantially up-regulated in the Δdgeo_0257 mutant strain. We revealed that the non-pigmented strain resulted from the genomic integration of ISDge7 family member IS elements, which were also highly up-regulated, particularly following oxidative stress. The transposition path for both transposases is a replicative mode. When exposed to oxidative stress in the absence of the putative DNA binding protein Dgeo_0257, a reddish D. geothermalis strain became non-pigmented. This transformation was facilitated by transposition of an ISDge7 family IS element into a gene encoding a key enzyme of carotenoid biosynthesis. Further, we present evidence of additional active transposition by the ISDge5 family IS elements, a gene that was up-regulated during the stationary phase regardless of the presence of oxidative stress.

Defensive Function of Transposable Elements in Bacteria

It has been widely debated whether transposable elements have a positive or a negative effect on their host cells. This study demonstrated that transposable elements, specifically insertion sequences (ISs), can adopt a defensive role in Escherichia coli. In three different E. coli strains (S17, DH5α, and Nissle 1917), IS1 and IS10 rapidly disrupted the I-CeuI gene (encoding I-CeuI endonuclease) on the plasmid pLO11-ICeuI as early as the first generation, despite the gene-circuit being under control of an arabinose promoter. Proteomics analysis showed that the protein abundance profile of E. coli DH5α with pLO11-ICeuI in the fifth generation was nearly opposite to that of control strain (E. coli with pLO11, no I-CeuI). The DNA damage caused by the leaky expression of I-CeuI was enough to trigger a SOS response and alter lipid synthesis, ribosomal activity, RNA/DNA metabolism, central dogma and cell cycle processes in E. coli DH5α. After the ISs disrupted the expression of I-CeuI, cells fully recovered by the 31st generation had a protein abundance profile similar to that of the control strain. This study showed that ISs readily mutated a harmful gene which subsequently restored host fitness. These observations have implications for the stability of designed gene circuits in synthetic biology.

DYNAMICS OF AN INSERTION SEQUENCE INFECTION IN A SPATIALLY STRUCTURED ENVIRONMENT

Bacterial insertion sequences (ISs), the simplest form of autonomous mobile DNA, depend on their prokaryote hosts to spread in a spatially structured environment. We use a spatially explicit metapopulation model to simulate the spread of an IS that can have both detrimental and beneficial effects on its host cell. We find that, on the one hand, the spatial structure of the metapopulation and cell dispersal between subpopulations have no strong effect on the time to full infection of the metapopulation. On the other hand, factors that influence the IS infection dynamics within a subpopulation have a strong effect on that time. These factors are mainly the fitness benefit of an IS and the rate of horizontal gene transfer. We also find that the infection process of a metapopulation is very erratic in its early phase. Finally, we show that the infection’s success depends critically on the initially infected subpopulation.

Characterization of a functional insertion sequence ISSau2 from Staphylococcus aureus

Background

ISSau2 has been suggested as a member of the IS150 f subgroup in the IS3 family. It encodes a fusion transposase OrfAB produced by programmed − 1 translational frameshifting with two overlapping reading frames orfA and orfB. To better characterize ISSau2, the binding and cleaving activities of the ISSau2 transposase and its transposition frequency were studied.

Results

The purified ISSau2 transposase OrfAB was a functional protein in vitro since it bound specifically to ISSau2 terminal inverted repeat sequences (IRs) and cleaved the transposon ends at the artificial mini-transposon pUC19-IRL-gfp-IRR. In addition, the transposition frequency of ISSau2 in vivo was approximately 1.76 ± 0.13 × 10^− 3, based on a GFP hop-on assay. Furthermore, OrfB cleaved IRs with the similar catalytic activity of OrfAB, while OrfA had no catalytic activity. Finally, either OrfA or OrfB significantly reduced the transposition of ISSau2 induced by OrfAB.

Conclusion

We have confirmed that ISSau2 is a member of IS150/IS3 family. The ISSau2 transposase OrfAB could bind to and cleave the specific fragments containing the terminal inverted repeat sequences and induce the transposition, suggesting that ISSau2 is at least partially functional. Meanwhile, both OrfA and OrfB inhibited the transposition by ISSau2. Our results will help understand biological roles of ISSau2 in its host S. aureus.

Disentangling the effects of selection and loss bias on gene dynamics

We combine mathematical modeling of genome evolution with comparative analysis of prokaryotic genomes to estimate the relative contributions of selection and intrinsic loss bias to the evolution of different functional classes of genes and mobile genetic elements (MGE). An exact solution for the dynamics of gene family size was obtained under a linear duplication-transfer-loss model with selection. With the exception of genes involved in information processing, particularly translation, which are maintained by strong selection, the average selection coefficient for most nonparasitic genes is low albeit positive, compatible with observed positive correlation between genome size and effective population size. Free-living microbes evolve under stronger selection for gene retention than parasites. Different classes of MGE show a broad range of fitness effects, from the nearly neutral transposons to prophages, which are actively eliminated by selection. Genes involved in antiparasite defense, on average, incur a fitness cost to the host that is at least as high as the cost of plasmids. This cost is probably due to the adverse effects of autoimmunity and curtailment of horizontal gene transfer caused by the defense systems and selfish behavior of some of these systems, such as toxin-antitoxin and restriction modification modules. Transposons follow a biphasic dynamics, with bursts of gene proliferation followed by decay in the copy number that is quantitatively captured by the model. The horizontal gene transfer to loss ratio, but not duplication to loss ratio, correlates with genome size, potentially explaining increased abundance of neutral and costly elements in larger genomes.

Transposable Elements Mediate Adaptive Debilitation of Flagella in Experimental Escherichia coli Populations

Although insertion sequence (IS) elements are generally considered genomic parasites, they can mediate adaptive genetic changes in bacterial genomes. We discovered that among 12 laboratory-evolved Escherichia coli populations, three had experienced at least six different IS1-mediated deletions of flagellar genes. These deletions all involved the master flagellar regulator flhDC, and as such completely incapacitate motility. Two lines of evidence strongly suggest that these deletions were adaptive in our evolution experiment: (1) parallel evolution in three independent populations is highly unlikely just by chance, and (2) one of these deletion mutations swept to fixation within ~1000 generations, which is over two million times faster than expected if this deletion was instead selectively neutral and thus evolving by genetic drift. Because flagella are energetically expensive to synthesize and operate, we suspect that debilitating their construction conferred a fitness advantage in our well-stirred evolution experiment. These findings underscore the important role that IS elements can play in mediating adaptive loss-of-function mutations in bacteria.

A shifting mutational landscape in 6 nutritional states: Stress-induced mutagenesis as a series of distinct stress input-mutation output relationships

Environmental stresses increase genetic variation in bacteria, plants, and human cancer cells. The linkage between various environments and mutational outcomes has not been systematically investigated, however. Here, we established the influence of nutritional stresses commonly found in the biosphere (carbon, phosphate, nitrogen, oxygen, or iron limitation) on both the rate and spectrum of mutations in Escherichia coli. We found that each limitation was associated with a remarkably distinct mutational profile. Overall mutation rates were not always elevated, and nitrogen, iron, and oxygen limitation resulted in major spectral changes but no net increase in rate. Our results thus suggest that stress-induced mutagenesis is a diverse series of stress input-mutation output linkages that is distinct in every condition. Environment-specific spectra resulted in the differential emergence of traits needing particular mutations in these settings. Mutations requiring transpositions were highest under iron and oxygen limitation, whereas base-pair substitutions and indels were highest under phosphate limitation. The unexpected diversity of input-output effects explains some important phenomena in the mutational biases of evolving genomes. The prevalence of bacterial insertion sequence transpositions in the mammalian gut or in anaerobically stored cultures is due to environmentally determined mutation availability. Likewise, the much-discussed genomic bias towards transition base substitutions in evolving genomes can now be explained as an environment-specific output. Altogether, our conclusion is that environments influence genetic variation as well as selection.

The impact of insertion sequences on bacterial genome plasticity and adaptability

Transposable elements (TE), small mobile genetic elements unable to exist independently of the host genome, were initially believed to be exclusively deleterious genomic parasites. However, it is now clear that they play an important role as bacterial mutagenic agents, enabling the host to adapt to new environmental challenges and to colonize new niches. This review focuses on the impact of insertion sequences (IS), arguably the smallest TE, on bacterial genome plasticity and concomitant adaptability of phenotypic traits, including resistance to antibacterial agents, virulence, pathogenicity and catabolism. The direct consequence of IS transposition is the insertion of one DNA sequence into another. This event can result in gene inactivation as well as in modulation of neighbouring gene expression. The latter is usually mediated by de-repression or by the introduction of a complete or partial promoter located within the element. Furthermore, transcription and transposition of IS are affected by host factors and in some cases by environmental signals offering the host an adaptive strategy and promoting genetic variability to withstand the environmental challenges.

Environment-directed activation of the Escherichia coliflhDC operon by transposons

The flagellar system in Escherichia coli K12 is expressed under the control of the flhDC-encoded master regulator FlhDC. Transposition of insertion sequence (IS) elements to the upstream flhDC promoter region up-regulates transcription of this operon, resulting in a more rapid motility. Wang and Wood (ISME J 2011;5:1517-1525) provided evidence that insertion of IS5 into upstream activating sites occurs at higher rates in semi-solid agar media in which swarming behaviour is allowed as compared with liquid or solid media where swarming cannot occur. We confirm this conclusion and show that three IS elements, IS1, IS3 and IS5, transpose to multiple upstream sites within a 370 bp region of the flhDC operon control region. Hot spots for IS insertion correlate with positions of stress-induced DNA duplex destabilization (SIDD). We show that IS insertion occurs at maximal rates in 0.24 % agar, with rates decreasing dramatically with increasing or decreasing agar concentrations. In mixed cultures, we show that these mutations preferentially arise from the wild-type parent at frequencies of up to 3×10-3 cell-1 day-1 when the inoculated parental and co-existing IS-activated mutant cells are entering the stationary growth phase. We rigorously show that the apparent increased mutation frequencies cannot be accounted for by increased swimming or by increased growth under the selective conditions used. Thus, our data are consistent with the possibility that appropriate environmental conditions, namely those that permit but hinder flagellar rotation, result in the activation of a mutational pathway that involves IS element insertion upstream of the flhDC operon.

Stability of Cross-Feeding Polymorphisms in Microbial Communities

Cross-feeding, a relationship wherein one organism consumes metabolites excreted by another, is a ubiquitous feature of natural and clinically-relevant microbial communities and could be a key factor promoting diversity in extreme and/or nutrient-poor environments. However, it remains unclear how readily cross-feeding interactions form, and therefore our ability to predict their emergence is limited. In this paper we developed a mathematical model parameterized using data from the biochemistry and ecology of an E. coli cross-feeding laboratory system. The model accurately captures short-term dynamics of the two competitors that have been observed empirically and we use it to systematically explore the stability of cross-feeding interactions for a range of environmental conditions. We find that our simple system can display complex dynamics including multi-stable behavior separated by a critical point. Therefore whether cross-feeding interactions form depends on the complex interplay between density and frequency of the competitors as well as on the concentration of resources in the environment. Moreover, we find that subtly different environmental conditions can lead to dramatically different results regarding the establishment of cross-feeding, which could explain the apparently unpredictable between-population differences in experimental outcomes. We argue that mathematical models are essential tools for disentangling the complexities of cross-feeding interactions.

A Mutational Hotspot and Strong Selection Contribute to the Order of Mutations Selected for during Escherichia coli Adaptation to the Gut

The relative role of drift versus selection underlying the evolution of bacterial species within the gut microbiota remains poorly understood. The large sizes of bacterial populations in this environment suggest that even adaptive mutations with weak effects, thought to be the most frequently occurring, could substantially contribute to a rapid pace of evolutionary change in the gut. We followed the emergence of intra-species diversity in a commensal Escherichia coli strain that previously acquired an adaptive mutation with strong effect during one week of colonization of the mouse gut. Following this first step, which consisted of inactivating a metabolic operon, one third of the subsequent adaptive mutations were found to have a selective effect as high as the first. Nevertheless, the order of the adaptive steps was strongly affected by a mutational hotspot with an exceptionally high mutation rate of 10-5. The pattern of polymorphism emerging in the populations evolving within different hosts was characterized by periodic selection, which reduced diversity, but also frequency-dependent selection, actively maintaining genetic diversity. Furthermore, the continuous emergence of similar phenotypes due to distinct mutations, known as clonal interference, was pervasive. Evolutionary change within the gut is therefore highly repeatable within and across hosts, with adaptive mutations of selection coefficients as strong as 12% accumulating without strong constraints on genetic background. In vivo competitive assays showed that one of the second steps (focA) exhibited positive epistasis with the first, while another (dcuB) exhibited negative epistasis. The data shows that strong effect adaptive mutations continuously recur in gut commensal bacterial species.

Transposon-mediated activation of the Escherichia coli glpFK operon is inhibited by specific DNA-binding proteins: Implications for stress-induced transposition events

Escherichia coli cells deleted for the cyclic AMP (cAMP) receptor protein (Crp) gene (Δcrp) cannot utilize glycerol because cAMP-Crp is a required activator of the glycerol utilization operon, glpFK. We have previously shown that a transposon, Insertion Sequence 5 (IS5), can insert into the upstream regulatory region of the operon to activate the glpFK promoter and enable glycerol utilization. GlpR, which represses glpFK transcription, binds to the glpFK upstream region near the site of IS5 insertion and inhibits insertion. By adding cAMP to the culture medium in ΔcyaA cells, we here show that the cAMP-Crp complex, which also binds to the glpFK upstream regulatory region, inhibits IS5 hopping into the activating site. Control experiments showed that the frequencies of mutations in response to cAMP were independent of parental cell growth rate and the selection procedure. These findings led to the prediction that glpFK-activating IS5 insertions can also occur in wild-type (Crp⁺) cells under conditions that limit cAMP production. Accordingly, we found that IS5 insertion into the activating site in wild-type cells is elevated in the presence of glycerol and a non-metabolizable sugar analogue that lowers cytoplasmic cAMP concentrations. The resultant IS5 insertion mutants arising in this minimal medium become dominant constituents of the population after prolonged periods of growth. The results show that DNA binding transcription factors can reversibly mask a favored transposon target site, rendering a hot spot for insertion less favored. Such mechanisms could have evolved by natural selection to overcome environmental adversity.

The dynamics of mobile promoters: Enhanced stability in promoter regions

Mobile promoters are emerging as a new class of mobile genetic elements, first identified by examining prokaryote genome sequences, and more recently confirmed by experimental observations in bacteria. Recent datasets have identified over 40,000 putative mobile promoters in sequenced prokaryote genomes, however only one-third of these are in regions of the genome directly upstream from coding sequences, that is, in promoter regions. The presence of many promoter sequences in non-promoter regions is unexplained. Here we develop a general mathematical model for the dynamics of mobile promoters, extending previous work to capture the dynamics both within and outside promoter regions. From this general model, we apply rigorous model selection techniques to identify which parameters are statistically justified in describing the available mobile promoter data, and find best-fit values of these parameters. Our results suggest that high rates of horizontal gene transfer maintain the population of mobile promoters in promoter regions, and that once established at these sites, mobile promoters are rarely lost, but are commonly copied to other genomic regions. In contrast, mobile promoter copies in non-promoter regions are more numerous and more volatile, experiencing substantially higher rates of duplication, loss and diversification.

Rapid compensatory evolution promotes the survival of conjugative plasmids

Conjugative plasmids play a vital role in bacterial adaptation through horizontal gene transfer. Explaining how plasmids persist in host populations however is difficult, given the high costs often associated with plasmid carriage. Compensatory evolution to ameliorate this cost can rescue plasmids from extinction. In a recently published study we showed that compensatory evolution repeatedly targeted the same bacterial regulatory system, GacA/GacS, in populations of plasmid-carrying bacteria evolving across a range of selective environments. Mutations in these genes arose rapidly and completely eliminated the cost of plasmid carriage. Here we extend our analysis using an individual based model to explore the dynamics of compensatory evolution in this system. We show that mutations which ameliorate the cost of plasmid carriage can prevent both the loss of plasmids from the population and the fixation of accessory traits on the bacterial chromosome. We discuss how dependent the outcome of compensatory evolution is on the strength and availability of such mutations and the rate at which beneficial accessory traits integrate on the host chromosome.

Adaptation to Parasites and Costs of Parasite Resistance in Mutator and Nonmutator Bacteria

Parasitism creates selection for resistance mechanisms in host populations and is hypothesized to promote increased host evolvability. However, the influence of these traits on host evolution when parasites are no longer present is unclear. We used experimental evolution and whole-genome sequencing of Escherichia coli to determine the effects of past and present exposure to parasitic viruses (phages) on the spread of mutator alleles, resistance, and bacterial competitive fitness. We found that mutator alleles spread rapidly during adaptation to any of four different phage species, and this pattern was even more pronounced with multiple phages present simultaneously. However, hypermutability did not detectably accelerate adaptation in the absence of phages and recovery of fitness costs associated with resistance. Several lineages evolved phage resistance through elevated mucoidy, and during subsequent evolution in phage-free conditions they rapidly reverted to nonmucoid, phage-susceptible phenotypes. Genome sequencing revealed that this phenotypic reversion was achieved by additional genetic changes rather than by genotypic reversion of the initial resistance mutations. Insertion sequence (IS) elements played a key role in both the acquisition of resistance and adaptation in the absence of parasites; unlike single nucleotide polymorphisms, IS insertions were not more frequent in mutator lineages. Our results provide a genetic explanation for rapid reversion of mucoidy, a phenotype observed in other bacterial species including human pathogens. Moreover, this demonstrates that the types of genetic change underlying adaptation to fitness costs, and consequently the impact of evolvability mechanisms such as increased point-mutation rates, depend critically on the mechanism of resistance.

Trade-Offs of Escherichia coli Adaptation to an Intracellular Lifestyle in Macrophages

The bacterium Escherichia coli exhibits remarkable genomic and phenotypic variation, with some pathogenic strains having evolved to survive and even replicate in the harsh intra-macrophage environment. The rate and effects of mutations that can cause pathoadaptation are key determinants of the pace at which E. coli can colonize such niches and become pathogenic. We used experimental evolution to determine the speed and evolutionary paths undertaken by a commensal strain of E. coli when adapting to intracellular life. We estimated the acquisition of pathoadaptive mutations at a rate of 10⁻⁶ per genome per generation, resulting in the fixation of more virulent strains in less than a hundred generations. Whole genome sequencing of independently evolved clones showed that the main targets of intracellular adaptation involved loss of function mutations in genes implicated in the assembly of the lipopolysaccharide core, iron metabolism and di- and tri-peptide transport, namely rfaI, fhuA and tppB, respectively. We found a substantial amount of antagonistic pleiotropy in evolved populations, as well as metabolic trade-offs, commonly found in intracellular bacteria with reduced genome sizes. Overall, the low levels of clonal interference detected indicate that the first steps of the transition of a commensal E. coli into intracellular pathogens are dominated by a few pathoadaptive mutations with very strong effects.

Enhanced production of recombinant proteins with Corynebacterium glutamicum by deletion of insertion sequences (IS elements)

BACKGROUND

In most bacteria, various jumping genetic elements including insertion sequences elements (IS elements) cause a variety of genetic rearrangements resulting in harmful effects such as genome and recombinant plasmid instability. The genetic stability of a plasmid in a host is critical for high-level production of recombinant proteins, and in this regard, the development of an IS element-free strain could be a useful strategy for the enhanced production of recombinant proteins. Corynebacterium glutamicum, which is a workhorse in the industrial-scale production of various biomolecules including recombinant proteins, also has several IS elements, and it is necessary to identify the critical IS elements and to develop IS element deleted strain.

RESULTS

From the cultivation of C. glutamicum harboring a plasmid for green fluorescent protein (GFP) gene expression, non-fluorescent clones were isolated by FACS (fluorescent activated cell sorting). All the isolated clones had insertions of IS elements in the GFP coding region, and two major IS elements (ISCg1 and ISCg2 families) were identified. By co-cultivating cells harboring either the isolated IS element-inserted plasmid or intact plasmid, it was clearly confirmed that cells harboring the IS element-inserted plasmids became dominant during the cultivation due to their growth advantage over cells containing intact plasmids, which can cause a significant reduction in recombinant protein production during cultivation. To minimize the harmful effects of IS elements on the expression of heterologous genes in C. glutamicum, two IS element free C. glutamicum strains were developed in which each major IS element was deleted, and enhanced productivity in the engineered C. glutamicum strain was successfully demonstrated with three models: GFP, poly(3-hydroxybutyrate) [P(3HB)] and γ-aminobutyrate (GABA).

CONCLUSIONS

Our findings clearly indicate that the hopping of IS elements could be detrimental to the production of recombinant proteins in C. glutamicum, emphasizing the importance of developing IS element free host strains.

Dynamics of bacterial insertion sequences: can transposition bursts help the elements persist?

Background

Currently there is no satisfactory explanation for why bacterial insertion sequences (ISs) widely occur across prokaryotes despite being mostly harmful to their host genomes. Rates of horizontal gene transfer are likely to be too low to maintain ISs within a population. IS-induced beneficial mutations may be important for both prevalence of ISs and microbial adaptation to changing environments but may be too rare to sustain IS elements in the long run. Environmental stress can induce elevated rates of IS transposition activities; such episodes are known as ‘transposition bursts’. By examining how selective forces and transposition events interact to influence IS dynamics, this study asks whether transposition bursts can lead to IS persistence.

Results

We show through a simulation model that ISs are gradually eliminated from a population even if IS transpositions occasionally cause advantageous mutations. With beneficial mutations, transposition bursts create variation in IS copy numbers and improve cell fitness on average. However, these benefits are not usually sufficient to overcome the negative selection against the elements, and transposition bursts amplify the mean fitness effect which, if negative, simply accelerates the extinction of ISs. If down regulation of transposition occurs, IS extinctions are reduced while ISs still generate variation amongst bacterial genomes.

Conclusions

Transposition bursts do not help ISs persist in a bacterial population in the long run because most burst-induced mutations are deleterious and therefore not favoured by natural selection. However, bursts do create more genetic variation through which occasional advantageous mutations can help organisms adapt. Regulation of IS transposition bursts and stronger positive selection of the elements interact to slow down the burst-induced extinction of ISs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0560-5) contains supplementary material, which is available to authorized users.

'Deadman' and 'Passcode' microbial kill switches for bacterial containment

Biocontainment systems that couple environmental sensing with circuit-based control of cell viability could be used to prevent escape of genetically modified microbes into the environment. Here we present two engineered safe-guard systems: the Deadman and Passcode kill switches. The Deadman kill switch uses unbalanced reciprocal transcriptional repression to couple a specific input signal with cell survival. The Passcode kill switch uses a similar two-layered transcription design and incorporates hybrid LacI/GalR family transcription factors to provide diverse and complex environmental inputs to control circuit function. These synthetic gene circuits efficiently kill Escherichia coli and can be readily reprogrammed to change their environmental inputs, regulatory architecture and killing mechanism.

Extinction probabilities and stationary distributions of mobile genetic elements in prokaryotes: The birth-death-diversification model

Theoretical approaches are essential to our understanding of the complex dynamics of mobile genetic elements (MGEs) within genomes. Recently, the birth-death-diversification model was developed to describe the dynamics of mobile promoters (MPs), a particular class of MGEs in prokaryotes. A unique feature of this model is that genetic diversification of elements was included. To explore the implications of diversification on the longterm fate of MGE lineages, in this contribution we analyze the extinction probabilities, extinction times and equilibrium solutions of the birth-death-diversification model. We find that diversification increases both the survival and growth rate of MGE families, but the strength of this effect depends on the rate of horizontal gene transfer (HGT). We also find that the distribution of MGE families per genome is not necessarily monotonically decreasing, as observed for MPs, but may have a peak in the distribution that is related to the HGT rate. For MPs specifically, we find that new families have a high extinction probability, and predict that the number of MPs is increasing, albeit at a very slow rate. Additionally, we develop an extension of the birth-death-diversification model which allows MGEs in different regions of the genome, for example coding and non-coding, to be described by different rates. This extension may offer a potential explanation as to why the majority of MPs are located in non-promoter regions of the genome.

Transposase interaction with the β sliding clamp: effects on insertion sequence proliferation and transposition rate

Insertion sequences (ISs) are ubiquitous and abundant mobile genetic elements in prokaryotic genomes. ISs often encode only one protein, the transposase, which catalyzes their transposition. Recent studies have shown that transposases of many different IS families interact with the β sliding clamp, a DNA replication factor of the host. However, it was unclear to what extent this interaction limits or favors the ability of ISs to colonize a chromosome from a phylogenetically-distant organism, or if the strength of this interaction affects the transposition rate. Here we describe the proliferation of a member of the IS1634 family in Acidiphilium over ~600 generations of cultured growth. We demonstrate that the purified transposase binds to the β sliding clamp of Acidiphilium, Leptospirillum and E. coli. Further, we also demonstrate that the Acidiphilium IS1634 transposase binds to the archaeal sliding clamp (PCNA) from Methanosarcina, and that the transposase encoded by Methanosarcina IS1634 binds to Acidiphilium β. Finally, we demonstrate that increasing the strength of the interaction between β and transposase results in a higher transposition rate in vivo. Our results suggest that the interaction could determine the potential of ISs to be mobilized in bacterial populations and also their ability to proliferate within chromosomes.

Predicting the Genetic Stability of Engineered DNA Sequences with the EFM Calculator

Unwanted evolution can rapidly degrade the performance of genetically engineered circuits and metabolic pathways installed in living organisms. We created the Evolutionary Failure Mode (EFM) Calculator to computationally detect common sources of genetic instability in an input DNA sequence. It predicts two types of mutational hotspots: deletions mediated by homologous recombination and indels caused by replication slippage on simple sequence repeats. We tested the performance of our algorithm on genetic circuits that were previously redesigned for greater evolutionary reliability and analyzed the stability of sequences in the iGEM Registry of Standard Biological Parts. More than half of the parts in the Registry are predicted to experience >100-fold elevated mutation rates due to the inclusion of unstable sequence configurations. We anticipate that the EFM Calculator will be a useful negative design tool for avoiding volatile DNA encodings, thereby increasing the evolutionary lifetimes of synthetic biology devices.

Quorum-sensing crosstalk-driven synthetic circuits: from unimodality to trimodality

Widespread quorum-sensing (QS) enables bacteria to communicate and plays a critical role in controlling bacterial virulence. However, effects of promiscuous QS crosstalk and its implications for gene regulation and cell decision-making remain largely unknown. Here we systematically studied the crosstalk between LuxR/I and LasR/I systems and found that QS crosstalk can be dissected into signal crosstalk and promoter crosstalk. Further investigations using synthetic positive feedback circuits revealed that signal crosstalk significantly decreases a circuit's bistable potential while maintaining unimodality. Promoter crosstalk, however, reproducibly generates complex trimodal responses resulting from noise-induced state transitions and host-circuit interactions. A mathematical model that integrates the circuit's nonlinearity, stochasticity, and host-circuit interactions was developed, and its predictions of conditions for trimodality were verified experimentally. Combining synthetic biology and mathematical modeling, this work sheds light on the complex behaviors emerging from QS crosstalk, which could be exploited for therapeutics and biotechnology.

Mutational signatures indicative of environmental stress in bacteria

Evolutionary innovations are dependent on mutations. Mutation rates are increased by adverse conditions in the laboratory, but there is no evidence that stressful environments that do not directly impact on DNA leave a mutational imprint on extant genomes. Mutational spectra in the laboratory are normally determined with unstressed cells but are unavailable with stressed bacteria. To by-pass problems with viability, selection effects, and growth rate differences due to stressful environments, in this study we used a set of genetically engineered strains to identify the mutational spectrum associated with nutritional stress. The strain set members each had a fixed level of the master regulator protein, RpoS, which controls the general stress response of Escherichia coli. By assessing mutations in cycA gene from 485 cycloserine resistant mutants collected from as many independent cultures with three distinct perceived stress (RpoS) levels, we were able establish a dose-dependent relationship between stress and mutational spectra. The altered mutational patterns included base pair substitutions, single base pair indels, longer indels, and transpositions of different insertion sequences. The mutational spectrum of low-RpoS cells closely matches the genome-wide spectrum previously generated in laboratory environments, while the spectra of high RpoS, high perceived stress cells more closely matches spectra found in comparisons of extant genomes. Our results offer an explanation of the uneven mutational profiles such as the transition-transversion biases observed in extant genomes and provide a framework for assessing the contribution of stress-induced mutagenesis to evolutionary transitions and the mutational emergence of antibiotic resistance and disease states.

Birth, death, and diversification of mobile promoters in prokaryotes

A previous study of prokaryotic genomes identified large reservoirs of putative mobile promoters (PMPs), that is, homologous promoter sequences associated with nonhomologous coding sequences. Here we extend this data set to identify the full complement of mobile promoters in sequenced prokaryotic genomes. The expanded search identifies nearly 40,000 PMP sequences, 90% of which occur in noncoding regions of the genome. To gain further insight from this data set, we develop a birth-death-diversification model for mobile genetic elements subject to sequence diversification; applying the model to PMPs we are able to quantify the relative importance of duplication, loss, horizontal gene transfer (HGT), and diversification to the maintenance of the PMP reservoir. The model predicts low rates of HGT relative to the duplication and loss of PMP copies, rapid dynamics of PMP families, and a pool of PMPs that exist as a single copy in a genome at any given time, despite their mobility. We report evidence of these "singletons" at high frequencies in prokaryotic genomes. We also demonstrate that including selection, either for or against PMPs, was not necessary to describe the observed data.