#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

The SNAP hypothesis: Chromosomal rearrangements could emerge from positive Selection during Niche Adaptation


Authors: Gerrit Brandis aff001;  Diarmaid Hughes aff001
Authors place of work: Department of Medical Biochemistry and Microbiology, Biomedical Center, Uppsala University, Uppsala, Sweden aff001
Published in the journal: The SNAP hypothesis: Chromosomal rearrangements could emerge from positive Selection during Niche Adaptation. PLoS Genet 16(3): e32767. doi:10.1371/journal.pgen.1008615
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1008615

Summary

The relative linear order of most genes on bacterial chromosomes is not conserved over evolutionary timescales. One explanation is that selection is weak, allowing recombination to randomize gene order by genetic drift. However, most chromosomal rearrangements are deleterious to fitness. In contrast, we propose the hypothesis that rearrangements in gene order are more likely the result of selection during niche adaptation (SNAP). Partial chromosomal duplications occur very frequently by recombination between direct repeat sequences. Duplicated regions may contain tens to hundreds of genes and segregate quickly unless maintained by selection. Bacteria exposed to non-lethal selections (for example, a requirement to grow on a poor nutrient) can adapt by maintaining a duplication that includes a gene that improves relative fitness. Further improvements in fitness result from the loss or inactivation of non-selected genes within each copy of the duplication. When genes that are essential in single copy are lost from different copies of the duplication, segregation is prevented even if the original selection is lifted. Functional gene loss continues until a new genetic equilibrium is reached. The outcome is a rearranged gene order. Mathematical modelling shows that this process of positive selection to adapt to a new niche can rapidly drive rearrangements in gene order to fixation. Signature features (duplication formation and divergence) of the SNAP model were identified in natural isolates from multiple species showing that the initial two steps in the SNAP process can occur with a remarkably high frequency. Further bioinformatic and experimental analyses are required to test if and to which extend the SNAP process acts on bacterial genomes.

Keywords:

Natural selection – Evolutionary genetics – Ecological niches – Bacterial genomics – Operons – Genetic drift – Frameshift mutation – Chromosomal duplications

Introduction

Genetic information is encoded in nucleic acid chromosomes organized as linear sequences of genes. Comparative genomic analyses support the hypothesis that life on earth has evolved from a universal common ancestor [16]. This genetic diversity of life reflects the interplay between selection for organisms to occupy and thrive in different environmental niches, and the operation of mechanisms that can change the existing nucleic acid sequence in a chromosome. The mechanisms of genetic change are errors in the accuracy of chromosome replication, and the recombination of sequences within and between chromosomes. The former mechanism can lead to sequence divergence between homologous genes in separate species, whereas the latter mechanism can create novel genes by fusion or splitting of existing genes, and can also move genes from one chromosomal location to another. Because organisms must maintain a high level of relative fitness to compete for resources to support survival, growth and replication, changes in individual gene sequences are often subject to selection to maintain or adapt their functionality in particular environments.

The relationship between selection, conservation of gene order on chromosomes, and relative bacterial fitness in different environments is less obvious. The requirement to integrate gene expression with chromosome replication is one force that shapes the linear organization of bacterial chromosomes. Bacterial genes are most often co-oriented with the direction of replication, and most of the highly conserved and highly expressed genes are located in the half of the chromosome closest to the origin of replication [7]. This replication-related selection can minimise transcription-translation collisions and takes advantage of gene dosage effects to increase expression of some genes but it is not clear that it explicitly selects for maintenance of an ancient linear gene order. A remarkable example of conservation of an ancient co-linear organization of gene is found for a large set of genes involved in transcription and translation [8]. This conserved cluster of operons includes: S10 (11 ribosomal proteins), spc (11 ribosomal proteins and SecY), alpha (4 ribosomal proteins and RpoA), rrnB (3 ribosomal rRNA and 2 tRNA genes), tufB (4 tRNA genes, EF-TuB), secE (SecE, NusG), rpoBC (4 ribosomal proteins, RpoB and RpoC) and str (2 ribosomal proteins, EF-G and EF-TuA). This gene/operon cluster was present in the last common ancestor of the bacteria and archaea [911]. Although in many species some of these operons have become separated by gene insertions, the ancient organization is conserved in many of the Enterobacteriaceae [1215]. The underlying selective mechanism has recently been linked to these operons being concatenated [16]. By experimentally manipulating one of the contiguous operon pairs, tufB-secE in Salmonella, it was shown that an inter-operon terminator-promoter overlap has a significant role in regulating gene expression and its interruption significantly reduces bacterial fitness. The other operons of the ancient cluster that remain contiguous in Salmonella (tufB-secE-rpoBC and S10-spc-alpha) are also each connected by an inter-operon terminator-promoter overlap. Accordingly, it was proposed that the concatenation of operons is an ancient feature of some operons that restricts the potential to rearrange particular regions of bacterial chromosomes and selects for the maintenance of a co-linear operon organization over billions of years [16].

The organization of many bacterial genes into multigene transcriptional units, operons, also suggests mechanisms that could act to conserve linear gene order [1719]. Within operons gene order might be maintained by selection for co-regulation, or for horizontal transfer of a fully functional unit. However, even the conservation of operon organization is generally low over evolutionary time spans for distantly related species [20, 21], although there are exceptions, for example, E. coli and S. enterica, where despite greater than 100 Myr of separation, co-linear gene order within operons, and throughout the chromosome is remarkably conserved [15, 22].

In spite of the examples of conservation above, the linear organization of homologous genes on bacterial chromosomes of different species is highly variable and for most homologous genes there is no long-range co-linearity in gene order [23, 24]. The standard interpretation for the low level of conservation is that selection to maintain linear gene order is weak and this allows changes in gene order to occur by genetic drift. In contrast, an in silico study of contiguous gene pairs across 126 bacterial genomes of different species found that the maintenance of contiguity was actually higher than predicted by experimental parameters, even for gene pairs not in operons, suggesting that many gene order rearrangements are deleterious and that purifying selection is operating [25]. This paradox could be resolved if gene order rearrangements during speciation did not arise primarily by genetic drift but were instead selected. We propose a radical alternative to the drift hypothesis: Selection for Niche Adaptation. The SNAP hypothesis, proposes that changes in relative gene order on bacterial chromosomes are driven by selection. During evolution the organisms that succeed are those that can best adapt to the available environmental niches (survival of the fittest). Such niches are not constant but can arise or change over time as a result of changes in environmental conditions, and because of changes wrought by the interactions of different organisms with both the organic and the physical environments. Our hypothesis is that rearrangements in chromosomal gene order can be selected indirectly as a result of selection acting on organisms (in particular microorganisms) to adapt to changing or novel environmental niches. On an evolutionary timescale the chromosomes of organisms adapting to a new niche would very rapidly ‘snap’ into a new gene order organization. The SNAP hypothesis is explained in words and figures in the text below, and modelled mathematically using reasonable experimentally-derived parameters.

