Insertion Sequence-Driven Evolution of Escherichia coli in Chemostats

Influence of insertion sequences on population structure of phytopathogenic bacteria in the Ralstonia solanacearum species complex

Ralstonia solanacearum species complex (RSSC) is a destructive group of plant pathogenic bacteria and the causative agent of bacterial wilt disease. Experimental studies have attributed RSSC virulence to insertion sequences (IS), transposable genetic elements which can both disrupt and activate host genes. Yet, the global diversity and distribution of RSSC IS are unknown. In this study, IS were bioinformatically identified in a diverse collection of 356 RSSC isolates representing five phylogenetic lineages and their diversity investigated based on genetic distance measures and comparisons with the ISFinder database. IS phylogenetic associations were determined based on their distribution across the RSSC phylogeny. Moreover, IS positions within genomes were characterised and their potential gene disruptions determined based on IS proximity to coding sequences. In total, we found 24732 IS belonging to eleven IS families and 26 IS subgroups with over half of the IS found in the megaplasmid. While IS families were generally widespread across the RSSC phylogeny, IS subgroups showed strong lineage-specific distributions and genetically similar bacterial isolates had similar IS contents. Similar associations with bacterial host genetic background were also observed with IS insertion positions which were highly conserved in closely related bacterial isolates. Finally, IS were found to disrupt genes with predicted functions in virulence, stress tolerance, and metabolism suggesting that they might be adaptive. This study highlights that RSSC insertion sequences track the evolution of their bacterial hosts potentially contributing to both intra- and inter-lineage genetic diversity.

Genomic landscapes of bacterial transposons and their applications in strain improvement

Transposons are mobile genetic elements that can give rise to gene mutation and genome rearrangement. Due to their mobility, transposons have been exploited as genetic tools for modification of plants, animals, and microbes. Although a plethora of reviews have summarized families of transposons, the transposons from fermentation bacteria have not been systematically documented, which thereby constrain the exploitation for metabolic engineering and synthetic biology purposes. In this review, we summarize the transposons from the most used fermentation bacteria including Escherichia coli, Bacillus subtilis, Lactococcus lactis, Corynebacterium glutamicum, Klebsiella pneumoniae, and Zymomonas mobilis by literature retrieval and data mining from GenBank and KEGG. We also outline the state-of-the-art advances in basic research and industrial applications especially when allied with other genetic tools. Overall, this review aims to provide valuable insights for transposon-mediated strain improvement. KEY POINTS: • The transposons from the most-used fermentation bacteria are systematically summarized. • The applications of transposons in strain improvement are comprehensively reviewed.

The dynamic network of IS30 transposition pathways

The E. coli element IS30 has adopted the copy-out-paste-in transposition mechanism that is prevalent in a number of IS-families. As an initial step, IS30 forms free circular transposition intermediates like IS minicircles or tandem IS-dimers by joining the inverted repeats of a single element or two, sometimes distantly positioned IS copies, respectively. Then, the active IR-IR junction of these intermediates reacts with the target DNA, which generates insertions, deletions, inversions or cointegrates. The element shows dual target specificity as it can insert into hot spot sequences or next to its inverted repeats. In this study the pathways of rearrangements of transposition-derived cointegrate-like structures were examined. The results showed that the probability of further rearrangements in these structures depends on whether the IS elements are flanked by hot spot sequences or take part in an IR-IR junction. The variability of the deriving products increases with the number of simultaneously available IRs and IR-IR joints in the cointegrates or the chromosome. Under certain conditions, the parental structures whose transposition formed the cointegrates are restored and persist among the rearranged products. Based on these findings, a novel dynamic model has been proposed for IS30, which possibly fits to other elements that have adopted the same transposition mechanism. The model integrates the known transposition pathways and the downstream rearrangements occurring after the formation of different cointegrate-like structures into a complex network. Important feature of this network is the presence of “feedback loops” and reversible transposition rearrangements that can explain how IS30 generates variability and preserves the original genetic constitution in the bacterial population, which contributes to the adaptability and evolution of host bacteria.

A synergistic arrangement of two unrelated IS elements facilitates adjacent deletion in Micrococcus luteus ATCC49732

Mutants of M. luteus strain ATCC49732 lacking the yellow pigment sarcinaxanthin were observed at an unexpectedly high frequency and the molecular basis was investigated. PCR probing revealed complete deletion of the crt biosynthetic operon in 11/14 mutants. Inverse PCR was used to identify a common breakpoint 35 kb downstream from crt precisely at the end of the right inverted repeat (IRR) of a partial ISMlu8 element that lies between two inversely oriented full-length ISMlu2. Three different breakpoints 5′ to crt were found with the sequence CTAG one bp 5′ to each novel junction. Analysis of 35 genomic sites with single ISMlu8 insertions showed that ISMlu8 transposase has high specificity for CTAG, implicating its key role in formation of the Δcrt deletions. No downstream deletion endpoints were observed at an immediately adjacent ISMlu8 with a nearly identical IRR in the same orientation and slightly closer to the crt operon, indicating that access of ISMlu8 transposase to the ISMlu2-flanked ISMlu8 IRR is greatly enhanced by the surrounding inverted repeat arrangement. The association of high frequency genomic rearrangement with this distinctive natural configuration of ISs from two different IS families offers a new insight into IS element evolutionary potential.

Beyond to the Stable: Role of the Insertion Sequences as Epidemiological Descriptors in Corynebacterium striatum

