Insertion sequence diversity in archaea

Hijackers, hitchhikers, or co-drivers? The mysteries of mobilizable genetic elements

Mobile genetic elements shape microbial gene repertoires and populations. Recent results reveal that many, possibly most, microbial mobile genetic elements require helpers to transfer between genomes, which we refer to as Hitcher Genetic Elements (hitchers or HGEs). They may be a large fraction of pathogenicity and resistance genomic islands, whose mechanisms of transfer have remained enigmatic for decades. Together with their helper elements and their bacterial hosts, hitchers form tripartite networks of interactions that evolve rapidly within a parasitism-mutualism continuum. In this emerging view of microbial genomes as communities of mobile genetic elements many questions arise. Which elements are being moved, by whom, and how? How often are hitchers costly hyper-parasites or beneficial mutualists? What is the evolutionary origin of hitchers? Are there key advantages associated with hitchers' lifestyle that justify their unexpected abundance? And why are hitchers systematically smaller than their helpers? In this essay, we start answering these questions and point ways ahead for understanding the principles, origin, mechanisms, and impact of hitchers in bacterial ecology and evolution.

Molecular identification of polymorphic transposable elements in populations of the invasive ant Cardiocondyla obscurior

Transposable elements (TEs) are found in virtually every eukaryotic genome and are important for generating de novo genetic variation. However, outside of costly and time-consuming whole-genome sequencing approaches, the set of available methods to study TE polymorphisms in non-model species is very limited. The Transposon Display (TD) is a simple yet effective technique to characterize polymorphisms across samples by identifying amplified fragment length polymorphisms using primers targeting specific TE families. So far, this technique has almost exclusively been used in plants. Here, we present an optimized TD protocol for insect species with small genomes such as ants (ca. 200-600 Mb). We characterized TE polymorphisms between two distinct genetic lineages of the invasive ant Cardiocondyla obscurior, as well as between neighboring populations of the New World lineage. We found active LTR/Ty3 retrotransposons, that contributed to the genetic diversification of populations in this species.

Evolutionary mining and functional characterization of TnpB nucleases identify efficient miniature genome editors

As the evolutionary ancestor of Cas12 nuclease, the transposon (IS200/IS605)-encoded TnpB proteins act as compact RNA-guided DNA endonucleases. To explore their evolutionary diversity and potential as genome editors, we screened TnpBs from 64 annotated IS605 members and identified 25 active in Escherichia coli, of which three are active in human cells. Further characterization of these 25 TnpBs enables prediction of the transposon-associated motif (TAM) and the right-end element RNA (reRNA) directly from genomic sequences. We established a framework for annotating TnpB systems in prokaryotic genomes and applied it to identify 14 additional candidates. Among these, ISAam1 (369 amino acids (aa)) and ISYmu1 (382 aa) TnpBs demonstrated robust editing activity across dozens of genomic loci in human cells. Both RNA-guided genome editors demonstrated similar editing efficiency as SaCas9 (1,053 aa) while being substantially smaller. The enormous diversity of TnpBs holds potential for the discovery of additional valuable genome editors.

'All About' Extremophiles

Despite common perception, most of Earth is what is often referred to as an 'extreme environment.' Yet to the organisms that call these places home, it is simply that (home). They have adapted to thrive in these environments and, in the process, have evolved many unique adaptations at the molecular- and 'omic-level. Scientists' interest in these organisms has typically been in how they and their products can be harnessed for biotechnological applications and the environments where they are found, while the general public's veers more toward a fascination with their deviation from the 'norm'. However, these organisms have so much more to tell us about Life and the myriad ways there are to perform 'simple' biological processes.

Membrane lipid and expression responses of Saccharolobus islandicus REY15A to acid and cold stress

Archaea adjust the number of cyclopentane rings in their glycerol dibiphytanyl glycerol tetraether (GDGT) membrane lipids as a homeostatic response to environmental stressors such as temperature, pH, and energy availability shifts. However, archaeal expression patterns that correspond with changes in GDGT composition are less understood. Here we characterize the acid and cold stress responses of the thermoacidophilic crenarchaeon Saccharolobus islandicus REY15A using growth rates, core GDGT lipid profiles, transcriptomics and proteomics. We show that both stressors result in impaired growth, lower average GDGT cyclization, and differences in gene and protein expression. Transcription data revealed differential expression of the GDGT ring synthase grsB in response to both acid stress and cold stress. Although the GDGT ring synthase encoded by grsB forms highly cyclized GDGTs with ≥5 ring moieties, S. islandicus grsB upregulation under acidic pH conditions did not correspond with increased abundances of highly cyclized GDGTs. Our observations highlight the inability to predict GDGT changes from transcription data alone. Broader analysis of transcriptomic data revealed that S. islandicus differentially expresses many of the same transcripts in response to both acid and cold stress. These included upregulation of several biosynthetic pathways and downregulation of oxidative phosphorylation and motility. Transcript responses specific to either of the two stressors tested here included upregulation of genes related to proton pumping and molecular turnover in acid stress conditions and upregulation of transposases in cold stress conditions. Overall, our study provides a comprehensive understanding of the GDGT modifications and differential expression characteristic of the acid stress and cold stress responses in S. islandicus.

G. de Lomana A, Gomes-Filho JV, Vêncio RZN, Moritz RL, Koide T, Baliga NS. A Genome-Scale Atlas Reveals Complex Interplay of Transcription and Translation in an Archaeon

