Causes of insertion sequences abundance in prokaryotic genomes

DYNAMICS OF AN INSERTION SEQUENCE INFECTION IN A SPATIALLY STRUCTURED ENVIRONMENT

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

Genomic repeats, misassembly and reannotation: a case study with long-read resequencing of Porphyromonas gingivalis reference strains

BACKGROUND

Without knowledge of their genomic sequences, it is impossible to make functional models of the bacteria that make up human and animal microbiota. Unfortunately, the vast majority of publicly available genomes are only working drafts, an incompleteness that causes numerous problems and constitutes a major obstacle to genotypic and phenotypic interpretation. In this work, we began with an example from the class Bacteroidia in the phylum Bacteroidetes, which is preponderant among human orodigestive microbiota. We successfully identify the genetic loci responsible for assembly breaks and misassemblies and demonstrate the importance and usefulness of long-read sequencing and curated reannotation.

RESULTS

We showed that the fragmentation in Bacteroidia draft genomes assembled from massively parallel sequencing linearly correlates with genomic repeats of the same or greater size than the reads. We also demonstrated that some of these repeats, especially the long ones, correspond to misassembled loci in three reference Porphyromonas gingivalis genomes marked as circularized (thus complete or finished). We prove that even at modest coverage (30X), long-read resequencing together with PCR contiguity verification (rrn operons and an integrative and conjugative element or ICE) can be used to identify and correct the wrongly combined or assembled regions. Finally, although time-consuming and labor-intensive, consistent manual biocuration of three P. gingivalis strains allowed us to compare and correct the existing genomic annotations, resulting in a more accurate interpretation of the genomic differences among these strains.

CONCLUSIONS

In this study, we demonstrate the usefulness and importance of long-read sequencing in verifying published genomes (even when complete) and generating assemblies for new bacterial strains/species with high genomic plasticity. We also show that when combined with biological validation processes and diligent biocurated annotation, this strategy helps reduce the propagation of errors in shared databases, thus limiting false conclusions based on incomplete or misleading information.

High abundance and expression of transposases in bacteria from the Baltic Sea

Transposases are mobile genetic elements suggested to have an important role in bacterial genome plasticity and host adaptation but their transcriptional activity in natural bacterial communities is largely unexplored. Here we analyzed metagenomes and -transcriptomes of size fractionated (0.1-0.8, 0.8-3.0 and 3.0-200 μm) bacterial communities from the brackish Baltic Sea, and adjacent marine waters. The Baltic Sea transposase levels, up to 1.7% of bacterial genes and 2% of bacterial transcripts, were considerably higher than in marine waters and similar to levels reported for extreme environments. Large variations in expression were found between transposase families and groups of bacteria, with a two-fold higher transcription in Cyanobacteria than in any other phylum. The community-level results were corroborated at the genus level by Synechococcus transposases reaching up to 5.2% of genes and 6.9% of transcripts, which is in contrast to marine Synechococcus that largely lack these genes. Levels peaked in Synechococcus from the largest size fraction, suggesting high frequencies of lateral gene transfer and high genome plasticity in colony-forming picocyanobacteria. Together, the results support an elevated rate of transposition-based genome change and adaptation in bacterial populations of the Baltic Sea, and possibly also of other highly dynamic estuarine waters.

Why Might Bacterial Pathogens Have Small Genomes?

Bacteria that cause serious disease often have smaller genomes, and fewer genes, than their nonpathogenic, or less pathogenic relatives. Here, we review evidence for the generality of this association, and summarise the various reasons why the association might hold. We focus on the population genetic processes that might lead to reductive genome evolution, and show how several of these could be connected to pathogenicity. We find some evidence for most of the processes having acted in bacterial pathogens, including several different modes of genome reduction acting in the same lineage. We argue that predictable processes of genome evolution might not reflect any common underlying process.

Genetic plasticity of the Shigella virulence plasmid is mediated by intra- and inter-molecular events between insertion sequences

Acquisition of a single copy, large virulence plasmid, pINV, led to the emergence of Shigella spp. from Escherichia coli. The plasmid encodes a Type III secretion system (T3SS) on a 30 kb pathogenicity island (PAI), and is maintained in a bacterial population through a series of toxin:antitoxin (TA) systems which mediate post-segregational killing (PSK). The T3SS imposes a significant cost on the bacterium, and strains which have lost the plasmid and/or genes encoding the T3SS grow faster than wild-type strains in the laboratory, and fail to bind the indicator dye Congo Red (CR). Our aim was to define the molecular events in Shigella flexneri that cause loss of Type III secretion (T3S), and to examine whether TA systems exert positional effects on pINV. During growth at 37°C, we found that deletions of regions of the plasmid including the PAI lead to the emergence of CR-negative colonies; deletions occur through intra-molecular recombination events between insertion sequences (ISs) flanking the PAI. Furthermore, by repositioning MvpAT (which belongs to the VapBC family of TA systems) near the PAI, we demonstrate that the location of this TA system alters the rearrangements that lead to loss of T3S, indicating that MvpAT acts both globally (by reducing loss of pINV through PSK) as well as locally (by preventing loss of adjacent sequences). During growth at environmental temperatures, we show for the first time that pINV spontaneously integrates into different sites in the chromosome, and this is mediated by inter-molecular events involving IS1294. Integration leads to reduced PAI gene expression and impaired secretion through the T3SS, while excision of pINV from the chromosome restores T3SS function. Therefore, pINV integration provides a reversible mechanism for Shigella to circumvent the metabolic burden imposed by pINV. Intra- and inter-molecular events between ISs, which are abundant in Shigella spp., mediate plasticity of S. flexneri pINV.

Community-like genome in single cells of the sulfur bacterium Achromatium oxaliferum

Polyploid bacteria are common, but the genetic and functional diversity resulting from polyploidy is unknown. Here we use single-cell genomics, metagenomics, single-cell amplicon sequencing, and fluorescence in situ hybridization, to show that individual cells of Achromatium oxaliferum, the world’s biggest known freshwater bacterium, harbor genetic diversity typical of whole bacterial communities. The cells contain tens of transposable elements, which likely cause the unprecedented diversity that we observe in the sequence and synteny of genes. Given the high within-cell diversity of the usually conserved 16S ribosomal RNA gene, we suggest that gene conversion occurs in multiple, separated genomic hotspots. The ribosomal RNA distribution inside the cells hints to spatially differential gene expression. We also suggest that intracellular gene transfer may lead to extensive gene reshuffling and increased diversity.

The cells of Achromatium bacteria are remarkably large and contain multiple chromosome copies. Here, Ionescu et al. show that chromosome copies within individual cells display high diversity, similar to that of bacterial communities, and contain tens of transposable elements.

How a Genetically Stable Extremophile Evolves: Modes of Genome Diversification in the Archaeon Sulfolobus acidocaldarius

