Distant horizontal gene transfer is rare for multiple families of prokaryotic insertion sequences

Pathogenicity and virulence of Shigella sonnei: A highly drug-resistant pathogen of increasing prevalence

Shigella spp. are the causative agent of shigellosis (or bacillary dysentery), a diarrhoeal disease characterized for the bacterial invasion of gut epithelial cells. Among the 4 species included in the genus, Shigella flexneri is principally responsible for the disease in the developing world while Shigella sonnei is the main causative agent in high-income countries. Remarkably, as more countries improve their socioeconomic conditions, we observe an increase in the relative prevalence of S. sonnei. To date, the reasons behind this change in aetiology depending on economic growth are not understood. S. flexneri has been widely used as a model to study the pathogenesis of the genus, but as more research data are collected, important discrepancies with S. sonnei have come to light. In comparison to S. flexneri, S. sonnei can be differentiated in numerous aspects; it presents a characteristic O-antigen identical to that of one serogroup of the environmental bacterium Plesiomonas shigelloides, a group 4 capsule, antibacterial mechanisms to outcompete and displace gut commensal bacteria, and a poorer adaptation to an intracellular lifestyle. In addition, the World Health Organization (WHO) have recognized the significant threat posed by antibiotic-resistant strains of S. sonnei, demanding new approaches. This review gathers knowledge on what is known about S. sonnei within the context of other Shigella spp. and aims to open the door for future research on understanding the increasing spread of this pathogen.

No Isolate, No Problem: Using a Novel Insertion Sequence PCR to Link Rats to Human Shigellosis Cases in an Underserved Urban Community

During an investigation into a cluster of Shigella flexneri serotype 2a cases in an underserved community, we assessed the relatedness of human and rat S. flexneri isolates utilizing a novel PCR targeting insertion sites (IS-PCR) of mobile elements in the Shigella genome characteristic of the cluster strain. Whole-genome sequences of S. flexneri (n = 50) associated with the cluster were analyzed. De novo genome assemblies were analyzed by a Geneious V10.2.6 motif search, and two unique IS were identified in all human Shigella sequences of the local cluster. Hydrolysis probe PCR assays were designed to detect these sequences consisting of forward and reverse primers to amplify across each insertion site and a hydrolysis probe spanning the insertion site. IS-PCR was performed for three Shigella PCR-positive culture-negative rat intestine specimens from this community. Both insertion sites were detected in the de novo genome assemblies of all clinical S. flexneri isolates (n = 50). Two of the three PCR-positive culture-negative rat samples were positive for both unique ISs identified in the human S. flexneri isolates, suggesting that the rat Shigella species strains were closely related to the human strains in the cluster. The cycle threshold (Ct) values were >35, indicating that the bacterial load was very low in the rat samples. Two unique IS were identified in clinical isolates from a community S. flexneri cluster. Both IS targets were identified in PCR-positive (Shigella spp.), culture-negative rat tissue and clinical isolates from humans, indicating relatedness. IMPORTANCE This article describes a novel molecular method to show relatedness between bacterial infections, which may not be able to grow in the laboratory due to treatment with antibiotics or for bacteria requiring unique conditions to grow well. Uniquely, we applied this technique to Shigella isolates from human cases associated with a local cluster in an underserved community, as well as rat samples from the same community. We believe that this novel approach can serve as a complementary method to support outbreak/cluster investigation for Shigella spp.

Insertion sequence contributes to the evolution and environmental adaptation of Acidithiobacillus

BACKGROUND

The genus Acidithiobacillus has been widely concerned due to its superior survival and oxidation ability in acid mine drainage (AMD). However, the contribution of insertion sequence (IS) to their biological evolution and environmental adaptation is very limited. ISs are the simplest kinds of mobile genetic elements (MGEs), capable of interrupting genes, operons, or regulating the expression of genes through transposition activity. ISs could be classified into different families with their own members, possessing different copies.

RESULTS

