Causes of insertion sequences abundance in prokaryotic genomes

A metagenomics pipeline reveals insertion sequence-driven evolution of the microbiota

Insertion sequence (IS) elements are mobile genetic elements in bacterial genomes that support adaptation. We developed a database of IS elements coupled to a computational pipeline that identifies IS element insertions in the microbiota. We discovered that diverse IS elements insert into the genomes of intestinal bacteria regardless of human host lifestyle. These insertions target bacterial accessory genes that aid in their adaptation to unique environmental conditions. Using IS expansion in Bacteroides, we show that IS activity leads to the insertion of "hot spots" in accessory genes. We show that IS insertions are stable and can be transferred between humans. Extreme environmental perturbations force IS elements to fall out of the microbiota, and many fail to rebound following homeostasis. Our work shows that IS elements drive bacterial genome diversification within the microbiota and establishes a framework for understanding how strain-level variation within the microbiota impacts human health.

Investigating the impact of insertion sequences and transposons in the genomes of the most significant phytopathogenic bacteria

Genetic variability in phytopathogens is one of the main problems encountered for effective plant disease control. This fact may be related to the presence of transposable elements (TEs), but little is known about their role in host genomes. Here, we performed the most comprehensive analysis of insertion sequences (ISs) and transposons (Tns) in the genomes of the most important bacterial plant pathogens. A total of 35 692 ISs and 71 transposons were identified in 270 complete genomes. The level of pathogen-host specialization was found to be a significant determinant of the element distribution among the species. Some Tns were identified as carrying virulence factors, such as genes encoding effector proteins of the type III secretion system and resistance genes for the antimicrobial streptomycin. Evidence for IS-mediated ectopic recombination was identified in Xanthomonas genomes. Moreover, we found that IS elements tend to be inserted in regions near virulence and fitness genes, such ISs disrupting avirulence genes in X. oryzae genomes. In addition, transcriptome analysis under different stress conditions revealed differences in the expression of genes encoding transposases in the Ralstonia solanacearum, X. oryzae, and P. syringae species. Lastly, we also investigated the role of Tns in regulation via small noncoding regulatory RNAs and found these elements may target plant-cell transcriptional activators. Taken together, the results indicate that TEs may have a fundamental role in variability and virulence in plant pathogenic bacteria.

Increasing transposase abundance with ocean depth correlates with a particle-associated lifestyle

Transposases are mobile genetic elements that move within and between genomes, promoting genomic plasticity in microorganisms. In marine microbial communities, the abundance of transposases increases with depth, but the reasons behind this trend remain unclear. Our analysis of metagenomes from the Tara Oceans and Malaspina Expeditions suggests that a particle-associated lifestyle is the main covariate for the high occurrence of transposases in the deep ocean, and this trend holds true for individual genomes as well as in a community-wide sense. We observed a strong and depth-independent correlation between transposase abundance and the presence of biofilm-associated genes, as well as the prevalence of secretory enzymes. This suggests that mobile genetic elements readily propagate among microbial communities within crowded biofilms. Furthermore, we show that particle association positively correlates with larger genome size, which is in turn associated with higher transposase abundance. Cassette sequences associated with transposons are enriched with genes related to defense mechanisms, which are more highly expressed in the deep sea. Thus, while transposons spread at the expense of their microbial hosts, they also introduce novel genes and potentially benefit the hosts in helping to compete for limited resources. Overall, our results suggest a new understanding of deep ocean particles as highways for gene sharing among defensively oriented microbial genomes.IMPORTANCEGenes can move within and between microbial genomes via mobile genetic elements, which include transposases and transposons. In the oceans, there is a puzzling increase in transposase abundance in microbial genomes as depth increases. To gain insight into this trend, we conducted an extensive analysis of marine microbial metagenomes and metatranscriptomes. We found a significant correlation between transposase abundance and a particle-associated lifestyle among marine microbes at both the metagenome and genome-resolved levels. We also observed a link between transposase abundance and genes related to defense mechanisms. These results suggest that as microbes become densely packed into crowded particles, mobile genes are more likely to spread and carry genetic material that provides a competitive advantage in crowded habitats. This may enable deep sea microbes to effectively compete in such environments.

Genomic insights into virulence, antimicrobial resistance, and adaptation acumen of Escherichia coli isolated from an urban environment

Populations of common commensal bacteria such as Escherichia coli undergo genetic changes by the acquisition of certain virulence and antimicrobial resistance (AMR) encoding genetic elements leading to the emergence of pathogenic strains capable of surviving in the previously uninhabited or protected niches. These bacteria are also reported to be prevalent in the environment where they survive by adopting various recombination strategies to counter microflora of the soil and water, under constant selection pressure(s). In this study, we performed molecular characterization, phenotypic AMR analysis, and whole genome sequencing (WGS) of E. coli (n = 37) isolated from soil and surface water representing the urban and peri-urban areas. The primary aim of this study was to understand the genetic architecture and pathogenic acumen exhibited by environmental E. coli. WGS-based analysis entailing resistome and virulome profiling indicated the presence of various virulence (adherence, iron uptake, and toxins) and AMR encoding genes, including blaNDM-5 in the environmental isolates. A majority of our isolates belonged to phylogroup B1 (73%). A few isolates in our collection were of sequence type(s) (ST) 58 and 224 that could have emerged recently as clonal lineages and might pose risk of infection/transmission. Mobile genetic elements (MGEs) such as plasmids (predominantly) of the IncF family, prophages, pipolins, and insertion elements such as IS1 and IS5 were also observed to exist, which may presumably aid in the propagation of genes encoding resistance against antimicrobial drugs. The observed high prevalence of MGEs associated with multidrug resistance in pathogenic E. coli isolates belonging to the phylogroup B1 underscores the need for extended surveillance to keep track of and prevent the transmission of the bacterium to certain vulnerable human and animal populations.

IMPORTANCE

Evolutionary patterns of E. coli bacteria convey that they evolve into highly pathogenic forms by acquiring fitness advantages, such as AMR, and various virulence factors through the horizontal gene transfer (HGT)-mediated acquisition of MGEs. However, limited research on the genetic profiles of environmental E. coli, particularly from India, hinders our understanding of their transition to pathogenic forms and impedes the adoption of a comprehensive approach to address the connection between environmentally dwelling E. coli populations and human and veterinary public health. This study focuses on high-resolution genomic analysis of the environmental E. coli isolates aiming to understand the genetic similarities and differences among isolates from different environmental niches and uncover the survival strategies employed by these bacteria to thrive in their surroundings. Our approach involved molecular characterization of environmental samples using PCR-based DNA fingerprinting and subsequent WGS analysis. This multidisciplinary approach is likely to provide valuable insights into the understanding of any potential spill-over to human and animal populations and locales. Investigating these environmental isolates has significant potential for developing epidemiological strategies against transmission and understanding niche-specific evolutionary patterns.

Phage-plasmids promote recombination and emergence of phages and plasmids

Phages and plasmids are regarded as distinct types of mobile genetic elements that drive bacterial evolution by horizontal gene transfer. However, the distinction between both types is blurred by the existence of elements known as prophage-plasmids or phage-plasmids, which transfer horizontally between cells as viruses and vertically within cellular lineages as plasmids. Here, we study gene flow between the three types of elements. We show that the gene repertoire of phage-plasmids overlaps with those of phages and plasmids. By tracking recent recombination events, we find that phage-plasmids exchange genes more frequently with plasmids than with phages, and that direct gene exchange between plasmids and phages is less frequent in comparison. The results suggest that phage-plasmids can mediate gene flow between plasmids and phages, including exchange of mobile element core functions, defense systems, and antibiotic resistance. Moreover, a combination of gene transfer and gene inactivation may result in the conversion of elements. For example, gene loss turns P1-like phage-plasmids into integrative prophages or into plasmids (that are no longer phages). Remarkably, some of the latter have acquired conjugation-related functions to became mobilisable by conjugation. Thus, our work indicates that phage-plasmids can play a key role in the transfer of genes across mobile elements within their hosts, and can act as intermediates in the conversion of one type of element into another.

From accurate genome sequence to biotechnological application: The thermophile Mycolicibacterium hassiacum as experimental model