The scale of post-transcriptional regulation and the implications of its interplay with other forms of regulation in environmental acclimation are underexplored for organisms of the domain Archaea. Here, we have investigated the scale of post-transcriptional regulation in the extremely halophilic archaeon Halobacterium salinarum NRC-1 by integrating the transcriptome-wide locations of transcript processing sites (TPSs) and SmAP1 binding, the genome-wide locations of antisense RNAs (asRNAs), and the consequences of RNase_2099C knockout on the differential expression of all genes. This integrated analysis has discovered that 54% of all protein-coding genes in the genome of this haloarchaeon are likely targeted by multiple mechanisms for putative post-transcriptional processing and regulation, with about 20% of genes likely being regulated by combinatorial schemes involving SmAP1, asRNAs, and RNase_2099C. Comparative analysis of mRNA levels (transcriptome sequencing [RNA-Seq]) and protein levels (sequential window acquisition of all theoretical fragment ion spectra mass spectrometry [SWATH-MS]) for 2,579 genes over four phases of batch culture growth in complex medium generated additional evidence for the conditional post-transcriptional regulation of 7% of all protein-coding genes. We demonstrate that post-transcriptional regulation may act to fine-tune specialized and rapid acclimation to stressful environments, e.g., as a switch to turn on gas vesicle biogenesis to promote vertical relocation under anoxic conditions and modulate the frequency of transposition by insertion sequence (IS) elements of the IS200/IS605, IS4, and ISH3 families. Findings from this study are provided as an atlas in a public Web resource (https://halodata.systemsbiology.net). IMPORTANCE While the transcriptional regulation landscape of archaea has been extensively investigated, we currently have limited knowledge about post-transcriptional regulation and its driving mechanisms in this domain of life. In this study, we collected and integrated omics data from multiple sources and technologies to infer post-transcriptionally regulated genes and the putative mechanisms modulating their expression at the protein level in Halobacterium salinarum NRC-1. The results suggest that post-transcriptional regulation may drive environmental acclimation by regulating hallmark biological processes. To foster discoveries by other research groups interested in the topic, we extended our integrated data to the public in the form of an interactive atlas (https://halodata.systemsbiology.net).

Hydrolytic enzyme screening and carotenoid production evaluation of halophilic archaea isolated from highly heavy metal-enriched solar saltern sediments

This paper aimed to screen the enzymatic activities and evaluate the carotenoid production level of twenty-two halophilic archaea isolated from Sfax solar saltern sediments. The molecular identification performed by sequencing the 16S rRNA genes showed that all strains have a high similarity degree (99.7-100%) with Halobacterium salinarum NRC-1. The strains were screened for the presence of eight hydrolase activities using agar plate-based assays. The most detected enzyme was gelatinase (77.27% of total strains), followed by protease (63.63%) and amylase activities (50%). The carotenoid production yields of the strains ranged between 2.027 and 14.880 mg/l. The UV-Visible spectroscopy of pigments revealed that it was a bacterioruberin type. When evaluated and compared to the standard β-carotene, the antioxidant capacities of these pigments showed a scavenging activity of more than 75% at a concentration of 5 μg/ml for three strains (AS16, AS17, and AS18). Then a sequence of one-step optimization processes was performed, using the one-factor-at-a-time approach, to define the optimum conditions for growth and carotenoid production of the highest carotenoid producing strain (AS17). Different environmental factors and nutritional conditions were tested. Variations in these factors were found to deeply influence growth and carotenoid production. A maximum carotenoid production (16.490 mg/l), higher than that of the control (14.880 mg/l), was observed at 37 °C, pH 7, 250 g/l of salinity, with 80% air phase in the flask at 110 rpm, in presence of light and in culture media containing (g/l) 10, yeast extract; 7.5, casamino acid; 20, MgSO4; 4, KCl; and 3, trisodium citrate.

The origin of RNA interference: Adaptive or neutral evolution?

The origin of RNA interference (RNAi) is usually explained by a defense-based hypothesis, in which RNAi evolved as a defense against transposable elements (TEs) and RNA viruses and was already present in the last eukaryotic common ancestor (LECA). However, since RNA antisense regulation and double-stranded RNAs (dsRNAs) are ancient and widespread phenomena, the origin of defensive RNAi should have occurred in parallel with its regulative functions to avoid imbalances in gene regulation. Thus, we propose a neutral evolutionary hypothesis for the origin of RNAi in which qualitative system drift from a prokaryotic antisense RNA gene regulation mechanism leads to the formation of RNAi through constructive neutral evolution (CNE). We argue that RNAi was already present in the ancestor of LECA before the need for a new defense system arose and that its presence helped to shape eukaryotic genomic architecture and stability.

Where does RNA interference come from? This Essay describes a new step-by-step evolutionary model of how RNA interference might have originated in early eukaryotes through neutral events from the molecular machinery present in prokaryotes.

What makes a megaplasmid?

Naturally occurring plasmids come in different sizes. The smallest are less than a kilobase of DNA, while the largest can be over three orders of magnitude larger. Historically, research has tended to focus on smaller plasmids that are usually easier to isolate, manipulate and sequence, but with improved genome assemblies made possible by long-read sequencing, there is increased appreciation that very large plasmids—known as megaplasmids—are widespread, diverse, complex, and often encode key traits in the biology of their host microorganisms. Why are megaplasmids so big? What other features come with large plasmid size that could affect bacterial ecology and evolution? Are megaplasmids 'just' big plasmids, or do they have distinct characteristics? In this perspective, we reflect on the distribution, diversity, biology, and gene content of megaplasmids, providing an overview to these large, yet often overlooked, mobile genetic elements.

This article is part of the theme issue ‘The secret lives of microbial mobile genetic elements’.

Assessment of the plasmidome of an extremophilic microbial community from the Diamante Lake, Argentina

Diamante Lake located at 4589 m.a.s.l. in the Andean Puna constitutes an extreme environment. It is exposed to multiple extreme conditions such as an unusually high concentration of arsenic (over 300 mg L^-1) and low oxygen pressure. Microorganisms thriving in the lake display specific genotypes that facilitate survival, which include at least a multitude of plasmid-encoded resistance traits. Hence, the genetic information provided by the plasmids essentially contributes to understand adaptation to different stressors. Though plasmids from cultivable organisms have already been analyzed to the sequence level, the impact of the entire plasmid-borne genetic information on such microbial ecosystem is not known. This study aims at assessing the plasmidome from Diamante Lake, which facilitates the identification of potential hosts and prediction of gene functions as well as the ecological impact of mobile genetic elements. The deep-sequencing analysis revealed a large fraction of previously unknown DNA sequences of which the majority encoded putative proteins of unknown function. Remarkably, functions related to the oxidative stress response, DNA repair, as well as arsenic- and antibiotic resistances were annotated. Additionally, all necessary capacities related to plasmid replication, mobilization and maintenance were detected. Sequences characteristic for megaplasmids and other already known plasmid-associated genes were identified as well. The study highlights the potential of the deep-sequencing approach specifically targeting plasmid populations as it allows to evaluate the ecological impact of plasmids from (cultivable and non-cultivable) microorganisms, thereby contributing to the understanding of the distribution of resistance factors within an extremophilic microbial community.

The lanthipeptide biosynthetic clusters of the domain Archaea

Research on Archaea's secondary metabolites is still lagging behind that of Bacteria and Eukarya. Our goal was to contribute to this knowledge gap by analyzing the lanthipeptide's clusters in Archaea. As previously proposed, Archaea encodes only class II synthetases (LanMs), which we found to be confined to the class Halobacteria (also known as haloarchaea). In total, we analyzed the phylogeny and the domains of 42 LanMs. Four types were identified, and the majority of them belong to the CCG group due to their cyclization domain, which includes LanMs of Cyanobacteria. Putative cognate peptides were predicted for most of LanMs and are a very diverse group of molecules that share a Kx(Y/F)(D/E)xx(F/Y) motif in their leader peptides. According to their homology, some of them were categorized into subfamilies, including Halolancins, Haladacins, Haloferaxcins and Halobiforcins. Many LanM genes were associated with mobile genetic elements, and their vicinities mainly encode ABC and MFS transporters, tailoring enzymes and uncharacterized proteins. Our results suggest that the biosynthesis of lanthipeptides in haloarchaea can entail distinct enzymology that must lead to the production of peptides with novel structures and unpredicted biological and ecological roles. Finally, an Haloferax mediterranei knockout, lacking its three lanM genes, was generated, and it was concluded that its antimicrobial activity is not primarily related to the production of lanthipeptides.

Archaeal tyrosine recombinases

The integration of mobile genetic elements into their host chromosome influences the immediate fate of cellular organisms and gradually shapes their evolution. Site-specific recombinases catalyzing this integration have been extensively characterized both in bacteria and eukarya. More recently, a number of reports provided the in-depth characterization of archaeal tyrosine recombinases and highlighted new particular features not observed in the other two domains. In addition to being active in extreme environments, archaeal integrases catalyze reactions beyond site-specific recombination. Some of these integrases can catalyze low-sequence specificity recombination reactions with the same outcome as homologous recombination events generating deep rearrangements of their host genome. A large proportion of archaeal integrases are termed suicidal due to the presence of a specific recombination target within their own gene. The paradoxical maintenance of integrases that disrupt their gene upon integration implies novel mechanisms for their evolution. In this review, we assess the diversity of the archaeal tyrosine recombinases using a phylogenomic analysis based on an exhaustive similarity network. We outline the biochemical, ecological and evolutionary properties of these enzymes in the context of the families we identified and emphasize similarities and differences between archaeal recombinases and their bacterial and eukaryal counterparts.

Integration of Complete Plasmids Containing Bont Genes into Chromosomes of Clostridium parabotulinum, Clostridium sporogenes, and Clostridium argentinense

At least 40 toxin subtypes of botulinum neurotoxins (BoNTs), a heterogenous group of bacterial proteins, are produced by seven different clostridial species. A key factor that drives the diversity of neurotoxigenic clostridia is the association of bont gene clusters with various genomic locations including plasmids, phages and the chromosome. Analysis of Clostridium sporogenes BoNT/B1 strain CDC 1632, C. argentinense BoNT/G strain CDC 2741, and Clostridium parabotulinum BoNT/B1 strain DFPST0006 genomes revealed bont gene clusters within plasmid-like sequences within the chromosome or nested in large contigs, with no evidence of extrachromosomal elements. A nucleotide sequence (255,474 bp) identified in CDC 1632 shared 99.5% identity (88% coverage) with bont/B1-containing plasmid pNPD7 of C. sporogenes CDC 67071; CDC 2741 contig AYSO01000020 (1.1 MB) contained a ~140 kb region which shared 99.99% identity (100% coverage) with plasmid pRSJ17_1 of C. argentinense BoNT/G strain 89G; and DFPST0006 contig JACBDK0100002 (573 kb) contained a region that shared 100% identity (99%) coverage with the bont/B1-containing plasmid pCLD of C. parabotulinum Okra. This is the first report of full-length plasmid DNA-carrying complete neurotoxin gene clusters integrated in three distinct neurotoxigenic species: C. parabotulinum, C. sporogenes and C. argentinense.

Comparative insights into genome signatures of ferric iron oxide- and anode-stimulated Desulfuromonas spp. strains

BACKGROUND

Halotolerant Fe (III) oxide reducers affiliated in the family Desulfuromonadaceae are ubiquitous and drive the carbon, nitrogen, sulfur and metal cycles in marine subsurface sediment. Due to their possible application in bioremediation and bioelectrochemical engineering, some of phylogenetically close Desulfuromonas spp. strains have been isolated through enrichment with crystalline Fe (III) oxide and anode. The strains isolated using electron acceptors with distinct redox potentials may have different abilities, for instance, of extracellular electron transport, surface recognition and colonization. The objective of this study was to identify the different genomic signatures between the crystalline Fe (III) oxide-stimulated strain AOP6 and the anode-stimulated strains WTL and DDH964 by comparative genome analysis.

RESULTS

The AOP6 genome possessed the flagellar biosynthesis gene cluster, as well as diverse and abundant genes involved in chemotaxis sensory systems and c-type cytochromes capable of reduction of electron acceptors with low redox potentials. The WTL and DDH964 genomes lacked the flagellar biosynthesis cluster and exhibited a massive expansion of transposable gene elements that might mediate genome rearrangement, while they were deficient in some of the chemotaxis and cytochrome genes and included the genes for oxygen resistance.

CONCLUSIONS

Our results revealed the genomic signatures distinctive for the ferric iron oxide- and anode-stimulated Desulfuromonas spp. strains. These findings highlighted the different metabolic abilities, such as extracellular electron transfer and environmental stress resistance, of these phylogenetically close bacterial strains, casting light on genome evolution of the subsurface Fe (III) oxide reducers.

The IS6 family, a clinically important group of insertion sequences including IS26

The IS6 family of bacterial and archaeal insertion sequences, first identified in the early 1980s, has proved to be instrumental in the rearrangement and spread of multiple antibiotic resistance. Two IS, IS26 (found in many enterobacterial clinical isolates as components of both chromosome and plasmids) and IS257 (identified in the plasmids and chromosomes of gram-positive bacteria), have received particular attention for their clinical impact. Although few biochemical data are available concerning the transposition mechanism of these elements, genetic studies have provided some interesting observations suggesting that members of the family might transpose using an unexpected mechanism. In this review, we present an overview of the family, the distribution and phylogenetic relationships of its members, their impact on their host genomes and analyse available data concerning the particular transposition pathways they may use. We also provide a mechanistic model that explains the recent observations on one of the IS6 family transposition pathways: targeted cointegrate formation between replicons.

Utilization of Phenol as Carbon Source by the Thermoacidophilic Archaeon Saccharolobus solfataricus P2 Is Limited by Oxygen Supply and the Cellular Stress Response

Present in many industrial effluents and as common degradation product of organic matter, phenol is a widespread compound which may cause serious environmental problems, due to its toxicity to animals and humans. Degradation of phenol from the environment by mesophilic bacteria has been studied extensively over the past decades, but only little is known about phenol biodegradation at high temperatures or low pH. In this work we studied phenol degradation in the thermoacidophilic archaeon Saccharolobus solfataricus P2 (basonym: Sulfolobus solfataricus) under extreme conditions (80°C, pH 3.5). We combined metabolomics and transcriptomics together with metabolic modeling to elucidate the organism’s response to growth with phenol as sole carbon source. Although S. solfataricus is able to utilize phenol for biomass production, the carbon source induces profound stress reactions, including genome rearrangement as well as a strong intracellular accumulation of polyamines. Furthermore, computational modeling revealed a 40% higher oxygen demand for substrate oxidation, compared to growth on glucose. However, only 16.5% of oxygen is used for oxidation of phenol to catechol, resulting in a less efficient integration of carbon into the biomass. Finally, our data underlines the importance of the phenol meta-degradation pathway in S. solfataricus and enables us to predict enzyme candidates involved in the degradation processes downstream of 2-hydroxymucconic acid.

Probing the Mobilome: Discoveries in the Dynamic Microbiome

There has been an explosion of metagenomic data representing human, animal, and environmental microbiomes. This provides an unprecedented opportunity for comparative and longitudinal studies of many functional aspects of the microbiome that go beyond taxonomic classification, such as profiling genetic determinants of antimicrobial resistance, interactions with the host, potentially clinically relevant functions, and the role of mobile genetic elements (MGEs). One of the most important but least studied of these aspects are the MGEs, collectively referred to as the 'mobilome'. Here we elaborate on the benefits and limitations of using different metagenomic protocols, discuss the relative merits of various sequencing technologies, and highlight relevant bioinformatics tools and pipelines to predict the presence of MGEs and their microbial hosts.

An analysis of the IS6/IS26 family of insertion sequences: is it a single family?

The relationships within a curated set of 112 insertion sequences (ISs) currently assigned to the IS6 family, here re-named the IS6/IS26 family, in the ISFinder database were examined. The encoded DDE transposases include a helix-helix-turn-helix (H-HTH) potential DNA binding domain N-terminal to the catalytic (DDE) domain, but 10 from Clostridia include one or two additional N-terminal domains. The transposase phylogeny clearly separated 75 derived from bacteria from 37 from archaea. The longer bacterial transposases also clustered separately. The 65 shorter bacterial transposases, including Tnp26 from IS26, formed six clades but share significant conservation in the H-HTH domain and in a short extension at the N-terminus, and several amino acids in the catalytic domain are completely or highly conserved. At the outer ends of these ISs, 14 bp were strongly conserved as terminal inverted repeats (TIRs) with the first two bases (GG) and the seventh base (G) present in all except one IS. The longer bacterial transposases are only distantly related to the short bacterial transposases, with only some amino acids conserved. The TIR consensus was longer and only one IS started with GG. The 37 archaeal transposases are only distantly related to either the short or the long bacterial transposases and different residues were conserved. Their TIRs are loosely related to the bacterial TIR consensus but are longer and many do not begin with GG. As they do not fit well with most bacterial ISs, the inclusion of the archaeal ISs and the longer bacterial ISs in the IS6/IS26 family is not appropriate.

Multiple origins of prokaryotic and eukaryotic single-stranded DNA viruses from bacterial and archaeal plasmids

Single-stranded (ss) DNA viruses are a major component of the earth virome. In particular, the circular, Rep-encoding ssDNA (CRESS-DNA) viruses show high diversity and abundance in various habitats. By combining sequence similarity network and phylogenetic analyses of the replication proteins (Rep) belonging to the HUH endonuclease superfamily, we show that the replication machinery of the CRESS-DNA viruses evolved, on three independent occasions, from the Reps of bacterial rolling circle-replicating plasmids. The CRESS-DNA viruses emerged via recombination between such plasmids and cDNA copies of capsid genes of eukaryotic positive-sense RNA viruses. Similarly, the rep genes of prokaryotic DNA viruses appear to have evolved from HUH endonuclease genes of various bacterial and archaeal plasmids. Our findings also suggest that eukaryotic polyomaviruses and papillomaviruses with dsDNA genomes have evolved via parvoviruses from CRESS-DNA viruses. Collectively, our results shed light on the complex evolutionary history of a major class of viruses revealing its polyphyletic origins.

Most single-stranded DNA viruses have small genomes replicated by rolling circle mechanism which is initiated by the Rep protein. Here, using sequence similarity network and phylogenetic analyses, Kazlauskas et al. show that viral Reps evolved from Reps of bacterial and archaeal plasmids on multiple independent occasions.

Methylation guide RNA evolution in archaea: structure, function and genomic organization of 110 C/D box sRNA families across six Pyrobaculum species

Archaeal homologs of eukaryotic C/D box small nucleolar RNAs (C/D box sRNAs) guide precise 2′-O-methyl modification of ribosomal and transfer RNAs. Although C/D box sRNA genes constitute one of the largest RNA gene families in archaeal thermophiles, most genomes have incomplete sRNA gene annotation because reliable, fully automated detection methods are not available. We expanded and curated a comprehensive gene set across six species of the crenarchaeal genus Pyrobaculum, particularly rich in C/D box sRNA genes. Using high-throughput small RNA sequencing, specialized computational searches and comparative genomics, we analyzed 526 Pyrobaculum C/D box sRNAs, organizing them into 110 families based on synteny and conservation of guide sequences which determine methylation targets. We examined gene duplications and rearrangements, including one family that has expanded in a pattern similar to retrotransposed repetitive elements in eukaryotes. New training data and inclusion of kink-turn secondary structural features enabled creation of an improved search model. Our analyses provide the most comprehensive, dynamic view of C/D box sRNA evolutionary history within a genus, in terms of modification function, feature plasticity, and gene mobility.

Integrated mobile genetic elements in Thaumarchaeota

To explore the diversity of mobile genetic elements (MGE) associated with archaea of the phylum Thaumarchaeota, we exploited the property of most MGE to integrate into the genomes of their hosts. Integrated MGE (iMGE) were identified in 20 thaumarchaeal genomes amounting to 2 Mbp of mobile thaumarchaeal DNA. These iMGE group into five major classes: (i) proviruses, (ii) casposons, (iii) insertion sequence-like transposons, (iv) integrative-conjugative elements and (v) cryptic integrated elements. The majority of the iMGE belong to the latter category and might represent novel families of viruses or plasmids. The identified proviruses are related to tailed viruses of the order Caudovirales and to tailless icosahedral viruses with the double jelly-roll capsid proteins. The thaumarchaeal iMGE are all connected within a gene sharing network, highlighting pervasive gene exchange between MGE occupying the same ecological niche. The thaumarchaeal mobilome carries multiple auxiliary metabolic genes, including multicopper oxidases and ammonia monooxygenase subunit C (AmoC), and stress response genes, such as those for universal stress response proteins (UspA). Thus, iMGE might make important contributions to the fitness and adaptation of their hosts. We identified several iMGE carrying type I-B CRISPR-Cas systems and spacers matching other thaumarchaeal iMGE, suggesting antagonistic interactions between coexisting MGE and symbiotic relationships with the ir archaeal hosts.

Into the Thermus Mobilome: Presence, Diversity and Recent Activities of Insertion Sequences Across Thermus spp

A high level of transposon-mediated genome rearrangement is a common trait among microorganisms isolated from thermal environments, probably contributing to the extraordinary genomic plasticity and horizontal gene transfer (HGT) observed in these habitats. In this work, active and inactive insertion sequences (ISs) spanning the sequenced members of the genus Thermus were characterized, with special emphasis on three T. thermophilus strains: HB27, HB8, and NAR1. A large number of full ISs and fragments derived from different IS families were found, concentrating within megaplasmids present in most isolates. Potentially active ISs were identified through analysis of transposase integrity, and domestication-related transposition events of ISTth7 were identified in laboratory-adapted HB27 derivatives. Many partial copies of ISs appeared throughout the genome, which may serve as specific targets for homologous recombination contributing to genome rearrangement. Moreover, recruitment of IS1000 32 bp segments as spacers for CRISPR sequence was identified, pointing to the adaptability of these elements in the biology of these thermophiles. Further knowledge about the activity and functional diversity of ISs in this genus may contribute to the generation of engineered transposons as new genetic tools, and enrich our understanding of the outstanding plasticity shown by these thermophiles.

Development of a high sensitivity TaqMan-based PCR assay for the specific detection of Mycobacterium tuberculosis complex in both pulmonary and extrapulmonary specimens

Tuberculosis (TB) causes a heavy health burden worldwide, especially in developing countries. The need for the rapid and accurate diagnosis of TB has not been satisfied, especially for extra-pulmonary specimens or specimens with acid fast stain (AFS)-negative condition. Development and validation of a novel, sensitive and specific assay for diagnosing TB is essential. We developed IS4 primer/probe based on insertion sequence 6110 (IS6110). A qPCR assay was designed for detecting a specific region in IS6110 by BLAST. The IS4 primer/probe concentration, qPCR efficiency and various of PCR additives were evaluated and optimized. Thirty-four species of commonly isolated microorganisms were used for evaluating the analytical specificity. Moreover, 130 clinical specimens were collected for evaluating the performance versus Cobas TaqMan MTB (CTM) assay kit and culture. The amplification efficiencies of IS4 were 99.61% and 102.61% without and with internal control DNA (Bacteriophage Lambda), respectively. Dimethyl sulfoxide outperformed glycerol or BSA for eliciting the most effective amplification and the lowest limit of detection. In evaluating the clinical performance, various specimen types were collected. IS4 demonstrated a high degree of agreement (kappa = 0.71) with CTM. The clinical sensitivity and specificity of IS4 and CTM were 92.11% (35/38), 82.61% (76/92), 84.21% (32/38) and 95.65% (88/92), respectively. The clinical sensitivity and specificity of IS4 were similar for both pulmonary [92.00% (23/25) and 76.92% (30/39), respectively] and extrapulmonary [92.31% (12/13) and 86.79% (46/53), respectively] specimens. Among AFS-negative cases, the clinical sensitivity and specificity remained 90.48% (19/21) and 83.91% (73/87), respectively, with culture as the gold standard. We concluded that IS4, a new primer/probe pair for TaqMan based qPCR assay, was developed, optimized, and validated for the sensitive and specific detection of TB among various specimen types. The performance was not compromised under AFS-negative conditions.

New Insights on the Evolution of Genome Content: Population Dynamics of Transposable Elements in Flies and Humans

Understanding the abundance, diversity, and distribution of TEs in genomes is crucial to understand genome structure, function, and evolution. Advances in whole-genome sequencing techniques, as well as in bioinformatics tools, have increased our ability to detect and analyze the transposable element content in genomes. In addition to reference genomes, we now have access to population datasets in which multiple individuals within a species are sequenced. In this chapter, we highlight the recent advances in the study of TE population dynamics focusing on fruit flies and humans, which represent two extremes in terms of TE abundance, diversity, and activity. We review the most recent methodological approaches applied to the study of TE dynamics as well as the new knowledge on host factors involved in the regulation of TE activity. In addition to transposition rates, we also focus on TE deletion rates and on the selective forces that affect the dynamics of TEs in genomes.

Phenotypic and genomic comparison of Photorhabdus luminescens subsp. laumondii TT01 and a widely used rifampicin-resistant Photorhabdus luminescens laboratory strain

BACKGROUND

Photorhabdus luminescens is an enteric bacterium, which lives in mutualistic association with soil nematodes and is highly pathogenic for a broad spectrum of insects. A complete genome sequence for the type strain P. luminescens subsp. laumondii TT01, which was originally isolated in Trinidad and Tobago, has been described earlier. Subsequently, a rifampicin resistant P. luminescens strain has been generated with superior possibilities for experimental characterization. This strain, which is widely used in research, was described as a spontaneous rifampicin resistant mutant of TT01 and is known as TT01-Rif^R.

RESULTS

Unexpectedly, upon phenotypic comparison between the rifampicin resistant strain and its presumed parent TT01, major differences were found with respect to bioluminescence, pigmentation, biofilm formation, haemolysis as well as growth. Therefore, we renamed the strain TT01-Rif^R to DJC. To unravel the genomic basis of the observed differences, we generated a complete genome sequence for strain DJC using the PacBio long read technology. As strain DJC was supposed to be a spontaneous mutant, only few sequence differences were expected. In order to distinguish these from potential sequencing errors in the published TT01 genome, we re-sequenced a derivative of strain TT01 in parallel, also using the PacBio technology. The two TT01 genomes differed at only 30 positions. In contrast, the genome of strain DJC varied extensively from TT01, showing 13,000 point mutations, 330 frameshifts, and 220 strain-specific regions with a total length of more than 300 kb in each of the compared genomes.

CONCLUSIONS

According to the major phenotypic and genotypic differences, the rifampicin resistant P. luminescens strain, now named strain DJC, has to be considered as an independent isolate rather than a derivative of strain TT01. Strains TT01 and DJC both belong to P. luminescens subsp. laumondii.

Genome size evolution in the Archaea

What determines variation in genome size, gene content and genetic diversity at the broadest scales across the tree of life? Much of the existing work contrasts eukaryotes with prokaryotes, the latter represented mainly by Bacteria. But any general theory of genome evolution must also account for the Archaea, a diverse and ecologically important group of prokaryotes that represent one of the primary domains of cellular life. Here, we survey the extant diversity of Bacteria and Archaea, and ask whether the general principles of genome evolution deduced from the study of Bacteria and eukaryotes also apply to the archaeal domain. Although Bacteria and Archaea share a common prokaryotic genome architecture, the extant diversity of Bacteria appears to be much higher than that of Archaea. Compared with Archaea, Bacteria also show much greater genome-level specialisation to specific ecological niches, including parasitism and endosymbiosis. The reasons for these differences in long-term diversification rates are unclear, but might be related to fundamental differences in informational processing machineries and cell biological features that may favour archaeal diversification in harsher or more energy-limited environments. Finally, phylogenomic analyses suggest that the first Archaea were anaerobic autotrophs that evolved on the early Earth.

Multiple serine transposase dimers assemble the transposon-end synaptic complex during IS607-family transposition

IS607-family transposons are unusual because they do not have terminal inverted repeats or generate target site duplications. They encode two protein-coding genes, but only tnpA is required for transposition. Our X-ray structures confirm that TnpA is a member of the serine recombinase (SR) family, but the chemically-inactive quaternary structure of the dimer, along with the N-terminal location of the DNA binding domain, are different from other SRs. TnpA dimers from IS1535 cooperatively associate with multiple subterminal repeats, which together with additional nonspecific binding, form a nucleoprotein filament on one transposon end that efficiently captures a second unbound end to generate the paired-end complex (PEC). Formation of the PEC does not require a change in the dimeric structure of the catalytic domain, but remodeling of the C-terminal α-helical region is involved. We posit that the PEC recruits a chemically-active conformer of TnpA to the transposon end to initiate DNA chemistry.

Transposons: the agents of antibiotic resistance in bacteria

Transposons are a group of mobile genetic elements that are defined as a DNA sequence. Transposons can jump into different places of the genome; for this reason, they are called jumping genes. However, some transposons are always kept at the insertion site in the genome. Most transposons are inactivated and as a result, cannot move. Transposons are divided into two main groups: retrotransposons (class І) and DNA transposons (class ІІ). Retrotransposons are often found in eukaryotes. DNA transposons can be found in both eukaryotes and prokaryotes. The bacterial transposons belong to the DNA transposons and the Tn family, which are usually the carrier of additional genes for antibiotic resistance. Transposons can transfer from a plasmid to other plasmids or from a DNA chromosome to plasmid and vice versa that cause the transmission of antibiotic resistance genes in bacteria. The treatment of bacterial infectious diseases is difficult because of existing antibiotic resistance that part of this antibiotic resistance is caused by transposons. Bacterial infectious diseases are responsible for the increasing rise in world mortality rate. In this review, transposons and their roles have been studied in bacterial antibiotic resistance, in detail.

Role of free DNA ends and protospacer adjacent motifs for CRISPR DNA uptake in Pyrococcus furiosus

To acquire CRISPR–Cas immunity against invasive mobile genetic elements, prokaryotes must first integrate fragments of foreign DNA into their genomic CRISPR arrays for use in future invader silencing. Here, we found that the hyperthermophilic archaeaon, Pyrococcus furiosus, actively incorporates DNA fragments (spacers) from both plasmid (foreign) and host genome (self) sequences into its seven CRISPR loci. The majority of new spacers were derived from DNA immediately downstream from a 5′-CCN-3′ protospacer adjacent motif (PAM) that is critical for invader targeting. Interestingly, spacers were preferentially acquired from genome or plasmid regions corresponding to active transposons, CRISPR loci, ribosomal RNA genes, rolling circle origins of replication, and areas where plasmids recombined with the host chromosome. A common feature of the highly sampled spacers is that they arise from DNA regions expected to undergo DNA nicking and/or double-strand breaks. Taken together with recent results from bacterial systems, our findings indicate that free DNA termini and PAMs are conserved features important for CRISPR spacer uptake in diverse prokaryotes and CRISPR–Cas systems. Moreover, lethal self-targeting by CRISPR systems may contribute to host genome stability by eliminating cells undergoing active transposon mobility or chromosomal uptake of autonomously replicating foreign mobile genetic elements.

Viruses of archaea: Structural, functional, environmental and evolutionary genomics

Viruses of archaea represent one of the most enigmatic parts of the virosphere. Most of the characterized archaeal viruses infect extremophilic hosts and display remarkable diversity of virion morphotypes, many of which have never been observed among viruses of bacteria or eukaryotes. The uniqueness of the virion morphologies is matched by the distinctiveness of the genomes of these viruses, with ∼75% of genes encoding unique proteins, refractory to functional annotation based on sequence analyses. In this review, we summarize the state-of-the-art knowledge on various aspects of archaeal virus genomics. First, we outline how structural and functional genomics efforts provided valuable insights into the functions of viral proteins and revealed intricate details of the archaeal virus-host interactions. We then highlight recent metagenomics studies, which provided a glimpse at the diversity of uncultivated viruses associated with the ubiquitous archaea in the oceans, including Thaumarchaeota, Marine Group II Euryarchaeota, and others. These findings, combined with the recent discovery that archaeal viruses mediate a rapid turnover of thaumarchaea in the deep sea ecosystems, illuminate the prominent role of these viruses in the biosphere. Finally, we discuss the origins and evolution of archaeal viruses and emphasize the evolutionary relationships between viruses and non-viral mobile genetic elements. Further exploration of the archaeal virus diversity as well as functional studies on diverse virus-host systems are bound to uncover novel, unexpected facets of the archaeal virome.

An Important Role for Purifying Selection in Archaeal Genome Evolution

As the null hypothesis of genome evolution, population genetic theory suggests that selection strength controls genome size. Through the process of genetic drift, this theory predicts that compact genomes are maintained by strong purifying selection while complex genomes are enabled by weak purifying selection. It offers a unifying framework that explains why prokaryotic genomes are much smaller than their eukaryotic counterparts. However, recent findings suggest that bigger prokaryotic genomes appear to experience stronger purifying selection, indicating that purifying selection may not dominate prokaryotic genome evolution. Since archaeal genomes were underrepresented in those studies, generalization of the conclusions to both archaeal and bacterial genomes may not be warranted. In this study, we revisited this matter by focusing on archaeal and bacterial genomes separately. We found that bigger bacterial genomes indeed experienced stronger purifying selection, but the opposite was observed in archaeal genomes. This new finding would predict an enrichment of noncoding sequences in large archaeal genomes, which was confirmed by an analysis of coding density. In contrast, coding density remained stable regardless of bacterial genome size. In conclusion, this study suggests that purifying selection may play a more important role in archaeal genome evolution than previously hypothesized, indicating that there could be a major difference between the evolutionary regimes of Archaea and Bacteria. IMPORTANCE The evolution of genome complexity is a fundamental question in biology. A hallmark of eukaryotic genome complexity is that larger genomes tend to have more noncoding sequences, which are believed to be minimal in archaeal and bacterial genomes. However, we found that archaeal genomes also possessed this eukaryotic feature while bacterial genomes did not. This could be predicted from our analysis on genetic drift, which showed a relaxation of purifying selection in larger archaeal genomes, also a eukaryotic feature. In contrast, the opposite was evident in bacterial genomes.

Metagenomic analysis reveals wastewater treatment plants as hotspots of antibiotic resistance genes and mobile genetic elements

The intensive use of antibiotics results in their continuous release into the environment and the subsequent widespread occurrence of antibiotic resistant bacteria (ARB), antibiotic resistance genes (ARGs) and mobile genetic elements (MGEs). This study used Illumina high-throughput sequencing to investigate the broad-spectrum profiles of both ARGs and MGEs in activated sludge and anaerobically digested sludge from a full-scale wastewater treatment plant. A pipeline for identifying antibiotic resistance determinants was developed that consisted of four categories: gene transfer potential, ARG potential, ARGs pathway and ARGs phylogenetic origin. The metagenomic analysis showed that the activated sludge and the digested sludge exhibited different microbial communities and changes in the types and occurrence of ARGs and MGEs. In total, 42 ARGs subtypes were identified in the activated sludge, while 51 ARG subtypes were detected in the digested sludge. Additionally, MGEs including plasmids, transposons, integrons (intI1) and insertion sequences (e.g. ISSsp4, ISMsa21 and ISMba16) were abundant in the two sludge samples. The co-occurrence pattern between ARGs and microbial taxa revealed by network analysis indicated that some environmental bacteria (e.g. Clostridium and Nitrosomonas) might be potential hosts of multiple ARGs. The findings increase our understanding of WWTPs as hotspots of ARGs and MGEs, and contribute towards preventing their release into the downstream environment.

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.

Metabolic and evolutionary patterns in the extremely acidophilic archaeon Ferroplasma acidiphilum YT

Ferroplasmaceae represent ubiquitous iron-oxidising extreme acidophiles with a number of unique physiological traits. In a genome-based study of Ferroplasma acidiphilum YT, the only species of the genus Ferroplasma with a validly published name, we assessed its central metabolism and genome stability during a long-term cultivation experiment. Consistently with physiology, the genome analysis points to F. acidiphilum YT having an obligate peptidolytic oligotrophic lifestyle alongside with anaplerotic carbon assimilation. This narrow trophic specialisation abridges the sugar uptake, although all genes for glycolysis and gluconeogenesis, including bifunctional unidirectional fructose 1,6-bisphosphate aldolase/phosphatase, have been identified. Pyruvate and 2-oxoglutarate dehydrogenases are substituted by 'ancient' CoA-dependent pyruvate and alpha-ketoglutarate ferredoxin oxidoreductases. In the lab culture, after ~550 generations, the strain exhibited the mutation rate of ≥1.3 × 10-8 single nucleotide substitutions per site per generation, which is among the highest values recorded for unicellular organisms. All but one base substitutions were G:C to A:T, their distribution between coding and non-coding regions and synonymous-to-non-synonymous mutation ratios suggest the neutral drift being a prevalent mode in genome evolution in the lab culture. Mutations in nature seem to occur with lower frequencies, as suggested by a remarkable genomic conservation in F. acidiphilum YT variants from geographically distant populations.

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.

DNA repeat sequences: diversity and versatility of functions

Although discovered decades ago, the molecular identification, the diversity and versatility of functions, and the evolutionary origin of repeat DNA sequences (REPs) containing palindromic units in prokaryotes are now bringing attention to a wide range of biological scientists. A brief account of the current state of the repeat DNA sequences is presented here.

The Tn3-family of Replicative Transposons

Transposons of the Tn3 family form a widespread and remarkably homogeneous group of bacterial transposable elements in terms of transposition functions and an extremely versatile system for mediating gene reassortment and genomic plasticity owing to their modular organization. They have made major contributions to antimicrobial drug resistance dissemination or to endowing environmental bacteria with novel catabolic capacities. Here, we discuss the dynamic aspects inherent to the diversity and mosaic structure of Tn3-family transposons and their derivatives. We also provide an overview of current knowledge of the replicative transposition mechanism of the family, emphasizing most recent work aimed at understanding this mechanism at the biochemical level. Previous and recent data are put in perspective with those obtained for other transposable elements to build up a tentative model linking the activities of the Tn3-family transposase protein with the cellular process of DNA replication, suggesting new lines for further investigation. Finally, we summarize our current view of the DNA site-specific recombination mechanisms responsible for converting replicative transposition intermediates into final products, comparing paradigm systems using a serine recombinase with more recently characterized systems that use a tyrosine recombinase.

Plasmids from Euryarchaeota

Many plasmids have been described in Euryarchaeota, one of the three major archaeal phyla, most of them in salt-loving haloarchaea and hyperthermophilic Thermococcales. These plasmids resemble bacterial plasmids in terms of size (from small plasmids encoding only one gene up to large megaplasmids) and replication mechanisms (rolling circle or theta). Some of them are related to viral genomes and form a more or less continuous sequence space including many integrated elements. Plasmids from Euryarchaeota have been useful for designing efficient genetic tools for these microorganisms. In addition, they have also been used to probe the topological state of plasmids in species with or without DNA gyrase and/or reverse gyrase. Plasmids from Euryarchaeota encode both DNA replication proteins recruited from their hosts and novel families of DNA replication proteins. Euryarchaeota form an interesting playground to test evolutionary hypotheses on the origin and evolution of viruses and plasmids, since a robust phylogeny is available for this phylum. Preliminary studies have shown that for different plasmid families, plasmids share a common gene pool and coevolve with their hosts. They are involved in gene transfer, mostly between plasmids and viruses present in closely related species, but rarely between cells from distantly related archaeal lineages. With few exceptions (e.g., plasmids carrying gas vesicle genes), most archaeal plasmids seem to be cryptic. Interestingly, plasmids and viral genomes have been detected in extracellular membrane vesicles produced by Thermococcales, suggesting that these vesicles could be involved in the transfer of viruses and plasmids between cells.

Sense overlapping transcripts in IS1341-type transposase genes are functional non-coding RNAs in archaea

The existence of sense overlapping transcripts that share regulatory and coding information in the same genomic sequence shows an additional level of prokaryotic gene expression complexity. Here we report the discovery of ncRNAs associated with IS1341-type transposase (tnpB) genes, at the 3'-end of such elements, with examples in archaea and bacteria. Focusing on the model haloarchaeon Halobacterium salinarum NRC-1, we show the existence of sense overlapping transcripts (sotRNAs) for all its IS1341-type transposases. Publicly available transcriptome compendium show condition-dependent differential regulation between sotRNAs and their cognate genes. These sotRNAs allowed us to find a UUCA tetraloop motif that is present in other archaea (ncRNA family HgcC) and in a H. salinarum intergenic ncRNA derived from a palindrome associated transposable elements (PATE). Overexpression of one sotRNA and the PATE-derived RNA harboring the tetraloop motif improved H. salinarum growth, indicating that these ncRNAs are functional.

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.

Evolutionary Trends of the Transposase-Encoding Open Reading Frames A and B (orfA and orfB) of the Mycobacterial IS6110 Insertion Sequence

Background

The IS6110 insertion sequence, a member of the IS3 family of insertion sequences, was found to be specific to the Mycobacterium tuberculosis complex (MTBC). Although IS6110 has been extensively characterized as a transposable genetic marker, the evolutionary history of its own transposase-encoding sequence has not, to the best of our knowledge, been investigated.

Methodology/Principal Findings

Here we explored the evolution of the IS6110 sequence by analysing the genetic variability and the selective forces acting on its transposase-encoding open reading frames (ORFs) A and B (orfA and orfB). For this purpose, we used a strain collection consisting of smooth tubercle bacilli (STB), an early branching lineage of the MTBC, and present-day M. tuberculosis strains representing the full breadth of genetic diversity in Tunisia. In each ORF, we found a major haplotype that dominated over a flat distribution of rare descendent haplotypes, consisting mainly of single- and double-nucleotide variant singletons. The predominant haplotypes consisted of both ancestral and present-day strains, suggesting that IS6110 acquisition predated the emergence of the MTBC. There was no evidence of recombination and both ORFs were subjected to strict purifying selection, as demonstrated by their dN/dS ratios (0.29 and 0.51, respectively), as well as their significantly negative Tajima’s D statistics. Strikingly, the purifying selection acting on orfA proved much more stringent, suggesting its critical role in regulating the transpositional process. Maximum likelihood analyses further excluded any possibility of positive selection acting on single amino acid residues.

Conclusions/Significance

Taken together our data fit with an evolutionary scenario according to which the observed variability pattern of the IS6110 transposase-encoding ORFs is generated mainly through random point mutations that accrued on a functionally optimal IS6110 copy, whose acquisition predated the emergence of the MTBC complex. Background selection acting against deleterious mutations led to an excess of low-frequency variants.

DNA transposons have colonized the genome of the giant virus Pandoravirus salinus

Background

Transposable elements are mobile DNA sequences that are widely distributed in prokaryotic and eukaryotic genomes, where they represent a major force in genome evolution. However, transposable elements have rarely been documented in viruses, and their contribution to viral genome evolution remains largely unexplored. Pandoraviruses are recently described DNA viruses with genome sizes that exceed those of some prokaryotes, rivaling parasitic eukaryotes. These large genomes appear to include substantial noncoding intergenic spaces, which provide potential locations for transposable element insertions. However, no mobile genetic elements have yet been reported in pandoravirus genomes.

Results

Here, we report a family of miniature inverted-repeat transposable elements (MITEs) in the Pandoravirus salinus genome, representing the first description of a virus populated with a canonical transposable element family that proliferated by transposition within the viral genome. The MITE family, which we name Submariner, includes 30 copies with all the hallmarks of MITEs: short length, terminal inverted repeats, TA target site duplication, and no coding capacity. Submariner elements show signs of transposition and are undetectable in the genome of Pandoravirus dulcis, the closest known relative Pandoravirus salinus. We identified a DNA transposon related to Submariner in the genome of Acanthamoeba castellanii, a species thought to host pandoraviruses, which contains remnants of coding sequence for a Tc1/mariner transposase. These observations suggest that the Submariner MITEs of P. salinus belong to the widespread Tc1/mariner superfamily and may have been mobilized by an amoebozoan host. Ten of the 30 MITEs in the P. salinus genome are located within coding regions of predicted genes, while others are close to genes, suggesting that these transposons may have contributed to viral genetic novelty.

Conclusions

Our discovery highlights the remarkable ability of DNA transposons to colonize and shape genomes from all domains of life, as well as giant viruses. Our findings continue to blur the division between viral and cellular genomes, adhering to the emerging view that the content, dynamics, and evolution of the genomes of giant viruses do not substantially differ from those of cellular organisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0145-1) contains supplementary material, which is available to authorized users.