Results

Genetic drift hypothesis

In the standard model, gene order on chromosomes is assumed to be under very weak selection and therefore subject to evolution by genetic drift associated with recombination. Several different types of recombinational event could be involved in rearranging the order of genes on a chromosome: inversion, transposition, deletion, and the acquisition of homologous genes by horizontal gene transfer (Fig 1A). In principle, the successive occurrence of one or more of these types of recombination event could ultimately lead to a significant rearrangement in the linear order of genes on a chromosome. However, in practice the relative fitness of intermediates, and the rates associated with each step in the process, will impose severe limitations on the drift hypothesis as a primary explanation for gene order rearrangements. For an environmentally well-adapted organism there will, in most cases, be no selective benefit associated with inverting, deleting, or transposing a chromosomal segment. Similarly, acquiring additional copies of existing genes by HGT and their insertion at a novel location is unlikely to increase fitness. Deletion or impairment of any essential gene will be lethal or will severely reduce fitness. For most non-lethal chromosomal rearrangements the expectation is that at best they will be neutral but are more likely to have a negative effect on relative fitness [25]. It is unlikely that chromosomal rearrangements, even when they are neutral with respect to fitness, will increase in frequency and reach fixation in a population. A second limitation on the drift hypothesis is the low frequency with which individual non-lethal recombination events, such as inversions, occur in bacterial populations [2629]. Significant gene order rearrangements between species would require a succession of non-lethal recombination events, each occurring with a low probability, and each reaching fixation in a population, to generate a significant shuffling of gene order as observed when comparing different species [2325]. In summary, while the recombinational mechanisms illustrated in Fig 1A could promote genome fluidity over successive cycles, if each event occurs at a low frequency, and without a positive selection, fixation would depend strongly on founder effects (small population bottlenecks). We do not rule out genetic drift as a contributing factor in gene order rearrangements but we think that our alternative hypothesis, SNAP, has some significant advantages in terms of the probability of occurring and being selected to fixation.

Fig. 1. Comparison of standard hypothesis of genomic rearrangements and the SNAP model.
Comparison of standard hypothesis of genomic rearrangements and the SNAP model.
(A) Changes in gene order caused by inversions, transpositions, deletions and re-acquisition. Genes in their original location are shown in dark blue, novel locations are indicated in light blue and genes acquired by horizontal gene transfer are shown in purple. (B) Selection under niche adaptation (SNAP). The gene under selection for duplication is shown in green, genes inactivated are marked with a red X.

The SNAP hypothesis

SNAP, Selection during Niche Adaptation is based on a sequential series of high frequency events and is driven by selection to fixation (Fig 1B). The SNAP hypothesis involves four sequential stages: Duplication, Selection, Inactivation, and Fixation.

(i) Duplication. Duplication of segments of a bacterial chromosome is a very frequent event, occurring spontaneously at rates of >10−2 to 10−5 by recombination between repetitive sequences [30]. The regions duplicated can vary in size from tens of kilobases up to more than a megabase [3032]. Duplications are intrinsically unstable and segregate unless maintained by selection [33].

(ii) Selection. Bacteria frequently live in sub-optimal environments, for example habitats that are nutrient-poor or mildly toxic. Under such conditions, duplications will be selectively maintained if they confer a fitness advantage, for example, if increased dosage of a nutrient transporter gene improves relative fitness [32]. Exposure to antibiotics is also known to select duplications, for example when the bacteria carry a gene encoding a sub-optimal antibiotic-degrading enzyme [34, 35]. In such cases the increase in gene dosage associated with a duplication or amplification provides a strong selective benefit in the particular environmental niche. In addition to having a gene dosage effect, a duplication could also confer a selective advantage by placing a gene under the control of an alternative potent promoter thus increasing its expression or altering its regulation [36]. Adaptive duplications could also be selected for fast growth in nutrient-rich environments. An example could be the occurrence of multiple rrn operons in many microbial species that may be a selected genetic mechanism contributing to fast growth [3741]. Also, the frequently observed duplication of the tuf gene, encoding elongation factor EF-Tu, may have been selected in different bacterial species because this duplication helps support faster growth rates than are supported by a single gene copy [4244].

(iii) Inactivation. A duplication is a double-edged sword. The regional duplication will be maintained by selection on the relevant gene(s) but the other genes in the duplicated region will not be under positive selection. Accordingly, most duplicated genes, even those that are essential as single copy genes, can accumulate mutations, either because they are not essential as duplicates, or because their duplication reduces fitness (resource wastage, interference with normal physiology) and there is a positive selection to remove their activity [45]. This process inevitably leads to the accumulation of inactivating mutations in the genes of the duplicated region that are not under positive selection. Gene-inactivating mutations (for example, frameshift, nonsense, deletion) occur with spontaneous rates of 10−5 to 10−6 per gene [46, 47]. Recombination between repeat sequences that lie within the duplicated region (IS elements for example, or other repeat sequences) could lead to a loss of parts of a duplication, including a copy of an essential gene, at much higher rates. We make the reasonable assumption that gene inactivation mutation will occur randomly with respect to each copy of the duplication.

(iv) Fixation. Inactivation of a different essential gene (or a gene required for high fitness) in each copy of the duplicated region will prevent segregation of the duplication. At this stage the duplication is fixed and the net outcome is a chromosome in which the remaining active genes have a rearranged order relative to the ancestral order (see Fig 1B). The remaining duplicated genes can continue to accumulate mutations (including deletions) in each copy of the duplicated region contributing further to rearrangements of the original gene order. In E. coli there are over 350 chromosomal genes that are essential for growth under rich medium conditions [48] but in general bacteria will have many other genes where inactivation would significantly reduce fitness, or be incompatible with growth under a variety of specific conditions [4951]. Accordingly, a duplicated region of 100 kb is likely to contain several essential genes providing mutational targets where inactivation will result in fixation of the duplication and a rearranged gene order on the chromosome.

The SNAP hypothesis does not rule out a role for genetic drift in causing gene order rearrangements. It is an alternative mechanism that has very significant advantages compared to genetic drift: it is associated with natural selection (bacteria adapting to a new environment), it is initiated at a very high frequency (spontaneous duplications), it is irreversible (once essential genes have been inactivated in each arm of the duplication), and it is driven to fixation by positive selection. Accordingly, we propose that positive selection might play a major role in driving change in the relative order of most genes on a bacterial chromosome.

Mathematical modelling of SNAP