In this study, the distribution and evolution of ISs, as well as the functions of the genes around ISs in 36 Acidithiobacillus genomes, were analyzed. The results showed that 248 members belonging to 23 IS families with a total of 10,652 copies were identified within the target genomes. The IS families and copy numbers among each species were significantly different, indicating that the IS distribution of Acidithiobacillus were not even. A. ferrooxidans had 166 IS members, which may develop more gene transposition strategies compared with other Acidithiobacillus spp. What's more, A. thiooxidans harbored the most IS copies, suggesting that their ISs were the most active and more likely to transpose. The ISs clustered in the phylogenetic tree approximately according to the family, which were mostly different from the evolutionary trends of their host genomes. Thus, it was suggested that the recent activity of ISs of Acidithiobacillus was not only determined by their genetic characteristics, but related with the environmental pressure. In addition, many ISs especially Tn3 and IS110 families were inserted around the regions whose functions were As/Hg/Cu/Co/Zn/Cd translocation and sulfur oxidation, implying that ISs could improve the adaptive capacities of Acidithiobacillus to the extremely acidic environment by enhancing their resistance to heavy metals and utilization of sulfur.

CONCLUSIONS

This study provided the genomic evidence for the contribution of IS to evolution and adaptation of Acidithiobacillus, opening novel sights into the genome plasticity of those acidophiles.

Differential Afa/Dr Fimbriae Expression in the Multidrug-Resistant Escherichia coli ST131 Clone

Many antibiotic resistant uropathogenic Escherichia coli (UPEC) strains belong to clones defined by their multilocus sequence type (ST), with ST131 being the most dominant. Although we have a good understanding of resistance development to fluoroquinolones and third-generation cephalosporins by ST131, our understanding of the virulence repertoire that has contributed to its global dissemination is limited. Here we show that the genes encoding Afa/Dr fimbriae, a group of adhesins strongly associated with UPEC that cause gestational pyelonephritis and recurrent cystitis, are found in approximately one third of all ST131 strains. Sequence comparison of the AfaE adhesin protein revealed a unique allelic variant carried by 82.9% of afa-positive ST131 strains. We identify the afa regulatory region as a hotspot for the integration of insertion sequence (IS) elements, all but one of which alter afa transcription. Close investigation demonstrated that the integration of an IS1 element in the afa regulatory region leads to increased expression of Afa/Dr fimbriae, promoting enhanced adhesion to kidney epithelial cells and suggesting a mechanism for altered virulence. Finally, we provide evidence for a more widespread impact of IS1 on ST131 genome evolution, suggesting that IS dynamics contribute to strain level microevolution that impacts ST131 fitness. IMPORTANCE E. coli ST131 is the most common antibiotic resistant UPEC clone associated with human urinary tract and bloodstream infections. Understanding the features of ST131 that have driven its global dissemination remains a critical priority if we are to counter its increasing antibiotic resistance. Here, we utilized a large collection of ST131 isolates to investigate the prevalence, regulation, and function of Afa/Dr fimbriae, a well-characterized UPEC colonization and virulence factor. We show that the afa genes are found frequently in ST131 and demonstrate how the integration of IS elements in the afa regulatory region modulates Afa expression, presenting an example of altered virulence capacity. We also exploit a curated set of ST131 genomes to map the integration of the antibiotic resistance-associated IS1 element in the ST131 pangenome, providing evidence for its widespread impact on ST131 genome evolution.

Impact of insertion sequences on convergent evolution of Shigella species

Shigella species are specialised lineages of Escherichia coli that have converged to become human-adapted and cause dysentery by invading human gut epithelial cells. Most studies of Shigella evolution have been restricted to comparisons of single representatives of each species; and population genomic studies of individual Shigella species have focused on genomic variation caused by single nucleotide variants and ignored the contribution of insertion sequences (IS) which are highly prevalent in Shigella genomes. Here, we investigate the distribution and evolutionary dynamics of IS within populations of Shigella dysenteriae Sd1, Shigella sonnei and Shigella flexneri. We find that five IS (IS1, IS2, IS4, IS600 and IS911) have undergone expansion in all Shigella species, creating substantial strain-to-strain variation within each population and contributing to convergent patterns of functional gene loss within and between species. We find that IS expansion and genome degradation are most advanced in S. dysenteriae and least advanced in S. sonnei; and using genome-scale models of metabolism we show that Shigella species display convergent loss of core E. coli metabolic capabilities, with S. sonnei and S. flexneri following a similar trajectory of metabolic streamlining to that of S. dysenteriae. This study highlights the importance of IS to the evolution of Shigella and provides a framework for the investigation of IS dynamics and metabolic reduction in other bacterial species.