In recent years, epidemiological studies of infectious agents have focused mainly on the pathogen and stable components of its genome. The use of these stable components makes it possible to know the evolutionary or epidemiological relationships of the isolates of a particular pathogen. Under this approach, focused on the pathogen, the identification of resistance genes is a complementary stage of a bacterial characterization process or an appendix of its epidemiological characterization, neglecting its genetic components’ acquisition or dispersal mechanisms. Today we know that a large part of antibiotic resistance is associated with mobile elements. Corynebacterium striatum, a bacterium from the normal skin microbiota, is also an opportunistic pathogen. In recent years, reports of infections and nosocomial outbreaks caused by antimicrobial multidrug-resistant C. striatum strains have been increasing worldwide. Despite the different existing mobile genomic elements, there is evidence that acquired resistance genes are coupled to insertion sequences in C. striatum. This perspective article reviews the insertion sequences linked to resistance genes, their relationship to evolutionary lineages, epidemiological characteristics, and the niches the strains inhabit. Finally, we evaluate the potential of the insertion sequences for their application as a descriptor of epidemiological scenarios, allowing us to anticipate the emergence of multidrug-resistant lineages.

Expanding Characterized Diversity and the Pool of Complete Genome Sequences of Methylococcus Species, the Bacteria of High Environmental and Biotechnological Relevance

The bacterial genus Methylococcus, which comprises aerobic thermotolerant methanotrophic cocci, was described half-a-century ago. Over the years, a member of this genus, Methylococcus capsulatus Bath, has become a major model organism to study genomic and metabolic basis of obligate methanotrophy. High biotechnological potential of fast-growing Methylococcus species, mainly as a promising source of feed protein, has also been recognized. Despite this big research attention, the currently cultured Methylococcus diversity is represented by members of the two species, M. capsulatus and M. geothermalis, while finished genome sequences are available only for two strains of these methanotrophs. This study extends the pool of phenotypically characterized Methylococcus strains with good-quality genome sequences by contributing four novel isolates of these bacteria from activated sludge, landfill cover soil, and freshwater sediments. The determined genome sizes of novel isolates varied between 3.2 and 4.0Mb. As revealed by the phylogenomic analysis, strains IO1, BH, and KN2 affiliate with M. capsulatus, while strain Mc7 may potentially represent a novel species. Highest temperature optima (45-50°C) and highest growth rates in bioreactor cultures (up to 0.3h-1) were recorded for strains obtained from activated sludge. The comparative analysis of all complete genomes of Methylococcus species revealed 4,485 gene clusters. Of these, pan-genome core comprised 2,331 genes (on average 51.9% of each genome), with the accessory genome containing 846 and 1,308 genes in the shell and the cloud, respectively. Independently of the isolation source, all strains of M. capsulatus displayed surprisingly high genome synteny and a striking similarity in gene content. Strain Mc7 from a landfill cover soil differed from other isolates by the high content of mobile genetic elements in the genome and a number of genome-encoded features missing in M. capsulatus, such as sucrose biosynthesis and the ability to scavenge phosphorus and sulfur from the environment.

Mobile genetic elements mediate the mixotrophic evolution of novel Alicyclobacillus species for acid mine drainage adaptation

Alicyclobacillus species inhabit diverse environments and have adapted to broad ranges of pH and temperature. However, their adaptive evolutions remain elusive, especially regarding the role of mobile genetic elements (MGEs). Here, we characterized the distributions and functions of MGEs in Alicyclobacillus species across five environments, including acid mine drainage (AMD), beverages, hot springs, sediments, and soils. Nine Alicyclobacillus strains were isolated from AMD and possessed larger genome sizes and more genes than those from other environments. Four AMD strains evolved to be mixotrophic and fell into distinctive clusters in phylogenetic tree. Four types of MGEs including genomic island (GI), insertion sequence (IS), prophage, and integrative and conjugative element (ICE) were widely distributed in Alicyclobacillus species. Further, AMD strains did not possess CRISPR-Cas systems, but had more GI, IS, and ICE, as well as more MGE-associated genes involved in the oxidation of iron and sulfide and the resistance of heavy metal and low temperature. These findings highlight the differences in phenotypes and genotypes between strains isolated from AMD and other environments and the important role of MGEs in rapid environment niche expansions.

A glimpse of potential transposable element impact on adaptation of the industrial yeast Saccharomyces cerevisiae

The adapted industrial yeast strain Saccharomyces cerevisiae NRRL Y-50049 is able to in situ detoxify major toxic aldehyde compounds derived from sugar conversion of lignocellulosic biomass while producing ethanol. Pathway-based studies on its mechanisms of tolerance have been reported previously, however, little is known about transposable element (TE) involvement in its adaptation to inhibitory compounds. This work presents a comparative dynamic transcription expression analysis in response to a toxic treatment between Y-50049 and its progenitor, an industrial type strain NRRL Y-12632, using a time-course study. At least 77 TEs from Y-50049 showed significantly increased expression compared with its progenitor, especially during the late lag phase. Sequence analysis revealed significant differences in TE sequences between the two strains. Y-50049 was also found to have a transposons of yeast 2 (Ty2) long terminal repeat-linked YAT1 gene showing significantly higher copy number changes than its progenitor. These results raise awareness of potential TE involvement in the adaptation of industrial yeast to the tolerance of toxic chemicals.