Mycobacteria constitute a large group of microorganisms belonging to the phylum Actinobacteria encompassing some of the most relevant pathogenic bacteria and many saprophytic isolates that share a unique and complex cell envelope. Also unique to this group is the extensive capability to use and synthesize sterols, a class of molecules that include active signalling compounds of pharmaceutical use. However, few mycobacterial species and strains have been established as laboratory models to date, Mycolicibacterium smegmatis mc2 155 being the most common one. In this work, we focus on the use of a thermophilic mycobacterium, Mycolicibacterium hassiacum, which grows optimally above 50°C, as an emerging experimental model valid to extend our general knowledge of mycobacterial biology as well as for application purposes. To that end, accurate genomic sequences are key for gene mining, the study of pathogenicity or lack thereof and the potential for gene transfer. The combination of long- and short-massive sequencing technologies is strictly necessary to remove biases caused by errors specific to long-reads technology. By doing so in M. hassiacum, we obtained from the curated genome clues regarding the genetic manipulation potential of this microorganism from the presence of insertion sequences, CRISPR-Cas, type VII ESX secretion systems, as well as lack of plasmids. Finally, as a proof of concept of the applicability of M. hassiacum as a laboratory and industrial model, we used this high-quality genome of M. hassiacum to successfully knockout a gene involved in the use of phytosterols as source of carbon and energy, using an improved gene cassette for thermostable selection and a transformation protocol at high temperature.

Identifying and tracking mobile elements in evolving compost communities yields insights into the nanobiome

Microbial evolution is driven by rapid changes in gene content mediated by horizontal gene transfer (HGT). While mobile genetic elements (MGEs) are important drivers of gene flux, the nanobiome-the zoo of Darwinian replicators that depend on microbial hosts-remains poorly characterised. New approaches are necessary to increase our understanding beyond MGEs shaping individual populations, towards their impacts on complex microbial communities. A bioinformatic pipeline (xenoseq) was developed to cross-compare metagenomic samples from microbial consortia evolving in parallel, aimed at identifying MGE dissemination, which was applied to compost communities which underwent periodic mixing of MGEs. We show that xenoseq can distinguish movement of MGEs from demographic changes in community composition that otherwise confounds identification, and furthermore demonstrate the discovery of various unexpected entities. Of particular interest was a nanobacterium of the candidate phylum radiation (CPR) which is closely related to a species identified in groundwater ecosystems (Candidatus Saccharibacterium), and appears to have a parasitic lifestyle. We also highlight another prolific mobile element, a 313 kb plasmid hosted by a Cellvibrio lineage. The host was predicted to be capable of nitrogen fixation, and acquisition of the plasmid coincides with increased ammonia production. Taken together, our data show that new experimental strategies combined with bioinformatic analyses of metagenomic data stand to provide insight into the nanobiome as a driver of microbial community evolution.

Gill-associated bacteria are homogeneously selected in amphibious mangrove crabs to sustain host intertidal adaptation

BACKGROUND

The transition from water to air is a key event in the evolution of many marine organisms to access new food sources, escape water hypoxia, and exploit the higher and temperature-independent oxygen concentration of air. Despite the importance of microorganisms in host adaptation, their contribution to overcoming the challenges posed by the lifestyle changes from water to land is not well understood. To address this, we examined how microbial association with a key multifunctional organ, the gill, is involved in the intertidal adaptation of fiddler crabs, a dual-breathing organism.

RESULTS

Electron microscopy revealed a rod-shaped bacterial layer tightly connected to the gill lamellae of the five crab species sampled across a latitudinal gradient from the central Red Sea to the southern Indian Ocean. The gill bacterial community diversity assessed with 16S rRNA gene amplicon sequencing was consistently low across crab species, and the same actinobacterial group, namely Ilumatobacter, was dominant regardless of the geographic location of the host. Using metagenomics and metatranscriptomics, we detected that these members of actinobacteria are potentially able to convert ammonia to amino acids and may help eliminate toxic sulphur compounds and carbon monoxide to which crabs are constantly exposed.

CONCLUSIONS

These results indicate that bacteria selected on gills can play a role in the adaptation of animals in dynamic intertidal ecosystems. Hence, this relationship is likely to be important in the ecological and evolutionary processes of the transition from water to air and deserves further attention, including the ontogenetic onset of this association. Video Abstract.

Persistence of antibiotic resistance from animal agricultural effluents to surface water revealed by genome-centric metagenomics

Concerns about antibiotic resistance genes (ARGs) released from wastewaters of livestock or fish farming into the natural environment are increasing, but studies on unculturable bacteria related to the dissemination of antibiotic resistance are limited. Here, we reconstructed 1100 metagenome-assembled genomes (MAGs) to assess the impact of microbial antibiotic resistome and mobilome in wastewaters discharged to Korean rivers. Our results indicate that ARGs harbored in the MAGs were disseminated from wastewater effluents into downstream rivers. Moreover, it was found that ARGs are more commonly co-localized with mobile genetic elements (MGEs) in agricultural wastewater than in river water. Among the effluent-derived phyla, uncultured members of the superphylum Patescibacteria possessed a high number of MGEs, along with co-localized ARGs. Our findings suggest that members of the Patesibacteria are a potential vector for propagating ARGs into the environmental community. Therefore, we propose that the dissemination of ARGs by uncultured bacteria should be further investigated in multiple environments.

Targeted IS-element sequencing uncovers transposition dynamics during selective pressure in enterococci

Insertion sequences (IS) are simple transposons implicated in the genome evolution of diverse pathogenic bacterial species. Enterococci have emerged as important human intestinal pathogens with newly adapted virulence potential and antibiotic resistance. These genetic features arose in tandem with large-scale genome evolution mediated by mobile elements. Pathoadaptation in enterococci is thought to be mediated in part by the IS element IS256 through gene inactivation and recombination events. However, the regulation of IS256 and the mechanisms controlling its activation are not well understood. Here, we adapt an IS256-specfic deep sequencing method to describe how chronic lytic phage infection drives widespread diversification of IS256 in E. faecalis and how antibiotic exposure is associated with IS256 diversification in E. faecium during a clinical human infection. We show through comparative genomics that IS256 is primarily found in hospital-adapted enterococcal isolates. Analyses of IS256 transposase gene levels reveal that IS256 mobility is regulated at the transcriptional level by multiple mechanisms in E. faecalis, indicating tight control of IS256 activation in the absence of selective pressure. Our findings reveal that stressors such as phages and antibiotic exposure drives rapid genome-scale transposition in the enterococci. IS256 diversification can therefore explain how selective pressures mediate evolution of the enterococcal genome, ultimately leading to the emergence of dominant nosocomial lineages that threaten human health.

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.

One to host them all: genomics of the diverse bacterial endosymbionts of the spider Oedothorax gibbosus

Bacterial endosymbionts of the groups Wolbachia, Cardinium and Rickettsiaceae are well known for their diverse effects on their arthropod hosts, ranging from mutualistic relationships to reproductive phenotypes. Here, we analysed a unique system in which the dwarf spider Oedothorax gibbosus is co-infected with up to five different endosymbionts affiliated with Wolbachia, 'Candidatus Tisiphia' (formerly Torix group Rickettsia), Cardinium and Rhabdochlamydia. Using short-read genome sequencing data, we show that the endosymbionts are heterogeneously distributed among O. gibbosus populations and are frequently found co-infecting spider individuals. To study this intricate host-endosymbiont system on a genome-resolved level, we used long-read sequencing to reconstruct closed genomes of the Wolbachia, 'Ca. Tisiphia' and Cardinium endosymbionts. We provide insights into the ecology and evolution of the endosymbionts and shed light on the interactions with their spider host. We detected high quantities of transposable elements in all endosymbiont genomes and provide evidence that ancestors of the Cardinium, 'Ca. Tisiphia' and Wolbachia endosymbionts have co-infected the same hosts in the past. Our findings contribute to broadening our knowledge about endosymbionts infecting one of the largest animal phyla on Earth and show the usefulness of transposable elements as an evolutionary 'contact-tracing' tool.