In order to analyze in molecular terms how Sulfolobus genomes diverge, damage-induced mutations and natural polymorphisms (PMs) were identified in laboratory constructs and wild-type isolates, respectively, of Sulfolobus acidocaldarius Among wild-type isolates drawn from one local population, pairwise nucleotide divergence averaged 4 × 10^-6, which is about 0.15% of the corresponding divergence reported for Sulfolobus islandicus The most variable features of wild-type S. acidocaldarius genomes were homopolymer (mononucleotide) tracts and longer tandem repeats, consistent with the spontaneous mutations that occur under laboratory conditions. Natural isolates, however, also revealed large insertions/deletions and inversions, which did not occur in any of the laboratory-manipulated strains. Several of the large insertions/deletions could be attributed to the integration or excision of mobile genetic elements (MGEs), and each MGE represented a distinct system of site-specific recombination. The mode of recombination associated with one MGE, a provirus related to Sulfolobus turreted icosahedral virus, was also seen in certain chromosomal inversions. Artificially induced mutations, non-MGE insertions/deletions, and small PMs exhibited different distributions over the genome, suggesting that large-scale patterning of Sulfolobus genomes begins early in the divergence process. Unlike induced mutations, natural base pair substitutions occurred in clusters, and one cluster exhibited properties expected of nonreciprocal recombination (gene conversion) between dispersed imperfect repeats. Taken together, the results identify simple replication errors, slipped-strand events promoted by tandem repeats, homologous recombination, and rearrangements promoted by MGEs as the primary sources of genetic variation for this extremely acidophilic archaeon in its geothermal environment.IMPORTANCE The optimal growth temperatures of hyperthermophilic archaea accelerate DNA decomposition, which is expected to make DNA repair especially important for their genetic stability, yet these archaea lack certain broadly conserved types of DNA repair proteins. In this study, the genome of the extreme thermoacidophile Sulfolobus acidocaldarius was found to be remarkably stable, accumulating few mutations in many (though not all) laboratory manipulations and in natural populations. Furthermore, all the genetic processes that were inferred to diversify these genomes also operate in mesophilic bacteria and eukaryotes. This suggests that a common set of mechanisms produces most of the genetic variation in all microorganisms, despite the fundamental differences in physiology, DNA repair systems, and genome structure represented in the three domains of life.

Disentangling the effects of selection and loss bias on gene dynamics

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

The Evolution of Bacterial Genome Architecture

The genome architecture of bacteria and eukaryotes evolves in opposite directions when subject to genetic drift, a difference that can be ascribed to the fact that bacteria exhibit a mutational bias that deletes superfluous sequences, whereas eukaryotes are biased toward large insertions. Expansion of eukaryotic genomes occurs through the addition of non-functional sequences, such as repetitive sequences and transposable elements, whereas variation in bacterial genome size is largely due to the acquisition and loss of functional accessory genes. These properties create the situation in which eukaryotes with very similar numbers of genes can have vastly different genome sizes, while in bacteria, gene number scales linearly with genome size. Some bacterial genomes, however, particularly those of species that undergo bottlenecks due to recent association with hosts, accumulate pseudogenes and mobile elements, conferring them a low gene content relative to their genome size. These non-functional sequences are gradually eroded and eliminated after long-term association with hosts, with the result that obligate symbionts have the smallest genomes of any cellular organism. The architecture of bacterial genomes is shaped by complex and diverse processes, but for most bacterial species, genome size is governed by a non-adaptive process, i.e., genetic drift coupled with a mutational bias toward deletions. Thus, bacteria with small effective population sizes typically have the smallest genomes. Some marine bacteria counter this near-universal trend: despite having immense population sizes, selection, not drift, acts to reduce genome size in response to metabolic constraints in their nutrient-limited environment.

Comparative genomics of Vibrio campbellii strains and core species of the Vibrio Harveyi clade

The core of the Vibrio Harveyi clade contains V. harveyi, V. campbellii, V. owensii, V. jasicida, and V. rotiferianus. They are well recognized aquatic animal pathogens, but misclassification has been common due to similarities in their rDNA sequences and phenotypes. To better understand their evolutionary relationships and functional features, we sequenced a shrimp pathogen strain V. harveyi 1114GL, reclassified it as V. campbellii and compared this and 47 other sequenced Vibrio genomes in the Harveryi clade. A phylogeny based on 1,775 genes revealed that both V. owensii and V. jasicida were closer to V. campbellii than to V. harveyi and that V. campbellii strains can be divided into two distinct groups. Species-specific genes such as intimin and iron acquisition genes were identified in V. campbellii. In particular, the 1114GL strain contains two bacterial immunoglobulin-like genes for cell adhesion with 22 Big_2 domains that have been extensively reshuffled and are by far the most expanded among all species surveyed in this study. The 1114GL strain differed from ATCC BAA-1116 by ~9% at the synonymous sites, indicating high diversity within V. campbellii. Our study revealed the characteristics of V. campbellii in the Harveyi clade and the genetic basis for their wide-spread pathogenicity.

Arsenic resistance strategy in Pantoea sp. IMH: Organization, function and evolution of ars genes

Pantoea sp. IMH is the only bacterium found in genus Pantoea with a high As resistance capacity, but its molecular mechanism is unknown. Herein, the organization, function, and evolution of ars genes in IMH are studied starting with analysis of the whole genome. Two ars systems - ars1 (arsR1B1C1H1) and ars2 (arsR2B2C2H2) - with low sequence homology and two arsC-like genes, were found in the IMH genome. Both ars1 and ars2 are involved in the As resistance, where ars1 is the major contributor at 15 °C and ars2 at 30 °C. The difference in the behavior of these two ars systems is attributed to the disparate activities of their arsR promoters at different temperatures. Sequence analysis based on concatenated ArsRBC indicates that ars1 and ars2 clusters may be acquired from Franconibacter helveticus LMG23732 and Serratia marcescens (plasmid R478), respectively, by horizontal gene transfer (HGT). Nevertheless, two arsC-like genes, probably arising from the duplication of arsC, do not contribute to the As resistance. Our results indicate that Pantoea sp. IMH acquired two different As resistance genetic systems by HGT, allowing the colonization of changing ecosystems, and highlighting the flexible adaptation of microorganisms to resist As.

Insertion sequences enrichment in extreme Red sea brine pool vent

Mobile genetic elements are major agents of genome diversification and evolution. Limited studies addressed their characteristics, including abundance, and role in extreme habitats. One of the rare natural habitats exposed to multiple-extreme conditions, including high temperature, salinity and concentration of heavy metals, are the Red Sea brine pools. We assessed the abundance and distribution of different mobile genetic elements in four Red Sea brine pools including the world's largest known multiple-extreme deep-sea environment, the Red Sea Atlantis II Deep. We report a gradient in the abundance of mobile genetic elements, dramatically increasing in the harshest environment of the pool. Additionally, we identified a strong association between the abundance of insertion sequences and extreme conditions, being highest in the harshest and deepest layer of the Red Sea Atlantis II Deep. Our comparative analyses of mobile genetic elements in secluded, extreme and relatively non-extreme environments, suggest that insertion sequences predominantly contribute to polyextremophiles genome plasticity.

Genetic and life-history traits associated with the distribution of prophages in bacteria

Nearly half of the sequenced bacteria are lysogens and many of their prophages encode adaptive traits. Yet, the variables driving prophage distribution remain undetermined. We identified 2246 prophages in complete bacterial genomes to study the genetic and life-history traits associated with lysogeny. While optimal growth temperatures and average cell volumes were not associated with lysogeny, prophages were more frequent in pathogens and in bacteria with small minimal doubling times. Their frequency also increased with genome size, but only for genomes smaller than 6 Mb. The number of spacers in CRISPR-Cas systems and the frequency of type III systems were anticorrelated with prophage frequency, but lysogens were more likely to encode type I and type II systems. The minimal doubling time was the trait most correlated with lysogeny, followed by genome size and pathogenicity. We propose that bacteria with highly variable growth rates often encounter lower opportunity costs for lysogeny relative to lysis. These results contribute to explain the paucity of temperate phages in certain bacterial clades and of bacterial lysogens in certain environments. They suggest that genetic and life-history traits affect the contributions of temperate phages to bacterial genomes.