Genome rearrangements and megaplasmid loss in the filamentous bacterium Kitasatospora viridifaciens are associated with protoplast formation and regeneration

Filamentous Actinobacteria are multicellular bacteria with linear replicons. Kitasatospora viridifaciens DSM 40239 contains a linear 7.8 Mb chromosome and an autonomously replicating plasmid KVP1 of 1.7 Mb. Here we show that lysozyme-induced protoplast formation of the multinucleated mycelium of K. viridifaciens drives morphological diversity. Characterisation and sequencing of an individual revertant colony that had lost the ability to differentiate revealed that the strain had not only lost most of KVP1 but also carried deletions in the right arm of the chromosome. Strikingly, the deletion sites were preceded by insertion sequence elements, suggesting that the rearrangements may have been caused by replicative transposition and homologous recombination between both replicons. These data indicate that protoplast formation is a stressful process that can lead to profound genetic changes.

Evolutionary Changes after Translational Challenges Imposed by Horizontal Gene Transfer

Genes acquired by horizontal gene transfer (HGT) may provide the recipient organism with potentially new functions, but proper expression level and integration of the transferred genes in the novel environment are not granted. Notably, transferred genes can differ from the receiving genome in codon usage preferences, leading to impaired translation and reduced functionality. Here, we characterize the genomic and proteomic changes undergone during experimental evolution of Escherichia coli after HGT of three synonymous versions, presenting very different codon usage preference, of an antibiotic resistance gene. The experimental evolution was conducted with and without the corresponding antibiotic and the mutational patterns and proteomic profiles after 1,000 generations largely depend on the experimental growth conditions (e.g., mutations in antibiotic off-target genes), and on the synonymous gene version transferred (e.g., mutations in genes responsive to translational stress). The transfer of an exogenous gene extensively modifies the whole proteome, and these proteomic changes are different for the different version of the transferred gene. Additionally, we identified conspicuous changes in global regulators and in intermediate metabolism, confirmed the evolutionary ratchet generated by mutations in DNA repair genes and highlighted the plasticity of bacterial genomes accumulating large and occasionally transient duplications. Our results support a central role of HGT in fuelling evolution as a powerful mechanism promoting rapid, often dramatic genotypic and phenotypic changes. The profound reshaping of the pre-existing geno/phenotype allows the recipient bacteria to explore new ways of functioning, far beyond the mere acquisition of a novel function.

Insertion sequences drive the emergence of a highly adapted human pathogen

Pseudomonas aeruginosa is a highly adaptive opportunistic pathogen that can have serious health consequences in patients with lung disorders. Taxonomic outliers of P. aeruginosa of environmental origin have recently emerged as infectious for humans. Here, we present the first genome-wide analysis of an isolate that caused fatal haemorrhagic pneumonia. In two clones, CLJ1 and CLJ3, sequentially recovered from a patient with chronic pulmonary disease, insertion of a mobile genetic element into the P. aeruginosa chromosome affected major virulence-associated phenotypes and led to increased resistance to the antibiotics used to combat the infection. Comparative genome, proteome and transcriptome analyses revealed that this ISL3-family insertion sequence disrupted the genes for flagellar components, type IV pili, O-specific antigens, translesion polymerase and enzymes producing hydrogen cyanide. Seven-fold more insertions were detected in the later isolate, CLJ3, than in CLJ1, some of which modified strain susceptibility to antibiotics by disrupting the genes for the outer-membrane porin OprD and the regulator of β-lactamase expression AmpD. In the Galleria mellonella larvae model, the two strains displayed different levels of virulence, with CLJ1 being highly pathogenic. This study revealed insertion sequences to be major players in enhancing the pathogenic potential of a P. aeruginosa taxonomic outlier by modulating both its virulence and its resistance to antimicrobials, and explains how this bacterium adapts from the environment to a human host.

Experimental Design, Population Dynamics, and Diversity in Microbial Experimental Evolution

In experimental evolution, laboratory-controlled conditions select for the adaptation of species, which can be monitored in real time. Despite the current popularity of such experiments, nature's most pervasive biological force was long believed to be observable only on time scales that transcend a researcher's life-span, and studying evolution by natural selection was therefore carried out solely by comparative means. Eventually, microorganisms' propensity for fast evolutionary changes proved us wrong, displaying strong evolutionary adaptations over a limited time, nowadays massively exploited in laboratory evolution experiments. Here, we formulate a guide to experimental evolution with microorganisms, explaining experimental design and discussing evolutionary dynamics and outcomes and how it is used to assess ecoevolutionary theories, improve industrially important traits, and untangle complex phenotypes. Specifically, we give a comprehensive overview of the setups used in experimental evolution. Additionally, we address population dynamics and genetic or phenotypic diversity during evolution experiments and expand upon contributing factors, such as epistasis and the consequences of (a)sexual reproduction. Dynamics and outcomes of evolution are most profoundly affected by the spatiotemporal nature of the selective environment, where changing environments might lead to generalists and structured environments could foster diversity, aided by, for example, clonal interference and negative frequency-dependent selection. We conclude with future perspectives, with an emphasis on possibilities offered by fast-paced technological progress. This work is meant to serve as an introduction to those new to the field of experimental evolution, as a guide to the budding experimentalist, and as a reference work to the seasoned expert.