A minimal mathematical model of SNAP is presented in Fig 2. The spontaneous rates of duplication and mutational gene inactivation used in modelling are taken from published literature [30, 46, 47]. The only variable parameter in the model is the range of potential effects on relative fitness of duplications and mutations within duplicated regions. The model makes the following assumptions: (i) regional duplications occur and can be maintained by selection for a phenotype that is satisfied by duplication of one or more genes encoded within the duplicated region; (ii) the duplicated regions contain at least two essential genes; (iii) gene inactivating mutations occur with normal rates and can inactivate different essential genes in each copy of the duplicated region; (iv) once two different essential genes have been inactivated in different copies of the duplicated region the duplication can no longer segregate to a single copy while maintaining the original gene order.

Fig. 2. Outline of the minimal population dynamics model.
Outline of the minimal population dynamics model.
The model allows the appearance of four types of cells: wild-type cells (Wt), cells with a duplication (Dup) of a region that includes two essential genes (green), and cells with the duplication and a single (Δ1) or double (Δ2) inactivation of essential genes. All possible directions and rates of evolution are displayed and the inactivation of two essential genes is assumed to stabilize the duplication.

In this model the wild-type spontaneously generates duplications that are stabilized by a selection for a phenotype (step 1). An essential gene within one copy of the duplicated region is mutationally inactivated (step 2). At this stage there are alternative paths. If the duplication is maintained there is the opportunity for an essential gene within the second copy of the duplicated region to be mutationally inactivated (step 3). Step 3 stabilizes the duplication with a novel linear gene order. Alternatively, if the duplication segregates (for example, because selection is relieved) the original gene order will be maintained. The minimal model is illustrated here with rates for each step that are conservative estimates based on experimentally determined values [30, 46, 47].

Using this minimal mathematical model, we have measured how changing the values assigned to the fitness parameters would influence the probability of fixing a rearranged gene order (Fig 3). In the absence of any selection or fitness costs, duplication and single gene inactivation occurs and reaches a steady state but does not go to fixation (Fig 3, panel A). Once selection and fitness costs are introduced (a novel environment where the duplication has a fitness advantage over the wild-type) the population carrying duplications increases dramatically and sub-populations carrying the single and double gene inactivation mutations increase in frequency (Fig 3, panels B and C). Adding the assumption that carrying duplicate genes confers a fitness cost leads to the rapid increase and subsequent fixation of the mutant population with the novel gene order (double gene inactivation) (Fig 3, panels D-F). This minimal model suggests that a novel gene order can be generated within a small number of generations if the initial duplication has a selective benefit over the wild-type and the inactivation of duplicate genes from either of the copies further improves fitness.

Fig. 3. Modelled population dynamics under varying selective conditions.
Modelled population dynamics under varying selective conditions.
Number of wild-type cells (Wt, black), cells with the duplication (Dup, green) and cells with the duplication that carry a single (Δ1, blue) or double (Δ2, red) essential gene inactivation mutation are shown as a function of time. (A-C) Strains carrying the duplication have, relative to the wild-type (A) equal fitness, (B) 5% fitness advantage, or (C) 25% fitness advantage. (D-F) Illustrate panel (C) with the added assumption that deleting unnecessary duplicate genes (Δ2) confers a fitness advantage of (D) 2%, (E) 5%, or (F) 10%. All models were run as serial transfers with a starting population of 106 wild-type cells, a total population size of 1010 cells and 108 cells transferred per bottleneck. The appearance and reversion of mutant populations was determined by a Monte Carlo procedure based on the frequencies displayed in Fig 2. The fitness parameters for the populations are shown in each panel. All graphs display the average of 100 independent runs. Models were run with Berkeley Madonna (version 9.1.14).

There are several additional features, that for simplicity, have been omitted from this minimal model but which may play either a restrictive or a positive role in this evolutionary process in different species, or under different selective conditions.

  • A feature of the model that potentially restricts its influence on genome rearrangements is the requirement that at least two essential genes be contained within the duplicated region. Essential genes are not expected to be evenly distributed throughout the genome, in which case for some duplications there might never be a transition from step 2 to step 3. This restriction will mostly affect smaller duplications in regions of the chromosome that are poor in essential genes but is less likely to affect large duplications. A counter argument is that under the actual conditions that are selective for maintenance of a duplication (e.g., growth in a challenging niche) many additional genes, even if not essential under all conditions, may be under strong selection to maintain fitness [44, 49, 50].

  • A feature of the model that potentially promotes gene order rearrangements is that many duplications will result in unbalanced chromosome replichores. These mutants will be under selection not only to maintain the duplication but also to rebalance their replichores so as to reduce associated fitness costs [5257]. An improvement in replichore balance could be achieved by a deletion or an inversion. Chromosomes that have undergone a process of duplication followed by inversion will be locked into a structure where the duplication can no longer easily be segregated. This sequence of events can help to promote genome rearrangements by effectively stabilizing a duplication even if the original selection is relieved.

  • The evolutionary process does not stop after an essential gene has been inactivated in each duplicated region. The fitness costs associated with having tens to hundreds of genes duplicated will act as a driving force for the continued selection and fixation of mutants that delete or otherwise inactivate all non-required extra copies of duplicated genes where such duplications have a negative impact on fitness.

  • Another feature that could promote rapid genetic change is the high prevalence of bacteria that are mutator clones with high mutation rates. Mutator bacteria are estimated to be up to 1% of natural isolates [5860], and even higher among some clinical isolates [61]. Mutator clones, including those caused by inactivation of the mismatch repair system, have not only a significantly increased rate of point mutation [62] but also a significantly higher rate of recombination that can cause chromosomal rearrangements including duplications, deletions and inversions [27, 63, 64]. Recombinational gene inactivation could also be caused by the movement of IS elements and transposons, the frequency of which will vary between species and potentially be influenced by the environment. With regard to mobile genetic elements (MGE) we note that care must be taken in estimating the number of duplications in genome sequences, to distinguish between those involving non-mobile sequences (the main focus of the SNAP hypothesis) and duplications arising from the movement of MGEs resulting in increased copy number.

Gene inactivation by point mutations occurring at a normal mutation rate (as modelled in Fig 3) leads to a very conservative estimate of gene inactivation rates, and if instead, deletion and insertional inactivation events dominate, and mutators play a significant role, then the rates of gene inactivation within a duplicated region of the chromosome could be much higher than in our simple model.

Identification of duplications in natural isolates

Available genome sequences from clinical and environmental isolates of Acinetobacter baumannii, Escherichia coli, Mycobacterium tuberculosis, and Pseudomonas aeruginosa were analysed to identify signature features (duplication formation and divergence) of the SNAP model. One hundred genome sequences for each species were downloaded from the Sequence Read Archive (SRA), assembled to a respective standard reference sequence, and duplications were identified based on increased sequence coverage. Duplications were present in 2–4% of the isolates of each species and ranged in length from 8 to 355 kb (Fig 4A). Further analysis of the duplicated sequences showed that two of the fourteen isolates (14%) contained diverging duplications, identified as having a mutation present in ~50% of the reads: A M. tuberculosis isolate had a frameshift mutation in one copy of MRA_RS09940, a glutamine synthetase gene (Fig 4B and 4C) and an E. coli isolate had a R276C mutation in one copy of the dacD gene encoding D-alanyl-D-alanine carboxypeptidase.