Hydrogen peroxide treatment induces the transposition of an insertion sequence in Deinococcus radiopugnans DY59

Deinococcus radiopugnans DY59 (formerly Deinococcus swuensis DY59) is a radiation-resistant bacterium isolated from soil. From the 3.5 Mb genomic DNA sequence of strain DY59 (December 2014), 31 insertion sequence (IS) elements of six IS families including IS1, IS4, IS5, IS66, IS630, and IS701 and five unclassified IS elements were detected. Upon induction of oxidative stress with 80 and 100 mM H₂O₂, the unique ISs of the IS4 family member were actively translocated into a carotenoid biosynthesis gene phytoene desaturase (QR90_10400), resulting in non-pigment phenotypic selection. Therefore, these active transpositions of a specific IS family member were induced by oxidative stress at 80 and 100 mM H₂O₂. Furthermore, D. radiopugnans DY59 exhibited extremely higher MIC values against H₂O₂ treatment. To explain this phenomenon, qRT-PCR was conducted to assess the expression levels of catalase and three LysR family regulators. Our findings indicated that the ISDrpg2 and ISDrpg3 elements of the IS4 family were actively transposed into the phytoene desaturase gene by H₂O₂ treatment via replicative transposition. However, high H₂O₂ resistance did not originate from H₂O₂-induced expression of catalase and LysR family regulators.

Repeat sequences limit the effectiveness of lateral gene transfer and favored the evolution of meiotic sex in early eukaryotes

The transition from prokaryotic lateral gene transfer to eukaryotic meiotic sex is poorly understood. Phylogenetic evidence suggests that it was tightly linked to eukaryogenesis, which involved an unprecedented rise in both genome size and the density of genetic repeats. Expansion of genome size raised the severity of Muller's ratchet, while limiting the effectiveness of lateral gene transfer (LGT) at purging deleterious mutations. In principle, an increase in recombination length combined with higher rates of LGT could solve this problem. Here, we show using a computational model that this solution fails in the presence of genetic repeats prevalent in early eukaryotes. The model demonstrates that dispersed repeat sequences allow ectopic recombination, which leads to the loss of genetic information and curtails the capacity of LGT to prevent mutation accumulation. Increasing recombination length in the presence of repeat sequences exacerbates the problem. Mutational decay can only be resisted with homology along extended sequences of DNA. We conclude that the transition to homologous pairing along linear chromosomes was a key innovation in meiotic sex, which was instrumental in the expansion of eukaryotic genomes and morphological complexity.

Comparative genome analysis of Vagococcus fluvialis reveals abundance of mobile genetic elements in sponge-isolated strains

BACKGROUND

Vagococcus fluvialis is a species of lactic acid bacteria found both free-living in river and seawater and associated to hosts, such as marine sponges. This species has been greatly understudied, with no complete genome assembly available to date, which is essential for the characterisation of the mobilome.

RESULTS

We sequenced and assembled de novo the complete genome sequences of five V. fluvialis isolates recovered from marine sponges. Pangenome analysis of the V. fluvialis species (total of 17 genomes) showed a high intraspecific diversity, with 45.5% of orthologous genes found to be strain specific. Despite this diversity, analyses of gene functions clustered all V. fluvialis species together and separated them from other sequenced Vagococcus species. V. fluvialis strains from different habitats were highly similar in terms of functional diversity but the sponge-isolated strains were enriched in several functions related to the marine environment. Furthermore, sponge-isolated strains carried a significantly higher number of mobile genetic elements (MGEs) compared to previously sequenced V. fluvialis strains from other environments. Sponge-isolated strains carried up to 4 circular plasmids each, including a 48-kb conjugative plasmid. Three of the five strains carried an additional circular extrachromosomal sequence, assumed to be an excised prophage as it contained mainly viral genes and lacked plasmid replication genes. Insertion sequences (ISs) were up to five times more abundant in the genomes of sponge-isolated strains compared to the others, including several IS families found exclusively in these genomes.

CONCLUSIONS

Our findings highlight the dynamics and plasticity of the V. fluvialis genome. The abundance of mobile genetic elements in the genomes of sponge-isolated V. fluvialis strains suggests that the mobilome might be key to understanding the genomic signatures of symbiosis in bacteria.

Genomic features of Mycoplasma bovis subtypes currently circulating in France

Background

Mycoplasma (M.) bovis is a major etiological agent of bovine respiratory disease, which is the most economically costly disease of cattle worldwide. Cattle disease surveillance on M. bovis is increasingly using gene-based techniques, such as multilocus sequence typing (MLST), or genome-based techniques such as core genome MLST that both require only partial genomic data. However, accurate up-to-date surveillance also demands complete, circular genomes that can be used as reference to track the evolution of the different lineages. Yet, in France, two of the main subtypes currently circulating still have no representing genome in public databases. Here, to address this gap, we provide and compare three new complete M. bovis genomes obtained from recent clinical isolates that represent major subtypes circulating in France and Europe.

Results

Genomes were obtained using a hybrid assembly strategy (Illumina and Nanopore) with fine-tuning of settings and inputs used in the Unicycler assembly pipeline, such as size selection of reads and quality trimming of the FASTQ files. The main characteristics and synteny of the genomes were compared. The three genomes mainly differed by their content in terms of mobile genetic elements, i.e. integrative conjugative elements (ICE) and insertion sequences (IS), a feature that impacts their structure. For instance, strain L15527, representing subtype3 (st3), harbours an exceptionally high number of ICEs, which results in a bigger-sized genome than all those previously described and could be associated with the propensity of st3 to gain and fix mutations through chromosomal transfer mechanisms. In contrast, strain F9160, of st1, is very close to the PG45 type strain isolated in 1961 in the USA, and harbours a huge number of IS. These features may be associated with an evolution towards a host-restricted state or in a “closed” host or environment reservoir until a recent re-emergence.

Conclusions

Whole-genome comparison of the three French M. bovis subtypes provides valuable resources for future studies combining epidemiology, phylogenetic data, and phylodynamic methods.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08818-9.

Metagenomic methylation patterns resolve bacterial genomes of unusual size and structural complexity

The plasticity of bacterial and archaeal genomes makes examining their ecological and evolutionary dynamics both exciting and challenging. The same mechanisms that enable rapid genomic change and adaptation confound current approaches for recovering complete genomes from metagenomes. Here, we use strain-specific patterns of DNA methylation to resolve complex bacterial genomes from long-read metagenomic data of a marine microbial consortium, the “pink berries” of the Sippewissett Marsh (USA). Unique combinations of restriction-modification (RM) systems encoded by the bacteria produced distinctive methylation profiles that were used to accurately bin and classify metagenomic sequences. Using this approach, we finished the largest and most complex circularized bacterial genome ever recovered from a metagenome (7.9 Mb with >600 transposons), the finished genome of Thiohalocapsa sp. PB-PSB1 the dominant bacteria in the consortia. From genomes binned by methylation patterns, we identified instances of horizontal gene transfer between sulfur-cycling symbionts (Thiohalocapsa sp. PB-PSB1 and Desulfofustis sp. PB-SRB1), phage infection, and strain-level structural variation. We also linked the methylation patterns of each metagenome-assembled genome with encoded DNA methyltransferases and discovered new RM defense systems, including novel associations of RM systems with RNase toxins.

Repeat dynamics across timescales: a perspective from sibling allotetraploid marsh orchids (Dactylorhiza majalis s.l.)

To provide insights into the fate of transposable elements (TEs) across timescales in a post-polyploidization context, we comparatively investigate five sibling Dactylorhiza allotetraploids (Orchidaceae) formed independently and sequentially between 500 and 100K generations ago by unidirectional hybridization between diploids D. fuchsii and D. incarnata. Our results first reveal that the paternal D. incarnata genome shows a marked increased content of LTR retrotransposons compared to the maternal species, reflected in its larger genome size and consistent with a previously hypothesized bottleneck. With regard to the allopolyploids, in the youngest D. purpurella both genome size and TE composition appear to be largely additive with respect to parents, whereas for polyploids of intermediate ages we uncover rampant genome expansion on a magnitude of multiple entire genomes of some plants such as Arabidopsis. The oldest allopolyploids in the series are not larger than the intermediate ones. A putative tandem repeat, potentially derived from a non-autonomous miniature inverted-repeat TE (MITE) drives much of the genome dynamics in the allopolyploids. The highly dynamic MITE-like element is found in higher proportions in the maternal diploid, D. fuchsii, but is observed to increase in copy number in both subgenomes of the allopolyploids. Altogether, the fate of repeats appears strongly regulated and therefore predictable across multiple independent allopolyploidization events in this system. Apart from the MITE-like element, we consistently document a mild genomic shock following the allopolyploidizations investigated here, which may be linked to their relatively large genome sizes, possibly associated with strong selection against further genome expansions.