Instability of the 16S rRNA methyltransferase-encoding npmA gene: why have bacterial cells possessing npmA not spread despite their high and broad resistance to aminoglycosides?

The NpmA bacterial 16S rRNA methyltransferase, which is identified from Escherichia coli strains, confers high resistance to many types of aminoglycoside upon its host cells. But despite its resistance-conferring ability, only two cases of its isolation from E. coli (14 years apart) have been reported to date. Here, we investigated the effect of the npmA gene on aminoglycoside resistance in Pseudomonas aeruginosa and Klebsiella pneumoniae and its stability in E. coli cells by comparing it with armA, another 16S rRNA methyltransferase gene currently spreading globally. As a result, we found that npmA conferred resistance to all types of aminoglycoside antibiotics we tested (except streptomycin) in both P. aeruginosa and K. pneumoniae, as well in E. coli. In addition, co-expression of armA and npmA resulted in an additive effect for the resistance. However, in return for the resistance, we also observed that the growth rates and the cell survivability of the strains transformed with the npmA-harboring plasmids were inferior than those of the control strains and that these plasmids were easily disrupted by IS10, IS1, and IS5 insertion sequences. We discuss these data in the context of the threat posed by pathogenic strains possessing npmA.

Co-evolution with recombination affects the stability of mobile genetic element insertions within gene families of Salmonella

Bacteria can have multiple copies of a gene at separate locations on the same chromosome. Some of these gene families, including tuf (translation elongation factor EF-Tu) and rrl (ribosomal RNA), encode functions critically important for bacterial fitness. Genes within these families are known to evolve in concert using homologous recombination to transfer genetic information from one gene to another. This mechanism can counteract the detrimental effects of nucleotide sequence divergence over time. Whether such mechanisms can also protect against the potentially lethal effects of mobile genetic element insertion is not well understood. To address this we constructed two different length insertion cassettes to mimic mobile genetic elements and inserted these into various positions of the tuf and rrl genes. We measured rates of recombinational repair that removed the inserted cassette and studied the underlying mechanism. Our results indicate that homologous recombination can protect the tuf and rrl genes from inactivation by mobile genetic elements, but for insertions within shorter gene sequences the efficiency of repair is very low. Intriguingly, we found that physical distance separating genes on the chromosome directly affects the rate of recombinational repair suggesting that relative location will influence the ability of homologous recombination to maintain homogeneity.

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.

Dynamics and genetic diversification of Escherichia coli during experimental adaptation to an anaerobic environment

Background

Many bacteria are facultative anaerobes, and can proliferate in both anoxic and oxic environments. Under anaerobic conditions, fermentation is the primary means of energy generation in contrast to respiration. Furthermore, the rates and spectra of spontaneous mutations that arise during anaerobic growth differ to those under aerobic growth. A long-term selection experiment was undertaken to investigate the genetic changes that underpin how the facultative anaerobe, Escherichia coli, adapts to anaerobic environments.

Methods

Twenty-one populations of E. coli REL4536, an aerobically evolved 10,000th generation descendent of the E. coli B strain, REL606, were established from a clonal ancestral culture. These were serially sub-cultured for 2,000 generations in a defined minimal glucose medium in strict aerobic and strict anaerobic environments, as well as in a treatment that fluctuated between the two environments. The competitive fitness of the evolving lineages was assessed at approximately 0, 1,000 and 2,000 generations, in both the environment of selection and the alternative environment. Whole genome re-sequencing was performed on random colonies from all lineages after 2,000-generations. Mutations were identified relative to the ancestral genome, and based on the extent of parallelism, traits that were likely to have contributed towards adaptation were inferred.

Results

There were increases in fitness relative to the ancestor among anaerobically evolved lineages when tested in the anaerobic environment, but no increases were found in the aerobic environment. For lineages that had evolved under the fluctuating regime, relative fitness increased significantly in the anaerobic environment, but did not increase in the aerobic environment. The aerobically-evolved lineages did not increase in fitness when tested in either the aerobic or anaerobic environments. The strictly anaerobic lineages adapted more rapidly to the anaerobic environment than did the fluctuating lineages. Two main strategies appeared to predominate during adaptation to the anaerobic environment: modification of energy generation pathways, and inactivation of non-essential functions. Fermentation pathways appeared to alter through selection for mutations in genes such as nadR, adhE, dcuS/R, and pflB. Mutations were frequently identified in genes for presumably dispensable functions such as toxin-antitoxin systems, prophages, virulence and amino acid transport. Adaptation of the fluctuating lineages to the anaerobic environments involved mutations affecting traits similar to those observed in the anaerobically evolved lineages.

Discussion

There appeared to be strong selective pressure for activities that conferred cell yield advantages during anaerobic growth, which include restoring activities that had previously been inactivated under long-term continuous aerobic evolution of the ancestor.

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.

Regulation of solvent tolerance in Pseudomonas putida S12 mediated by mobile elements

Organic solvent‐tolerant bacteria are outstanding and versatile hosts for the bio‐based production of a broad range of generally toxic aromatic compounds. The energetically costly solvent tolerance mechanisms are subject to multiple levels of regulation, involving among other mobile genetic elements. The genome of the solvent‐tolerant Pseudomonas putida S12 contains many such mobile elements that play a major role in the regulation and adaptation to various stress conditions, including the regulation of expression of the solvent efflux pump SrpABC. We recently sequenced the genome of P. putida S12. Detailed annotation identified a threefold higher copy number of the mobile element ISS12 in contrast to earlier observations. In this study, we describe the mobile genetic elements and elaborate on the role of ISS12 in the establishment and maintenance of solvent tolerance in P. putida. We identified three different variants of ISS12 of which a single variant exhibits a high translocation rate. One copy of this variant caused a loss of solvent tolerance in the sequenced strain by disruption of srpA. Solvent tolerance could be restored by applying selective pressure, leading to a clean excision of the mobile element.