Fig. 4. Duplications identified in natural isolates of A. baumannii, E. coli, M. tuberculosis, and P. aeruginosa.
Duplications identified in natural isolates of <i>A</i>. <i>baumannii</i>, <i>E</i>. <i>coli</i>, <i>M</i>. <i>tuberculosis</i>, and <i>P</i>. <i>aeruginosa</i>.
(A) One hundred whole genome sequences per species were downloaded from the SRA and analysed for regions with increased coverage. Duplicated regions are indicated with green bars and represent unique segments of the chromosomes. (B) Read coverage analysis of a chromosomal section within the M. tuberculosis isolate with a 21 kb duplication. The blue shades (top to bottom) represent the maximum, average and minimum read coverage on a sliding 1 kb window. Genes within the chromosomal segment are indicated below. The duplicated region contains 21 genes and the frameshift mutation that is present in one copy of the glutamine synthetase gene is indicated with a dotted red line. (C) Sequence analysis of the frameshift insertion within the glutamine synthetase gene (~25% of reads shown). The consensus sequence is shown as sequence logo on the top with the reads below. Residues in the reads that match the reference are shown as dots. The insertion of a thymine is indicated in red. The site of the insertion has a 155-fold coverage and the frameshift present in 49% of reads.

The number of identified duplications in this dataset is most likely an underestimate. Culturing isolates under laboratory conditions to obtain pure cultures will remove the conditions that selected the duplication and will lead to segregation unless the duplication is stabilized. The fact that multiple isolates with duplications were identified for every species shows that duplications of chromosomal regions are very common among natural isolates. These duplications were stable enough to be present after laboratory culture conditions and to acquire mutations in one of the duplicated copies. The M. tuberculosis isolate that had an inactivating frameshift mutation in one copy of a glutamine synthetase gene (Fig 4C) represents in principle the Δ1 mutant class predicted in the SNAP model (Fig 2).

Discussion

Understanding drivers and mechanisms of genetic change is fundamental to understanding the diversity of life on earth. This diversity of lifeforms has evolved from a common ancestor by mutation and recombination of existing genetic material. Most research in this area has focused on the causes, and selection, of changes in gene sequences, and there has been much less research into the causes, and selection, of changes at the level of the chromosome [7]. Current theory interprets the widespread diversity in chromosomal gene order as evidence of very weak selection, with rearrangements occurring by genetic drift. Accordingly, rearrangements in gene order that are not counter-selected can accumulate by successive recombinational events (inversions, transpositions, deletions, and re-acquisitions by HGT) leading ultimately to a shuffled set of genes [25]. However, experimental evidence shows that most individual chromosomal rearrangements reduce fitness, creating a barrier to their fixation [28, 29]. The major advantages of the SNAP hypothesis over the genetic drift hypothesis are: (i) it is associated with an important lifestyle event (entry into a new ecological niche); (ii) it is initiated by a high-frequency event (partial chromosome duplication); (iii) it is driven by positive selection (adaptation to the new niche by increased gene dosage); (iv) selection to reduce the dosage of non-selected genes drives the loss of function or deletion of many duplicated genes; (v) the loss of essential genes in each copy of the duplicated region traps the rearrangement; (vi) a rearranged gene order becomes fixed in the niche-adapted bacterial variant. An additional consequence is that bacteria with a novel gene order will be genetically more isolated, contributing to the process of species separation in bacteria.

Most bacterial genes are organized into multigene transcriptional units, operons, that can be physiologically advantageous in terms of transcriptional co-regulation of genes with intersecting functionalities [1719]. The organization of genes into operons is likely to act as a selective force resisting disruptive rearrangements in linear gene order within the operon if that reduces relative fitness. In this regard, finely regulated operons may be under stronger positive selection and able to resist disruptive rearrangements more than poorly regulated operons. However, even for the tryptophan operon, a classic whole-pathway operon with an ancient history (present in the common ancestor of Bacteria and Archaea), phylogenetic analysis has revealed many differences in gene order in different bacterial lineages [65]. Operons can also be advantageous for their member genes on an evolutionary timescale, by increasing the likelihood that the genes contained within the operon can benefit from horizontal gene transfer events by being transferred as part of a fully functional unit [20, 21]. Re-ordering linear gene order is however, not just a potential disrupter of operons. Rearrangements in linear gene order can act to create novel transcriptional units with potential selective value if they increase fitness of the organism [66, 67]. Accordingly, the pathway to fixation of a new gene order during the process of SNAP could involve a series of different selection processes: selection to maintain the initially selected gene dosage benefit, selection to reduce the negative effect of costly duplications, and selection to maintain fortuitously created novel regulatory units arising during the fixation process.

The SNAP hypothesis as outlined here is a dynamic process that begins with high-frequency spontaneous duplications of chromosome segments [30] that are maintained by selection for increased gene dosage [3032, 68], and ultimately, through a process of mutation and recombination, driven by selection for high fitness, results in the fixation a new linear gene order (Fig 1B). The high frequency of chromosome segment duplications predicts that occasionally the duplication should be retained, either by selection for gene dosage or as a result of mutational fixation. Genome analyses provide evidence for some bacterial genes arising by duplication [6974]. One interesting example is that duplicated segments have been found in the genomes of Mycobacterial species, ranging in size from 30 to 350 kb [7577] suggesting that they are maintained, or very frequently generated, by selection. The frequent presence of multiple copies of ribosomal RNA operons in bacterial genomes is a classic example of duplicated genes that are stably maintained on evolutionary timescales. It is assumed that these operons have a common evolutionary origin and that the presence of multiple copies in many [37] but not all [78, 79] bacterial species is most probably explained by duplication of chromosomal regions. The selection for different copy numbers correlates closely with growth rate [37] but there is evidence that selection for adaptation to different ecological niches and for the ability to respond efficiently to the availability of resources also plays a significant role [38].

To search for genomic evidence relevant to the SNAP hypothesis we examined recent genome sequence data deposited at the Sequence Read Archive. We chose, without any pre-screening, one hundred genome sequences from each of four clinically important bacterial species: E. coli, P. aeruginosa, A. baumannii, and M. tuberculosis (SI, Table). We searched the raw sequence reads for evidence of partial chromosomal duplications (step 1 in the model), and mutations within one copy of a duplicated region (step 2 in the model). We found duplicated regions in the genome sequences of all four species at frequencies of 2 to 4%, and we also observed mutations at 50% frequency in 2 of the 14 duplicated regions (Fig 4). These mutations included one frameshift mutation in a duplicated region of M. tuberculosis that is expected to inactivate the gene (glutamine synthetase) and this represents a good example of the second step in the model (Fig 2). Given that these clinical samples do not represent bacteria encountering a novel environment, and that the genomic DNA was prepared for sequencing without special selection to maintain unstable duplications, these data show that the initial two steps in the SNAP process can occur with a remarkably high frequency.