How Do Transposable Elements Activate Expression of Transcriptionally Silent Antibiotic Resistance Genes?

The rapidly emerging phenomenon of antibiotic resistance threatens to substantially reduce the efficacy of available antibacterial therapies. Dissemination of resistance, even between phylogenetically distant bacterial species, is mediated mainly by mobile genetic elements, considered to be natural vectors of horizontal gene transfer. Transposable elements (TEs) play a major role in this process—due to their highly recombinogenic nature they can mobilize adjacent genes and can introduce them into the pool of mobile DNA. Studies investigating this phenomenon usually focus on the genetic load of transposons and the molecular basis of their mobility. However, genes introduced into evolutionarily distant hosts are not necessarily expressed. As a result, bacterial genomes contain a reservoir of transcriptionally silent genetic information that can be activated by various transposon-related recombination events. The TEs themselves along with processes associated with their transposition can introduce promoters into random genomic locations. Thus, similarly to integrons, they have the potential to convert dormant genes into fully functional antibiotic resistance determinants. In this review, we describe the genetic basis of such events and by extension the mechanisms promoting the emergence of new drug-resistant bacterial strains.

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

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

Genetic Characteristics and Enzymatic Activities of Bacillus velezensis KS04AU as a Stable Biocontrol Agent against Phytopathogens

Bacillus velezensis has a broad application in the agricultural and industrial sectors for its biocontrol properties and its potential active secondary metabolites. The defined phenotypic characteristics of a strain vary according to its ecosystem. We report the complete genomic analysis of B. velezensis KS04AU compared to four strains of B. velezensis (SRCM102752, ONU-553, FZB42, and JS25R) and two closely related Bacillus amyloliquefaciens (LL3 and IT-45). A total of 4771 protein coding genes comprises the KS04AU genome, in comparison with 3334 genes core genes found in the six other strains and the remaining 1437 shell genes. Average nucleotide identity of the target strain to the six other strains showed 99.65% to B. velezensis ONU-553, sharing 60 orthologous genes. Secondary metabolite gene cluster analysis of all strains showed that KS04AU has a mersacidin cluster gene, which is absent in the genome of the other strains. PHASTER analysis also showed KS04AU harboring two phages (Aeribacllus AP45 NC_048651 and Paenibacillus_Tripp NC_028930), which were also unique in comparison with the other strains. Analysis on anti-microbial resistance genes showed no difference in the genome of KS04AU to any of the other genomes, with the exception of B. amyloliquefaciens IT-45 which had one unique small multidrug-resistance antibiotic efflux-pump gene (qacJ). The CRISPR-Cas systems in the strains were also compared showing one CRISPR gene found only in KS04AU. Hydrolytic activity, antagonistic activity against phytopathogens (Fusarium oxysporum, Fusarium graminearum, Alternaria alternata and Pseudomonas syringae) and biocontrol against tomato foot and root rot experiments were carried out. B. velezensis KS04AU inhibits the growth of all phytopathogens tested, produces hydrolytic activity, and reduces Fusarium oxysporum f.sp. radicis-lycopersici (Forl) ZUM2407 lesions up to 46.02 ± 0.12%. The obtained results confirm B. velezensis KS04AU as a potential biocontrol strain for plant protection.

From a large-scale genomic analysis of insertion sequences to insights into their regulatory roles in prokaryotes

Background

Insertion sequences (ISs) are mobile repeat sequences and most of them can copy themselves to new host genome locations, leading to genome plasticity and gene regulation in prokaryotes. In this study, we present functional and evolutionary relationships between IS and neighboring genes in a large-scale comparative genomic analysis.

Results

IS families were located in all prokaryotic phyla, with preferential occurrence of IS3, IS4, IS481, and IS5 families in Alpha-, Beta-, and Gammaproteobacteria, Actinobacteria and Firmicutes as well as in eukaryote host-associated organisms and autotrophic opportunistic pathogens. We defined the concept of the IS-Gene couple (IG), which allowed to highlight the functional and regulatory impacts of an IS on the closest gene. Genes involved in transcriptional regulation and transport activities were found overrepresented in IG. In particular, major facilitator superfamily (MFS) transporters, ATP-binding proteins and transposases raised as favorite neighboring gene functions of IS hotspots. Then, evolutionary conserved IS-Gene sets across taxonomic lineages enabled the classification of IS-gene couples into phylum, class-to-genus, and species syntenic IS-Gene couples. The IS5, IS21, IS4, IS607, IS91, ISL3 and IS200 families displayed two to four times more ISs in the phylum and/or class-to-genus syntenic IGs compared to other IS families. This indicates that those families were probably inserted earlier than others and then subjected to horizontal transfer, transposition and deletion events over time. In phylum syntenic IG category, Betaproteobacteria, Crenarchaeota, Calditrichae, Planctomycetes, Acidithiobacillia and Cyanobacteria phyla act as IS reservoirs for other phyla, and neighboring gene functions are mostly related to transcriptional regulators. Comparison of IS occurrences with predicted regulatory motifs led to ~ 26.5% of motif-containing ISs with 2 motifs per IS in average. These results, concomitantly with short IS-Gene distances, suggest that those ISs would interfere with the expression of neighboring genes and thus form strong candidates for an adaptive pairing.

Conclusions

All together, our large-scale study provide new insights into the IS genetic context and strongly suggest their regulatory roles.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08678-3.

Intragenic codon usage in proteobacteria: Translational selection, IS expansion and genomic shrinkage

This manuscript presents a method to systematically study intragenic variations in codon usage using correspondence analysis and the effective number of codons. The method is applied to >1100 proteobacteria. Codon usage biases (measured as inertia) increases with genome size, the same is true for the percentage of inertia explained by the first axis. It is shown that there is often a relaxed or more uniform codon usage near the gene termini. Ithis is not seen n small genomes, notably those of intracellular organisms like Buchnera aphidicola or Rickettsia prowazekii where translational selection plays less of a role. When genes from E. coli, for which translational selection is well described, are split into low, intermediate and high expression, respectively, it is shown that the intragenic codon usage pattern with more uniform usage at termini exist across all three expression groups. Furthermore, the correspondence analysis reveals a unique pattern in Bordetella pertussis due to IS expansion. This study thus shows that translational selection, genome shrinkage and IS expansion result in characteristic patterns in intragenic codon usage.

inPOSE: A Flexible Toolbox for Chromosomal Cloning and Amplification of Bacterial Transgenes

Cloning the genes and operons encoding heterologous functions in bacterial hosts is now almost exclusively carried out using plasmid vectors. This has multiple drawbacks, including the need for constant selection and variation in copy numbers. The chromosomal integration of transgenes has always offered a viable alternative; however, to date, it has been of limited use due to its tedious nature and often being limited to a single copy. We introduce here a strategy that uses bacterial insertion sequences, which are the simplest autonomous transposable elements to insert and amplify genetic cargo into a bacterial chromosome. Transgene insertion can take place either as transposition or homologous recombination, and copy number amplification is achieved using controlled copy-paste transposition. We display the successful use of IS1 and IS3 for this purpose in Escherichia coli cells using various selection markers. We demonstrate the insertion of selectable genes, an unselectable gene and a five-gene operon in up to two copies in a single step. We continue with the amplification of the inserted cassette to double-digit copy numbers within two rounds of transposase induction and selection. Finally, we analyze the stability of the cloned genetic constructs in the lack of selection and find it to be superior to all investigated plasmid-based systems. Due to the ubiquitous nature of transposable elements, we believe that with proper design, this strategy can be adapted to numerous other bacterial species.

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

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