CHO Cells adapted to inorganic phosphate limitation show higher growth and higher pyruvate carboxylase flux in phosphate replete conditions

Inorganic phosphate (P_i ) is an essential ion involved in diverse cellular processes including metabolism. Changes in cellular metabolism upon long term adaptation to P_i limitation have been reported in E. coli. Given the essential role of P_i , adaptation to P_i limitation may also result in metabolic changes in animal cells. In this study, we have adapted CHO cells producing recombinant IgG to limiting P_i conditions for 75 days. Not surprisingly, adapted cells showed better survival under P_i limitation. Here, we report the finding that such cells also showed better growth characteristics compared to control in batch culture replete with P_i (higher peak density and integral viable cell density), accompanied by a lower specific oxygen uptake rate and cytochrome oxidase activity towards the end of exponential phase. Surprisingly, the adapted cells grew to a lower peak density under glucose limitation. This suggests long term P_i limitation may lead to selection for an altered metabolism with higher dependence on glucose availability for biomass assimilation compared to control. Steady state U-¹³ C glucose labeling experiments suggest that adapted cells have a higher pyruvate carboxylase flux. Consistent with this observation, supplementation with aspartate abolished the peak density difference whereas supplementation with serine did not abolish the difference. This supports the hypothesis that cell growth in the adapted culture might be higher due to a higher pyruvate carboxylase flux. Decreased fitness under carbon limitation and mutations in the sucABCD operon has been previously reported in E. coli upon long term adaptation to P_i limitation, suggestive of a similarity in cellular response among such diverse species. © 2017 American Institute of Chemical Engineers Biotechnol. Prog., 33:749-758, 2017.

Quantitative assessment of insertion sequence impact on bacterial genome architecture

Insertion sequence (IS) elements are important mediators of genome plasticity and can lead to phenotypic changes with evolutionary significance. In multidrug-resistant Acinetobacter baumannii and Klebsiella pneumoniae, IS elements have contributed significantly to the mobilization of genes that encode resistance to antimicrobial drugs. A systematic analysis of IS elements is needed for a more comprehensive understanding of their evolutionary impact. We developed a computational approach (ISseeker) to annotate IS elements in draft genome assemblies and applied the method to analysis of IS elements in all publicly available A. baumannii(>1000) and K. pneumoniae(>800) genome sequences, in a phylogenetic context. Most IS elements in A. baumanniigenomes are species-specific ISAba elements, whereas K. pneumoniaegenomes contain significant numbers of both ISKpn elements and elements that are found throughout the Enterobacteriaceae. A. baumanniigenomes have a higher density of IS elements than K. pneumoniae, averaging ~33 vs ~27 copies per genome. In K. pneumoniae, several insertion sites are shared by most genomes in the ST258 clade, whereas in A. baumannii, different IS elements are abundant in different phylogenetic groups, even among closely related Global Clone 2 strains. IS elements differ in the distribution of insertion locations relative to genes, with some more likely to disrupt genes and others predominantly in intergenic regions. Several genes and intergenic regions had multiple independent insertion events, suggesting that those events may confer a selective advantage. Genome- and taxon-wide characterization of insertion locations revealed that IS elements have been active contributors to genome diversity in both species.

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.

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.

Simple phenotypic sweeps hide complex genetic changes in populations

Changes in allele frequencies and the fixation of beneficial mutations are central to evolution. The precise relationship between mutational and phenotypic sweeps is poorly described however, especially when multiple alleles are involved. Here, we investigate these relationships in a bacterial population over 60 days in a glucose-limited chemostat in a large population. High coverage metagenomic analysis revealed a disconnection between smooth phenotypic sweeps and the complexity of genetic changes in the population. Phenotypic adaptation was due to convergent evolution and involved soft sweeps by 7–26 highly represented alleles of several genes in different combinations. Allele combinations spread from undetectably low baselines, indicating that minor subpopulations provide the basis of most innovations. A hard sweep was also observed, involving a single combination of rpoS, mglD, malE, sdhC, and malT mutations sweeping to greater than 95% of the population. Other mutant genes persisted but at lower abundance, including hfq, consistent with its demonstrated frequency-dependent fitness under glucose limitation. Other persistent, newly identified low-frequency mutations were in the aceF, galF, ribD and asm genes, in noncoding regulatory regions, three large indels and a tandem duplication; these were less affected by fluctuations involving more dominant mutations indicating separate evolutionary paths. Our results indicate a dynamic subpopulation structure with a minimum of 42 detectable mutations maintained over 60 days. We also conclude that the massive population-level mutation supply in combination with clonal interference leads to the soft sweeps observed, but not to the exclusion of an occasional hard sweep.

The functional basis of adaptive evolution in chemostats