The computational model and the genome-level analysis of natural isolates sequences indicate that the SNAP process can act on bacterial genomes. Nevertheless, so far there is no direct empirical evidence that genome rearrangements in modern bacterial species have been caused by SNAP. A complicating factor is that once the SNAP process is completed there is no genome feature left that is unique to the model. A possible bioinformatic approach to test the hypothesis would be a high-throughput analysis of modern bacterial chromosomes to search for intermediate steps of the SNAP process. For example, a larger than expected number of duplicate genes and/or pseudogenes with matching active copies could be the remains of old duplications. Alternatively, a long-term adaptation experiment of a bacterial clone to a novel environment (e.g. growth on a poor carbon source) could be analysed to experimentally identify and validate each of the proposed steps of the SNAP hypothesis.

In summary, the SNAP hypothesis is based on a sequential series of high-frequency events (ecological and genetic), that can selectively drive a process leading with a high probability to rearrangements in chromosomal gene order, and possibly also contributing to creating species barriers between bacteria.

Methods

Mathematical model

The mathematical model was designed to model 1000 generations of growth of a wild-type population (Wt). The model allows the appearance of cells with a small duplication (Dup) that includes two essential genes, and cells with the duplication and a single (Δ1) or double (Δ2) inactivation of essential genes. Rates of duplication formation and mutational gene inactivation were estimated based on previous experimental data [30, 46, 47]. All possible directions and rates of evolution are displayed in Fig 2 and the inactivation of two essential genes is assumed to stabilize the duplication. Fitness effects of duplications and gene inactivations were the variable parameter of the model and are displayed in Fig 3.

The bacterial growth rate is a monotonically increasing function of the concentration of a limiting resource, R (mg L-1) [80]

where Vi is the relative fitness of the ith strain of bacteria and k is the concentration of the resource at which Vi is at half its maximum value. With these definitions the change in densities of bacterial populations and the concentration of resources are given by the following two coupled differential equations:
where ni is the density of strain i (cfu mL-1) and e is the conversion efficiency parameter (μg cell-1). The standard parameters Rt = 0 = 100 mg L-1, k = 1 mg L-1, and e = 10−9 μg cell-1 result in a growth cycle that leads to a final density of approximately 1010 cfu mL-1. After every cycle the culture is 100-fold diluted (108 cells per bottleneck) into fresh media and grown to full density. Serial passaging was repeated until a total growth of 1000 generations. A Monte Carlo procedure was used to determine the appearance of Wt, Dup, Δ1 and Δ2 cells. The probability pi>j(t) that a cell j is generated from a cell i at time point t is
where gi is the number of generations of growth of the strain i at time point t, and μi>j is the mutation/recombination rate to generate cell j from cell i. A random number x (0 < x < 1) is generated. A single cell of strain j generated at time point t if x < pi>j(t). The simulation was programmed in Berkeley Madonna (Version 9.1.14) and run with varying fitness values. All results are averages of 100 independent simulations.

Analysis of natural isolates

Genome analyses were performed using the CLC Genomics Workbench version 11.0.1 (Qiagen). Whole genome sequencing reads were downloaded from the Sequence Read Archive for one hundred natural isolates per species (S1 Table). The downloaded reads were trimmed and mapped to a respective standard reference sequence (Trim settings: Quality limit: 0.05; Ambiguous limit: 2. Mapping settings: Match score: 1; Mismatch score: 2; Cost of insertions and deletions: Linear gap cost; Insertion cost: 3; Deletion cost: 3; Insertion open cost: 6; Insertion extend cost: 1; Deletion open cost: 6; Deletion extend cost: 1; Length fraction: 0.5; Similarity fraction 0.8; Auto-detect paired distances; Non-specific match handling: Map randomly. Reference sequences: A. baumannii str. ACICU: NC_010611; E. coli K-12 str. MG1655: NC_000913; M. tuberculosis str. H37Ra: NC_009525; P. aeruginosa str. PAO1: NC_002516). Duplications were identified based on visual assessment of the CLC sequence coverage tracks. See Fig 4B for an example of an identified duplication. All isolates containing duplications are highlighted yellow in S1 Table.

Supporting information

S1 Table [xlsx]
SRA metadata tables.


Zdroje

1. Woese CR. Interpreting the universal phylogenetic tree. Proc Natl Acad Sci USA. 2000;97(15):8392–6. doi: 10.1073/pnas.97.15.8392 WOS:000088273900039. 10900003

2. Koonin EV. Comparative genomics, minimal gene-sets and the last universal common ancestor. Nat Rev Microbiol. 2003;1(2):127–36. doi: 10.1038/nrmicro751 WOS:000220402500014. 15035042

3. Koonin EV. Carl Woese's vision of cellular evolution and the domains of life. RNA Biol. 2014;11(3):197–204. doi: 10.4161/rna.27673 WOS:000334999500006. 24572480

4. Forterre P. The universal tree of life: an update. Front Microbiol. 2015;6. doi: 10.3389/fmicb.2015.00717 WOS:000358717800001. 26257711

5. Booth A, Mariscal C, Doolittle WF. The Modern Synthesis in the Light of Microbial Genomics. Annu Rev Microbiol. 2016;70:279–97. doi: 10.1146/annurev-micro-102215-095456 WOS:000383052200016. 27482743

6. Weiss MC, Preiner M, Xavier JC, Zimorski V, Martin WF. The last universal common ancestor between ancient Earth chemistry and the onset of genetics. PLoS Genet. 2018;14(8). doi: 10.1371/journal.pgen.1007518 WOS:000443389100009. 30114187

7. Touchon M, Rocha EPC. Coevolution of the Organization and Structure of Prokaryotic Genomes. CSH Perspect Biol. 2016;8(1). doi: 10.1101/cshperspect.a018168 WOS:000371181000007. 26729648

8. Tamames J. Evolution of gene order conservation in prokaryotes. Genome Biol. 2001;2(6). WOS:000207584100012.

9. Wachtershauser G. Towards a reconstruction of ancestral genomes by gene cluster alignment. Syst Appl Microbiol. 1998;21(4):473–7. WOS:000078011000001.

10. Coenye T, Vandamme P. Organisation of the S10, spc and alpha ribosomal protein gene clusters in prokaryotic genomes. Fems Microbiol Lett. 2005;242(1):117–26. doi: 10.1016/j.femsle.2004.10.050 WOS:000226264100016. 15621428