Transposable elements promote the evolution of genome streamlining

Eukaryotes and prokaryotes have distinct genome architectures, with marked differences in genome size, the ratio of coding/non-coding DNA, and the abundance of transposable elements (TEs). As TEs replicate independently of their hosts, the proliferation of TEs is thought to have driven genome expansion in eukaryotes. However, prokaryotes also have TEs in intergenic spaces, so why do prokaryotes have small, streamlined genomes? Using an in silico model describing the genomes of single-celled asexual organisms that coevolve with TEs, we show that TEs acquired from the environment by horizontal gene transfer can promote the evolution of genome streamlining. The process depends on local interactions and is underpinned by rock-paper-scissors dynamics in which populations of cells with streamlined genomes beat TEs, which beat non-streamlined genomes, which beat streamlined genomes, in continuous and repeating cycles. Streamlining is maladaptive to individual cells, but improves lineage viability by hindering the proliferation of TEs. Streamlining does not evolve in sexually reproducing populations because recombination partially frees TEs from the deleterious effects they cause. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Microbial defenses against mobile genetic elements and viruses: Who defends whom from what?

Prokaryotes have numerous mobile genetic elements (MGEs) that mediate horizontal gene transfer (HGT) between cells. These elements can be costly, even deadly, and cells use numerous defense systems to filter, control, or inactivate them. Recent studies have shown that prophages, conjugative elements, their parasites (phage satellites and mobilizable elements), and other poorly described MGEs encode defense systems homologous to those of bacteria. These constitute a significant fraction of the repertoire of cellular defense genes. As components of MGEs, these defense systems have presumably evolved to provide them, not the cell, adaptive functions. While the interests of the host and MGEs are aligned when they face a common threat such as an infection by a virulent phage, defensive functions carried by MGEs might also play more selfish roles to fend off other antagonistic MGEs or to ensure their maintenance in the cell. MGEs are eventually lost from the surviving host genomes by mutational processes and their defense systems can be co-opted when they provide an advantage to the cell. The abundance of defense systems in MGEs thus sheds new light on the role, effect, and fate of the so-called "cellular defense systems," whereby they are not only merely microbial defensive weapons in a 2-partner arms race, but also tools of intragenomic conflict between multiple genetic elements with divergent interests that shape cell fate and gene flow at the population level.

Lanza V, Rodríguez-Beltrán J, Galán JC, San Millán A, Cantón R, Coque TM. Evolutionary Pathways and Trajectories in Antibiotic Resistance

Evolution is the hallmark of life. Descriptions of the evolution of microorganisms have provided a wealth of information, but knowledge regarding "what happened" has precluded a deeper understanding of "how" evolution has proceeded, as in the case of antimicrobial resistance. The difficulty in answering the "how" question lies in the multihierarchical dimensions of evolutionary processes, nested in complex networks, encompassing all units of selection, from genes to communities and ecosystems. At the simplest ontological level (as resistance genes), evolution proceeds by random (mutation and drift) and directional (natural selection) processes; however, sequential pathways of adaptive variation can occasionally be observed, and under fixed circumstances (particular fitness landscapes), evolution is predictable. At the highest level (such as that of plasmids, clones, species, microbiotas), the systems' degrees of freedom increase dramatically, related to the variable dispersal, fragmentation, relatedness, or coalescence of bacterial populations, depending on heterogeneous and changing niches and selective gradients in complex environments. Evolutionary trajectories of antibiotic resistance find their way in these changing landscapes subjected to random variations, becoming highly entropic and therefore unpredictable. However, experimental, phylogenetic, and ecogenetic analyses reveal preferential frequented paths (highways) where antibiotic resistance flows and propagates, allowing some understanding of evolutionary dynamics, modeling and designing interventions. Studies on antibiotic resistance have an applied aspect in improving individual health, One Health, and Global Health, as well as an academic value for understanding evolution. Most importantly, they have a heuristic significance as a model to reduce the negative influence of anthropogenic effects on the environment.

Frequency, composition and mobility of Escherichia coli-derived transposable elements in holdings of plasmid repositories

By providing the scientific community with uniform and standardized resources of consistent quality, plasmid repositories play an important role in enabling scientific reproducibility. Plasmids containing insertion sequence elements (IS elements) represent a challenge from this perspective, as they can change the plasmid structure and function. In this study, we conducted a systematic analysis of a subset of plasmid stocks distributed by plasmid repositories (The Arabidopsis Biological Resource Center and Addgene) which carry unintended integrations of bacterial mobile genetic elements. The integration of insertion sequences was most often found in, but not limited to, pBR322‐derived vectors, and did not affect the function of the specific plasmids. In certain cases, the entire stock was affected, but the majority of the stocks tested contained a mixture of the wild‐type and the mutated plasmids, suggesting that the acquisition of IS elements likely occurred after the plasmids were acquired by the repositories. However, comparison of the sequencing results of the original samples revealed that some plasmids already carried insertion mutations at the time of donation. While an extensive BLAST analysis of 47 877 plasmids sequenced from the Addgene repository uncovered IS elements in only 1.12%, suggesting that IS contamination is not widespread, further tests showed that plasmid integration of IS elements can propagate in conventional Escherichia coli hosts over a few tens of generations. Use of IS‐free E. coli hosts prevented the emergence of IS insertions as well as that of small indels, suggesting that the use of IS‐free hosts by donors and repositories could help limit unexpected and unwanted IS integrations into plasmids.

Insights on mauritiana-like Elements Diversity in Mayetiola destructor and M. hordei (Diptera: Cecidomyiidae)

Mariner-like elements (MLEs) are class II transposons belonging to the Tc1-mariner family, that have successfully invaded many insect genomes. In the current study, the availability of the Hessian fly Mayetiola destructor genome has enabled us to perform in silico analysis of MLEs using as query the previously described mariner element (Desmar1) belonging to mauritiana subfamily. Eighteen mauritiana-like elements were detected and were clustered into three main groups named Desmar1-like, MauCons1 and MauCons2. Subsequently, in vitro analysis was carried out to investigate mauritiana-like elements in M. destructor as well as in Mayetiola hordei using primers designed from TIRs of the previously identified MLEs. PCR amplifications were successful and a total of 12 and 17 mauritiana-like elements were discovered in M. destructor and M. hordei, respectively. Sequence analyses of mauritiana-like elements obtained in silico and in vitro have showed that MauCons1 and MauCons2 elements share low similarity with Desmar1 ranging from 50% to 55% suggesting different groups under mauritiana subfamily have invaded the genomes of M. destructor and M. hordei. These groups are likely inherited by vertical transmission that subsequently underwent different evolutionary histories. This work describes new mauritiana-like elements in M. destructor that are distinct from the previouslydiscovered Desmar1 and provides the first evidence of MLEs belonging to mauritiana subfamily in M. hordei.

Influence of Redox Imbalances on the Transposition of Insertion Sequences in Deinococcus geothermalis

The transposition of insertion sequence elements was evaluated among different Deinococcus geothermalis lineages, including the wild-type, a cystine importer-disrupted mutant, a complemented strain, and a cystine importer-overexpressed strain. Cellular growth reached early exponential growth at OD₆₀₀ 2.0 and late exponential growth at OD₆₀₀ 4.0. Exposing the cells to hydrogen peroxide (80–100 mM) resulted in the transposition of insertion sequences (ISs) in genes associated with the carotenoid biosynthesis pathway. Particularly, ISDge7 (an IS5 family member) and ISDge5 (an IS701 family member) from the cystine importer-disrupted mutant were transposed into phytoene desaturase (dgeo_0524) via replicative transposition. Further, the cystine importer-overexpressed strain Δdgeo_1985R showed transposition of both ISDge2 and ISDge5 elements. In contrast, IS transposition was not detected in the complementary strain. Interestingly, a cystine importer-overexpressing strain exhibited streptomycin resistance, indicating that point mutation occurred in the rpsL (dgeo_1873) gene encoding ribosomal protein S12. qRT-PCR analyses were then conducted to evaluate the expression of oxidative stress response genes, IS elements, and low-molecular-weight thiol compounds such as mycothiol and bacillithiol. Nevertheless, the mechanisms that trigger IS transposition in redox imbalance conditions remain unclear. Here, we report that the active transposition of different IS elements was affected by intracellular redox imbalances caused by cystine importer deficiencies or overexpression.