Two of the central problems in biology are determining the molecular basis of adaptive evolution and understanding how cells regulate their growth. The chemostat is a device for culturing cells that provides great utility in tackling both of these problems: it enables precise control of the selective pressure under which organisms evolve and it facilitates experimental control of cell growth rate. The aim of this review is to synthesize results from studies of the functional basis of adaptive evolution in long-term chemostat selections using Escherichia coli and Saccharomyces cerevisiae. We describe the principle of the chemostat, provide a summary of studies of experimental evolution in chemostats, and use these studies to assess our current understanding of selection in the chemostat. Functional studies of adaptive evolution in chemostats provide a unique means of interrogating the genetic networks that control cell growth, which complements functional genomic approaches and quantitative trait loci (QTL) mapping in natural populations. An integrated approach to the study of adaptive evolution that accounts for both molecular function and evolutionary processes is critical to advancing our understanding of evolution. By renewing efforts to integrate these two research programs, experimental evolution in chemostats is ideally suited to extending the functional synthesis to the study of genetic networks.

Large chromosomal rearrangements during a long-term evolution experiment with Escherichia coli

Large-scale rearrangements may be important in evolution because they can alter chromosome organization and gene expression in ways not possible through point mutations. In a long-term evolution experiment, twelve Escherichia coli populations have been propagated in a glucose-limited environment for over 25 years. We used whole-genome mapping (optical mapping) combined with genome sequencing and PCR analysis to identify the large-scale chromosomal rearrangements in clones from each population after 40,000 generations. A total of 110 rearrangement events were detected, including 82 deletions, 19 inversions, and 9 duplications, with lineages having between 5 and 20 events. In three populations, successive rearrangements impacted particular regions. In five populations, rearrangements affected over a third of the chromosome. Most rearrangements involved recombination between insertion sequence (IS) elements, illustrating their importance in mediating genome plasticity. Two lines of evidence suggest that at least some of these rearrangements conferred higher fitness. First, parallel changes were observed across the independent populations, with ~65% of the rearrangements affecting the same loci in at least two populations. For example, the ribose-utilization operon and the manB-cpsG region were deleted in 12 and 10 populations, respectively, suggesting positive selection, and this inference was previously confirmed for the former case. Second, optical maps from clones sampled over time from one population showed that most rearrangements occurred early in the experiment, when fitness was increasing most rapidly. However, some rearrangements likely occur at high frequency and may have simply hitchhiked to fixation. In any case, large-scale rearrangements clearly influenced genomic evolution in these populations.

Bacterial chromosomes are dynamic structures shaped by long histories of evolution. Among genomic changes, large-scale DNA rearrangements can have important effects on the presence, order, and expression of genes. Whole-genome sequencing that relies on short DNA reads cannot identify all large-scale rearrangements. Therefore, deciphering changes in the overall organization of genomes requires alternative methods, such as optical mapping. We analyzed the longest-running microbial evolution experiment (more than 25 years of evolution in the laboratory) by optical mapping, genome sequencing, and PCR analyses. We found multiple large genome rearrangements in all 12 independently evolving populations. In most cases, it is unclear whether these changes were beneficial themselves or, alternatively, hitchhiked to fixation with other beneficial mutations. In any case, many genome rearrangements accumulated over decades of evolution, providing these populations with genetic plasticity reminiscent of that observed in some pathogenic bacteria.

Evidence that the insertion events of IS2 transposition are biased towards abrupt compositional shifts in target DNA and modulated by a diverse set of culture parameters

Insertion specificity of mobile genetic elements is a rather complex aspect of DNA transposition, which, despite much progress towards its elucidation, still remains incompletely understood. We report here the results of a meta-analysis of IS2 target sites from genomic, phage, and plasmid DNA and find that newly acquired IS2 elements are consistently inserted around abrupt DNA compositional shifts, particularly in the form of switch sites of GC skew. The results presented in this study not only corroborate our previous observations that both the insertion sequence (IS) minicircle junction and target region adopt intrinsically bent conformations in IS2, but most interestingly, extend this requirement to other families of IS elements. Using this information, we were able to pinpoint regions with high propensity for transposition and to predict and detect, de novo, a novel IS2 insertion event in the 3' region of the gfp gene of a reporter plasmid. We also found that during amplification of this plasmid, process parameters such as scale, culture growth phase, and medium composition exacerbate IS2 transposition, leading to contamination levels with potentially detrimental clinical effects. Overall, our findings provide new insights into the role of target DNA structure in the mechanism of transposition of IS elements and extend our understanding of how culture conditions are a relevant factor in the induction of genetic instability.

DnaA dynamics could be linked with fitness cost in bacteria

Fitness cost is the measure of the metabolic burden of unneeded gene expression. It is defined as the lag in bacterial cells growth harboring unneeded genes relative to unburdened cells. Separate cells can concurrently adapt to the burden, demonstrating a decrease in or even a disappearance of the lag. The precise mechanisms of this adaptation are not clearly understood. One possibility is that an increased amount of free ribosomes "absorb" the unnecessary burden. In this work, the mechanism by which an increased concentration of ribosomes could result in faster growth and mask the unneeded gene expression burden is discussed. The initiation time of chromosome replication by the initiator protein DnaA, for which the accumulation speed depends on the ribosomes amount, is taken into account.

The first steps of adaptation of Escherichia coli to the gut are dominated by soft sweeps