A cis-encoded sRNA, Hfq and mRNA secondary structure act independently to suppress IS200 transposition

IS200 is found throughout Enterobacteriaceae and transposes at a notoriously low frequency. In addition to the transposase protein (TnpA), IS200 encodes an uncharacterized Hfq-binding sRNA that is encoded opposite to the tnpA 5'UTR. In the current work we asked if this sRNA represses tnpA expression. We show here that the IS200 sRNA (named art200 for antisense regulator of transposase IS200) basepairs with tnpA to inhibit translation initiation. Unexpectedly, art200-tnpA pairing is limited to 40 bp, despite 90 nt of perfect complementarity. Additionally, we show that Hfq and RNA secondary structure in the tnpA 5'UTR each repress tnpA expression in an art200-independent manner. Finally, we show that disrupting translational control of tnpA expression leads to increased IS200 transposition in E. coli. The current work provides new mechanistic insight into why IS200 transposition is so strongly suppressed. The possibility of art200 acting in trans to regulate a yet-unidentified target is discussed as well as potential applications of the IS200 system for designing novel riboregulators.

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.

pTC Plasmids from Sulfolobus Species in the Geothermal Area of Tengchong, China: Genomic Conservation and Naturally-Occurring Variations as a Result of Transposition by Mobile Genetic Elements