Horizontal Gene Transfer and Endophytes: An Implication for the Acquisition of Novel Traits

Horizontal gene transfer (HGT), an important evolutionary mechanism observed in prokaryotes, is the transmission of genetic material across phylogenetically distant species. In recent years, the availability of complete genomes has facilitated the comprehensive analysis of HGT and highlighted its emerging role in the adaptation and evolution of eukaryotes. Endophytes represent an ecologically favored association, which highlights its beneficial attributes to the environment, in agriculture and in healthcare. The HGT phenomenon in endophytes, which features an important biological mechanism for their evolutionary adaptation within the host plant and simultaneously confers "novel traits" to the associated microbes, is not yet completely understood. With a focus on the emerging implications of HGT events in the evolution of biological species, the present review discusses the occurrence of HGT in endophytes and its socio-economic importance in the current perspective. To our knowledge, this review is the first report that provides a comprehensive insight into the impact of HGT in the adaptation and evolution of endophytes.

Beneath the surface: Evolution of methane activity in the bacterial multicomponent monooxygenases

The bacterial multicomponent monooxygenase (BMM) family has evolved to oxidise a wide array of hydrocarbon substrates of importance to environmental emissions and biotechnology: foremost amongst these is methane, which requires among the most powerful oxidant in biology to activate. To understand how the BMM evolved methane oxidation activity, we investigated the changes in the enzyme family at different levels: operonic, phylogenetic analysis of the catalytic hydroxylase, subunit or folding factor presence, and sequence-function analysis across the entirety of the BMM phylogeny. Our results show that the BMM evolution of new activities was enabled by incremental increases in oxidative power of the active site, and these occur in multiple branches of the hydroxylase phylogenetic tree. While the hydroxylase primary sequence changes that resulted in increased oxidative power of the enzyme appear to be minor, the principle evolutionary advances enabling methane activity occurred in the other components of the BMM complex and in the recruitment of stability proteins. We propose that enzyme assembly and stabilization factors have independently-evolved multiple times in the BMM family to support enzymes that oxidise increasingly difficult substrates. Herein, we show an important example of evolution of catalytic function where modifications to the active site and substrate accessibility, which are the usual focus of enzyme evolution, are overshadowed by broader scale changes to structural stabilization and non-catalytic unit development. Retracing macroscale changes during enzyme evolution, as demonstrated here, should find ready application to other enzyme systems and in protein design.

Barcode index number, taxonomic rank and modes of speciation: examples from fish

Species delimitation by DNA sequence data or DNA barcoding is successful, as confirmed by the vast BOLD data base. However, the theory that would explain this fact has not been developed yet. An approach based on Barcoding Index Number (BIN), suggested in the assignment, allows delimiting of taxa of three ranks (species, genera, and families) and statistical validation with a high precision of delimiting (over 80%), as well as shows for majority of Co-1-based single gene trees good correspondence between their topology and conventional taxa content for analyzed fish species (R² ≈ 0.84-0.98). Knowledge of deviations from these data can help to find out new taxa and improve biodiversity description. It is concluded that delimiting is successful for bulk of cases because the geographic mode of speciation prevails in nature. It takes a long time for new taxa to form in isolation, which allows accumulation of random mutations and many different nucleotide substitutions between them that can be detected by molecular markers and give unique DNA barcodes. The use of BIN approach, described here, can aid greatly in making this important question clearer especially under wider examination of other organisms.

Massive horizontal transfer of transposable elements in insects