Prokaryotic and Eukaryotic Horizontal Transfer of Sailor (DD82E), a New Superfamily of IS630-Tc1-Mariner DNA Transposons

Simple Summary

Transposable elements, including DNA transposons, play a significant role in genetic material exchanges between prokaryotes and eukaryotes. Comparative profiling of the evolution pattern of DNA transposons between prokaryotes and eukaryotes may identify potential genetic material exchanges between them and provide insights into the evolutionary history of prokaryotic and eukaryotic genomes. The members of the IS630-Tc1-mariner (ITm) group may represent the most diverse and widely distributed DNA transposons in nature, and the discovery of new members of this group is highly expected based on the increasing availability of genome sequencing data. We discovered a new superfamily (termed Sailor) belonging to the ITm hyperfamily, which differed from the known superfamilies of Tc1/mariner, DDxD/pogo and DD34E/Gambol, regarding phylogenetic position and catalytic domain. Our data revealed that Sailor was distributed in both prokaryotes and eukaryotes and suggested that horizontal transfer (HT) events of Sailor may occur from prokaryotic to eukaryotic genomes. Finally, internal transmissions of Sailor in prokaryotes and eukaryotes were also detected.

Abstract

Here, a new superfamily of IS630-Tc1-mariner (ITm) DNA transposons, termed Sailor, is identified, that is characterized by a DD82E catalytic domain and is distinct from all previously known superfamilies of the ITm group. Phylogenetic analyses revealed that Sailor forms a monophyletic clade with a more intimate link to the clades of Tc1/mariner and DD34E/Gambol. Sailor was detected in both prokaryotes and eukaryotes and invaded a total of 256 species across six kingdoms. Sailor is present in nine species of bacteria, two species of plantae, four species of protozoa, 23 species of Chromista, 12 species of Fungi and 206 species of animals. Moreover, Sailor is extensively distributed in invertebrates (a total of 206 species from six phyla) but is absent in vertebrates. Sailor transposons are 1.38–6.98 kb in total length and encoded transposases of ~676 aa flanked by TIRs with lengths between 18, 1362 and 4 bp (TATA) target-site duplications. Furthermore, our analysis provided strong evidence of Sailor transmissions from prokaryotes to eukaryotes and internal transmissions in both. These data update the classification of the ITm group and will contribute to the understanding of the evolution of ITm transposons and that of their hosts.

Genome Plasticity by Insertion Sequences Learned From a Case of Radiation-Resistant Bacterium Deinococcus geothermalis

The genome of the radiation-resistant bacterium Deinococcus geothermalis contains 19 types of insertion sequences (ISs), including 93 total transposases (Tpases) in 73 full-length ISs from the main chromosome and 2 mega plasmids. In this study, 68 ISs from the D. geothermalis genome were extracted to implicate the earlier genome before its mutation by transposition of ISs. The total size of eliminated ISs from genome was 78.85 kb. From these in silico corrections of mutation by the ISs, we have become aware of some bioinformatics factualness as follows: (1) can reassemble the disrupted genes if the exact IS region was eliminated, (2) can configure the schematic clustering of major DDE type Tpases, (3) can determine IS integration order across multiple hot spots, and (4) can compare genetic relativeness by the partial synteny analysis between D. geothermalis and Deinococcus strain S9. Recently, we found that several IS elements actively transferred to other genomic sites under hydrogen peroxide-induced oxidative stress conditions, resulting in the inactivation of functional genes. Therefore, the single species genome’s mobilome study provides significant support to define bacterial genome plasticity and molecular evolution from past and present progressive transposition events.

The Evolution of Interdependence in a Four-Way Mealybug Symbiosis

Mealybugs are insects that maintain intracellular bacterial symbionts to supplement their nutrient-poor plant sap diets. Some mealybugs have a single betaproteobacterial endosymbiont, a Candidatus Tremblaya species (hereafter Tremblaya) that alone provides the insect with its required nutrients. Other mealybugs have two nutritional endosymbionts that together provision these same nutrients, where Tremblaya has gained a gammaproteobacterial partner that resides in its cytoplasm. Previous work had established that Pseudococcus longispinus mealybugs maintain not one but two species of gammaproteobacterial endosymbionts along with Tremblaya. Preliminary genomic analyses suggested that these two gammaproteobacterial endosymbionts have large genomes with features consistent with a relatively recent origin as insect endosymbionts, but the patterns of genomic complementarity between members of the symbiosis and their relative cellular locations were unknown. Here, using long-read sequencing and various types of microscopy, we show that the two gammaproteobacterial symbionts of P. longispinus are mixed together within Tremblaya cells, and that their genomes are somewhat reduced in size compared with their closest nonendosymbiotic relatives. Both gammaproteobacterial genomes contain thousands of pseudogenes, consistent with a relatively recent shift from a free-living to an endosymbiotic lifestyle. Biosynthetic pathways of key metabolites are partitioned in complex interdependent patterns among the two gammaproteobacterial genomes, the Tremblaya genome, and horizontally acquired bacterial genes that are encoded on the mealybug nuclear genome. Although these two gammaproteobacterial endosymbionts have been acquired recently in evolutionary time, they have already evolved codependencies with each other, Tremblaya, and their insect host.

How sequence populations persist inside bacterial genomes

Compared to their eukaryotic counterparts, bacterial genomes are small and contain extremely tightly packed genes. Repetitive sequences are rare but not completely absent. One of the most common repeat families is REPINs. REPINs can replicate in the host genome and form populations that persist for millions of years. Here, we model the interactions of these intragenomic sequence populations with the bacterial host. We first confirm well-established results, in the presence and absence of horizontal gene transfer (hgt) sequence populations either expand until they drive the host to extinction or the sequence population gets purged from the genome. We then show that a sequence population can be stably maintained, when each individual sequence provides a benefit that decreases with increasing sequence population size. Maintaining a sequence population of stable size also requires the replication of the sequence population to be costly to the host, otherwise the sequence population size will increase indefinitely. Surprisingly, in regimes with high hgt rates, the benefit conferred by the sequence population does not have to exceed the damage it causes to its host. Our analyses provide a plausible scenario for the persistence of sequence populations in bacterial genomes. We also hypothesize a limited biologically relevant parameter range for the provided benefit, which can be tested in future experiments.

Traveler, a New DD35E Family of Tc1/Mariner Transposons, Invaded Vertebrates Very Recently

The discovery of new members of the Tc1/mariner superfamily of transposons is expected based on the increasing availability of genome sequencing data. Here, we identified a new DD35E family termed Traveler (TR). Phylogenetic analyses of its DDE domain and full-length transposase showed that, although TR formed a monophyletic clade, it exhibited the highest sequence identity and closest phylogenetic relationship with DD34E/Tc1. This family displayed a very restricted taxonomic distribution in the animal kingdom and was only detected in ray-finned fish, anura, and squamata, including 91 vertebrate species. The structural organization of TRs was highly conserved across different classes of animals. Most intact TR transposons had a length of ∼1.5 kb (range 1,072-2,191 bp) and harbored a single open reading frame encoding a transposase of ∼340 aa (range 304-350 aa) flanked by two short-terminal inverted repeats (13-68 bp). Several conserved motifs, including two helix-turn-helix motifs, a GRPR motif, a nuclear localization sequence, and a DDE domain, were also identified in TR transposases. This study also demonstrated the presence of horizontal transfer events of TRs in vertebrates, whereas the average sequence identities and the evolutionary dynamics of TR elements across species and clusters strongly indicated that the TR family invaded the vertebrate lineage very recently and that some of these elements may be currently active, combining the intact TR copies in multiple lineages of vertebrates. These data will contribute to the understanding of the evolutionary history of Tc1/mariner transposons and that of their hosts.

Insertion-sequence-mediated mutations both promote and constrain evolvability during a long-term experiment with bacteria