Plasmids occur frequently in Archaea. A novel plasmid (denoted pTC1) containing typical conjugation functions has been isolated from Sulfolobus tengchongensis RT8-4, a strain obtained from a hot spring in Tengchong, China, and characterized. The plasmid is a circular double-stranded DNA molecule of 20,417 bp. Among a total of 26 predicted pTC1 ORFs, 23 have homologues in other known Sulfolobus conjugative plasmids (CPs). pTC1 resembles other Sulfolobus CPs in genome architecture, and is most highly conserved in the genomic region encoding conjugation functions. However, attempts to demonstrate experimentally the capacity of the plasmid for conjugational transfer were unsuccessful. A survey revealed that pTC1 and its closely related plasmid variants were widespread in the geothermal area of Tengchong. Variations of the plasmids at the target sites for transposition by an insertion sequence (IS) and a miniature inverted-repeat transposable element (MITE) were readily detected. The IS was efficiently inserted into the pTC1 genome, and the inserted sequence was inactivated and degraded more frequently in an imprecise manner than in a precise manner. These results suggest that the host organism has evolved a strategy to maintain a balance between the insertion and elimination of mobile genetic elements to permit genomic plasticity while inhibiting their fast spreading.

The diversity of prokaryotic DDE transposases of the mutator superfamily, insertion specificity, and association with conjugation machineries

Transposable elements (TEs) are major components of both prokaryotic and eukaryotic genomes and play a significant role in their evolution. In this study, we have identified new prokaryotic DDE transposase families related to the eukaryotic Mutator-like transposases. These genes were retrieved by cascade PSI-Blast using as initial query the transposase of the streptococcal integrative and conjugative element (ICE) TnGBS2. By combining secondary structure predictions and protein sequence alignments, we predicted the DDE catalytic triad and the DNA-binding domain recognizing the terminal inverted repeats. Furthermore, we systematically characterized the organization and the insertion specificity of the TEs relying on these prokaryotic Mutator-like transposases (p-MULT) for their mobility. Strikingly, two distant TE families target their integration upstream σ_A dependent promoters. This allowed us to identify a transposase sequence signature associated with this unique insertion specificity and to show that the dissymmetry between the two inverted repeats is responsible for the orientation of the insertion. Surprisingly, while DDE transposases are generally associated with small and simple transposons such as insertion sequences (ISs), p-MULT encoding TEs show an unprecedented diversity with several families of IS, transposons, and ICEs ranging in size from 1.1 to 52 kb.