Recombination frequency in plasmid DNA containing direct repeats—predictive correlation with repeat and intervening sequence length

Laboratory evolution of the bacterial genome structure through insertion sequence activation

The genome structure fundamentally shapes bacterial physiology, ecology, and evolution. Though insertion sequences (IS) are known drivers of drastic evolutionary changes in the genome structure, the process is typically slow and challenging to observe in the laboratory. Here, we developed a system to accelerate IS-mediated genome structure evolution by introducing multiple copies of a high-activity IS in Escherichia coli. We evolved the bacteria under relaxed neutral conditions, simulating those leading to IS expansion in host-restricted endosymbionts and pathogens. Strains accumulated a median of 24.5 IS insertions and underwent over 5% genome size changes within ten weeks, comparable to decades-long evolution in wild-type strains. The detected interplay of frequent small deletions and rare large duplications updates the view of genome reduction under relaxed selection from a simple consequence of the deletion bias to a nuanced picture including transient expansions. The high IS activity resulted in structural variants of IS and the emergence of composite transposons, illuminating potential evolutionary pathways for ISs and composite transposons. The extensive genome rearrangements we observed establish a baseline for assessing the fitness effects of IS insertions, genome size changes, and rearrangements, advancing our understanding of how mobile elements shape bacterial genomes.

The Causes for Genomic Instability and How to Try and Reduce Them Through Rational Design of Synthetic DNA

Genetic engineering has revolutionized our ability to manipulate DNA and engineer organisms for various applications. However, this approach can lead to genomic instability, which can result in unwanted effects such as toxicity, mutagenesis, and reduced productivity. To overcome these challenges, smart design of synthetic DNA has emerged as a promising solution. By taking into consideration the intricate relationships between gene expression and cellular metabolism, researchers can design synthetic constructs that minimize metabolic stress on the host cell, reduce mutagenesis, and increase protein yield. In this chapter, we summarize the main challenges of genomic instability in genetic engineering and address the dangers of unknowingly incorporating genomically unstable sequences in synthetic DNA. We also demonstrate the instability of those sequences by the fact that they are selected against conserved sequences in nature. We highlight the benefits of using ESO, a tool for the rational design of DNA for avoiding genetically unstable sequences, and also summarize the main principles and working parameters of the software that allow maximizing its benefits and impact.

pHAPE: a plasmid for production of DNA size marker ladders for gel electrophoresis