11. Barloy-Hubler F, Lelaure V, Galibert F. Ribosomal protein gene cluster analysis in eubacterium genomics: homology between Sinorhizobium meliloti strain 1021 and Bacillus subtilis. Nucleic Acids Res. 2001;29(13):2747–56. doi: 10.1093/nar/29.13.2747 PubMed Central PMCID: PMC55768. 11433019

12. Brocks JJ, Schaeffer P. Okenane, a biomarker for purple sulfur bacteria (Chromatiaceae), and other new carotenoid derivatives from the 1640 Ma Barney Creek Formation. Geochim Cosmochim Ac. 2008;72(5):1396–414. doi: 10.1016/j.gca.2007.12.006 WOS:000254198000010.

13. Marin J, Battistuzzi FU, Brown AC, Hedges SB. The Timetree of Prokaryotes: New Insights into Their Evolution and Speciation. Mol Biol Evol. 2017;34(2):437–46. doi: 10.1093/molbev/msw245 WOS:000396511300012. 27965376

14. Brocks JJ, Love GD, Summons RE, Knoll AH, Logan GA, Bowden SA. Biomarker evidence for green and purple sulphur bacteria in a stratified Palaeoproterozoic sea. Nature. 2005;437(7060):866–70. doi: 10.1038/nature04068 16208367.

15. McClelland M, Sanderson KE, Spieth J, Clifton SW, Latreille P, Courtney L, et al. Complete genome sequence of Salmonella enterica serovar Typhimurium LT2. Nature. 2001;413(6858):852–6. doi: 10.1038/35101614 11677609.

16. Brandis G, Cao S, Hughes D. Operon concatenation is an ancient feature that restricts the potential to rearrange bacterial chromosomes. Mol Biol Evol. 2019;36(9):1990–2000. doi: 10.1093/molbev/msz129 31132113.

17. Tamames J, Casari G, Ouzounis C, Valencia A. Conserved clusters of functionally related genes in two bacterial genomes. J Mol Evol. 1997;44(1):66–73. doi: 10.1007/pl00006122 WOS:A1997WD35900007. 9010137

18. Ermolaeva MD, White O, Salzberg SL. Prediction of operons in microbial genomes. Nucleic Acids Res. 2001;29(5):1216–21. doi: 10.1093/nar/29.5.1216 WOS:000167240500024. 11222772

19. Moreno-Hagelsieb G, Trevino V, Perez-Rueda E, Smith TF, Collado-Vides J. Transcription unit conservation in the three domains of life: a perspective from Escherichia coli. Trends Genet. 2001;17(4):175–7. doi: 10.1016/s0168-9525(01)02241-7 WOS:000168718300004. 11275307

20. Lawrence JG, Roth JR. Selfish operons: Horizontal transfer may drive the evolution of gene clusters. Genetics. 1996;143(4):1843–60. WOS:A1996VA24400030. 8844169

21. Itoh T, Takemoto K, Mori H, Gojobori T. Evolutionary instability of operon structures disclosed by sequence comparisons of complete microbial genomes. Mol Biol Evol. 1999;16(3):332–46. doi: 10.1093/oxfordjournals.molbev.a026114 WOS:000079160500003. 10331260

22. Ochman H, Groisman EA. The origin and evolution of species differences in Escherichia coli and Salmonella typhimurium. EXS. 1994;69:479–93. doi: 10.1007/978-3-0348-7527-1_27 7994120.

23. Tatusov RL, Mushegian AR, Bork P, Brown NP, Hayes WS, Borodovsky M, et al. Metabolism and evolution of Haemophilus influenzae deduced from a whole-genome comparison with Escherichia coli. Curr Biol. 1996;6(3):279–91. doi: 10.1016/s0960-9822(02)00478-5 WOS:A1996UC44000022. 8805245

24. Koonin EV, Mushegian AR, Rudd KE. Sequencing and analysis of bacterial genomes. Curr Biol. 1996;6(4):404–16. doi: 10.1016/s0960-9822(02)00508-0 WOS:A1996UH68400022. 8723345

25. Rocha EPC. Inference and analysis of the relative stability of bacterial chromosomes. Mol Biol Evol. 2006;23(3):513–22. doi: 10.1093/molbev/msj052 WOS:000235610300005. 16280545

26. Brandis G, Cao S, Hughes D. Co-evolution with recombination affects the stability of mobile genetic element insertions within gene families of Salmonella. Mol Microbiol. 2018: doi: 10.1111/mmi.13959 29603442.

27. Hughes D. Co-evolution of the tuf genes links gene conversion with the generation of chromosomal inversions. J Mol Biol. 2000;297(2):355–64. doi: 10.1006/jmbi.2000.3587 10715206.

28. Hughes D. Evaluating genome dynamics: the constraints on rearrangements within bacterial genomes. Genome Biol. 2000;1(6). WOS:000207583400002.

29. Hughes D. Impact of homologous recombination on genome organization and stability. In: Charlebois RL, editor. Organization of the prokaryotic genome. Washington DC, USA: ASM Press; 1999. p. 109–28.

30. Anderson P, Roth J. Spontaneous tandem genetic duplications in Salmonella typhimurium arise by unequal recombination between rRNA (rrn) cistrons. Proc Natl Acad Sci U S A. 1981;78(5):3113–7. doi: 10.1073/pnas.78.5.3113 6789329.

31. Straus DS, Hoffmann GR. Selection for a large genetic duplication in Salmonella typhimurium. Genetics. 1975;80(2):227–37. 1093939.

32. Sonti RV, Roth JR. Role of gene duplications in the adaptation of Salmonella typhimurium to growth on limiting carbon sources. Genetics. 1989;123(1):19–28. 2680755.

33. Andersson DI, Hughes D. Gene amplification and adaptive evolution in bacteria. Annu Rev Genet. 2009;43:167–95. doi: 10.1146/annurev-genet-102108-134805 19686082.

34. Sandegren L, Andersson DI. Bacterial gene amplification: implications for the evolution of antibiotic resistance. Nat Rev Microbiol. 2009;7(8):578–88. doi: 10.1038/nrmicro2174 19609259.

35. Sun S, Berg OG, Roth JR, Andersson DI. Contribution of gene amplification to evolution of increased antibiotic resistance in Salmonella typhimurium. Genetics. 2009;182(4):1183–95. doi: 10.1534/genetics.109.103028 19474201; PubMed Central PMCID: PMC2728858.

36. Anderson RP, Roth JR. Tandem genetic duplications in phage and bacteria. Annu Rev Microbiol. 1977;31:473–505. doi: 10.1146/annurev.mi.31.100177.002353 334045.

37. Roller BR, Stoddard SF, Schmidt TM. Exploiting rRNA operon copy number to investigate bacterial reproductive strategies. Nat Microbiol. 2016;1(11):16160. doi: 10.1038/nmicrobiol.2016.160 27617693; PubMed Central PMCID: PMC5061577.