Horizontal transfer (HT) of genetic material is central to the architecture and evolution of prokaryote genomes. Within eukaryotes, the majority of HTs reported so far are transfers of transposable elements (TEs). These reports essentially come from studies focusing on specific lineages or types of TEs. Because of the lack of large-scale survey, the amount and impact of HT of TEs (HTT) in eukaryote evolution, as well as the trends and factors shaping these transfers, are poorly known. Here, we report a comprehensive analysis of HTT in 195 insect genomes, representing 123 genera and 13 of the 28 insect orders. We found that these insects were involved in at least 2,248 HTT events that essentially occurred during the last 10 My. We show that DNA transposons transfer horizontally more often than retrotransposons, and unveil phylogenetic relatedness and geographical proximity as major factors facilitating HTT in insects. Even though our study is restricted to a small fraction of insect biodiversity and to a recent evolutionary timeframe, the TEs we found to be horizontally transferred generated up to 24% (2.08% on average) of all nucleotides of insect genomes. Together, our results establish HTT as a major force shaping insect genome evolution.

Four genes essential for recombination define GInts, a new type of mobile genomic island widespread in bacteria

Integrases are a family of tyrosine recombinases that are highly abundant in bacterial genomes, actively disseminating adaptive characters such as pathogenicity determinants and antibiotics resistance. Using comparative genomics and functional assays, we identified a novel type of mobile genetic element, the GInt, in many diverse bacterial groups but not in archaea. Integrated as genomic islands, GInts show a tripartite structure consisting of the ginABCD operon, a cargo DNA region from 2.5 to at least 70 kb, and a short AT-rich 3' end. The gin operon is characteristic of GInts and codes for three putative integrases and a small putative helix-loop-helix protein, all of which are essential for integration and excision of the element. Genes in the cargo DNA are acquired mostly from phylogenetically related bacteria and often code for traits that might increase fitness, such as resistance to antimicrobials or virulence. GInts also tend to capture clusters of genes involved in complex processes, such as the biosynthesis of phaseolotoxin by Pseudomonas syringae. GInts integrate site-specifically, generating two flanking direct imperfect repeats, and excise forming circular molecules. The excision process generates sequence variants at the element attachment site, which can increase frequency of integration and drive target specificity.

A Comprehensive Study of Cyanobacterial Morphological and Ecological Evolutionary Dynamics through Deep Geologic Time

Cyanobacteria have exerted a profound influence on the progressive oxygenation of Earth. As a complementary approach to examining the geologic record—phylogenomic and trait evolutionary analyses of extant species can lead to new insights. We constructed new phylogenomic trees and analyzed phenotypic trait data using novel phylogenetic comparative methods. We elucidated the dynamics of trait evolution in Cyanobacteria over billion-year timescales, and provide evidence that major geologic events in early Earth’s history have shaped—and been shaped by—evolution in Cyanobacteria. We identify a robust core cyanobacterial phylogeny and a smaller set of taxa that exhibit long-branch attraction artifacts. We estimated the age of nodes and reconstruct the ancestral character states of 43 phenotypic characters. We find high levels of phylogenetic signal for nearly all traits, indicating the phylogeny carries substantial predictive power. The earliest cyanobacterial lineages likely lived in freshwater habitats, had small cell diameters, were benthic or sessile, and possibly epilithic/endolithic with a sheath. We jointly analyzed a subset of 25 binary traits to determine whether rates of trait evolution have shifted over time in conjunction with major geologic events. Phylogenetic comparative analysis reveal an overriding signal of decreasing rates of trait evolution through time. Furthermore, the data suggest two major rate shifts in trait evolution associated with bursts of evolutionary innovation. The first rate shift occurs in the aftermath of the Great Oxidation Event and “Snowball Earth” glaciations and is associated with decrease in the evolutionary rates around 1.8–1.6 Ga. This rate shift seems to indicate the end of a major diversification of cyanobacterial phenotypes–particularly related to traits associated with filamentous morphology, heterocysts and motility in freshwater ecosystems. Another burst appears around the time of the Neoproterozoic Oxidation Event in the Neoproterozoic, and is associated with the acquisition of traits involved in planktonic growth in marine habitats. Our results demonstrate how uniting genomic and phenotypic datasets in extant bacterial species can shed light on billion-year old events in Earth’s history.