DNA size markers (also known as 'molecular weight markers' or 'DNA ladders') are an essential tool when using gel electrophoresis to identify and purify nucleic acids. However, the cost of these DNA ladders is not insignificant and, over time, impinges on the funds available for research and training in molecular biology. Here, we describe a method for the generation of 'pHAPE', a plasmid from which a variety of DNA ladders can be generated via simple restriction enzyme digestions. The pHAPE plasmid can be generated by mutagenesis of the commonly used pBluescript II SK+ phagemid followed by insertion of a 7141 bp sequence (comprised of three smaller, synthetic fragments). Our use of pHAPE allows us some small relief from the ever-rising costs of performing molecular biology experiments ('Don't worry, pHAPE').

Selfish, promiscuous and sometimes useful: how mobile genetic elements drive horizontal gene transfer in microbial populations

Horizontal gene transfer (HGT) drives microbial adaptation but is often under the control of mobile genetic elements (MGEs) whose interests are not necessarily aligned with those of their hosts. In general, transfer is costly to the donor cell while potentially beneficial to the recipients. The diversity and plasticity of cell–MGEs interactions, and those among MGEs, result in complex evolutionary processes where the source, or even the existence of selection for maintaining a function in the genome, is often unclear. For example, MGE-driven HGT depends on cell envelope structures and defense systems, but many of these are transferred by MGEs themselves. MGEs can spur periods of intense gene transfer by increasing their own rates of horizontal transmission upon communicating, eavesdropping, or sensing the environment and the host physiology. This may result in high-frequency transfer of host genes unrelated to the MGE. Here, we review how MGEs drive HGT and how their transfer mechanisms, selective pressures and genomic traits affect gene flow, and therefore adaptation, in microbial populations. The encoding of many adaptive niche-defining microbial traits in MGEs means that intragenomic conflicts and alliances between cells and their MGEs are key to microbial functional diversification.

This article is part of a discussion meeting issue ‘Genomic population structures of microbial pathogens’.

Evolutionary Stability Optimizer (ESO): A Novel Approach to Identify and Avoid Mutational Hotspots in DNA Sequences While Maintaining High Expression Levels

Modern synthetic biology procedures rely on the ability to generate stable genetic constructs that keep their functionality over long periods of time. However, maintenance of these constructs requires energy from the cell and thus reduces the host’s fitness. Natural selection results in loss-of-functionality mutations that negate the expression of the construct in the population. Current approaches for the prevention of this phenomenon focus on either small-scale, manual design of evolutionary stable constructs or the detection of mutational sites with unstable tendencies. We designed the Evolutionary Stability Optimizer (ESO), a software tool that enables the large-scale automatic design of evolutionarily stable constructs with respect to both mutational and epigenetic hotspots and allows users to define custom hotspots to avoid. Furthermore, our tool takes the expression of the input constructs into account by considering the guanine-cytosine (GC) content and codon usage of the host organism, balancing the trade-off between stability and gene expression, allowing to increase evolutionary stability while maintaining the high expression. In this study, we present the many features of the ESO and show that it accurately predicts the evolutionary stability of endogenous genes. The ESO was created as an easy-to-use, flexible platform based on the notion that directed genetic stability research will continue to evolve and revolutionize current applications of synthetic biology. The ESO is available at the following link: https://www.cs.tau.ac.il/~tamirtul/ESO/.

Characterization of Proteus vulgaris Strain P3M, a Foodborne Multidrug-Resistant Bacterium Isolated from Penaeus vannamei in China

Proteus vulgaris is an important foodborne opportunistic pathogen, both environmentally and clinically. The use of appropriate antibiotics has significant therapeutic effects, but has led to the emergence and spread of drug-resistant strains. In this study, a P. vulgaris strain, designated "P3M," was isolated from Penaeus vannamei in Tianjin, China. The whole genome of P3M was sequenced, generating detailed information, including the key genes involved in important metabolic pathways and their physiological functions. A total of 218 antibiotic resistance genes (ARGs) were predicted in the genome. The determination of various minimum inhibitory concentrations indicated that P3M is a multidrug-resistant (MDR) bacterium, with significant resistance to 16 antibiotics in seven categories. Determination of fractional inhibitory concentration index showed that the combination of ciprofloxacin plus tetracycline exhibited synergistic antimicrobial activity. Bioinformatics and phylogenetic analyses detected the presence of two two-component systems that mediate multidrug resistance and several mobile genetic elements involved in the horizontal transfer of ARGs in P3M. P. vulgaris strains represent a serious challenge to clinicians and infection control teams for its ubiquity worldwide and close relevance with human life. To the best of our knowledge, we report the first isolation and characterization of an important foodborne MDR P. vulgaris strain, and this study will provide necessary theoretical basis for the selection and clinical use of the appropriate antibiotics.

Functional Identification and Evolutionary Analysis of Two Novel Plasmids Mediating Quinolone Resistance in Proteus vulgaris

Plasmid-mediated quinolone resistance (PMQR) remains one of the main mechanisms of bacterial quinolone resistance and plays an important role in the transmission of antibiotic resistance genes (ARGs). In this study, two novel plasmids, p3M-2A and p3M-2B, which mediate quinolone resistance in Proteus vulgaris strain 3M (P3M) were identified. Of these, only p3M-2B appeared to be a qnrD-carrying plasmid. Both p3M-2A and p3M-2B could be transferred into Escherichia coli, and the latter caused a twofold change in ciprofloxacin resistance, according to the measured minimum inhibitory concentration (MIC). Plasmid curing/complementation and qRT-PCR results showed that p3M-2A can directly regulate the expression of qnrD in p3M-2B under treatment with ciprofloxacin, in which process, ORF1 was found to play an important role. Sequence alignments and phylogenetic analysis revealed the evolutionary relationships of all reported qnrD-carrying plasmids and showed that ORF1–4 in p3M-2B is the most conserved backbone for the normal function of qnrD-carrying plasmids. The identified direct repeats (DR) suggested that, from an evolutionary perspective, p3M-2B may have originated from the 2683-bp qnrD-carrying plasmid and may increase the possibility of plasmid recombination and then of qnrD transfer. To the best of our knowledge, this is the first identification of a novel qnrD-carrying plasmid isolated from a P. vulgaris strain of shrimp origin and a plasmid that plays a regulatory role in qnrD expression. This study also sheds new light on plasmid evolution and on the mechanism of horizontal transfer of ARGs encoded by plasmids.

Mapping Protein-Protein Interaction Interface Peptides with Jun-Fos Assisted Phage Display and Deep Sequencing

Protein-protein interactions govern many cellular processes, and identifying binding interaction sites on proteins can facilitate the discovery of inhibitors to block such interactions. Here we identify peptides from a randomly fragmented plasmid encoding the β-lactamase inhibitory protein (BLIP) and the Lac repressor (LacI) that represent regions of protein-protein interactions. We utilized a Jun-Fos-assisted phage display system that has previously been used to screen cDNA and genomic libraries to identify antibody antigens. Affinity selection with polyclonal antibodies against LacI or BLIP resulted in the rapid enrichment of in-frame peptides from various regions of the proteins. Further, affinity selection with β-lactamase enriched peptides that encompass regions of BLIP previously shown to contribute strongly to the binding energy of the BLIP/β-lactamase interaction, i.e., hotspot residues. Further, one of the regions enriched by affinity selection encompassed a disulfide-constrained region of BLIP that forms part of the BLIP interaction surface in the native complex that we show also binds to β-lactamase as a disulfide-constrained macrocycle peptide with a K_D of ∼1 μM. Fragmented open reading frame (ORF) libraries may efficiently identify such naturally constrained peptides at protein-protein interaction interfaces. With sufficiently deep coverage of ORFs by peptide-coding inserts, phage display and deep sequencing can provide detailed information on the domains or peptides that contribute to an interaction. Such information should enable the design of potentially therapeutic macrocycles or peptidomimetics that block the interaction.

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.

Advances in host and vector development for the production of plasmid DNA vaccines

Recent developments in DNA vaccine research provide a new momentum for this rather young and potentially disruptive technology. Gene-based vaccines are capable of eliciting protective immunity in humans to persistent intracellular pathogens, such as HIV, malaria, and tuberculosis, for which the conventional vaccine technologies have failed so far. The recent identification and characterization of genes coding for tumor antigens has stimulated the development of DNA-based antigen-specific cancer vaccines. Although most academic researchers consider the production of reasonable amounts of plasmid DNA (pDNA) for immunological studies relatively easy to solve, problems often arise during this first phase of production. In this chapter we review the current state of the art of pDNA production at small (shake flasks) and mid-scales (lab-scale bioreactor fermentations) and address new trends in vector design and strain engineering. We will guide the reader through the different stages of process design starting from choosing the most appropriate plasmid backbone, choosing the right Escherichia coli (E. coli) strain for production, and cultivation media and scale-up issues. In addition, we will address some points concerning the safety and potency of the produced plasmids, with special focus on producing antibiotic resistance-free plasmids. The main goal of this chapter is to make immunologists aware of the fact that production of the pDNA vaccine has to be performed with as much as attention and care as the rest of their research.

Mutation detection in plasmid-based biopharmaceuticals

As the number of applications involving therapeutic plasmid DNA (pDNA) increases worldwide, there is a growing concern over maintaining rigorous quality control through a panel of high-quality assays. For this reason, efficient, cost-effective and sensitive technologies enabling the identification of genetic variants and unwanted side products are needed to successfully establish the identity and stability of a plasmid-based biopharmaceutical. This review highlights several bioinformatic tools for ab initio detection of potentially unstable DNA regions, as well as techniques used for mutation detection in nucleic acids, with particular emphasis on pDNA.

Analysis of DNA repeats in bacterial plasmids reveals the potential for recurrent instability events

Structural instability has been frequently observed in natural plasmids and vectors used for protein expression or DNA vaccine development. However, there is a lack of information concerning hotspot mapping, namely, DNA repeats or sequences identical to the host genome. This led us to evaluate the abundance and distribution of direct, inverted, and tandem repeats with high recombination potential in 36 natural plasmids from ten bacterial genera, as well as in several widely used bacterial and mammalian expression vectors. In natural plasmids, we observed an overrepresentation of close direct repeats in comparison to inverted ones and a preferential location of repeats with high recombination potential in intergenic regions, suggesting a highly plastic and dynamic behavior. In plasmid vectors, we found a high density of repeats within eukaryotic promoters and non-coding sequences. As a result of this in silico analysis, we detected a spontaneous recombination between two 21-bp direct repeats present in the human cytomegalovirus early enhancer/promoter (huCMV EEP) of the pCIneo plasmid. This finding is of particular importance, as the huCMV EEP is one of the most frequently used regulatory elements in plasmid vectors. Because pDNA integration into host gDNA can have adverse consequences in terms of plasmid processing and host safety, we also mapped several regions with high probability to mediate integration into the Escherichia coli or human genomes. Like repeated regions, some of these were located in non-coding regions of the plasmids, thus being preferential targets to be removed.

Effects of environmental stress on stability of tandem repeats in Escherichia coli O157:H7

Multilocus variable-number tandem-repeat analysis (MLVA) is used for source tracking Escherichia coli O157:H7 in agricultural environments. Tandem repeats were stable after limited replication but changed after exposure to irradiation, elevated temperatures, and starvation conditions. The pO157 plasmid was frequently lost under these stress conditions. Environmental stresses may increase phylogenetic diversity as measured by MLVA.

Genesis, effects and fates of repeats in prokaryotic genomes

DNA repeats are causes and consequences of genome plasticity. Repeats are created by intrachromosomal recombination or horizontal transfer. They are targeted by recombination processes leading to amplifications, deletions and rearrangements of genetic material. The identification and analysis of repeats in nearly 700 genomes of bacteria and archaea is facilitated by the existence of sequence data and adequate bioinformatic tools. These have revealed the immense diversity of repeats in genomes, from those created by selfish elements to the ones used for protection against selfish elements, from those arising from transient gene amplifications to the ones leading to stable duplications. Experimental works have shown that some repeats do not carry any adaptive value, while others allow functional diversification and increased expression. All repeats carry some potential to disorganize and destabilize genomes. Because recombination and selection for repeats vary between genomes, the number and types of repeats are also quite diverse and in line with ecological variables, such as host-dependent associations or population sizes, and with genetic variables, such as the recombination machinery. From an evolutionary point of view, repeats represent both opportunities and problems. We describe how repeats are created and how they can be found in genomes. We then focus on the functional and genomic consequences of repeats that dictate their fate.

High frequency plasmid recombination mediated by 28 bp direct repeats

The stability in Escherichia coli of a mammalian expression vector (pCIneo) and its derivative candidate DNA vaccine (pGPV-PV) is described. These multicopy pMB1-type plasmids were found to recombine in several recA E. coli strains due to the presence of two 28 bp direct repeats flanking intervening sequences of 1.6 kb (pCIneo) and 3.2 kb (pGPV-PV). In this recombination event, one of the direct repeats and the intervening sequence were deleted or duplicated, originating monomeric or/and hetero-dimeric plasmid forms, respectively. Additionally, the plasmid rearrangement led to the acquisition of a kanamycin resistance phenotype. Recombination frequencies between 7.8 x 10(-7) and 3.1 x 10(-5) were determined for DH5alpha and JM109(DE3) strains, respectively. Higher recombination frequencies were found in cells previously grown up to stationary growth phase being the monomeric plasmid form the prevalent one. Real-time PCR quantification revealed the presence of approximately 1.5 x 10(4) recombined molecules per 2 x 10(5 )cells pre-kanamycin exposure. Under selective pressure of this antibiotic, the number of recombined molecules increased approximately 2,000-fold in a 48-h period replacing the original plasmid form. The high frequency, at which deletion-duplication occurred in the absence of kanamycin selective pressure, should be regarded as a safety concern. This work highlights the impact of mutational hot spots on expression and cloning plasmid vectors and the need to carefully design plasmid vectors.