38. Klappenbach JA, Dunbar JM, Schmidt TM. rRNA operon copy number reflects ecological strategies of bacteria. Appl Environ Microbiol. 2000;66(4):1328–33. doi: 10.1128/aem.66.4.1328-1333.2000 10742207; PubMed Central PMCID: PMC91988.

39. Gyorfy Z, Draskovits G, Vernyik V, Blattner FF, Gaal T, Posfai G. Engineered ribosomal RNA operon copy-number variants of E. coli reveal the evolutionary trade-offs shaping rRNA operon number. Nucleic Acids Res. 2015;43(3):1783–94. doi: 10.1093/nar/gkv040 25618851; PubMed Central PMCID: PMC4330394.

40. Valdivia-Anistro JA, Eguiarte-Fruns LE, Delgado-Sapien G, Marquez-Zacarias P, Gasca-Pineda J, Learned J, et al. Variability of rRNA Operon Copy Number and Growth Rate Dynamics of Bacillus Isolated from an Extremely Oligotrophic Aquatic Ecosystem. Front Microbiol. 2015;6:1486. doi: 10.3389/fmicb.2015.01486 26779143; PubMed Central PMCID: PMC4700252.

41. Yano K, Masuda K, Akanuma G, Wada T, Matsumoto T, Shiwa Y, et al. Growth and sporulation defects in Bacillus subtilis mutants with a single rrn operon can be suppressed by amplification of the rrn operon. Microbiol. 2016;162(1):35–45. doi: 10.1099/mic.0.000207 26518335.

42. Kacar B, Garmendia E, Tuncbag N, Andersson DI, Hughes D. Functional Constraints on Replacing an Essential Gene with Its Ancient and Modern Homologs. mBio. 2017;8(4):e01276–17. ARTN e01276-17 doi: 10.1128/mBio.01276-17 WOS:000409384300045. 28851849

43. Garmendia E, Brandis G, Hughes D. Transcriptional Regulation Buffers Gene Dosage Effects on a Highly Expressed Operon in Salmonella. mBio. 2018;9(5). doi: 10.1128/mBio.01446-18 30206172; PubMed Central PMCID: PMC6134099.

44. Tubulekas I, Hughes D. Growth and translation elongation rate are sensitive to the concentration of EF-Tu. Mol Microbiol. 1993a;8(4):761–70. doi: 10.1111/j.1365-2958.1993.tb01619.x 8332067.

45. Adler M, Anjum M, Berg OG, Andersson DI, Sandegren L. High fitness costs and instability of gene duplications reduce rates of evolution of new genes by duplication-divergence mechanisms. Mol Biol Evol. 2014;31(6):1526–35. doi: 10.1093/molbev/msu111 24659815.

46. Andersson DI, Hughes D, Roth JR. The origin of mutants under selection: Interactions of mutation, growth and selection. 2011. In: EcoSal-Escherichia coli and Salmonella: Cellular and Molecular Biology [Internet]. Washington, DC.: ASM Press. Available from: http://www.ecosal.org.

47. Praski Alzrigat L, Huseby DL, Brandis G, Hughes D. Fitness cost constrains the spectrum of marR mutations in ciprofloxacin-resistant Escherichia coli. J Antimicrob Chemother. 2017;72(11):3016–24. doi: 10.1093/jac/dkx270 28962020; PubMed Central PMCID: PMC5890708.

48. Goodall ECA, Robinson A, Johnston IG, Jabbari S, Turner KA, Cunningham AF, et al. The Essential Genome of Escherichia coli K-12. mBio. 2018;9(1). doi: 10.1128/mBio.02096-17 29463657; PubMed Central PMCID: PMC5821084.

49. Chaudhuri RR, Morgan E, Peters SE, Pleasance SJ, Hudson DL, Davies HM, et al. Comprehensive assignment of roles for Salmonella typhimurium genes in intestinal colonization of food-producing animals. PLoS Genet. 2013;9(4):e1003456. doi: 10.1371/journal.pgen.1003456 23637626; PubMed Central PMCID: PMC3630085.

50. Vohra P, Chaudhuri RR, Mayho M, Vrettou C, Chintoan-Uta C, Thomson NR, et al. Retrospective application of transposon-directed insertion-site sequencing to investigate niche-specific virulence of Salmonella Typhimurium in cattle. BMC Genomics. 2019;20(1):20. doi: 10.1186/s12864-018-5319-0 30621582; PubMed Central PMCID: PMC6325888.

51. Lawley TD, Chan K, Thompson LJ, Kim CC, Govoni GR, Monack DM. Genome-wide screen for Salmonella genes required for long-term systemic infection of the mouse. PLoS Pathog. 2006;2(2):e11. doi: 10.1371/journal.ppat.0020011 16518469; PubMed Central PMCID: PMC1383486.

52. Campo N, Dias MJ, Daveran-Mingot ML, Ritzenthaler P, Le Bourgeois P. Chromosomal constraints in Gram-positive bacteria revealed by artificial inversions. Mol Microbiol. 2004;51(2):511–22. doi: 10.1046/j.1365-2958.2003.03847.x 14756790.

53. Liu GR, Liu WQ, Johnston RN, Sanderson KE, Li SX, Liu SL. Genome plasticity and ori-ter rebalancing in Salmonella typhi. Mol Biol Evol. 2006;23(2):365–71. doi: 10.1093/molbev/msj042 16237205.

54. Savic DJ, Nguyen SV, McCullor K, McShan WM. Biological impact of a large-scale genomic inversion that grossly disrupts the relative positions of the origin and terminus loci of the Streptococcus pyogenes chromosome. J Bacteriol. 2019;201(17). doi: 10.1128/JB.00090-19 31235514; PubMed Central PMCID: PMC6689312.

55. Lesterlin C, Pages C, Dubarry N, Dasgupta S, Cornet F. Asymmetry of chromosome replichores renders the DNA translocase activity of FtsK essential for cell division and cell shape maintenance in Escherichia coli. PLoS Genet. 2008;4(12):e1000288. doi: 10.1371/journal.pgen.1000288 19057667; PubMed Central PMCID: PMC2585057.

56. Darling AE, Miklos I, Ragan MA. Dynamics of genome rearrangement in bacterial populations. PLoS Genet. 2008;4(7):e1000128. doi: 10.1371/journal.pgen.1000128 18650965; PubMed Central PMCID: PMC2483231.

57. Esnault E, Valens M, Espeli O, Boccard F. Chromosome structuring limits genome plasticity in Escherichia coli. PLoS Genet. 2007;3(12):e226. doi: 10.1371/journal.pgen.0030226 18085828; PubMed Central PMCID: PMC2134941.

58. Taddei F, Radman M, Maynard-Smith J, Toupance B, Gouyon PH, Godelle B. Role of mutator alleles in adaptive evolution. Nature. 1997;387(6634):700–2. doi: 10.1038/42696 9192893.