Insertion sequences (IS) are ubiquitous bacterial mobile genetic elements, and the mutations they cause can be deleterious, neutral, or beneficial. The long-term dynamics of IS elements and their effects on bacteria are poorly understood, including whether they are primarily genomic parasites or important drivers of adaptation by natural selection. Here, we investigate the dynamics of IS elements and their contribution to genomic evolution and fitness during a long-term experiment with Escherichia coli. IS elements account for ~35% of the mutations that reached high frequency through 50,000 generations in those populations that retained the ancestral point-mutation rate. In mutator populations, IS-mediated mutations are only half as frequent in absolute numbers. In one population, an exceptionally high ~8-fold increase in IS150 copy number is associated with the beneficial effects of early insertion mutations; however, this expansion later slowed down owing to reduced IS150 activity. This population also achieves the lowest fitness, suggesting that some avenues for further adaptation are precluded by the IS150-mediated mutations. More generally, across all populations, we find that higher IS activity becomes detrimental to adaptation over evolutionary time. Therefore, IS-mediated mutations can both promote and constrain evolvability.

Insertion sequences (IS) are common mobile genetic elements in bacteria, but their effects on bacterial evolution are not well understood. Here, Consuegra and colleagues investigate the dynamics and fitness consequences of IS elements in E. coli over 50,000 generations.

Assessing the Genomic Variability of Gardnerella vaginalis through Comparative Genomic Analyses: Evolutionary and Ecological Implications

Gardnerella vaginalis is described as a common anaerobic vaginal bacterium whose presence may correlate with vaginal dysbiotic conditions. In the current study, we performed phylogenomic analyses of 72 G. vaginalis genome sequences, revealing noteworthy genome differences underlying a polyphyletic organization of this taxon. Particularly, the genomic survey revealed that this species may actually include nine distinct genotypes (GGtype1 to GGtype9). Furthermore, the observed link between sialidase and phylogenomic grouping provided clues of a connection between virulence potential and the evolutionary history of this microbial taxon. Specifically, based on the outcomes of these in silico analyses, GGtype3, GGtype7, GGtype8, and GGtype9 appear to have virulence potential since they exhibited the sialidase gene in their genomes. Notably, the analysis of 34 publicly available vaginal metagenomic samples allowed us to trace the distribution of the nine G. vaginalis genotypes identified in this study among the human population, highlighting how differences in genetic makeup could be related to specific ecological properties. Furthermore, comparative genomic analyses provided details about the G. vaginalis pan- and core genome contents, including putative genetic elements involved in the adaptation to the ecological niche as well as many putative virulence factors. Among these putative virulence factors, particularly noteworthy genes identified were the gene encoding cholesterol-dependent cytolysin (CDC) toxin vaginolysin and genes related to microbial biofilm formation, iron uptake, adhesion to the vaginal epithelium, as well as macrolide antibiotic resistance.IMPORTANCE The identification of nine different genotypes among members of G. vaginalis allowed us to distinguish an uneven distribution of virulence-associated genetic traits within this taxon and thus suggest the potential occurrence of putative pathogen and commensal G. vaginalis strains. These findings, coupled with metagenomics microbial profiling of human vaginal microbiota, permitted us to get insights into the distribution of the genotypes among the human population, highlighting the presence of different structural communities in terms of G. vaginalis genotypes.

Comparative Pathogenomics of Aeromonas veronii from Pigs in South Africa: Dominance of the Novel ST657 Clone

The pathogenomics of carbapenem-resistant Aeromonas veronii (A. veronii) isolates recovered from pigs in KwaZulu-Natal, South Africa, was explored by whole genome sequencing on the Illumina MiSeq platform. Genomic functional annotation revealed a vast array of similar central networks (metabolic, cellular, and biochemical). The pan-genome analysis showed that the isolates formed a total of 4349 orthologous gene clusters, 4296 of which were shared; no unique clusters were observed. All the isolates had similar resistance phenotypes, which corroborated their chromosomally mediated resistome (bla_CPHA3 and bla_OXA-12) and belonged to a novel sequence type, ST657 (a satellite clone). Isolates in the same sub-clades clustered according to their clonal lineages and host. Mobilome analysis revealed the presence of chromosome-borne insertion sequence families. The estimated pathogenicity score (P_score ≈ 0.60) indicated their potential pathogenicity in humans. Furthermore, these isolates carried several virulence factors (adherence factors, toxins, and immune evasion), in different permutations and combinations, indicating a differential ability to establish infection. Phylogenomic and metadata analyses revealed a predilection for water environments and aquatic animals, with more recent reports in humans and food animals across geographies, making A. veronii a potential One Health indicator bacterium.

Molecular analyses of the gill symbiosis of the bathymodiolin mussel Gigantidas platifrons

Although the deep-sea bathymodiolin mussels have been intensively studied as a model of animal-bacteria symbiosis, it remains challenging to assess the host-symbiont interactions due to the complexity of the symbiotic tissue-the gill. Using cold-seep mussel Gigantidas platifrons as a model, we isolated the symbiont harboring bacteriocytes and profiled the transcriptomes of the three major parts of the symbiosis-the gill, the bacteriocyte, and the symbiont. This breakdown of the complex symbiotic tissue allowed us to characterize the host-symbiont interactions further. Our data showed that the gill's non-symbiotic parts play crucial roles in maintaining and protecting the symbiosis; the bacteriocytes supply the symbiont with metabolites, control symbiont population, and shelter the symbiont from phage infection; the symbiont dedicates to the methane oxidation and energy production. This study demonstrates that the bathymodiolin symbiosis interacts at the tissue, cellular, and molecular level, maintaining high efficiency and harmonic chemosynthetic micro niche.

Phylogenetic background and habitat drive the genetic diversification of Escherichia coli

Escherichia coli is mostly a commensal of birds and mammals, including humans, where it can act as an opportunistic pathogen. It is also found in water and sediments. We investigated the phylogeny, genetic diversification, and habitat-association of 1,294 isolates representative of the phylogenetic diversity of more than 5,000 isolates from the Australian continent. Since many previous studies focused on clinical isolates, we investigated mostly other isolates originating from humans, poultry, wild animals and water. These strains represent the species genetic diversity and reveal widespread associations between phylogroups and isolation sources. The analysis of strains from the same sequence types revealed very rapid change of gene repertoires in the very early stages of divergence, driven by the acquisition of many different types of mobile genetic elements. These elements also lead to rapid variations in genome size, even if few of their genes rise to high frequency in the species. Variations in genome size are associated with phylogroup and isolation sources, but the latter determine the number of MGEs, a marker of recent transfer, suggesting that gene flow reinforces the association of certain genetic backgrounds with specific habitats. After a while, the divergence of gene repertoires becomes linear with phylogenetic distance, presumably reflecting the continuous turnover of mobile element and the occasional acquisition of adaptive genes. Surprisingly, the phylogroups with smallest genomes have the highest rates of gene repertoire diversification and fewer but more diverse mobile genetic elements. This suggests that smaller genomes are associated with higher, not lower, turnover of genetic information. Many of these genomes are from freshwater isolates and have peculiar traits, including a specific capsule, suggesting adaptation to this environment. Altogether, these data contribute to explain why epidemiological clones tend to emerge from specific phylogenetic groups in the presence of pervasive horizontal gene transfer across the species.

Previous large scale studies on the evolution of E. coli focused on clinical isolates emphasizing virulence and antibiotic resistance in medically important lineages. Yet, most E. coli strains are either human commensals or not associated with humans at all. Here, we analyzed a large collection of non-clinical isolates of the species to assess the mechanisms of gene repertoire diversification in the light of isolation sources and phylogeny. We show that gene repertoires evolve so rapidly by the high turnover of mobile genetic elements that epidemiologically indistinguishable strains can be phenotypically extremely heterogeneous, illustrating the velocity of bacterial adaptation and the importance of accounting for the information on the whole genome at the epidemiological scale. Phylogeny and habitat shape the genetic diversification of E. coli to similar extents. Surprisingly, freshwater strains seem specifically adapted to this environment, breaking the paradigm that E. coli environmental isolates are systematically fecal contaminations. As a consequence, the evolution of this species is also shaped by environmental habitats, and it may diversify by acquiring genes and mobile elements from environmental bacteria (and not just from gut bacteria). This may facilitate the acquisition of virulence factors and antibiotic resistance in the strains that become pathogenic.