Regulation of genetic flux between bacteria by restriction-modification systems

Restriction-modification (R-M) systems are often regarded as bacteria's innate immune systems, protecting cells from infection by mobile genetic elements (MGEs). Their diversification has been recently associated with the emergence of particularly virulent lineages. However, we have previously found more R-M systems in genomes carrying more MGEs. Furthermore, it has been suggested that R-M systems might favor genetic transfer by producing recombinogenic double-stranded DNA ends. To test whether R-M systems favor or disfavor genetic exchanges, we analyzed their frequency with respect to the inferred events of homologous recombination and horizontal gene transfer within 79 bacterial species. Genetic exchanges were more frequent in bacteria with larger genomes and in those encoding more R-M systems. We created a recognition target motif predictor for Type II R-M systems that identifies genomes encoding systems with similar restriction sites. We found more genetic exchanges between these genomes, independently of their evolutionary distance. Our results reconcile previous studies by showing that R-M systems are more abundant in promiscuous species, wherein they establish preferential paths of genetic exchange within and between lineages with cognate R-M systems. Because the repertoire and/or specificity of R-M systems in bacterial lineages vary quickly, the preferential fluxes of genetic transfer within species are expected to constantly change, producing time-dependent networks of gene transfer.

Identification and analysis of integrons and cassette arrays in bacterial genomes

Integrons recombine gene arrays and favor the spread of antibiotic resistance. Their broader roles in bacterial adaptation remain mysterious, partly due to lack of computational tools. We made a program – IntegronFinder – to identify integrons with high accuracy and sensitivity. IntegronFinder is available as a standalone program and as a web application. It searches for attC sites using covariance models, for integron-integrases using HMM profiles, and for other features (promoters, attI site) using pattern matching. We searched for integrons, integron-integrases lacking attC sites, and clusters of attC sites lacking a neighboring integron-integrase in bacterial genomes. All these elements are especially frequent in genomes of intermediate size. They are missing in some key phyla, such as α-Proteobacteria, which might reflect selection against cell lineages that acquire integrons. The similarity between attC sites is proportional to the number of cassettes in the integron, and is particularly low in clusters of attC sites lacking integron-integrases. The latter are unexpectedly abundant in genomes lacking integron-integrases or their remains, and have a large novel pool of cassettes lacking homologs in the databases. They might represent an evolutionary step between the acquisition of genes within integrons and their stabilization in the new genome.

Microvirga massiliensis sp. nov., the human commensal with the largest genome

Microvirga massiliensis sp. nov. strain JC119(T) is a bacteria isolated in Marseille from a stool sample collected in Senegal. The 16S rRNA (JF824802) of M. massiliensis JC119(T) revealed 95% sequence identity with Microvirga lotononidis WSM3557(T) (HM362432). This bacterium is aerobic, gram negative, catalase positive, and oxidase negative. The draft genome of M. massiliensis JC119(T) comprises a 9,207,211-bp-long genome that is the largest bacterial genome of an isolate in humans. The genome exhibits a G+C content of 63.28% and contains 8685 protein-coding genes and 77 RNA genes, including 21 rRNA genes. Here, we describe the features of M. massiliensis JC119(T), together with the genome sequence information and its annotation.

Coevolution Theory of the Genetic Code at Age Forty: Pathway to Translation and Synthetic Life

The origins of the components of genetic coding are examined in the present study. Genetic information arose from replicator induction by metabolite in accordance with the metabolic expansion law. Messenger RNA and transfer RNA stemmed from a template for binding the aminoacyl-RNA synthetase ribozymes employed to synthesize peptide prosthetic groups on RNAs in the Peptidated RNA World. Coevolution of the genetic code with amino acid biosynthesis generated tRNA paralogs that identify a last universal common ancestor (LUCA) of extant life close to Methanopyrus, which in turn points to archaeal tRNA introns as the most primitive introns and the anticodon usage of Methanopyrus as an ancient mode of wobble. The prediction of the coevolution theory of the genetic code that the code should be a mutable code has led to the isolation of optional and mandatory synthetic life forms with altered protein alphabets.

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.

Coevolution of the Organization and Structure of Prokaryotic Genomes

The cytoplasm of prokaryotes contains many molecular machines interacting directly with the chromosome. These vital interactions depend on the chromosome structure, as a molecule, and on the genome organization, as a unit of genetic information. Strong selection for the organization of the genetic elements implicated in these interactions drives replicon ploidy, gene distribution, operon conservation, and the formation of replication-associated traits. The genomes of prokaryotes are also very plastic with high rates of horizontal gene transfer and gene loss. The evolutionary conflicts between plasticity and organization lead to the formation of regions with high genetic diversity whose impact on chromosome structure is poorly understood. Prokaryotic genomes are remarkable documents of natural history because they carry the imprint of all of these selective and mutational forces. Their study allows a better understanding of molecular mechanisms, their impact on microbial evolution, and how they can be tinkered in synthetic biology.

Comparative Genomics of Pathogenic and Nonpathogenic Strains of Xanthomonas arboricola Unveil Molecular and Evolutionary Events Linked to Pathoadaptation

The bacterial species Xanthomonas arboricola contains plant pathogenic and nonpathogenic strains. It includes the pathogen X. arboricola pv. juglandis, causing the bacterial blight of Juglans regia. The emergence of a new bacterial disease of J. regia in France called vertical oozing canker (VOC) was previously described and the causal agent was identified as a distinct genetic lineage within the pathovar juglandis. Symptoms on walnut leaves and fruits are similar to those of a bacterial blight but VOC includes also cankers on trunk and branches. In this work, we used comparative genomics and physiological tests to detect differences between four X. arboricola strains isolated from walnut tree: strain CFBP 2528 causing walnut blight (WB), strain CFBP 7179 causing VOC and two nonpathogenic strains, CFBP 7634 and CFBP 7651, isolated from healthy walnut buds. Whole genome sequence comparisons revealed that pathogenic strains possess a larger and wider range of mobile genetic elements than nonpathogenic strains. One pathogenic strain, CFBP 7179, possessed a specific integrative and conjugative element (ICE) of 95 kb encoding genes involved in copper resistance, transport and regulation. The type three effector repertoire was larger in pathogenic strains than in nonpathogenic strains. Moreover, CFBP 7634 strain lacked the type three secretion system encoding genes. The flagellar system appeared incomplete and nonfunctional in the pathogenic strain CFBP 2528. Differential sets of chemoreceptor and different repertoires of genes coding adhesins were identified between pathogenic and nonpathogenic strains. Besides these differences, some strain-specific differences were also observed. Altogether, this study provides valuable insights to highlight the mechanisms involved in ecology, environment perception, plant adhesion and interaction, leading to the emergence of new strains in a dynamic environment.

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.

The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes

Whole-genome sequences are now available for many microbial species and clades, however existing whole-genome alignment methods are limited in their ability to perform sequence comparisons of multiple sequences simultaneously. Here we present the Harvest suite of core-genome alignment and visualization tools for the rapid and simultaneous analysis of thousands of intraspecific microbial strains. Harvest includes Parsnp, a fast core-genome multi-aligner, and Gingr, a dynamic visual platform. Together they provide interactive core-genome alignments, variant calls, recombination detection, and phylogenetic trees. Using simulated and real data we demonstrate that our approach exhibits unrivaled speed while maintaining the accuracy of existing methods. The Harvest suite is open-source and freely available from: http://github.com/marbl/harvest.

ISQuest: finding insertion sequences in prokaryotic sequence fragment data

MOTIVATION

Insertion sequences (ISs) are transposable elements present in most bacterial and archaeal genomes that play an important role in genomic evolution. The increasing availability of sequenced prokaryotic genomes offers the opportunity to study ISs comprehensively, but development of efficient and accurate tools is required for discovery and annotation. Additionally, prokaryotic genomes are frequently deposited as incomplete, or draft stage because of the substantial cost and effort required to finish genome assembly projects. Development of methods to identify IS directly from raw sequence reads or draft genomes are therefore desirable. Software tools such as Optimized Annotation System for Insertion Sequences and IScan currently identify IS elements in completely assembled and annotated genomes; however, to our knowledge no methods have been developed to identify ISs from raw fragment data or partially assembled genomes. We have developed novel methods to solve this computationally challenging problem, and implemented these methods in the software package ISQuest. This software identifies bacterial ISs and their sequence elements-inverted and direct repeats-in raw read data or contigs using flexible search parameters. ISQuest is capable of finding ISs in hundreds of partially assembled genomes within hours, making it a valuable high-throughput tool for a global search of IS elements. We tested ISQuest on simulated read libraries of 3810 complete bacterial genomes and plasmids in GenBank and were capable of detecting 82% of the ISs and transposases annotated in GenBank with 80% sequence identity.

CONTACT

abiswas@cs.odu.edu.

Everyman's Guide to Bacterial Insertion Sequences

The number and diversity of known prokaryotic insertion sequences (IS) have increased enormously since their discovery in the late 1960s. At present the sequences of more than 4000 different IS have been deposited in the specialized ISfinder database. Over time it has become increasingly apparent that they are important actors in the evolution of their host genomes and are involved in sequestering, transmitting, mutating and activating genes, and in the rearrangement of both plasmids and chromosomes. This review presents an overview of our current understanding of these transposable elements (TE), their organization and their transposition mechanism as well as their distribution and genomic impact. In spite of their diversity, they share only a very limited number of transposition mechanisms which we outline here. Prokaryotic IS are but one example of a variety of diverse TE which are being revealed due to the advent of extensive genome sequencing projects. A major conclusion from sequence comparisons of various TE is that frontiers between the different types are becoming less clear. We detail these receding frontiers between different IS-related TE. Several, more specialized chapters in this volume include additional detailed information concerning a number of these.

In a second section of the review, we provide a detailed description of the expanding variety of IS, which we have divided into families for convenience. Our perception of these families continues to evolve and families emerge regularly as more IS are identified. This section is designed as an aid and a source of information for consultation by interested specialist readers.

Escherichia coli O-Antigen Gene Clusters of Serogroups O62, O68, O131, O140, O142, and O163: DNA Sequences and Similarity between O62 and O68, and PCR-Based Serogrouping

The DNA sequence of the O-antigen gene clusters of Escherichia coli serogroups O62, O68, O131, O140, O142, and O163 was determined, and primers based on the wzx (O-antigen flippase) and/or wzy (O-antigen polymerase) genes within the O-antigen gene clusters were designed and used in PCR assays to identify each serogroup. Specificity was tested with E. coli reference strains, field isolates belonging to the target serogroups, and non-E. coli bacteria. The PCR assays were highly specific for the respective serogroups; however, the PCR assay targeting the O62 wzx gene reacted positively with strains belonging to E. coli O68, which was determined by serotyping. Analysis of the O-antigen gene cluster sequences of serogroups O62 and O68 reference strains showed that they were 94% identical at the nucleotide level, although O62 contained an insertion sequence (IS) element located between the rmlA and rmlC genes within the O-antigen gene cluster. A PCR assay targeting the rmlA and rmlC genes flanking the IS element was used to differentiate O62 and O68 serogroups. The PCR assays developed in this study can be used for the detection and identification of E. coli O62/O68, O131, O140, O142, and O163 strains isolated from different sources.

Comparative Genomics of Pathogenic and Nonpathogenic Strains of Xanthomonas arboricola Unveil Molecular and Evolutionary Events Linked to Pathoadaptation

The bacterial species Xanthomonas arboricola contains plant pathogenic and nonpathogenic strains. It includes the pathogen X. arboricola pv. juglandis, causing the bacterial blight of Juglans regia. The emergence of a new bacterial disease of J. regia in France called vertical oozing canker (VOC) was previously described and the causal agent was identified as a distinct genetic lineage within the pathovar juglandis. Symptoms on walnut leaves and fruits are similar to those of a bacterial blight but VOC includes also cankers on trunk and branches. In this work, we used comparative genomics and physiological tests to detect differences between four X. arboricola strains isolated from walnut tree: strain CFBP 2528 causing walnut blight (WB), strain CFBP 7179 causing VOC and two nonpathogenic strains, CFBP 7634 and CFBP 7651, isolated from healthy walnut buds. Whole genome sequence comparisons revealed that pathogenic strains possess a larger and wider range of mobile genetic elements than nonpathogenic strains. One pathogenic strain, CFBP 7179, possessed a specific integrative and conjugative element (ICE) of 95 kb encoding genes involved in copper resistance, transport and regulation. The type three effector repertoire was larger in pathogenic strains than in nonpathogenic strains. Moreover, CFBP 7634 strain lacked the type three secretion system encoding genes. The flagellar system appeared incomplete and nonfunctional in the pathogenic strain CFBP 2528. Differential sets of chemoreceptor and different repertoires of genes coding adhesins were identified between pathogenic and nonpathogenic strains. Besides these differences, some strain-specific differences were also observed. Altogether, this study provides valuable insights to highlight the mechanisms involved in ecology, environment perception, plant adhesion and interaction, leading to the emergence of new strains in a dynamic environment.

The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes

Whole-genome sequences are now available for many microbial species and clades, however existing whole-genome alignment methods are limited in their ability to perform sequence comparisons of multiple sequences simultaneously. Here we present the Harvest suite of core-genome alignment and visualization tools for the rapid and simultaneous analysis of thousands of intraspecific microbial strains. Harvest includes Parsnp, a fast core-genome multi-aligner, and Gingr, a dynamic visual platform. Together they provide interactive core-genome alignments, variant calls, recombination detection, and phylogenetic trees. Using simulated and real data we demonstrate that our approach exhibits unrivaled speed while maintaining the accuracy of existing methods. The Harvest suite is open-source and freely available from: http://github.com/marbl/harvest.

Insertion sequence elements in Lactococcus garvieae