The genomic signatures of Shigella evolution, adaptation and geographical spread

Shigella spp. are some of the key pathogens responsible for the global burden of diarrhoeal disease. These facultative intracellular bacteria belong to the family Enterobacteriaceae, together with other intestinal pathogens, such as Escherichia coli and Salmonella spp. The genus Shigella comprises four different species, each consisting of several serogroups, all of which show phenotypic similarity, including invasive pathogenicity. DNA sequencing suggests that this similarity results from the convergent evolution of different Shigella spp. founders. Here, we review the evolutionary relationships between Shigella spp. and E . coli, and we highlight how the genomic plasticity of these bacteria and their acquisition of a distinctive virulence plasmid have enabled the development of such highly specialized pathogens. Furthermore, we discuss the insights that genotyping and whole-genome sequencing have provided into the phylogenetics and intercontinental spread of Shigella spp.

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

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

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.

Phylogeny and Taxonomy of Archaea: A Comparison of the Whole-Genome-Based CVTree Approach with 16S rRNA Sequence Analysis

A tripartite comparison of Archaea phylogeny and taxonomy at and above the rank order is reported: (1) the whole-genome-based and alignment-free CVTree using 179 genomes; (2) the 16S rRNA analysis exemplified by the All-Species Living Tree with 366 archaeal sequences; and (3) the Second Edition of Bergey's Manual of Systematic Bacteriology complemented by some current literature. A high degree of agreement is reached at these ranks. From the newly proposed archaeal phyla, Korarchaeota, Thaumarchaeota, Nanoarchaeota and Aigarchaeota, to the recent suggestion to divide the class Halobacteria into three orders, all gain substantial support from CVTree. In addition, the CVTree helped to determine the taxonomic position of some newly sequenced genomes without proper lineage information. A few discrepancies between the CVTree and the 16S rRNA approaches call for further investigation.

Spaces of the possible: universal Darwinism and the wall between technological and biological innovation

Innovations in biological evolution and in technology have many common features. Some of them involve similar processes, such as trial and error and horizontal information transfer. Others describe analogous outcomes such as multiple independent origins of similar innovations. Yet others display similar temporal patterns such as episodic bursts of change separated by periods of stasis. We review nine such commonalities, and propose that the mathematical concept of a space of innovations, discoveries or designs can help explain them. This concept can also help demolish a persistent conceptual wall between technological and biological innovation.

Cargo capacity of phages and plasmids and other factors influencing horizontal transfers of prokaryote transposable elements

Horizontal transfer of transposable elements (TEs) plays a key role in prokaryote genome evolution. Most TEs do not encode the enzymatic machinery allowing them to transfer between host cells and it is widely assumed in the literature that horizontal transfer of prokaryote TEs is mediated by other mobile genetic elements such as phages and plasmids. In a recent study, we have shown that phages are less tolerant to insertion sequences (IS, the most frequent class of prokaryote TEs) and therefore have a lower cargo capacity than plasmids. Consequently, while our analysis confirmed the crucial role of plasmids as efficient vehicles of IS horizontal transfer, we concluded that phages are unlikely to efficiently shuttle IS elements between prokaryotes. Here, we discuss whether or not the distribution pattern observed for IS elements in phages and plasmids also holds for other TEs, such as transposons and mobile introns. We also further explore various factors that may impact the relative capacity of phages and plasmids to mediate TE horizontal transfer among prokaryotes.

The influence of the accessory genome on bacterial pathogen evolution

Bacterial pathogens exhibit significant variation in their genomic content of virulence factors. This reflects the abundance of strategies pathogens evolved to infect host organisms by suppressing host immunity. Molecular arms-races have been a strong driving force for the evolution of pathogenicity, with pathogens often encoding overlapping or redundant functions, such as type III protein secretion effectors and hosts encoding ever more sophisticated immune systems. The pathogens' frequent exposure to other microbes, either in their host or in the environment, provides opportunities for the acquisition or interchange of mobile genetic elements. These DNA elements accessorize the core genome and can play major roles in shaping genome structure and altering the complement of virulence factors. Here, we review the different mobile genetic elements focusing on the more recent discoveries and highlighting their role in shaping bacterial pathogen evolution.

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).