The accumulation of adaptive mutations is essential for survival in novel environments. However, in clonal populations with a high mutational supply, the power of natural selection is expected to be limited. This is due to clonal interference - the competition of clones carrying different beneficial mutations - which leads to the loss of many small effect mutations and fixation of large effect ones. If interference is abundant, then mechanisms for horizontal transfer of genes, which allow the immediate combination of beneficial alleles in a single background, are expected to evolve. However, the relevance of interference in natural complex environments, such as the gut, is poorly known. To address this issue, we have developed an experimental system which allows to uncover the nature of the adaptive process as Escherichia coli adapts to the mouse gut. This system shows the invasion of beneficial mutations in the bacterial populations and demonstrates the pervasiveness of clonal interference. The observed dynamics of change in frequency of beneficial mutations are consistent with soft sweeps, where different adaptive mutations with similar phenotypes, arise repeatedly on different haplotypes without reaching fixation. Despite the complexity of this ecosystem, the genetic basis of the adaptive mutations revealed a striking parallelism in independently evolving populations. This was mainly characterized by the insertion of transposable elements in both coding and regulatory regions of a few genes. Interestingly, in most populations we observed a complete phenotypic sweep without loss of genetic variation. The intense clonal interference during adaptation to the gut environment, here demonstrated, may be important for our understanding of the levels of strain diversity of E. coli inhabiting the human gut microbiota and of its recombination rate.

Adaptation to novel environments involves the accumulation of beneficial mutations. If these are rare the process will proceed slowly with each one sweeping to fixation on its own. On the contrary if they are common in clonal populations, individuals carrying different beneficial alleles will experience intense competition and only those clones carrying the stronger effect mutations will leave a future line of descent. This phenomenon is known as clonal interference and the extent to which it occurs in natural environments is unknown. One of the most complex natural environments for E. coli is the mammalian intestine, where it evolves in the presence of many species comprising the gut microbiota. We have studied the dynamics of adaptation of E. coli populations evolving in this relevant ecosystem. We show that clonal interference is pervasive in the mouse gut and that the targets of natural selection are similar in independently E. coli evolving populations. These results illustrate how experimental evolution in natural environments allows us to dissect the mechanisms underlying adaptation and its complex dynamics and further reveal the importance of mobile genetic elements in contributing to the adaptive diversification of bacterial populations in the gut.

A case of adaptation through a mutation in a tandem duplication during experimental evolution in Escherichia coli

Background

DNA duplications constitute important precursors for genome variation. Here we analyzed an unequal duplication harboring a beneficial mutation that may provide alternative evolutionary outcomes.

Results

We characterized this evolutionary event during experimental evolution for only 100 generations of an Escherichia coli strain under glucose limitation within chemostats. By combining Insertion Sequence based Restriction Length Polymorphism experiments, pulsed field gel electrophoresis and two independent genome re-sequencing experiments, we identified an evolved lineage carrying a 180 kb duplication of the 46’ region of the E. coli chromosome. This evolved duplication revealed a heterozygous state, with one copy harboring a 2668 bp deletion that included part of the ogrK gene and both the yegR and yegS genes. By genetically manipulating ancestral and evolved strains, we showed that the single yegS inactivation was sufficient to confer a frequency dependent fitness increase under the chemostat selective conditions in both the ancestor and evolved genetic contexts, implying that the duplication itself was not a direct fitness contributor. Nonetheless, the heterozygous duplicated state was relatively stable in the conditions prevailing during evolution in chemostats, in striking contrast to non selective conditions in which the duplication resolved at high frequency into either its ancestral or deleted copy.

Conclusions

Our results suggest that the duplication state may constitute a second order selection process providing higher evolutionary potential. Moreover, its heterozygous nature may provide differential evolutionary opportunities in alternating environments. Our results also highlighted how careful analyses of whole genome data are needed to identify such complex rearrangements.

Mutation accumulation and fitness in mutator subpopulations of Escherichia coli

Bacterial populations in clinical and laboratory settings contain a significant proportion of mutants with elevated mutation rates (mutators). Mutators have a particular advantage when multiple beneficial mutations are needed for fitness, as in antibiotic resistance. Nevertheless, high mutation rates potentially lead to increasing numbers of deleterious mutations and subsequently to the decreased fitness of mutators. To test how fitness changed with mutation accumulation, genome sequencing and fitness assays of nine Escherichia coli mutY mutators were undertaken in an evolving chemostat population at three time points. Unexpectedly, the fitness in members of the mutator subpopulation became constant despite a growing number of mutations over time. To test if the accumulated mutations affected fitness, we replaced each of the known beneficial mutations with wild-type alleles in a mutator isolate. We found that the other 25 accumulated mutations were not deleterious. Our results suggest that isolates with deleterious mutations are eliminated by competition in a continuous culture, leaving mutators with mostly neutral mutations. Interestingly, the mutator-non-mutator balance in the population reversed after the fitness plateau of mutators was reached, suggesting that the mutator-non-mutator ratio in populations has more to do with competition between members of the population than the accumulation of deleterious mutations.

Analysis of the Lactobacillus casei supragenome and its influence in species evolution and lifestyle adaptation

Background

The broad ecological distribution of L. casei makes it an insightful subject for research on genome evolution and lifestyle adaptation. To explore evolutionary mechanisms that determine genomic diversity of L. casei, we performed comparative analysis of 17 L. casei genomes representing strains collected from dairy, plant, and human sources.

Results