Early fate of exogenous promoters in E. coli

Gene gain by horizontal gene transfer is a major pathway of genome innovation in bacteria. The current view posits that acquired genes initially need to be silenced and that a bacterial chromatin protein, H-NS, plays a role in this silencing. However, we lack direct observation of the early fate of a horizontally transferred gene to prove this theory. We combine sequencing, flow cytometry and sorting, followed by microscopy to monitor gene expression and its variability after large-scale random insertions of a reporter gene in a population of Escherichia coli bacteria. We find that inserted promoters have a wide range of gene-expression variability related to their location. We find that high-expression clones carry insertions that are not correlated with H-NS binding. Conversely, binding of H-NS correlates with silencing. Finally, while most promoters show a common level of extrinsic noise, some insertions show higher noise levels. Analysis of these high-noise clones supports a scenario of switching due to transcriptional interference from divergent ribosomal promoters. Altogether, our findings point to evolutionary pathways where newly-acquired genes are not necessarily silenced, but may immediately explore a wide range of expression levels to probe the optimal ones.

Comparative genomics of Aeromonas veronii: Identification of a pathotype impacting aquaculture globally

Aeromonas veronii is a gram-negative species abundant in aquatic environments that causes disease in humans as well as terrestrial and aquatic animals. In the current study, 41 publicly available A. veronii genomes were compared to investigate distribution of putative virulence genes, global dissemination of pathotypes, and potential mechanisms of virulence. The complete genome of A. veronii strain ML09-123 from an outbreak of motile aeromonas septicemia in farm-raised catfish in the southeastern United States was included. Dissemination of A. veronii strain types was discovered in dispersed geographical locations. Isolate ML09-123 is highly similar to Chinese isolate TH0426, suggesting the two strains have a common origin and may represent a pathotype impacting aquaculture in both countries. Virulence of strain ML09-123 in catfish in a dose-dependent manner was confirmed experimentally. Subsystem category disposition showed the majority of genomes exhibit similar distribution of genomic elements. The type I secretion system (T1SS), type II secretion system (T2SS), type 4 pilus (T4P), and flagellum core elements are conserved in all A. veronii genomes, whereas the type III secretion system (T3SS), type V secretion system (T5SS), type VI secretion system (T6SS), and tight adherence (TAD) system demonstrate variable dispersal. Distribution of mobile elements is dependent on host and geographic origin, suggesting this species has undergone considerable genetic exchange. The data presented here lends insight into the genomic variation of A. veronii and identifies a pathotype impacting aquaculture globally.

Synthetic Genome Defenses against Selfish DNA Elements Stabilize Engineered Bacteria against Evolutionary Failure

Mobile genetic elements drive evolution by disrupting genes and rearranging genomes. Eukaryotes have evolved epigenetic mechanisms, including DNA methylation and RNA interference, that silence mobile elements and thereby preserve the integrity of their genomes. We created an artificial reprogrammable epigenetic system based on CRISPR interference to give engineered bacteria a similar line of defense against transposons and other selfish elements in their genomes. We demonstrate that this CRISPR interference against mobile elements (CRISPRi-ME) approach can be used to simultaneously repress two different transposon families in Escherichia coli, thereby increasing the evolutionary stability of costly protein expression. We further show that silencing a transposon in Acinetobacter baylyi ADP1 reduces mutation rates by a factor of 5, nearly as much as deleting all copies of this element from its genome. By deploying CRISPRi-ME on a broad-host-range vector, we have created a generalizable platform for stabilizing the genomes of engineered bacterial cells for applications in metabolic engineering and synthetic biology.

Genomic and physiological characterization and description of Marinobacter gelidimuriae sp. nov., a psychrophilic, moderate halophile from Blood Falls, an antarctic subglacial brine

Antarctic subice environments are diverse, underexplored microbial habitats. Here, we describe the ecophysiology and annotated genome of a Marinobacter strain isolated from a cold, saline, iron-rich subglacial outflow of the Taylor Glacier, Antarctica. This strain (BF04_CF4) grows fastest at neutral pH (range 6-10), is psychrophilic (range: 0°C-20°C), moderately halophilic (range: 0.8%-15% NaCl) and hosts genes encoding potential low temperature and high salt adaptations. The predicted proteome suggests it utilizes fewer charged amino acids than a mesophilic Marinobacter strain. BF04_CF4 has increased concentrations of membrane unsaturated fatty acids including palmitoleic (33%) and oleic (27.5%) acids that may help maintain cell membrane fluidity at low temperatures. The genome encodes proteins for compatible solute biosynthesis and transport, which are known to be important for growth in saline environments. Physiological verification of predicted metabolic functions demonstrate BF04_CF4 is capable of denitrification and may facilitate iron oxidation. Our data indicate that strain BF04_CF4 represents a new Marinobacter species, Marinobacter gelidimuriae sp. nov., that appears well suited for the subglacial environment it was isolated from. Marinobacter species have been isolated from other cold, saline environments in the McMurdo Dry Valleys and permanently cold environments globally suggesting that this lineage is cosmopolitan and ecologically relevant in icy brines.

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.

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.

Exploration of Survival Traits, Probiotic Determinants, Host Interactions, and Functional Evolution of Bifidobacterial Genomes Using Comparative Genomics

Members of the genus Bifidobacterium are found in a wide-range of habitats and are used as important probiotics. Thus, exploration of their functional traits at the genus level is of utmost significance. Besides, this genus has been demonstrated to exhibit an open pan-genome based on the limited number of genomes used in earlier studies. However, the number of genomes is a crucial factor for pan-genome calculations. We have analyzed the pan-genome of a comparatively larger dataset of 215 members of the genus Bifidobacterium belonging to different habitats, which revealed an open nature. The pan-genome for the 56 probiotic and human-gut strains of this genus, was also found to be open. The accessory- and unique-components of this pan-genome were found to be under the operation of Darwinian selection pressure. Further, their genome-size variation was predicted to be attributed to the abundance of certain functions carried by genomic islands, which are facilitated by insertion elements and prophages. In silico functional and host-microbe interaction analyses of their core-genome revealed significant genomic factors for niche-specific adaptations and probiotic traits. The core survival traits include stress tolerance, biofilm formation, nutrient transport, and Sec-secretion system, whereas the core probiotic traits are imparted by the factors involved in carbohydrate- and protein-metabolism and host-immunomodulations.

Genomes of ubiquitous marine and hypersaline Hydrogenovibrio, Thiomicrorhabdus and Thiomicrospira spp. encode a diversity of mechanisms to sustain chemolithoautotrophy in heterogeneous environments

Chemolithoautotrophic bacteria from the genera Hydrogenovibrio, Thiomicrorhabdus and Thiomicrospira are common, sometimes dominant, isolates from sulfidic habitats including hydrothermal vents, soda and salt lakes and marine sediments. Their genome sequences confirm their membership in a deeply branching clade of the Gammaproteobacteria. Several adaptations to heterogeneous habitats are apparent. Their genomes include large numbers of genes for sensing and responding to their environment (EAL- and GGDEF-domain proteins and methyl-accepting chemotaxis proteins) despite their small sizes (2.1-3.1 Mbp). An array of sulfur-oxidizing complexes are encoded, likely to facilitate these organisms' use of multiple forms of reduced sulfur as electron donors. Hydrogenase genes are present in some taxa, including group 1d and 2b hydrogenases in Hydrogenovibrio marinus and H. thermophilus MA2-6, acquired via horizontal gene transfer. In addition to high-affinity cbb₃ cytochrome c oxidase, some also encode cytochrome bd-type quinol oxidase or ba₃ -type cytochrome c oxidase, which could facilitate growth under different oxygen tensions, or maintain redox balance. Carboxysome operons are present in most, with genes downstream encoding transporters from four evolutionarily distinct families, which may act with the carboxysomes to form CO₂ concentrating mechanisms. These adaptations to habitat variability likely contribute to the cosmopolitan distribution of these organisms.