Lateral transfers of insertion sequences between Wolbachia, Cardinium and Rickettsia bacterial endosymbionts

Various bacteria live exclusively within arthropod cells and collectively act as an important driver of arthropod evolutionary ecology. Whereas rampant intra-generic DNA transfers were recently shown to have a pivotal role in the evolution of the most common of these endosymbionts, Wolbachia, the present study show that inter-generic DNA transfers also commonly take place, constituting a potent source of rapid genomic change. Bioinformatic, molecular and phylogenetic data provide evidence that a selfish genetic element, the insertion sequence ISRpe1, is widespread in the Wolbachia, Cardinium and Rickettsia endosymbionts and experiences recent (and likely ongoing) transfers over long evolutionary distances. Although many ISRpe1 copies were clearly expanding and leading to rapid endosymbiont diversification, degraded copies are also frequently found, constituting an unusual genomic fossil record suggestive of ancient ISRpe1 expansions. Overall, the present data highlight how ecological connections within the arthropod intracellular environment facilitate lateral DNA transfers between distantly related bacterial lineages.

Co-phylog: an assembly-free phylogenomic approach for closely related organisms

With the advent of high-throughput sequencing technologies, the rapid generation and accumulation of large amounts of sequencing data pose an insurmountable demand for efficient algorithms for constructing whole-genome phylogenies. The existing phylogenomic methods all use assembled sequences, which are often not available owing to the difficulty of assembling short-reads; this obstructs phylogenetic investigations on species without a reference genome. In this report, we present co-phylog, an assembly-free phylogenomic approach that creates a ‘micro-alignment’ at each ‘object’ in the sequence using the ‘context’ of the object and calculates pairwise distances before reconstructing the phylogenetic tree based on those distances. We explored the parameters’ usages and the optimal working range of co-phylog, assessed co-phylog using the simulated next-generation sequencing (NGS) data and the real NGS raw data. We also compared co-phylog method with traditional alignment and alignment-free methods and illustrated the advantages and limitations of co-phylog method. In conclusion, we demonstrated that co-phylog is efficient algorithm and that it delivers high resolution and accurate phylogenies using whole-genome unassembled sequencing data, especially in the case of closely related organisms, thereby significantly alleviating the computational burden in the genomic era.

Multiple inter-kingdom horizontal gene transfers in the evolution of the phosphoenolpyruvate carboxylase gene family

Pepcase is a gene encoding phosphoenolpyruvate carboxylase that exists in bacteria, archaea and plants,playing an important role in plant metabolism and development. Most plants have two or more pepcase genes belonging to two gene sub-families, while only one gene exists in other organisms. Previous research categorized one plant pepcase gene as plant-type pepcase (PTPC) while the other as bacteria-type pepcase (BTPC) because of its similarity with the pepcase gene found in bacteria. Phylogenetic reconstruction showed that PTPC is the ancestral lineage of plant pepcase, and that all bacteria, protistpepcase and BTPC in plants are derived from a lineage of pepcase closely related with PTPC in algae. However, their phylogeny contradicts the species tree and traditional chronology of organism evolution. Because the diversification of bacteria occurred much earlier than the origin of plants, presumably all bacterialpepcase derived from the ancestral PTPC of algal plants after divergingfrom the ancestor of vascular plant PTPC. To solve this contradiction, we reconstructed the phylogeny of pepcase gene family. Our result showed that both PTPC and BTPC are derived from an ancestral lineage of gamma-proteobacteriapepcases, possibly via an ancient inter-kingdom horizontal gene transfer (HGT) from bacteria to the eukaryotic common ancestor of plants, protists and cellular slime mold. Our phylogenetic analysis also found 48other pepcase genes originated from inter-kingdom HGTs. These results imply that inter-kingdom HGTs played important roles in the evolution of the pepcase gene family and furthermore that HGTsare a more frequent evolutionary event than previouslythought.