Insertion sequences are the simplest intracellular Mobile Genetic Elements which can occur in very high numbers in prokaryotic genomes, where they play an important evolutionary role by promoting genome plasticity. As such, the studies on the diversity and distribution of insertion sequences in genomes not yet investigated can contribute to improve the knowledge on a bacterial species and to identify new transposable elements. The present work describes the occurrence of insertion sequences in Lactococcus garvieae, an opportunistic emerging zoonotic and human pathogen, also associated with different food matrices. To date, no insertion elements have been described for L. garvieae in the IS element database. The analysis of the twelve published L. garvieae genomes identified 15 distinct insertion sequences that are members of the IS3, IS982, IS6, IS21 and IS256 families, including five new elements. Most of the insertion sequences in L. garvieae show substantial homology to the Lactococcus lactis elements, suggesting the movement of IS between these two species phylogenetically closely related. ISLL6 elements belonging to IS3 family were most abundant, with several copies distributed in 9 of the 12 genomes analyzed. An alignment analysis of two complete genomes carrying multi-copies of this insertion sequence indicates a possible involvement of ISLL6 in chromosomal rearrangement.

Mobile genetic elements in the bacterial phylum Acidobacteria

Analysis of the genome of Candidatus Solibacter usitatus Ellin6076, a member of the phylum Acidobacteria, revealed a large number of genes associated with mobile genetic elements. These genes encoded transposases, insertion sequence elements and phage integrases. When the amino acid sequences of the mobile element-associated genes were compared, many of them had high (90–100%) amino acid sequence identities, suggesting that these genes may have recently duplicated and dispersed throughout the genome. Although phage integrase encoding genes were prevalent in the Can. S. usitatus Ellin6076 genome, no intact prophage regions were found. This suggests that the Can. S. usitatus Ellin6076 large genome arose by horizontal gene transfer via ancient bacteriophage and/or plasmid-mediated transduction, followed by widespread small-scale gene duplications, resulting in an increased number of paralogs encoding traits that could provide selective metabolic, defensive and regulatory advantages in the soil environment. Here we examine the mobile element repertoire of Can. S. usitatus Ellin6076 in comparison to other genomes from the Acidobacteria phylum, reviewing published studies and contributing some new analyses. We also discuss the presence and potential roles of mobile elements in members of this phylum that inhabit a variety of environments.

From passengers to drivers: Impact of bacterial transposable elements on evolvability

Microbes have several mechanisms that promote evolutionary adaptation in stressful environments. The corresponding molecular pathways promote diversity through modulating rates of recombination, mutation or influence the activity of transposable genetic elements. Recent experimental studies suggest an evolutionary conflict between these mechanisms. Specifically, presence of mismatch repair mutator alleles in a bacterial population dramatically reduced fixation of bacterial insertion sequence elements. When rare, these elements had only a limited impact on adaptive evolution compared with other mutation-generating pathways. IS elements may initially spread like molecular parasites, but once present in many copies in a given genome, they might become generators of novelty during bacterial evolution.

Activity spectrum of colicins produced by Shigella sonnei and genetic mechanism of colicin resistance in conspecific S. sonnei strains and Escherichia coli

Colicin-mediated killing is an example of allelopathy, which has been found among several bacteria. Screening of 42 strains of Shigella sonnei isolated from diarrheal patients revealed that 39 (93%) S. sonnei strains were positive for colicin production against Escherichia coli DH5α. In the PCR-based detection of the colicin types, 36 (92.3%) were identified as E3, 2 (5.1%) as E3 and E8, and 1 (2.6%) as E3 and E2. Representative S. sonnei strains producing heterologous colicins exhibited antagonism against diarrheagenic Escherichia coli (DEC) groups. Although it is known that mutation in the colicin receptor renders the host resistant to colicin, there is a dearth of information on the genetic characterization of such mutants. In the fluctuation test, colicin-resistant E. coli mutants were found to occur spontaneously at the rates of 2.51 × 10(-8) and 5.52 × 10(-8) per generation when exposed to colicins E3 and E8 and colicins E3 and E2, respectively. Genotypic characterization of colicin-resistant E. coli (EC(Cr)) and S. sonnei (SS(Cr)) strains displayed mutations in the btuB gene, which encodes the receptor for vitamin B12 uptake. This gene was interrupted by various insertion sequences, such as IS1, IS2, and IS911. Complementation of EC(Cr) and SS(Cr) with plasmid-borne btuB (pbtuB) accomplished restoration of the colicin-susceptible phenotype. The vitamin B12 uptake assay gave an insight into the physiological relevance of the btuB mutation. Our studies provide insights into the latent influence of S. sonnei colicins in governing the existence of some of the shigellae and all of the DEC and the genetic mechanism underlying the emergence of resistance.

Natural selection for operons depends on genome size

In prokaryotes, genome size is associated with metabolic versatility, regulatory complexity, effective population size, and horizontal transfer rates. We therefore analyzed the covariation of genome size and operon conservation to assess the evolutionary models of operon formation and maintenance. In agreement with previous results, intraoperonic pairs of essential and of highly expressed genes are more conserved. Interestingly, intraoperonic pairs of genes are also more conserved when they encode proteins at similar cell concentrations, suggesting a role of cotranscription in diminishing the cost of waste and shortfall in gene expression. Larger genomes have fewer and smaller operons that are also less conserved. Importantly, lower conservation in larger genomes was observed for all classes of operons in terms of gene expression, essentiality, and balanced protein concentration. We reached very similar conclusions in independent analyses of three major bacterial clades (α- and β-Proteobacteria and Firmicutes). Operon conservation is inversely correlated to the abundance of transcription factors in the genome when controlled for genome size. This suggests a negative association between the complexity of genetic networks and operon conservation. These results show that genome size and/or its proxies are key determinants of the intensity of natural selection for operon organization. Our data fit better the evolutionary models based on the advantage of coregulation than those based on genetic linkage or stochastic gene expression. We suggest that larger genomes with highly complex genetic networks and many transcription factors endure weaker selection for operons than smaller genomes with fewer alternative tools for genetic regulation.

Population genomics of transposable elements in Drosophila

Studies of the population dynamics of transposable elements (TEs) in Drosophila melanogaster indicate that consistent forces are affecting TEs independently of their modes of transposition and regulation. New sequencing technologies enable biologists to sample genomes at an unprecedented scale in order to quantify genome-wide polymorphism for annotated and novel TE insertions. In this review, we first present new insights gleaned from high-throughput data for population genomics studies of D. melanogaster. We then consider the latest population genomics models for TE evolution and present examples of functional evidence revealed by genome-wide studies of TE population dynamics in D. melanogaster. Although most of the TE insertions are deleterious or neutral, some TE insertions increase the fitness of the individual that carries them and play a role in genome adaptation.

The genome of Paenibacillus sabinae T27 provides insight into evolution, organization and functional elucidation of nif and nif-like genes

Background

Most biological nitrogen fixation is catalyzed by the molybdenum nitrogenase. This enzyme is a complex which contains the MoFe protein encoded by nifDK and the Fe protein encoded by nifH. In addition to nifHDK, nifHDK-like genes were found in some Archaea and Firmicutes, but their function is unclear.

Results