59. LeClerc JE, Li B, Payne WL, Cebula TA. High mutation frequencies among Escherichia coli and Salmonella pathogens. Science. 1996;274(5290):1208–11. doi: 10.1126/science.274.5290.1208 8895473.

60. Gross MD, Siegel EC. Incidence of mutator strains in Escherichia coli and coliforms in nature. Mutat Res. 1981;91(2):107–10. doi: 10.1016/0165-7992(81)90081-6 7019693.

61. Ellington MJ, Livermore DM, Pitt TL, Hall LM, Woodford N. Mutators among CTX-M beta-lactamase-producing Escherichia coli and risk for the emergence of fosfomycin resistance. J Antimicrob Chemother. 2006;58(4):848–52. doi: 10.1093/jac/dkl315 16891630.

62. Marinus MG. DNA Mismatch Repair. EcoSal Plus. 2012;5(1). doi: 10.1128/ecosalplus.7.2.5 26442827; PubMed Central PMCID: PMC4231543.

63. Petit MA, Dimpfl J, Radman M, Echols H. Control of large chromosomal duplications in Escherichia coli by the mismatch repair system. Genetics. 1991;129(2):327–32. 1743481; PubMed Central PMCID: PMC1204626.

64. Bzymek M, Saveson CJ, Feschenko VV, Lovett ST. Slipped misalignment mechanisms of deletion formation: in vivo susceptibility to nucleases. J Bacteriol. 1999;181(2):477–82. 9882661; PubMed Central PMCID: PMC93401.

65. Xie G, Keyhani NO, Bonner CA, Jensen RA. Ancient origin of the tryptophan operon and the dynamics of evolutionary change. Microbiol Mol Biol Rev. 2003;67(3):303–42. doi: 10.1128/MMBR.67.3.303-342.2003 12966138.

66. Fondi M, Emiliani G, Fani R. Origin and evolution of operons and metabolic pathways. Res Microbiol. 2009;160(7):502–12. doi: 10.1016/j.resmic.2009.05.001 19465116.

67. Reams AB, Neidle EL. Selection for gene clustering by tandem duplication. Annu Rev Microbiol. 2004;58:119–42. doi: 10.1146/annurev.micro.58.030603.123806 15487932.

68. Romero D, Palacios R. Gene amplification and genomic plasticity in prokaryotes. Annu Rev Genet. 1997;31:91–111. doi: 10.1146/annurev.genet.31.1.91 9442891.

69. Hooper SD, Berg OG. On the nature of gene innovation: duplication patterns in microbial genomes. Mol Biol Evol. 2003;20(6):945–54. doi: 10.1093/molbev/msg101 12716994.

70. Jordan IK, Makarova KS, Spouge JL, Wolf YI, Koonin EV. Lineage-specific gene expansions in bacterial and archaeal genomes. Genome Res. 2001;11(4):555–65. doi: 10.1101/gr.166001 11282971; PubMed Central PMCID: PMC311027.

71. Bosserman RE, Thompson CR, Nicholson KR, Champion PA. Esx paralogs are functionally equivalent to ESX-1 proteins but are dispensable for virulence in Mycobacterium marinum. J Bacteriol. 2018;200(11). doi: 10.1128/JB.00726-17 29555701; PubMed Central PMCID: PMC5952400.

72. Saier MH Jr., Paulsen IT, Sliwinski MK, Pao SS, Skurray RA, Nikaido H. Evolutionary origins of multidrug and drug-specific efflux pumps in bacteria. FASEB J. 1998;12(3):265–74. doi: 10.1096/fasebj.12.3.265 9506471.

73. Perrin E, Fondi M, Bosi E, Mengoni A, Buroni S, Scoffone VC, et al. Subfunctionalization influences the expansion of bacterial multidrug antibiotic resistance. BMC Genomics. 2017;18(1):834. doi: 10.1186/s12864-017-4222-4 29084524; PubMed Central PMCID: PMC5663151.

74. Shah S, Cannon JR, Fenselau C, Briken V. A duplicated ESAT-6 region of ESX-5 is involved in protein export and virulence of Mycobacteria. Infect Immun. 2015;83(11):4349–61. doi: 10.1128/IAI.00827-15 26303392; PubMed Central PMCID: PMC4598393.

75. Brosch R, Gordon SV, Buchrieser C, Pym AS, Garnier T, Cole ST. Comparative genomics uncovers large tandem chromosomal duplications in Mycobacterium bovis BCG Pasteur. Yeast. 2000;17(2):111–23. doi: 10.1002/1097-0061(20000630)17:2<111::AID-YEA17>3.0.CO;2-G 10900457; PubMed Central PMCID: PMC2448323.

76. Galamba A, Soetaert K, Wang XM, De Bruyn J, Jacobs P, Content J. Disruption of adhC reveals a large duplication in the Mycobacterium smegmatis mc(2)155 genome. Microbiol. 2001;147(Pt 12):3281–94. doi: 10.1099/00221287-147-12-3281 11739760.

77. Domenech P, Rog A, Moolji JU, Radomski N, Fallow A, Leon-Solis L, et al. Origins of a 350-kilobase genomic duplication in Mycobacterium tuberculosis and its impact on virulence. Infect Immun. 2014;82(7):2902–12. doi: 10.1128/IAI.01791-14 24778110; PubMed Central PMCID: PMC4097636.

78. Andersson SG, Zomorodipour A, Winkler HH, Kurland CG. Unusual organization of the rRNA genes in Rickettsia prowazekii. J Bacteriol. 1995;177(14):4171–5. doi: 10.1128/jb.177.14.4171-4175.1995 7608097; PubMed Central PMCID: PMC177156.

79. Bercovier H, Kafri O, Sela S. Mycobacteria possess a surprisingly small number of ribosomal RNA genes in relation to the size of their genome. Biochem Biophys Res Commun. 1986;136(3):1136–41. doi: 10.1016/0006-291x(86)90452-3 3013168.

80. Monod J. The Growth of Bacterial Cultures. Annu Rev Microbiol. 1949;3:371–94. doi: 10.1146/annurev.mi.03.100149.002103 WOS:A1949XS94800016.


Článek vyšel v časopise

PLOS Genetics


2020 Číslo 3
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Důležitost adherence při depresivním onemocnění
nový kurz
Autoři: MUDr. Eliška Bartečková, Ph.D.

Koncepce osteologické péče pro gynekology a praktické lékaře
Autoři: MUDr. František Šenk

Sekvenční léčba schizofrenie
Autoři: MUDr. Jana Hořínková, Ph.D.

Hypertenze a hypercholesterolémie – synergický efekt léčby
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Multidisciplinární zkušenosti u pacientů s diabetem
Autoři: Prof. MUDr. Martin Haluzík, DrSc., prof. MUDr. Vojtěch Melenovský, CSc., prof. MUDr. Vladimír Tesař, DrSc.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#