OASIS: an automated program for global investigation of bacterial and archaeal insertion sequences

Insertion sequences (ISs) are simple transposable elements present in most bacterial and archaeal genomes and play an important role in genomic evolution. The recent expansion of sequenced genomes offers the opportunity to study ISs comprehensively, but this requires efficient and accurate tools for IS annotation. We have developed an open-source program called OASIS, or Optimized Annotation System for Insertion Sequences, which automatically annotates ISs within sequenced genomes. OASIS annotations of 1737 bacterial and archaeal genomes offered an unprecedented opportunity to examine IS evolution. At a broad scale, we found that most IS families are quite widespread; however, they are not present randomly across taxa. This may indicate differential loss, barriers to exchange and/or insufficient time to equilibrate across clades. The number of ISs increases with genome length, but there is both tremendous variation and no increase in IS density for genomes >2 Mb. At the finer scale of recently diverged genomes, the proportion of shared IS content falls sharply, suggesting loss and/or emergence of barriers to successful cross-infection occurs rapidly. Surprisingly, even after controlling for 16S rRNA sequence divergence, the same ISs were more likely to be shared between genomes labeled as the same species rather than as different species.

Structural and functional consequences of mutating a proteobacteria-specific surface residue in the catalytic domain of Escherichia coli GluRS

Nucleotides whose mutations seriously affect glutamylation efficiency are experimentally known for Escherichia coli tRNA(Glu). However, not much is known about functional hotspots on the complementary enzyme, glutamyl-tRNA synthetase (GluRS). From structural and functional studies on an Arg266Leu mutant of E. coli GluRS, we demonstrate that Arg266 is essential for efficient glutamylation of tRNA(Glu). Consistent with this result, we found that Arg266 is a conserved signature of proteobacterial GluRS. In contrast, most non-proteobacterial GluRS contain Leu, and never Arg, at this position. Our results imply a unique strategy of glutamylation of tRNA(Glu) in proteobacteria under phylum-specific evolutionary compulsions.

Acinetobacter insertion sequence ISAba11 belongs to a novel family that encodes transposases with a signature HHEK motif

Experimental and in silico PCR analysis targeting ISAba11 and TnAbaR islands in 196 epidemiologically unrelated Acinetobacter strains representative of ≥19 species were performed. The first two Acinetobacter baumannii ISAba11 elements identified had been found to map to the same site on TnAbaR transposons. However, no further evidence of physical linkage between the two elements was demonstrated. Indeed, examination of 25 definite or putative insertion sites suggested limited sequence specificity. Importantly, an aacC1-tagged version of ISAba11 was shown to actively transpose in A. baumannii. Similarity searches identified nine iso-ISAba11 elements in Acinetobacter and one in Enhydrobacter and single representatives of four distant homologs in bacteria belonging to the phyla "Cyanobacteria" and Proteobacteria. Phylogenetic, sequence, and structural analyses of ISAba11 and/or its associated transposase (Tnp(ISAba11)) suggested that these elements be assigned to a new family. All five homologs encode transposases with a shared extended signature comprising 16 invariant residues within the N2, N3, and C1 regions, four of which constituted the cardinal ISAba11 family HHEK motif that is substituted for the YREK DNA binding motif conserved in the IS4 family. Additionally, ISAba11 family members were associated with either no flanking direct repeat (DR) or an ISAba11-typical 5-bp DR and possessed variable-length terminal inverted repeats that exhibited extensive intrafamily sequence identity. Given the limited pairwise identity among Tnp(ISAba11) homologs and the observed restricted distribution of ISAba11, we propose that substantial gaps persist in the evolutionary record of ISAba11 and that this element represents a recent though potentially highly significant entrant into the A. baumannii gene pool.

A bacterial genome in transition--an exceptional enrichment of IS elements but lack of evidence for recent transposition in the symbiont Amoebophilus asiaticus

Background