We sequenced the genome of Paenibacillus sabinae T27. A total of 4,793 open reading frames were predicted from its 5.27 Mb genome. The genome of P. sabinae T27 contains fifteen nitrogen fixation (nif) genes, including three nifH, one nifD, one nifK, four nifB, two nifE, two nifN, one nifX and one nifV. Of the 15 nif genes, eight nif genes (nifB, nifH, nifD, nifK, nifE, nifN, nifX and nifV) and two non-nif genes (orf1 and hesA) form a complete nif gene cluster. In addition to the nif genes, there are nitrogenase-like genes, including two nifH-like genes and five pairs of nifDK-like genes. IS elements on the flanking regions of nif and nif-like genes imply that these genes might have been obtained by horizontal gene transfer. Phylogenies of the concatenated 8 nif gene (nifB, nifH, nifD, nifK, nifE, nifN, nifX and nifV) products suggest that P. sabinae T27 is closely related to Frankia. RT-PCR analysis showed that the complete nif gene cluster is organized as an operon. We demonstrated that the complete nif gene cluster under the control of σ⁷⁰-dependent promoter enabled Escherichia coli JM109 to fix nitrogen. Also, here for the first time we demonstrated that unlike nif genes, the transcriptions of nifHDK-like genes were not regulated by ammonium and oxygen, and nifH-like or nifD-like gene could not restore the nitrogenase activity of Klebsiella pneumonia nifH⁻ and nifD⁻ mutant strains, respectively, suggesting that nifHDK-like genes were not involved in nitrogen fixation.

Conclusions

Our data and analysis reveal the contents and distribution of nif and nif-like genes and contribute to the study of evolutionary history of nitrogen fixation in Paenibacillus. For the first time we demonstrated that the transcriptions of nifHDK-like genes were not regulated by ammonium and oxygen and nifHDK-like genes were not involved in nitrogen fixation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-723) contains supplementary material, which is available to authorized users.

The primary transcriptome of the marine diazotroph Trichodesmium erythraeum IMS101

Blooms of the dinitrogen-fixing marine cyanobacterium Trichodesmium considerably contribute to new nitrogen inputs into tropical oceans. Intriguingly, only 60% of the Trichodesmium erythraeum IMS101 genome sequence codes for protein, compared with ~85% in other sequenced cyanobacterial genomes. The extensive non-coding genome fraction suggests space for an unusually high number of unidentified, potentially regulatory non-protein-coding RNAs (ncRNAs). To identify the transcribed fraction of the genome, here we present a genome-wide map of transcriptional start sites (TSS) at single nucleotide resolution, revealing the activity of 6,080 promoters. We demonstrate that T. erythraeum has the highest number of actively splicing group II introns and the highest percentage of TSS yielding ncRNAs of any bacterium examined to date. We identified a highly transcribed retroelement that serves as template repeat for the targeted mutation of at least 12 different genes by mutagenic homing. Our findings explain the non-coding portion of the T. erythraeum genome by the transcription of an unusually high number of non-coding transcripts in addition to the known high incidence of transposable elements. We conclude that riboregulation and RNA maturation-dependent processes constitute a major part of the Trichodesmium regulatory apparatus.

Genome Sequence of Lactobacillus delbrueckii subsp. lactis CNRZ327, a Dairy Bacterium with Anti-Inflammatory Properties

Lactobacillus delbrueckii subsp. lactis CNRZ327 is a dairy bacterium with anti-inflammatory properties both in vitro and in vivo. Here, we report the genome sequence of this bacterium, which appears to contain no less than 215 insertion sequence (IS) elements, an exceptionally high number regarding the small genome size of the strain.

Large-scale genomic analysis suggests a neutral punctuated dynamics of transposable elements in bacterial genomes

Insertion sequences (IS) are the simplest and most abundant form of transposable DNA found in bacterial genomes. When present in multiple copies, it is thought that they can promote genomic plasticity and genetic exchange, thus being a major force of evolutionary change. The main processes that determine IS content in genomes are, though, a matter of debate. In this work, we take advantage of the large amount of genomic data currently available and study the abundance distributions of 33 IS families in 1811 bacterial chromosomes. This allows us to test simple models of IS dynamics and estimate their key parameters by means of a maximum likelihood approach. We evaluate the roles played by duplication, lateral gene transfer, deletion and purifying selection. We find that the observed IS abundances are compatible with a neutral scenario where IS proliferation is controlled by deletions instead of purifying selection. Even if there may be some cases driven by selection, neutral behavior dominates over large evolutionary scales. According to this view, IS and hosts tend to coexist in a dynamic equilibrium state for most of the time. Our approach also allows for a detection of recent IS expansions, and supports the hypothesis that rapid expansions constitute transient events—punctuations—during which the state of coexistence of IS and host becomes perturbated.

Insertion sequences (IS) are mobile genetic elements found in most prokaryotic genomes. They are able to autonomously change position and proliferate in chromosomes. The nature of the coevolutionary dynamics of IS with the genome that hosts them is a matter of debate: Do IS proliferate to the point of causing the extinction of the host? Is it possible that IS and hosts stably coexist? Can environmental perturbations cause IS expansions? What is the role of selection in controlling IS copy number? In this study, we have analysed abundance patterns of IS families to test two different evolutionary hypotheses: in the first one IS evolve neutrally, while in the second case they are affected by selection. Our results indicate that, most of the time, IS and their hosts coexist stably in a neutral scenario where the proliferation of IS through duplications and lateral gene transfer is balanced by regular deletions. Occasionally, though, this balance may be disrupted, causing temporary explosions of IS abundance.

Implications of streamlining theory for microbial ecology

Whether a small cell, a small genome or a minimal set of chemical reactions with self-replicating properties, simplicity is beguiling. As Leonardo da Vinci reportedly said, 'simplicity is the ultimate sophistication'. Two diverging views of simplicity have emerged in accounts of symbiotic and commensal bacteria and cosmopolitan free-living bacteria with small genomes. The small genomes of obligate insect endosymbionts have been attributed to genetic drift caused by small effective population sizes (Ne). In contrast, streamlining theory attributes small cells and genomes to selection for efficient use of nutrients in populations where Ne is large and nutrients limit growth. Regardless of the cause of genome reduction, lost coding potential eventually dictates loss of function. Consequences of reductive evolution in streamlined organisms include atypical patterns of prototrophy and the absence of common regulatory systems, which have been linked to difficulty in culturing these cells. Recent evidence from metagenomics suggests that streamlining is commonplace, may broadly explain the phenomenon of the uncultured microbial majority, and might also explain the highly interdependent (connected) behavior of many microbial ecosystems. Streamlining theory is belied by the observation that many successful bacteria are large cells with complex genomes. To fully appreciate streamlining, we must look to the life histories and adaptive strategies of cells, which impose minimum requirements for complexity that vary with niche.

Birth, death, and diversification of mobile promoters in prokaryotes

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

Bacterial insertion sequences: their genomic impact and diversity

Insertion sequences (ISs), arguably the smallest and most numerous autonomous transposable elements (TEs), are important players in shaping their host genomes. This review focuses on prokaryotic ISs. We discuss IS distribution and impact on genome evolution. We also examine their effects on gene expression, especially their role in activating neighbouring genes, a phenomenon of particular importance in the recent upsurge of bacterial antibiotic resistance. We explain how ISs are identified and classified into families by a combination of characteristics including their transposases (Tpases), their overall genetic organisation and the accessory genes which some ISs carry. We then describe the organisation of autonomous and nonautonomous IS‐related elements. This is used to illustrate the growing recognition that the boundaries between different types of mobile element are becoming increasingly difficult to define as more are being identified. We review the known Tpase types, their different catalytic activities used in cleaving and rejoining DNA strands during transposition, their organisation into functional domains and the role of this in regulation. Finally, we consider examples of prokaryotic IS domestication. In a more speculative section, we discuss the necessity of constructing more quantitative dynamic models to fully appreciate the continuing impact of TEs on prokaryotic populations.