Differences in L. casei genome inventory revealed an open pan-genome comprised of 1,715 core and 4,220 accessory genes. Extrapolation of pan-genome data indicates L. casei has a supragenome approximately 3.2 times larger than the average genome of individual strains. Evidence suggests horizontal gene transfer from other bacterial species, particularly lactobacilli, has been important in adaptation of L. casei to new habitats and lifestyles, but evolution of dairy niche specialists also appears to involve gene decay.

Conclusions

Genome diversity in L. casei has evolved through gene acquisition and decay. Acquisition of foreign genomic islands likely confers a fitness benefit in specific habitats, notably plant-associated niches. Loss of unnecessary ancestral traits in strains collected from bacterial-ripened cheeses supports the hypothesis that gene decay contributes to enhanced fitness in that niche. This study gives the first evidence for a L. casei supragenome and provides valuable insights into mechanisms for genome evolution and lifestyle adaptation of this ecologically flexible and industrially important lactic acid bacterium. Additionally, our data confirm the Distributed Genome Hypothesis extends to non-pathogenic, ecologically flexible species like L. casei.

Culture history and population heterogeneity as determinants of bacterial adaptation: the adaptomics of a single environmental transition

Diversity in adaptive responses is common within species and populations, especially when the heterogeneity of the frequently large populations found in environments is considered. By focusing on events in a single clonal population undergoing a single transition, we discuss how environmental cues and changes in growth rate initiate a multiplicity of adaptive pathways. Adaptation is a comprehensive process, and stochastic, regulatory, epigenetic, and mutational changes can contribute to fitness and overlap in timing and frequency. We identify culture history as a major determinant of both regulatory adaptations and microevolutionary change. Population history before a transition determines heterogeneities due to errors in translation, stochastic differences in regulation, the presence of aged, damaged, cheating, or dormant cells, and variations in intracellular metabolite or regulator concentrations. It matters whether bacteria come from dense, slow-growing, stressed, or structured states. Genotypic adaptations are history dependent due to variations in mutation supply, contingency gene changes, phase variation, lateral gene transfer, and genome amplifications. Phenotypic adaptations underpin genotypic changes in situations such as stress-induced mutagenesis or prophage induction or in biofilms to give a continuum of adaptive possibilities. Evolutionary selection additionally provides diverse adaptive outcomes in a single transition and generally does not result in single fitter types. The totality of heterogeneities in an adapting population increases the chance that at least some individuals meet immediate or future challenges. However, heterogeneity complicates the adaptomics of single transitions, and we propose that subpopulations will need to be integrated into future population biology and systems biology predictions of bacterial behavior.

The multiplicity of divergence mechanisms in a single evolving population

Background

Evolutionary divergence is common within bacterial species and populations, even during a single bacterial infection. We use large-scale genomic and phenotypic analysis to identify the extent of diversification in controlled experimental populations and apply these data to differentiate between several potential mechanisms of evolutionary divergence.

Results

We defined testable differences between five proposed mechanisms and used experimental evolution studies to follow eight glucose-limited Escherichia coli chemostat populations at two growth rates. Simple phenotypic tests identified 11 phenotype combinations evolving under glucose limitation. Each evolved population exhibited 3 to 5 different combinations of the 11 phenotypic clusters. Genome sequencing of a representative of each phenotypic cluster from each population identified 193 mutations in 48 isolates. Only two of the 48 strains had evolved identically. Convergent paths to the same phenotype occurred, but two pleiotropic mutations were unique to slow-growing bacteria, permitting them greater phenotypic variance. Indeed, greater diversity arose in slower-growing, more stressed cultures. Mutation accumulation, hypermutator presence and fitness mechanisms varied between and within populations, with the evolved fitness considerably more uniform with fast growth cultures. Negative frequency-dependent fitness was shown by a subset of isolates.

Conclusions

Evolutionary diversity is unlikely to be explained by any one of the available mechanisms. For a large population as used in this study, our results suggest that multiple mechanisms contribute to the mix of phenotypes and evolved fitness types in a diversifying population. Another major conclusion is that the capacity of a population to diversify is a function of growth rate.

Deletion and acquisition of genomic content during early stage adaptation of Pseudomonas aeruginosa to a human host environment

Adaptation of bacterial pathogens to a permanently host-associated lifestyle by means of deletion or acquisition of genetic material is usually examined through comparison of present-day isolates to a distant theoretical ancestor. This limits the resolution of the adaptation process. We conducted a retrospective study of the dissemination of the P.aeruginosa DK2 clone type among patients suffering from cystic fibrosis, sequencing the genomes of 45 isolates collected from 16 individuals over 35 years. Analysis of the genomes provides a high-resolution examination of the dynamics and mechanisms of the change in genetic content during the early stage of host adaptation by this P.aeruginosa strain as it adapts to the cystic fibrosis (CF) lung of several patients. Considerable genome reduction is detected predominantly through the deletion of large genomic regions, and up to 8% of the genome is deleted in one isolate. Compared with in vitro estimates the resulting average deletion rates are 12- to 36-fold higher. Deletions occur through both illegitimate and homologous recombination, but they are not IS element mediated as previously reported for early stage host adaptation. Uptake of novel DNA sequences during infection is limited as only one prophage region was putatively inserted in one isolate, demonstrating that early host adaptation is characterized by the reduction of genomic repertoire rather than acquisition of novel functions. Finally, we also describe the complete genome of this highly adapted pathogenic strain of P.aeruginosa to strengthen the genetic basis, which serves to help our understanding of microbial evolution in a natural environment.