Insertion sequence (IS) elements are important mediators of genome plasticity and are widespread among bacterial and archaeal genomes. The 1.88 Mbp genome of the obligate intracellular amoeba symbiont Amoebophilus asiaticus contains an unusually large number of transposase genes (n = 354; 23% of all genes).

Results

The transposase genes in the A. asiaticus genome can be assigned to 16 different IS elements termed ISCaa1 to ISCaa16, which are represented by 2 to 24 full-length copies, respectively. Despite this high IS element load, the A. asiaticus genome displays a GC skew pattern typical for most bacterial genomes, indicating that no major rearrangements have occurred recently. Additionally, the high sequence divergence of some IS elements, the high number of truncated IS element copies (n = 143), as well as the absence of direct repeats in most IS elements suggest that the IS elements of A. asiaticus are transpositionally inactive. Although we could show transcription of 13 IS elements, we did not find experimental evidence for transpositional activity, corroborating our results from sequence analyses. However, we detected contiguous transcripts between IS elements and their downstream genes at nine loci in the A. asiaticus genome, indicating that some IS elements influence the transcription of downstream genes, some of which might be important for host cell interaction.

Conclusions

Taken together, the IS elements in the A. asiaticus genome are currently in the process of degradation and largely represent reflections of the evolutionary past of A. asiaticus in which its genome was shaped by their activity.

Do phages efficiently shuttle transposable elements among prokaryotes?

The distribution, dynamics, and evolution of insertion sequences (IS), the most frequent class of prokaryotic transposable elements, are conditioned by their ability to horizontally transfer between cells. IS horizontal transfer (HT) requires shuttling by other mobile genetic elements. It is widely assumed in the literature that these vectors are phages and plasmids. By examining the relative abundance of IS in 454 plasmid and 446 phage genomes, we found that IS are very frequent in plasmids but, surprisingly, very rare in phages. Our results indicate that IS rarity in phages reflects very strong and efficient postinsertional purifying selection, mainly caused by a higher density of deleterious insertion sites in phages compared to plasmids. As they do not tolerate IS insertions, we conclude that phages may be rather poor vectors of IS HT in prokaryotes, in sharp contrast with the conventional view.

Mobile DNA can drive lineage extinction in prokaryotic populations

Natural selection ultimately acts on genes and other DNA sequences. Adaptations that are good for the gene can have adverse effects at higher levels of organization, including the individual or the population. Mobile genetic elements illustrate this principle well, because they can self-replicate within a genome at a cost to their host. As they are costly and can be transmitted horizontally, mobile elements can be seen as genomic parasites. It has been suggested that mobile elements may cause the extinction of their host populations. In organisms with very large populations, such as most bacteria, individual selection is highly effective in purging genomes of deleterious elements, suggesting that extinction is unlikely. Here we investigate the conditions under which mobile DNA can drive bacterial lineages to extinction. We use a range of epidemiological and ecological models to show that harmful mobile DNA can invade, and drive populations to extinction, provided their transmission rate is high and that mobile element-induced mortality is not too high. Population extinction becomes more likely when there are more elements in the population. Even if elements are costly, extinction can still occur because of the combined effect of horizontal gene transfer, a mortality induced by mobile elements. Our study highlights the potential of mobile DNA to be selected at the population level, as well as at the individual level.

Horizontal gene transfer in evolution: facts and challenges

The contribution of horizontal gene transfer to evolution has been controversial since it was suggested to be a force driving evolution in the microbial world. In this paper, I review the current standpoint on horizontal gene transfer in evolutionary thinking and discuss how important horizontal gene transfer is in evolution in the broad sense, and particularly in prokaryotic evolution. I review recent literature, asking, first, which processes are involved in the evolutionary success of transferred genes and, secondly, about the extent of horizontal gene transfer towards different evolutionary times. Moreover, I discuss the feasibility of reconstructing ancient phylogenetic relationships in the face of horizontal gene transfer. Finally, I discuss how horizontal gene transfer fits in the current neo-Darwinian evolutionary paradigm and conclude there is a need for a new evolutionary paradigm that includes horizontal gene transfer as well as other mechanisms in the explanation of evolution.