Insertion sequence distribution bias in Archaea

Insertion sequences (IS) are common transposable elements in Archaea. Intergenic IS elements are usually less harmful than intragenic ISs, simply because they are less likely to disrupt host gene function. However, because regulatory sequences are intergenic and upstream of genes, we hypothesized that not all intergenic regions are selectively equivalent for IS insertion. We tested this hypothesis by analyzing the distributions of intergenic IS elements within 155 fully sequenced archaeal genomes. Of the 22 genomes with enough IS elements for statistical analysis, five have significantly fewer ISs between divergently oriented neighboring genes than expected by chance, and seven have significantly more ISs between convergently oriented genes. Furthermore, of the 85 genomes with at least one expected IS within each of the three possible neighboring gene orientations (i.e., divergent, convergent, and tandem), 73 genomes have fewer ISs between divergently oriented genes than expected, and 60 have more ISs between convergently oriented genes than expected (both values deviate significantly from binomial probabilities of random distribution). We suspect that these non-random IS distributions are molded by natural selection resulting from differential disruption of neighboring gene regulation, and that this selective pressure has affected transposable element distributions in prokaryotes for billions of years.

Comparison of assembled Clostridium botulinum A1 genomes revealed their evolutionary relationship

Clostridium botulinum encompasses bacteria that produce at least one of the seven serotypes of botulinum neurotoxin (BoNT/A-G). The availability of genome sequences of four closely related Type A1 or A1(B) strains, as well as the A1-specific microarray, allowed the analysis of their genomic organizations and evolutionary relationship. The four genomes share >90% core genes and >96% functional groups. Phylogenetic analysis based on COG shows closer relations of the A1(B) strain, NCTC 2916, to B1 and F1 than A1 strains. Alignment of the genomes of the three A1 strains revealed a highly similar chromosomal structure with three small gaps in the genome of ATCC 19397 and one additional gap in the genome of Hall A, suggesting ATCC 19379 as an evolutionary intermediate between Hall A and ATCC 3502. Analyses of the four gap regions indicated potential horizontal gene transfer and recombination events important for the evolution of A1 strains.

Growth temperature and genome size in bacteria are negatively correlated, suggesting genomic streamlining during thermal adaptation

Prokaryotic genomes are small and compact. Either this feature is caused by neutral evolution or by natural selection favoring small genomes—genome streamlining. Three separate prior lines of evidence argue against streamlining for most prokaryotes. We find that the same three lines of evidence argue for streamlining in the genomes of thermophile bacteria. Specifically, with increasing habitat temperature and decreasing genome size, the proportion of genomic DNA in intergenic regions decreases. Furthermore, with increasing habitat temperature, generation time decreases. Genome-wide selective constraints do not decrease as in the reduced genomes of host-associated species. Reduced habitat variability is not a likely explanation for the smaller genomes of thermophiles. Genome size may be an indirect target of selection due to its association with cell volume. We use metabolic modeling to demonstrate that known changes in cell structure and physiology at high temperature can provide a selective advantage to reduce cell volume at high temperatures.

Horizontal transfer and evolution of prokaryote transposable elements in eukaryotes

Horizontal transfer (HT) of transposable elements (TEs) plays a key role in prokaryotic evolution, and mounting evidence suggests that it has also had an important impact on eukaryotic evolution. Although many prokaryote-to-prokaryote and eukaryote-to-eukaryote HTs of TEs have been characterized, only few cases have been reported between prokaryotes and eukaryotes. Here, we carried out a comprehensive search for all major groups of prokaryotic insertion sequences (ISs) in 430 eukaryote genomes. We uncovered a total of 80 sequences, all deriving from the IS607 family, integrated in the genomes of 14 eukaryote species belonging to four distinct phyla (Amoebozoa, Ascomycetes, Basidiomycetes, and Stramenopiles). Given that eukaryote IS607-like sequences are most closely related to cyanobacterial IS607 and that their phylogeny is incongruent with that of their hosts, we conclude that the presence of IS607-like sequences in eukaryotic genomes is the result of several HT events. Selection analyses further suggest that our ability to detect these prokaryote TEs today in eukaryotes is because HT of these sequences occurred recently and/or some IS607 elements were domesticated after HT, giving rise to new eukaryote genes. Supporting the recent age of some of these HTs, we uncovered intact full-length, potentially active IS607 copies in the amoeba Acanthamoeba castellani. Overall, our study shows that prokaryote-to-eukaryote HT of TEs occurred at relatively low frequency during recent eukaryote evolution and it sets IS607 as the most widespread TE (being present in prokaryotes, eukaryotes, and viruses).

Rates of transposition in Escherichia coli

The evolutionary role of transposable elements (TEs) is still highly controversial. Two key parameters, the transposition rate (u and w, for replicative and non-replicative transposition) and the excision rate (e) are fundamental to understanding their evolution and maintenance in populations. We have estimated u, w and e for six families of TEs (including eight members: IS1, IS2, IS3, IS4, IS5, IS30, IS150 and IS186) in Escherichia coli, using a mutation accumulation (MA) experiment. In this experiment, mutations accumulate essentially at the rate at which they appear, during a period of 80 500 (1610 generations × 50 lines) generations, and spontaneous transposition events can be detected. This differs from other experiments in which insertions accumulated under strong selective pressure or over a limited genomic target. We therefore provide new estimates for the spontaneous rates of transposition and excision in E. coli. We observed 25 transposition and three excision events in 50 MA lines, leading to overall rate estimates of u ∼ 1.15 × 10(-5), w ∼ 4 × 10(-8) and e ∼ 1.08 × 10(-6) (per element, per generation). Furthermore, extensive variation between elements was found, consistent with previous knowledge of the mechanisms and regulation of transposition for the different elements.

Improving prokaryotic transposable elements identification using a combination of de novo and profile HMM methods

BACKGROUND

Insertion Sequences (ISs) and their non-autonomous derivatives (MITEs) are important components of prokaryotic genomes inducing duplication, deletion, rearrangement or lateral gene transfers. Although ISs and MITEs are relatively simple and basic genetic elements, their detection remains a difficult task due to their remarkable sequence diversity. With the advent of high-throughput genome and metagenome sequencing technologies, the development of fast, reliable and sensitive methods of ISs and MITEs detection become an important challenge. So far, almost all studies dealing with prokaryotic transposons have used classical BLAST-based detection methods against reference libraries. Here we introduce alternative methods of detection either taking advantages of the structural properties of the elements (de novo methods) or using an additional library-based method using profile HMM searches.

RESULTS

In this study, we have developed three different work flows dedicated to ISs and MITEs detection: the first two use de novo methods detecting either repeated sequences or presence of Inverted Repeats; the third one use 28 in-house transposase alignment profiles with HMM search methods. We have compared the respective performances of each method using a reference dataset of 30 archaeal and 30 bacterial genomes in addition to simulated and real metagenomes. Compared to a BLAST-based method using ISFinder as library, de novo methods significantly improve ISs and MITEs detection. For example, in the 30 archaeal genomes, we discovered 30 new elements (+20%) in addition to the 141 multi-copies elements already detected by the BLAST approach. Many of the new elements correspond to ISs belonging to unknown or highly divergent families. The total number of MITEs has even doubled with the discovery of elements displaying very limited sequence similarities with their respective autonomous partners (mainly in the Inverted Repeats of the elements). Concerning metagenomes, with the exception of short reads data (<300 bp) for which both techniques seem equally limited, profile HMM searches considerably ameliorate the detection of transposase encoding genes (up to +50%) generating low level of false positives compare to BLAST-based methods.

CONCLUSION

Compared to classical BLAST-based methods, the sensitivity of de novo and profile HMM methods developed in this study allow a better and more reliable detection of transposons in prokaryotic genomes and metagenomes. We believed that future studies implying ISs and MITEs identification in genomic data should combine at least one de novo and one library-based method, with optimal results obtained by running the two de novo methods in addition to a library-based search. For metagenomic data, profile HMM search should be favored, a BLAST-based step is only useful to the final annotation into groups and families.

Manipulating or superseding host recombination functions: a dilemma that shapes phage evolvability

Phages, like many parasites, tend to have small genomes and may encode autonomous functions or manipulate those of their hosts'. Recombination functions are essential for phage replication and diversification. They are also nearly ubiquitous in bacteria. The E. coli genome encodes many copies of an octamer (Chi) motif that upon recognition by RecBCD favors repair of double strand breaks by homologous recombination. This might allow self from non-self discrimination because RecBCD degrades DNA lacking Chi. Bacteriophage Lambda, an E. coli parasite, lacks Chi motifs, but escapes degradation by inhibiting RecBCD and encoding its own autonomous recombination machinery. We found that only half of 275 lambdoid genomes encode recombinases, the remaining relying on the host's machinery. Unexpectedly, we found that some lambdoid phages contain extremely high numbers of Chi motifs concentrated between the phage origin of replication and the packaging site. This suggests a tight association between replication, packaging and RecBCD-mediated recombination in these phages. Indeed, phages lacking recombinases strongly over-represent Chi motifs. Conversely, phages encoding recombinases and inhibiting host recombination machinery select for the absence of Chi motifs. Host and phage recombinases use different mechanisms and the latter are more tolerant to sequence divergence. Accordingly, we show that phages encoding their own recombination machinery have more mosaic genomes resulting from recent recombination events and have more diverse gene repertoires, i.e. larger pan genomes. We discuss the costs and benefits of superseding or manipulating host recombination functions and how this decision shapes phage genome structure and evolvability.

Bacterial viruses, called bacteriophages, are extremely abundant in the biosphere. They have key roles in the regulation of bacterial populations and in the diversification of bacterial genomes. Among these viruses, lambdoid phages are very abundant in enterobacteria and exchange genetic material very frequently. This latter process is thought to increase phage diversity and therefore facilitate adaptation to hosts. Recombination is also essential for the replication of many lambdoid phages. Lambdoids have been described to encode their own recombination genes and inhibit their hosts'. In this study, we show that lambdoids are split regarding their capacity to encode autonomous recombination functions and that this affects the abundance of recombination-related sequence motifs. Half of the phages encode an autonomous system and inhibit their hosts'. The trade-off between superseding and manipulating the hosts' recombination functions has important consequences. The phages encoding autonomous recombination functions have more diverse gene repertoires and recombine more frequently. Viruses, as many other parasites, have small genomes and depend on their hosts for several housekeeping functions. Hence, they often face trade-offs between supersession and manipulation of molecular machineries. Our results suggest these trade-offs may shape viral gene repertoires, their sequence composition and even influence their evolvability.

GTAG- and CGTC-tagged palindromic DNA repeats in prokaryotes

BACKGROUND

REPs (Repetitive Extragenic Palindromes) are small (20-40 bp) palindromic repeats found in high copies in some prokaryotic genomes, hypothesized to play a role in DNA supercoiling, transcription termination, mRNA stabilization.

RESULTS

We have monitored a large number of REP elements in prokaryotic genomes, and found that most can be sorted into two large DNA super-families, as they feature at one end unpaired motifs fitting either the GTAG or the CGTC consensus. Tagged REPs have been identified in >80 species in 8 different phyla. GTAG and CGTC repeats reside predominantly in microorganisms of the gamma and alpha division of Proteobacteria, respectively. However, the identification of members of both super- families in deeper branching phyla such Cyanobacteria and Planctomycetes supports the notion that REPs are old components of the bacterial chromosome. On the basis of sequence content and overall structure, GTAG and CGTC repeats have been assigned to 24 and 4 families, respectively. Of these, some are species-specific, others reside in multiple species, and several organisms contain different REP types. In many families, most units are close to each other in opposite orientation, and may potentially fold into larger secondary structures. In different REP-rich genomes the repeats are predominantly located between unidirectionally and convergently transcribed ORFs. REPs are predominantly located downstream from coding regions, and many are plausibly transcribed and function as RNA elements. REPs located inside genes have been identified in several species. Many lie within replication and global genome repair genes. It has been hypothesized that GTAG REPs are miniature transposons mobilized by specific transposases known as RAYTs (REP associated tyrosine transposases). RAYT genes are flanked either by GTAG repeats or by long terminal inverted repeats (TIRs) unrelated to GTAG repeats. Moderately abundant families of TIRs have been identified in multiple species.

CONCLUSIONS

CGTC REPs apparently lack a dedicated transposase. Future work will clarify whether these elements may be mobilized by RAYTs or other transposases, and assess if de-novo formation of either GTAG or CGTC repeats type still occurs.

A new subgroup of the IS3 family and properties of its representative member ISPpy1

Recently, we described a novel insertion element, ISPpy1, isolated from a permafrost strain of Psychrobacter maritimus. In this work, we demonstrated that ISPpy1 is a member of a novel subgroup of the IS3 family of insertion sequences (ISs) that was not identified and characterized previously. IS elements of this subgroup termed the ISPpy1 subgroup are broadly distributed among different taxa of Eubacteria, including Geobacteraceae, Chlorobiaceae, Desulfobacteraceae, Methylobacteriaceae, Nitrosomonadaceae and Cyanobacteria. While displaying characteristic features of the IS3-family elements, ISPpy1 subgroup elements exhibit some unusual features. In particular, most of them have longer terminal repeats with unconventional ends and frameshifting box with an atypical organization, and, unlike many other IS3-family elements, do not exhibit any distinct IS specificity. We studied the transposition and mutagenic properties of a representative member of this subgroup, ISPpy1 and showed that in contrast to the original P. maritimus host, in a heterologous host, Escherichia coli K-12, it is able to translocate with extremely high efficiency into the chromosome, either by itself or as a part of a composite transposon containing two ISPpy1 copies. The majority of transposants carry multiple chromosomal copies (up to 12) of ISPpy1. It was discovered that ISPpy1 is characterized by a marked mutagenic activity in E. coli: its chromosomal insertions generate various types of mutations, including auxotrophic, pleiotropic and rifampicin-resistance mutations. The distribution of IS elements of the novel subgroup among different bacteria, their role in the formation of composite transposons and the horizontal transfer of genes are examined and discussed.