The impact of insertion sequences on bacterial genome plasticity and adaptability

Insertion sequence IS6770 modulates potassium symporter kup transcription in Enterococcus faecalis JH2-2 under low pH conditions

Enterococcus faecalis is a phylogenetically and industrially relevant microorganism associated with Lactic Acid Bacteria. Some strains of this bacterium are employed as probiotics in commercial applications, while others serve as the principal component in starter cultures for artisanal regional cheese production. However, over the last decade, this species has emerged as an opportunistic multiresistant pathogen, raising concerns about its impact on human health. Recently, we identified multiple potassium transporter systems in E. faecalis, including the Ktr systems (KtrAB and KtrAD), Kup, KimA and Kdp complex (KdpFABC). Nevertheless, the physiological significance of these proteins remains not fully understood. In this study, we observed that the kup gene promoter region in the JH2-2 strain was modified due to the insertion of a complete copy of the IS6770 insertion sequence. Consequently, we investigated the influence of IS6770 on the expression of the kup gene. To achieve this, we conducted a mapping of the promoter region of this gene in the E. faecalis JH2-2 strain, employing fluorescence gene reporters. In addition, a transcriptional analysis of the kup gene was executed in a strain derived from E. faecalis V583 that lacks the IS30-related insertion element, facilitating the identification of the transcriptional start site. Next, the expression of the kup gene was evaluated via RT-qPCR under different pH stressful conditions. A strong upregulation of the kup gene was observed at an initial pH of 5.0 in the strain derived from E. faecalis V583. However, the activation of transcription was not observed in the E. faecalis JH2-2 strain due to the hindrance caused by the presence of IS6770. Besides that, our computational analysis of E. faecalis genomes elucidates a plausible association between transposition and the regulation of the kup gene. Remarkably, the ubiquitous presence of IS6770 throughout the phylogenetic tree implies its ancient existence within E. faecalis. Moreover, the recurrent co-occurrence of IS6770 with the kup gene, observed in 30 % of IS6770-positive strains, alludes to the potential involvement of this genomic arrangement in the adaptive strategies of E. faecalis across diverse niches.

A study at the wildlife-livestock interface unveils the potential of feral swine as a reservoir for extended-spectrum β-lactamase-producing Escherichia coli

Wildlife is known to serve as carriers and sources of antimicrobial resistance (AMR). Due to their unrestricted movements and behaviors, they can spread antimicrobial resistant bacteria among livestock, humans, and the environment, thereby accelerating the dissemination of AMR. Extended-spectrum β-lactamase (ESBL)-producing Enterobacteriaceae is one of major concerns threatening human and animal health, yet transmission mechanisms at the wildlife-livestock interface are not well understood. Here, we investigated the mechanisms of ESBL-producing bacteria spreading across various hosts, including cattle, feral swine, and coyotes in the same habitat range, as well as from environmental samples over a two-year period. We report a notable prevalence and clonal dissemination of ESBL-producing E. coli in feral swine and coyotes, suggesting their persistence and adaptation within wildlife hosts. In addition, in silico studies showed that horizontal gene transfer, mediated by conjugative plasmids and insertion sequences elements, may play a key role in spreading the ESBL genes among these bacteria. Furthermore, the shared gut resistome of cattle and feral swine suggests the dissemination of antibiotic resistance genes at the wildlife-livestock interface. Taken together, our results suggest that feral swine may serve as a reservoir of ESBL-producing E. coli.

Metagenomic profiling uncovers microbiota and antibiotic resistance patterns across human, chicken, pig fecal, and soil environments

The ongoing and progressive evolution of antibiotic resistance presents escalating challenges for the clinical management and prevention of bacterial infections. Understanding the makeup of resistance genomes and accurately quantifying the current abundance of antibiotic resistance genes (ARGs) are crucial for assessing the threat of antimicrobial resistance (AMR) to public health. This comprehensive study investigated the distribution and diversity of bacterial community composition, ARGs, and virulence factors (VFs) across human, chicken, pig fecal, and soil microbiomes in various provinces of China. As a result, multidrug resistance was identified across all samples. Core ARGs primarily related to multidrug, MLS (Macrolides-Lincosamide-Streptogramins), and tetracycline resistance were characterized. A significant correlation between ARGs and bacterial taxa was observed, especially in soil samples. Probiotic strains such as Lactobacillus harbored ARGs, potentially contributing to the dissemination of antibiotic resistance. We screened subsets of ARGs from samples from different sources as indicators to assess the level of ARGs contamination in samples, with high accuracy. These results underline the complex relationship between microbial communities, resistance mechanisms, and environmental factors, emphasizing the importance of continued research and monitoring to better understand these dynamics.

Multidrug resistance in Salmonella isolates of swine origin: mobile genetic elements and plasmids associated with cephalosporin resistance with potential transmission to humans

The emergence of foodborne Salmonella strains carrying antimicrobial resistance (AMR) in mobile genetic elements (MGE) is a significant public health threat in a One Health context requiring continuous surveillance. Resistance to ciprofloxacin and cephalosporins is of particular concern. Since pigs are a relevant source of foodborne Salmonella for human beings, we studied transmissible AMR genes and MGE in a collection of 83 strains showing 9 different serovars and 15 patterns of multidrug resistant (MDR) previously isolated from pigs raised in the conventional breeding system of Northern Spain. All isolates were susceptible to ciprofloxacin and three isolates carried blaCMY-2 or blaCTX-M-9 genes responsible for cefotaxime resistance. Filter mating experiments showed that the two plasmids carrying blaCTX-M-9 were conjugative while that carrying blaCMY-2 was self-transmissible by transformation. Whole-genome sequencing and comparative analyses were performed on the isolates and plasmids. The IncC plasmid pSB109, carrying blaCMY-2, was similar to one found in S. Reading from cattle, indicating potential horizontal transfer between serovars and animal sources. The IncHI2 plasmids pSH102 in S. Heidelberg and pSTM45 in S. Typhimurium ST34, carrying blaCTX-M-9, shared similar backbones and two novel "complex class 1 integrons" containing different AMR and heavy metal genes. Our findings emphasize the importance of sequencing techniques to identify emerging AMR regions in conjugative and stable plasmids from livestock production. The presence of MGE carrying clinically relevant AMR genes raises public health concerns, requiring monitoring to mitigate the emergence of bacteria carrying AMR genes and subsequent spread through animals and food.IMPORTANCEThe emergence of foodborne Salmonella strains carrying antimicrobial resistance (AMR) in mobile genetic elements (MGE) is a significant public health threat in a One Health context. Since pigs are a relevant source of foodborne Salmonella for humans, in this study, we investigate different aspects of AMR in a collection of 83 Salmonella showing nine different serovars and 15 patterns of multidrug resistant (MDR) isolated from pigs raised in the conventional breeding system. Our findings emphasize the importance of sequencing techniques to identify emerging AMR regions in conjugative and stable plasmids from livestock production. The presence of MGE carrying clinically relevant AMR genes raises public health concerns, requiring monitoring to mitigate the emergence of bacteria carrying AMR genes and subsequent spread through animals and food.

Promoter recruitment drives the emergence of proto-genes in a long-term evolution experiment with Escherichia coli

The phenomenon of de novo gene birth-the emergence of genes from non-genic sequences-has received considerable attention due to the widespread occurrence of genes that are unique to particular species or genomes. Most instances of de novo gene birth have been recognized through comparative analyses of genome sequences in eukaryotes, despite the abundance of novel, lineage-specific genes in bacteria and the relative ease with which bacteria can be studied in an experimental context. Here, we explore the genetic record of the Escherichia coli long-term evolution experiment (LTEE) for changes indicative of "proto-genic" phases of new gene birth in which non-genic sequences evolve stable transcription and/or translation. Over the time span of the LTEE, non-genic regions are frequently transcribed, translated and differentially expressed, with levels of transcription across low-expressed regions increasing in later generations of the experiment. Proto-genes formed downstream of new mutations result either from insertion element activity or chromosomal translocations that fused preexisting regulatory sequences to regions that were not expressed in the LTEE ancestor. Additionally, we identified instances of proto-gene emergence in which a previously unexpressed sequence was transcribed after formation of an upstream promoter, although such cases were rare compared to those caused by recruitment of preexisting promoters. Tracing the origin of the causative mutations, we discovered that most occurred early in the history of the LTEE, often within the first 20,000 generations, and became fixed soon after emergence. Our findings show that proto-genes emerge frequently within evolving populations, can persist stably, and can serve as potential substrates for new gene formation.

Horizontal gene transfer in plant microbiomes: integrons as hotspots for cross-species gene exchange

Plant microbiomes play important roles in plant health and fitness. Bacterial horizontal gene transfer (HGT) can influence plant health outcomes, driving the spread of both plant growth-promoting and phytopathogenic traits. However, community dynamics, including the range of genetic elements and bacteria involved in this process are still poorly understood. Integrons are genetic elements recently shown to be abundant in plant microbiomes, and are associated with HGT across broad phylogenetic boundaries. They facilitate the spread of gene cassettes, small mobile elements that collectively confer a diverse suite of adaptive functions. Here, we analysed 5,565 plant-associated bacterial genomes to investigate the prevalence and functional diversity of integrons in this niche. We found that integrons are particularly abundant in the genomes of Pseudomonadales, Burkholderiales, and Xanthomonadales. In total, we detected nearly 9,000 gene cassettes, and found that many could be involved in plant growth promotion or phytopathogenicity, suggesting that integrons might play a role in bacterial mutualistic or pathogenic lifestyles. The rhizosphere was enriched in cassettes involved in the transport and metabolism of diverse substrates, suggesting that they may aid in adaptation to this environment, which is rich in root exudates. We also found that integrons facilitate cross-species HGT, which is particularly enhanced in the phyllosphere. This finding may provide an ideal opportunity to promote plant growth by fostering the spread of genes cassettes relevant to leaf health. Together, our findings suggest that integrons are important elements in plant microbiomes that drive HGT, and have the potential to facilitate plant host adaptation.

A pan-genomic assessment: Delving into the genome of the marine epiphyte Bacillus altitudinis strain 19_A and other very close Bacillus strains from multiple environments

Marine macroalgae are the habitat of epiphytic bacteria and provide several conditions for a beneficial biological interaction to thrive. Although Bacillus is one of the most abundant epiphytic genera, genomic information on marine macroalgae-associated Bacillus species remains scarce. In this study, we further investigated our previously published genome of the epiphytic strain Bacillus altitudinis 19_A to find features that could be translated to potential metabolites produced by this microorganism, as well as genes that play a role in its interaction with its macroalgal host. To achieve this goal, we performed a pan-genome analysis of Bacillus sp. and a codon bias assessment, including the genome of the strain Bacillus altitudinis 19_A and 29 complete genome sequences of closely related Bacillus strains isolated from soil, marine environments, plants, extreme environments, air, and food. This genomic analysis revealed that Bacillus altitudinis 19_A possessed unique genes encoding proteins involved in horizontal gene transfer, DNA repair, transcriptional regulation, and bacteriocin biosynthesis. In this comparative analysis, codon bias was not associated with the habitat of the strains studied. Some accessory genes were identified in the Bacillus altitudinis 19_A genome that could be related to its epiphytic lifestyle, as well as gene clusters for the biosynthesis of a sporulation-killing factor and a bacteriocin, showing their potential as a source of antimicrobial peptides. Our results provide a comprehensive view of the Bacillus altitudinis 19_A genome to understand its adaptation to the marine environment and its potential as a producer of bioactive compounds.

Incidence of carbapenem-resistant Escherichia coli ST2437 of clinical origin harbouring blaOXA-144 gene: a report from India

AIMS

Carbapenem-resistant Escherichia coli has been categorized as a pathogen of critical priority by the World Health Organization as it is highly infectious with high mortality and morbidity rates and widespread transmission potential. Carbapenem resistance is primarily mediated by carbapenemase-encoding genes and, additionally, through intrinsic factors. In India, over the years, carbapenemase-encoding genes have been reported from diverse clinically significant pathogens. The present study identifies E. coli of clinical origin that harbours blaOXA-144.

METHODS AND RESULTS

The study isolate was obtained from a tertiary referral hospital in northeast India. Carbapenemase production was investigated through culture on chromogenic agar and Rapidec Carba NP test as per manufacturer's instructions. Susceptibility of the isolate was performed by the Kirby-Bauer disc diffusion method and agar dilution method following CLSI guidelines. PCR targeting carbapenemase-encoding genes was performed, followed by transformation and conjugation experiments. Whole-genome sequencing of the isolate was done through the Illumina sequencing platform and the data were analysed using the Centre for Genomic Epidemiology database. BJD_EC180 is 6 919 180 bp in length and consists of six rRNA operons, 111 tRNA, and 6849 predicted protein-coding sequences. BJD_EC180 belonged to ST2437 and harboured the carbapenemase-encoding gene blaOXA-144 with ISAba1 upstream, along with multiple antibiotic resistance genes conferring clinical resistance towards beta-lactams, aminoglycosides, amphenicols, sulphonamides, tetracyclines, trimethoprim, and rifampin.

CONCLUSIONS

Carbapenem-resistant E. coli harbouring blaOXA-144 associated with insertion sequence pose a serious health threat as their mobilization into carbapenem non-susceptible strains that will contribute to the resistance burden and therefore, needs urgent monitoring.

Lineage-specific evolution of Aquibium, a close relative of Mesorhizobium, during habitat adaptation

The novel genus Aquibium that lacks nitrogenase was recently reclassified from the Mesorhizobium genus. The genomes of Aquibium species isolated from water were smaller and had higher GC contents than those of Mesorhizobium species. Six Mesorhizobium species lacking nitrogenase were found to exhibit low similarity in the average nucleotide identity values to the other 24 Mesorhizobium species. Therefore, they were classified as the non-N2-fixing Mesorhizobium lineage (N-ML), an evolutionary intermediate species. The results of our phylogenomic analyses and the loss of Rhizobiales-specific fur/mur indicated that Mesorhizobium species may have evolved from Aquibium species through an ecological transition. Halotolerant and alkali-resistant Aquibium and Mesorhizobium microcysteis belonging to N-ML possessed many tripartite ATP-independent periplasmic transporter and sodium/proton antiporter subunits composed of seven genes (mrpABCDEFG). These genes were not present in the N2-fixing Mesorhizobium lineage (ML), suggesting that genes acquired for adaptation to highly saline and alkaline environments were lost during the evolution of ML as the habitat changed to soil. Land-to-water habitat changes in Aquibium species, close relatives of Mesorhizobium species, could have influenced their genomic evolution by the gain and loss of genes. Our study indicated that lineage-specific evolution could have played a significant role in shaping their genome architecture and conferring their ability to thrive in different habitats.IMPORTANCEPhylogenetic analyses revealed that the Aquibium lineage (AL) and non-N2-fixing Mesorhizobium lineage (N-ML) were monophyletically grouped into distinct clusters separate from the N2-fixing Mesorhizobium lineage (ML). The N-ML, an evolutionary intermediate species having characteristics of both ancestral and descendant species, could provide a genomic snapshot of the genetic changes that occur during adaptation. Genomic analyses of AL, N-ML, and ML revealed that changes in the levels of genes related to transporters, chemotaxis, and nitrogen fixation likely reflect adaptations to different environmental conditions. Our study sheds light on the complex and dynamic nature of the evolution of rhizobia in response to changes in their environment and highlights the crucial role of genomic analysis in understanding these processes.

The first vanE-type vancomycin resistant Enterococcus faecalis isolates in Norway - phenotypic and molecular characteristics

OBJECTIVES

We aimed to characterize the vanE cluster and its genetic support in the first Norwegian vanE-type isolates and assess genetic relatedness to other vanE isolates.

METHODS

Two vanE-type vancomycin resistant Enterococcus faecalis (vanE-VREfs) isolates (E1 and E2) recovered from the same patient 30 months apart were examined for antimicrobial susceptibility, genome sequence, vancomycin resistance induction, vanE transferability, genome mutation rate, and phylogenetic relationship to E. faecalis closed genomes and two vanE-VREfs from North America.

RESULTS

The ST34 E1 and E2 strains expressed low-level vancomycin resistance and susceptibility to teicoplanin. Their vanE gene clusters were part of a non-transferable Tn6202. The histidine kinase part of vanSE was expressed although a premature stop codon (E1) and insertion of a transposase (E2) truncated their vanSE gene. The vancomycin resistance phenotype in E1 was inducible while constitutive in E2. E1 showed a 125-fold higher mutation rate than E2. Variant calling showed 60 variants but nearly identical chromosomal gene content and synteny between the isolates. Their genomes also showed high similarity to another ST34 vanE-VREfs from Canada.

CONCLUSION

In-depth genomic analyses of the first two vanE-VREfs found in Europe identified a single chromosomal insertion site of two variants of vanE-conferring Tn6202. Single nucleotide polymorphism (SNP) and core genome multilocus sequence type (cgMLST) analyses show the genomes are different. This can be explained by the high mutation rate of E1 and acquisition of different mobile genetic elements; thus, we believe the two isolates from the same patient are genetically related. Genome similarities also suggest relatedness between the Canadian and Norwegian vanE-VREfs.

Characterisation of IS1311 in Mycobacterium avium subspecies paratuberculosis genomes: Typing, continental clustering, microbial evolution and host adaptation

Johne's disease (JD), caused by Mycobacterium avium subspecies paratuberculosis (MAP) is a global burden for livestock producers and has an association with Crohn's disease in humans. Within MAP there are two major lineages, S/Type I/TypeIII and C/Type II, that vary in phenotype including culturability, host preference and virulence. These lineages have been identified using the IS1311 element, which contains a conserved, single nucleotide polymorphism. IS1311 and the closely related IS1245 element belong to the IS256 family of insertion sequences, are dispersed throughout M. avium taxa but remain poorly characterised. To investigate the distribution and diversity of IS1311 in MAP, 805 MAP genomes were collated from public databases. IS1245 was absent, while IS1311 sequence, copy number and insertion loci were conserved between MAP S lineages and varied within the MAP C lineage. One locus was specific to the S strains, which contained nine IS1311 copies. In contrast, C strains contained either seven or eight IS1311 loci. Most insertion loci were associated with the boundaries of homologous regions that had undergone genome rearrangement between the MAP lineages, suggesting that this sequence may be a driver of recombination. Phylogenomic geographic clustering of MAP subtypes was demonstrated for the first time, at continental scale, and indicated that there may have been recent MAP transmission between Europe and North America, in contrast to Australia where importation of live ruminants is generally prohibited. This investigation confirmed the utility of IS1311 typing in epidemiological studies and resolved anomalies in past studies. The results shed light on potential mechanisms of niche/host adaptation, virulence of MAP and global transmission dynamics.

Emerging Issues and Initial Insights into Bacterial Biofilms: From Orthopedic Infection to Metabolomics

Bacterial biofilms, enigmatic communities of microorganisms enclosed in an extracellular matrix, still represent an open challenge in many clinical contexts, including orthopedics, where biofilm-associated bone and joint infections remain the main cause of implant failure. This study explores the scenario of biofilm infections, with a focus on those related to orthopedic implants, highlighting recently emerged substantial aspects of the pathogenesis and their potential repercussions on the clinic, as well as the progress and gaps that still exist in the diagnostics and management of these infections. The classic mechanisms through which biofilms form and the more recently proposed new ones are depicted. The ways in which bacteria hide, become impenetrable to antibiotics, and evade the immune defenses, creating reservoirs of bacteria difficult to detect and reach, are delineated, such as bacterial dormancy within biofilms, entry into host cells, and penetration into bone canaliculi. New findings on biofilm formation with host components are presented. The article also delves into the emerging and critical concept of immunometabolism, a key function of immune cells that biofilm interferes with. The growing potential of biofilm metabolomics in the diagnosis and therapy of biofilm infections is highlighted, referring to the latest research.

The Transposition of Insertion Sequences in Sigma-Factor- and LysR-Deficient Mutants of Deinococcus geothermalis

Some insertion sequence (IS) elements were actively transposed using oxidative stress conditions, including gamma irradiation and hydrogen peroxide treatment, in Deinococcus geothermalis, a radiation-resistant bacterium. D. geothermalis wild-type (WT), sigma factor gene-disrupted (∆dgeo_0606), and LysR gene-disrupted (∆dgeo_1692) mutants were examined for IS induction that resulted in non-pigmented colonies after gamma irradiation (5 kGy) exposure. The loss of pigmentation occurred because dgeo_0524, which encodes a phytoene desaturase in the carotenoid pathway, was disrupted by the transposition of IS elements. The types and loci of the IS elements were identified as ISDge2 and ISDge6 in the ∆dgeo_0606 mutant and ISDge5 and ISDge7 in the ∆dgeo_1692 mutant, but were not identified in the WT strain. Furthermore, 80 and 100 mM H2O2 treatments induced different transpositions of IS elements in ∆dgeo_0606 (ISDge5, ISDge6, and ISDge7) and WT (ISDge6). However, no IS transposition was observed in the ∆dgeo_1692 mutant. The complementary strain of the ∆dgeo_0606 mutation showed recovery effects in the viability assay; however, the growth-delayed curve did not return because the neighboring gene dgeo_0607 was overexpressed, probably acting as an anti-sigma factor. The expression levels of certain transposases, recognized as pivotal contributors to IS transposition, did not precisely correlate with active transposition in varying oxidation environments. Nevertheless, these findings suggest that specific IS elements integrated into dgeo_0524 in a target-gene-deficient and oxidation-source-dependent manner.

Insertion sequence excision is enhanced by a protein that catalyzes branch migration and promotes microhomology-mediated end joining

Insertion sequence (IS)-excision enhancer (IEE) promotes the excision of ISs in the genome of enterohemorrhagic Escherichia coli O157. Because IEE-dependent IS excision occurs in the presence of transposase, the process of IS transposition may be involved in IS excision; however, little is understood about the molecular mechanisms of IS excision. Our in vitro analysis revealed that IEE exhibits DNA-dependent ATPase activity, which is activated by branched DNA. IEE also catalyzes the branch migration of fork-structured DNA. These results suggest that IEE remodels branched structures of the IS transposition intermediate. Sequence analysis of recombination sites in IS-excision products suggested that microhomologous sequences near the ends of the IS are involved in IS excision. IEE promoted microhomology-mediated end joining (MMEJ), in which base pairing between 6-nucleotides complementary ends of two 3'-protruding DNAs and subsequent elongation of the paired DNA strand occurred. IS-excision frequencies were significantly decreased in cells producing IEE mutants that had lost either branch migration or MMEJ activity, which suggests that these activities of IEE are required for IS excision. Based on our results, we propose a model for IS excision triggered by IEE and transposase.

PlasmidEC and gplas2: an optimized short-read approach to predict and reconstruct antibiotic resistance plasmids in Escherichia coli

Accurate reconstruction of Escherichia coli antibiotic resistance gene (ARG) plasmids from Illumina sequencing data has proven to be a challenge with current bioinformatic tools. In this work, we present an improved method to reconstruct E. coli plasmids using short reads. We developed plasmidEC, an ensemble classifier that identifies plasmid-derived contigs by combining the output of three different binary classification tools. We showed that plasmidEC is especially suited to classify contigs derived from ARG plasmids with a high recall of 0.941. Additionally, we optimized gplas, a graph-based tool that bins plasmid-predicted contigs into distinct plasmid predictions. Gplas2 is more effective at recovering plasmids with large sequencing coverage variations and can be combined with the output of any binary classifier. The combination of plasmidEC with gplas2 showed a high completeness (median=0.818) and F1-Score (median=0.812) when reconstructing ARG plasmids and exceeded the binning capacity of the reference-based method MOB-suite. In the absence of long-read data, our method offers an excellent alternative to reconstruct ARG plasmids in E. coli.

Genomic analysis of a cAmpC (CMY-41)-producing Citrobacter freundii ST64 isolated from patient

Antibiotic resistance in Citrobacter freundii is a public health concern. This study evaluated the closed genome of a C. freundii isolated from the stool of a hospitalized patient initially related to a Salmonella outbreak. Confirmation of the isolate was determined by whole-genome sequencing. Nanopore sequencing was performed using a MinION with a Flongle flow cell. Assembly using SPAdes and Unicycler yielded a closed genome annotated by National Center for Biotechnology Information Prokaryotic Genome Annotation Pipeline. Genomic analyses employed MLST 2.0, ResFinder4.1, PlasmidFinder2.1, and VFanalyzer. Phylogenetic comparison utilized the Center for Food Safety and Applied Nutrition (CFSAN)-single nucleotide polymorphism pipeline and Genetic Algorithm for Rapid Likelihood Inference. Antimicrobial susceptibility was tested by broth microdilution following Clinical and Laboratory Standards Institute criteria. Multi-locus sequence type in silico analysis assigned the C. freundii as sequence type 64 and the blaCMY-41 gene was detected in resistome investigation. The susceptibility to antibiotics, determined using Sensititre® plates, revealed resistance to aztreonam, colistin, cefoxitin, amoxicillin/clavulanic acid, sulfisoxazole, ampicillin, and streptomycin. The genetic relatedness of the C. freundii CFSAN077772 with publicly available C. freundii genomes revealed a close relationship to a C. freundii SRR1186659, isolated in 2009 from human stool in Tanzania. In addition, C. freundii CFSAN077772 is nested in the same cluster with C. freundii clinical strains isolated in Denmark, Mexico, Myanmar, and Canada, suggesting a successful intercontinental spread.

Within-Host Evolution of Bacterial Pathogens in Acute and Chronic Infection

Bacterial pathogens undergo remarkable adaptive change in response to the selective forces they encounter during host colonization and infection. Studies performed over the past few decades have demonstrated that many general evolutionary processes can be discerned during the course of host adaptation, including genetic diversification of lineages, clonal succession events, convergent evolution, and balanced fitness trade-offs. In some cases, elevated mutation rates resulting from mismatch repair or proofreading deficiencies accelerate evolution, and active mobile genetic elements or phages may facilitate genome plasticity. The host immune response provides another critical component of the fitness landscapes guiding adaptation, and selection operating on pathogens at this level may lead to immune evasion and the establishment of chronic infection. This review summarizes recent advances in this field, with a special focus on different forms of bacterial genome plasticity in the context of infection, and considers clinical consequences of adaptive changes for the host.

Annotation and Comparative Genomics of Prokaryotic Transposable Elements

The data generated in nearly 30 years of bacterial genome sequencing has revealed the abundance of transposable elements (TE) and their importance in genome and transcript remodeling through the mediation of DNA insertions and deletions, structural rearrangements, and regulation of gene expression. Furthermore, what we have learned from studying transposition mechanisms and their regulation in bacterial TE is fundamental to our current understanding of TE in other organisms because much of what has been observed in bacteria is conserved in all domains of life. However, unlike eukaryotic TE, prokaryotic TE sequester and transmit important classes of genes that impact host fitness, such as resistance to antibiotics and heavy metals and virulence factors affecting animals and plants, among other acquired traits. This provides dynamism and plasticity to bacteria, which would otherwise be propagated clonally. The insertion sequences (IS), the simplest form of prokaryotic TE, are autonomous and compact mobile genetic elements. These can be organized into compound transposons, in which two similar IS can flank any DNA segment and render it transposable. Other more complex structures, called unit transposons, can be grouped into four major families (Tn3, Tn7, Tn402, Tn554) with specific genetic characteristics. This chapter will revisit the prominent structural features of these elements, focusing on a genomic annotation framework and comparative analysis. Relevant aspects of TE will also be presented, stressing their key position in genome impact and evolution, especially in the emergence of antimicrobial resistance and other adaptive traits.

Comparative genomics reveals the correlations of stress response genes and bacteriophages in developing antibiotic resistance of Staphylococcus saprophyticus

IMPORTANCE

Staphylococcus saprophyticus is the second most common bacteria associated with urinary tract infections (UTIs) in women. The antimicrobial treatment regimen for uncomplicated UTI is normally nitrofurantoin, trimethoprim-sulfamethoxazole (TMP-SMX), or a fluoroquinolone without routine susceptibility testing of S. saprophyticus recovered from urine specimens. However, TMP-SMX-resistant S. saprophyticus has been detected recently in UTI patients, as well as in our cohort. Herein, we investigated the understudied resistance patterns of this pathogenic species by linking genomic antibiotic resistance gene (ARG) content to susceptibility phenotypes. We describe ARG associations with known and novel SCCmec configurations as well as phage elements in S. saprophyticus, which may serve as intervention or diagnostic targets to limit resistance transmission. Our analyses yielded a comprehensive database of phenotypic data associated with the ARG sequence in clinical S. saprophyticus isolates, which will be crucial for resistance surveillance and prediction to enable precise diagnosis and effective treatment of S. saprophyticus UTIs.

High genomic plasticity and unique features of Xanthomonas translucens pv. graminis revealed through comparative analysis of complete genome sequences

BACKGROUND

Xanthomonas translucens pv. graminis (Xtg) is a major bacterial pathogen of economically important forage grasses, causing severe yield losses. So far, genomic resources for this pathovar consisted mostly of draft genome sequences, and only one complete genome sequence was available, preventing comprehensive comparative genomic analyses. Such comparative analyses are essential in understanding the mechanisms involved in the virulence of pathogens and to identify virulence factors involved in pathogenicity.

RESULTS

In this study, we produced high-quality, complete genome sequences of four strains of Xtg, complementing the recently obtained complete genome sequence of the Xtg pathotype strain. These genomic resources allowed for a comprehensive comparative analysis, which revealed a high genomic plasticity with many chromosomal rearrangements, although the strains were highly related. A high number of transposases were exclusively found in Xtg and corresponded to 413 to 457 insertion/excision transposable elements per strain. These mobile genetic elements are likely to be involved in the observed genomic plasticity and may play an important role in the adaptation of Xtg. The pathovar was found to lack a type IV secretion system, and it possessed the smallest set of type III effectors in the species. However, three XopE and XopX family effectors were found, while in the other pathovars of the species two or less were present. Additional genes that were specific to the pathovar were identified, including a unique set of minor pilins of the type IV pilus, 17 TonB-dependent receptors (TBDRs), and 11 plant cell wall degradative enzymes.

CONCLUSION

These results suggest a high adaptability of Xtg, conferred by the abundance of mobile genetic elements, which could play a crucial role in pathogen adaptation. The large amount of such elements in Xtg compared to other pathovars of the species could, at least partially, explain its high virulence and broad host range. Conserved features that were specific to Xtg were identified, and further investigation will help to determine genes that are essential to pathogenicity and host adaptation of Xtg.

Short-term evolution of antibiotic responses in highly dynamic environments favors loss of regulation

Microbes inhabit natural environments that are remarkably dynamic, with sudden environmental shifts that require immediate action by the cell. To cope with changing environments, microbes are equipped with regulated response mechanisms that are only activated when needed. However, when exposed to extreme environments such as clinical antibiotic treatments, complete loss of regulation is frequently observed. Although recent studies suggest that the initial evolution of microbes in new environments tends to favor mutations in regulatory pathways, it is not clear how this evolution is affected by how quickly conditions change (i.e. dynamics), or which mechanisms are commonly used to implement new regulation. Here, we perform experimental evolution on continuous cultures of E. coli carrying the tetracycline resistance tet operon to identify specific types of mutations that adapt drug responses to different dynamical regimens of drug administration. When cultures are evolved under gradually increasing tetracycline concentrations, we observe no mutations in the tet operon, but a predominance of fine-tuning mutations increasing the affinity of alternative efflux pump AcrB to tetracycline. When cultures are instead periodically exposed to large drug doses, all populations developed transposon insertions in repressor TetR, resulting in loss of regulation of efflux pump TetA. We use a mathematical model of the dynamics of antibiotic responses to show that sudden exposure to large drug concentrations can overwhelm regulated responses, which cannot induce resistance fast enough, resulting in fitness advantage for constitutive expression of resistance. These results help explain the loss of regulation of antibiotic resistance by opportunistic pathogens evolving in clinical environments. Our experiment supports the notion that initial evolution in new ecological niches proceeds largely through regulatory mutations and suggests that transposon insertions are a main mechanism driving this process.

Promoter capture drives the emergence of proto-genes in Escherichia coli

The phenomenon of de novo gene birth-the emergence of genes from non-genic sequences-has received considerable attention due to the widespread occurrence of genes that are unique to particular species or genomes. Most instances of de novo gene birth have been recognized through comparative analyses of genome sequences in eukaryotes, despite the abundance of novel, lineage-specific genes in bacteria and the relative ease with which bacteria can be studied in an experimental context. Here, we explore the genetic record of the Escherichia coli Long-Term Evolution Experiment (LTEE) for changes indicative of "proto-genic" phases of new gene birth in which non-genic sequences evolve stable transcription and/or translation. Over the time-span of the LTEE, non-genic regions are frequently transcribed, translated and differentially expressed, thereby serving as raw material for new gene emergence. Most proto-genes result either from insertion element activity or chromosomal translocations that fused pre-existing regulatory sequences to regions that were not expressed in the LTEE ancestor. Additionally, we identified instances of proto-gene emergence in which a previously unexpressed sequence was transcribed after formation of an upstream promoter. Tracing the origin of the causative mutations, we discovered that most occurred early in the history of the LTEE, often within the first 20,000 generations, and became fixed soon after emergence. Our findings show that proto-genes emerge frequently within evolving populations, persist stably, and can serve as potential substrates for new gene formation.

Emergence of Raoultella ornithinolytica in human infections from different hospitals in Ecuador with OXA-48-producing resistance

Purpose

The purpose of this study was to highlight the clinical and molecular features of 13 Raoultella ornithinolytica strains isolated from clinical environments in Ecuador, and to perform comparative genomics with previously published genomes of Raoultella spp. As Raoultella is primarily found in environmental, clinical settings, we focused our work on identifying mechanisms of resistance that can provide this bacterium an advantage to establish and persist in hospital environments.

Methods

We analyzed 13 strains of Raoultella ornithinolytica isolated from patients with healthcare associated infections (HAI) in three hospitals in Quito and one in Santo Domingo de Los Tsáchilas, Ecuador, between November 2017 and April 2018. These isolates were subjected to phenotypic antimicrobial susceptibility testing, end-point polymerase chain reaction (PCR) to detect the presence of carbapenemases and whole-genome sequencing.

Results

Polymerase chain reaction revealed that seven isolates were positive isolates for blaOXA-48 and one for blaKPC-2 gene. Of the seven strains that presented the blaOXA-48 gene, six harbored it on an IncFII plasmid, one was inserted into the bacterial chromosome. The blaKPC gene was detected in an IncM2/IncR plasmid. From the bioinformatics analysis, nine genomes had the gene blaOXA-48, originating from Ecuador. Moreover, all R. ornithinolytica strains contained the ORN-1 gene, which confers resistance for β-lactams, such as penicillins and cephalosporins. Comparative genome analysis of the strains showed that the pangenome of R. ornithinolytica is considered an open pangenome, with 27.77% of core genes, which could be explained by the fact that the antibiotic resistance genes in the ancestral reconstruction are relatively new, suggesting that this genome is constantly incorporating new genes.

Conclusion

These results reveal the genome plasticity of R. ornithinolytica, particularly in acquiring antibiotic-resistance genes. The genomic surveillance and infectious control of these uncommon species are important since they may contribute to the burden of antimicrobial resistance and human health.

Antimicrobial resistome and mobilome in the urban river affected by combined sewer overflows and wastewater treatment effluent

The dissemination of antimicrobial resistance in the environment is an emerging global health problem. Wastewater treatment effluent and combined sewer overflows (CSOs) are major sources of antimicrobial resistance in urban rivers. This study aimed to clarify the effect of municipal wastewater treatment effluent and CSO on antimicrobial resistance genes (ARGs), mobile gene elements, and the microbial community in an urban river. The ARG abundance per 16S-based microbial population in the target river was 0.37-0.54 and 0.030-0.097 during the CSO event and dry weather, respectively. During the CSO event, the antimicrobial resistome in the river shifted toward a higher abundance of ARGs to clinically important drug classes, including macrolide, fluoroquinolone, and β-lactam, whereas ARGs to sulfonamide and multidrug by efflux pump were relatively abundant in dry weather. The abundance of intI1 and tnpA genes were highly associated with the total ARG abundance, suggesting their potential application as an indicator for estimating resistome contamination. Increase of prophage during the CSO event suggested that impact of CSO has a greater potential for horizontal gene transfer (HGT) via transduction. Consequently, CSO not only increases the abundance of ARGs to clinically important antimicrobials but also possibly enhances potential of HGT in urban rivers.

Hypermutator emergence in experimental Escherichia coli populations is stress-type dependent

Genotypes exhibiting an increased mutation rate, called hypermutators, can propagate in microbial populations because they can have an advantage due to the higher supply of beneficial mutations needed for adaptation. Although this is a frequently observed phenomenon in natural and laboratory populations, little is known about the influence of parameters such as the degree of maladaptation, stress intensity, and the genetic architecture for adaptation on the emergence of hypermutators. To address this knowledge gap, we measured the emergence of hypermutators over ~1,000 generations in experimental Escherichia coli populations exposed to different levels of osmotic or antibiotic stress. Our stress types were chosen based on the assumption that the genetic architecture for adaptation differs between them. Indeed, we show that the size of the genetic basis for adaptation is larger for osmotic stress compared to antibiotic stress. During our experiment, we observed an increased emergence of hypermutators in populations exposed to osmotic stress but not in those exposed to antibiotic stress, indicating that hypermutator emergence rates are stress type dependent. These results support our hypothesis that hypermutator emergence is linked to the size of the genetic basis for adaptation. In addition, we identified other parameters that covaried with stress type (stress level and IS transposition rates) that might have contributed to an increased hypermutator provision and selection. Our results provide a first comparison of hypermutator emergence rates under varying stress conditions and point towards complex interactions of multiple stress-related factors on the evolution of mutation rates.

Insertion sequence transposition inactivates CRISPR-Cas immunity

CRISPR-Cas immunity systems safeguard prokaryotic genomes by inhibiting the invasion of mobile genetic elements. Here, we screened prokaryotic genomic sequences and identified multiple natural transpositions of insertion sequences (ISs) into cas genes, thus inactivating CRISPR-Cas defenses. We then generated an IS-trapping system, using Escherichia coli strains with various ISs and an inducible cas nuclease, to monitor IS insertions into cas genes following the induction of double-strand DNA breakage as a physiological host stress. We identified multiple events mediated by different ISs, especially IS1 and IS10, displaying substantial relaxed target specificity. IS transposition into cas was maintained in the presence of DNA repair machinery, and transposition into other host defense systems was also detected. Our findings highlight the potential of ISs to counter CRISPR activity, thus increasing bacterial susceptibility to foreign DNA invasion.

A subclass of the IS1202 family of bacterial insertion sequences targets XerCD recombination sites

We describe here a new family of IS which are related to IS1202, originally isolated from Streptococcus pneumoniae in the mid-1990s and previously tagged as an emerging IS family in the ISfinder database. Members of this family have impacted some important properties of their hosts. We describe here another potentially important property of certain family members: specific targeting of xrs recombination sites. The family could be divided into three subgroups based on their transposase sequences and the length on the target repeats (DR) they generate on insertion: subgroup IS1202 (24-29 bp); ISTde1 (15-18 bp); and ISAba32 (5-6 bp). Members of the ISAba32 subgroup were repeatedly found abutting Xer recombinase recombination sites (xrs), separated by an intervening copy of a DR. These xrs sites, present in multiple copies in a number of Acinetobacter plasmids flanking antibiotic resistance genes, were proposed to form a new type of mobile genetic element using the chromosomally-encoded XerCD recombinase for mobility. Transposase alignments identified subgroup-specific indels which may be responsible for the differences in the transposition properties of the three subgroups (i.e. DR length and target specificity). We propose that this collection of IS be classed as a new insertion sequence family: the IS1202 family composed of three subgroups, only one of which specifically targets plasmid-borne xrs. We discuss the implications of xrs targeting for gene mobility.

Influence of insertion sequences on population structure of phytopathogenic bacteria in the Ralstonia solanacearum species complex

Ralstonia solanacearum species complex (RSSC) is a destructive group of plant pathogenic bacteria and the causative agent of bacterial wilt disease. Experimental studies have attributed RSSC virulence to insertion sequences (IS), transposable genetic elements which can both disrupt and activate host genes. Yet, the global diversity and distribution of RSSC IS are unknown. In this study, IS were bioinformatically identified in a diverse collection of 356 RSSC isolates representing five phylogenetic lineages and their diversity investigated based on genetic distance measures and comparisons with the ISFinder database. IS phylogenetic associations were determined based on their distribution across the RSSC phylogeny. Moreover, IS positions within genomes were characterised and their potential gene disruptions determined based on IS proximity to coding sequences. In total, we found 24732 IS belonging to eleven IS families and 26 IS subgroups with over half of the IS found in the megaplasmid. While IS families were generally widespread across the RSSC phylogeny, IS subgroups showed strong lineage-specific distributions and genetically similar bacterial isolates had similar IS contents. Similar associations with bacterial host genetic background were also observed with IS insertion positions which were highly conserved in closely related bacterial isolates. Finally, IS were found to disrupt genes with predicted functions in virulence, stress tolerance, and metabolism suggesting that they might be adaptive. This study highlights that RSSC insertion sequences track the evolution of their bacterial hosts potentially contributing to both intra- and inter-lineage genetic diversity.

A Genomic Snapshot of Antibiotic-ResistantEnterococcus faecalis within Public Hospital Environments in South Africa

Enterococci are among the most common opportunistic hospital pathogens. This study used whole-genome sequencing (WGS) and bioinformatics to determine the antibiotic resistome, mobile genetic elements, clone and phylogenetic relationship of Enterococcus faecalis isolated from hospital environments in South Africa. This study was carried out from September to November 2017. Isolates were recovered from 11 frequently touched sites by patients and healthcare workers in different wards at 4 levels of healthcare (A, B, C, and D) in Durban, South Africa. Out of the 245 identified E. faecalis isolates, 38 isolates underwent whole-genome sequencing (WGS) on the Illumina MiSeq platform, following microbial identification and antibiotic susceptibility tests. The tet(M) (31/38, 82%) and erm(C) (16/38, 42%) genes were the most common antibiotic-resistant genes found in isolates originating from different hospital environments which corroborated with their antibiotic resistance phenotypes. The isolates harboured mobile genetic elements consisting of plasmids (n = 11) and prophages (n = 14) that were mostly clone-specific. Of note, a large number of insertion sequence (IS) families were found on the IS3 (55%), IS5 (42%), IS1595 (40%), and Tn3 transposons the most predominant. Microbial typing using WGS data revealed 15 clones with 6 major sequence types (ST) belonging to ST16 (n = 7), ST40 (n = 6), ST21 (n = 5), ST126 (n = 3), ST23 (n = 3), and ST386 (n = 3). Phylogenomic analysis showed that the major clones were mostly conserved within specific hospital environments. However, further metadata insights revealed the complex intraclonal spread of these E. faecalis major clones between the sampling sites within each specific hospital setting. The results of these genomic analyses will offer insights into antibiotic-resistantE. faecalis in hospital environments relevant to the design of optimal infection prevention strategies in hospital settings.

Woronichinia naegeliana: A common nontoxigenic component of temperate freshwater cyanobacterial blooms with 30% of its genome in transposons

Monitoring in the U.S. state of Washington across the period 2007-2019 showed that Woronichinia has been present in many lakes state-wide. This cyanobacterium was commonly dominant or sub-dominant in cyanobacterial blooms in the wet temperate region west of the Cascade Mountains. In these lakes, Woronichinia often co-existed with Microcystis, Dolichospermum and Aphanizomenon flos-aquae and the cyanotoxin microcystin has often been present in those blooms, although it has not been known whether Woronichinia is a toxin producer. We report the first complete genome of Woronichinia naegeliana WA131, assembled from the metagenome of a sample collected from Wiser Lake, Washington, in 2018. The genome contains no genes for cyanotoxin biosynthesis or taste-and-odor compounds, but there are biosynthetic gene clusters for other bioactive peptides, including anabaenopeptins, cyanopeptolins, microginins and ribosomally produced, post-translationally modified peptides. Genes for photosynthesis, nutrient acquisition, vitamin synthesis and buoyancy that are typical of bloom-forming cyanobacteria are present, although nitrate and nitrite reductase genes are conspicuously absent. However, the 7.9 Mbp genome is 3-4 Mbp larger than those of the above-mentioned frequently co-existing cyanobacteria. The increased genome size is largely due to an extraordinary number of insertion sequence elements (transposons), which account for 30.3% of the genome and many of which are present in multiple copies. The genome contains a relatively large number of pseudogenes, 97% of which are transposase genes. W. naegeliana WA131 thus seems to be able to limit the potentially deleterious effects of high rates of recombination and transposition to the mobilome fraction of its genome.

Recent Trends in Antimicrobial Resistance among Anaerobic Clinical Isolates

Anaerobic bacteria are normal inhabitants of the human commensal microbiota and play an important role in various human infections. Tedious and time-consuming, antibiotic susceptibility testing is not routinely performed in all clinical microbiology laboratories, despite the increase in antibiotic resistance among clinically relevant anaerobes since the 1990s. β-lactam and metronidazole are the key molecules in the management of anaerobic infections, to the detriment of clindamycin. β-lactam resistance is usually mediated by the production of β-lactamases. Metronidazole resistance remains uncommon, complex, and not fully elucidated, while metronidazole inactivation appears to be a key mechanism. The use of clindamycin, a broad-spectrum anti-anaerobic agent, is becoming problematic due to the increase in resistance rate in all anaerobic bacteria, mainly mediated by Erm-type rRNA methylases. Second-line anti-anaerobes are fluoroquinolones, tetracyclines, chloramphenicol, and linezolid. This review aims to describe the up-to-date evolution of antibiotic resistance, give an overview, and understand the main mechanisms of resistance in a wide range of anaerobes.

Magnetic biochar promotes the risk of mobile genetic elements propagation in sludge anaerobic digestion

Mobile genetic elements (MGEs) mediated horizontal gene transfer is the primary reason for the propagation of antibiotic resistance genes in environment. The behavior of MGEs under magnetic biochar pressure in sludge anaerobic digestion (AD) is still unknown. This study evaluated the effects of different dosage magnetic biochar on the MGEs in AD reactors. The results showed that the biogas yield was highest (106.68 ± 1.16 mL g^-1 VS_added) with adding optimal dosage of magnetic biochar (25 mg g^-1 TS_added), due to it increased the microorganism's abundance involved in hydrolysis and methanogenesis. While, the total absolute abundance of MGEs in the reactors with magnetic biochar addition increased by 11.58%-77.37% compared with the blank reactor. When the dosage of magnetic biochar was 12.5 mg g^-1 TS_added, the relative abundance of most MGEs was the highest. The enrichment effect on ISCR1 was the most significant, and the enrichment rate reached 158.90-214.16%. Only the intI1 abundance was reduced and the removal rates yield 14.38-40.00%, which was inversely proportional to the dosage of magnetic biochar. Co-occurrence network explored that Proteobacteria (35.64%), Firmicutes (19.80%) and Actinobacteriota (15.84%) were the main potential host of MGEs. Magnetic biochar changed MGEs abundance by affecting the potential MGEs-host community structure and abundance. Redundancy analysis and variation partitioning analysis showed that the combined effect of polysaccharides, protein and sCOD exhibited the greatest contribution (accounted for 34.08%) on MGEs variation. These findings demonstrated that magnetic biochar increases the risk of MGEs proliferation in AD system.

Genomic distribution of the insertion sequence IS711 reveal a potential role in Brucella genome plasticity and host preference

The Insertion Sequence 711 (IS711) is linked to the Brucella genus. Mapping the genomic distribution of IS711 can help understand this insertion element's biological and evolutionary role. This work aimed to delineate the genomic distribution of the IS711 element and to study its association with Brucella evolution. A total of 124 genomes representing 9 Brucella species were searched using BLASTn sequence alignment tool to identify complete and truncated copies of IS711. Based on the genomic context, each IS711 locus was assigned a code using the initial letters of its neighboring genes. Various tools were used to annotate the neighboring genes and determine the shared synteny around orthologous IS711 loci. The tool Islandviewer 4 was used to scan for genomic islands. The Codon Tree method was used to build phylogenetic trees of B. melitensis, B. abortus, and B. suis genomes. The phylogenetic trees of the three species were analyzed, taking into account the genomic distribution patterns of IS711. The result of IS711 frequency analysis showed a relatively conserved number of copies/genome for the different species and for some biovars. The analysis showed that Brucella species with a relatively low IS711 copy number (4-8 copies/genome) are linked to domestic animals as primary hosts and have potential for zoonotic transmission. However, species with a relatively higher copy number (12-30 copies/genome) are less zoonotic and tend to be linked with wild animals as primary hosts. Analyzing the genomic distribution map of IS711 loci showed several unique patterns of IS711 distribution that are correlated with the evolution of Brucella species and biovars. The results also showed that 46.2% of the conserved IS711 elements are located within genomic islands. Based on our results and previous data, we postulate a model explaining the IS711 role in Brucella evolution. We assume that during the transition from a free-living to an intracellular lifestyle, a descendant of the Brucella genus had acquired a progenitor sequence of the IS711. Subsequently, a burst in IS711 transposition occurred. This parasitic expansion can be deleterious and has to be counteracted by evolutionary forces to prevent lineage extension and to promote adaptation to host. Similar to other plasmid-free pathogenic α-Proteobacteria bacteria, the balance of expansion and reduction of insertion elements could be one of the mechanisms to control genome reduction and streamlining. We hypothesize that the IS711-mediated genomic changes and other small sequence nucleotide changes in specific orthologous genes could significantly contribute to Brucella's evolution and adaptation to different animal hosts.

Variable resistance to zinc intoxication among Streptococcus agalactiae reveals a novel IS1381 insertion element within the zinc efflux transporter gene czcD

Streptococcus agalactiae, also known as group B Streptococcus, is an important human and animal pathogen. Zinc (Zn) is an essential trace element for normal bacterial physiology but intoxicates bacteria at high concentrations. Molecular systems for Zn detoxification exist in S. agalactiae, however the degree to which Zn detoxification may vary among different S. agalactiae isolates is not clear. We measured resistance to Zn intoxication in a diverse collection of clinical isolates of S. agalactiae by comparing the growth of the bacteria in defined conditions of Zn stress. We found significant differences in the ability of different S. agalactiae isolates to resist Zn intoxication; some strains such as S. agalactiae 18RS21 were able to survive and grow at 3.8-fold higher levels of Zn stress compared to other reference strains such as BM110 (6.4mM vs 1.68mM Zn as inhibitory, respectively). We performed in silico analysis of the available genomes of the S. agalactiae isolates used in this study to examine the sequence of czcD, which encodes an efflux protein for Zn that supports resistance in S. agalactiae. Interestingly, this revealed the presence of a mobile insertion sequence (IS) element, termed IS1381, in the 5' region of czcD in S. agalactiae strain 834, which was hyper-resistant to Zn intoxication. Interrogating a wider collection of S. agalactiae genomes revealed identical placement of IS1381 in czcD in other isolates from the clonal-complex-19 (CC19) 19 lineage. Collectively, these results show a resistance spectrum among S. agalactiae isolates enables survival in varying degrees of Zn stress, and this phenotypic variability has implications for understanding bacterial survival in metal stress.

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.

Mutation bias and the predictability of evolution

Predicting evolutionary outcomes is an important research goal in a diversity of contexts. The focus of evolutionary forecasting is usually on adaptive processes, and efforts to improve prediction typically focus on selection. However, adaptive processes often rely on new mutations, which can be strongly influenced by predictable biases in mutation. Here, we provide an overview of existing theory and evidence for such mutation-biased adaptation and consider the implications of these results for the problem of prediction, in regard to topics such as the evolution of infectious diseases, resistance to biochemical agents, as well as cancer and other kinds of somatic evolution. We argue that empirical knowledge of mutational biases is likely to improve in the near future, and that this knowledge is readily applicable to the challenges of short-term prediction. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

IS21 family transposase cleaved donor complex traps two right-handed superhelical crossings

Transposases are ubiquitous enzymes that catalyze DNA rearrangement events with broad impacts on gene expression, genome evolution, and the spread of drug-resistance in bacteria. Here, we use biochemical and structural approaches to define the molecular determinants by which IstA, a transposase present in the widespread IS21 family of mobile elements, catalyzes efficient DNA transposition. Solution studies show that IstA engages the transposon terminal sequences to form a high-molecular weight complex and promote DNA integration. A 3.4 Å resolution structure of the transposase bound to transposon ends corroborates our biochemical findings and reveals that IstA self-assembles into a highly intertwined tetramer that synapses two supercoiled terminal inverted repeats. The three-dimensional organization of the IstA•DNA cleaved donor complex reveals remarkable similarities with retroviral integrases and classic transposase systems, such as Tn7 and bacteriophage Mu, and provides insights into IS21 transposition.

Metagenomic profiling of antibiotic resistance genes in Red Sea brine pools

Antibiotic resistance (AR) is an alarming global health concern, causing an annual death rate of more than 35,000 deaths in the US. AR is a natural phenomenon, reported in several pristine environments. In this study, we report AR in pristine Red Sea deep brine pools. Antimicrobial resistance genes (ARGs) were detected for several drug classes with tetracycline and macrolide resistance being the most abundant. As expected, ARGs abundance increased in accordance with the level of human impact with pristine Red Sea samples having the lowest mean ARG level followed by estuary samples, while activated sludge samples showed a significantly higher ARG level. ARG hierarchical clustering grouped drug classes for which resistance was detected in Atlantis II Deep brine pool independent of the rest of the samples. ARG abundance was significantly lower in the Discovery Deep brine pool. A correlation between integrons and ARGs abundance in brine pristine samples could be detected, while insertion sequences and plasmids showed a correlation with ARGs abundance in human-impacted samples not seen in brine pristine samples. This suggests different roles of distinct mobile genetic elements (MGEs) in ARG distribution in pristine versus human-impacted sites. Additionally, we showed the presence of mobile antibiotic resistance genes in the Atlantis II brine pool as evidenced by the co-existence of integrases and plasmid replication proteins on the same contigs harboring predicted multidrug-resistant efflux pumps. This study addresses the role of non-pathogenic environmental bacteria as a silent reservoir for ARGs, and the possible horizontal gene transfer mechanism mediating ARG acquisition.

Genomic plasticity and adaptive capacity of the quaternary alkyl-ammonium compound and copper tolerant Acinetobacter bohemicus strain QAC-21b isolated from pig manure

Here, we present the genomic characterization of an Acinetobacter bohemicus strain QAC-21b which was isolated in the presence of a quaternary alky-ammonium compound (QAAC) from manure of a conventional German pig farm. The genetic determinants for QAAC, heavy metal and antibiotic resistances are reported based of the whole genome shotgun sequence and physiological growth tests. A. bohemicus QAC-21b grew in a species typical manner well at environmental temperatures but not at 37 °C. The strain showed tolerance to QAACs and copper but was susceptible to antibiotics relevant for Acinetobacter treatments. The genome of QAC-21b contained several Acinetobacter typical QAAC and heavy metal transporting efflux pumps coding genes, but no key genes for acquired antimicrobial resistances. The high genomic content of transferable genetic elements indicates that this bacterium can be involved in the transmission of antimicrobial resistances, if it is released with manure as organic fertilizer on agricultural fields. The genetic content of the strain was compared to that of two other A. bohemicus strains, the type strain ANC 3994^T, isolated from forest soil, and KCTC 42081, originally described as A. pakistanensis, a metal resistant strain isolated from a wastewater treatment pond. In contrast to the forest soil strain, both strains from anthropogenically impacted sources showed genetic features indicating their evolutionary adaptation to the anthropogenically impacted environments. Strain QAC-21b will be used as model strain to study the transmission of antimicrobial resistance to environmentally adapted Acinetobacter in agricultural environments receiving high content of pollutants with organic fertilizers from livestock husbandry.

Metaproteomic and gene expression analysis of interspecies interactions in a PAH-degrading synthetic microbial consortium constructed with the key microbes of a natural consortium

Polycyclic Aromatic Hydrocarbons (PAHs) impose adverse effects on the environment and human life. The use of synthetic microbial consortia is promising in bioremediation of contaminated sites with these pollutants. However, the design of consortia taking advantage of natural interactions has been poorly explored. In this study, a dual synthetic bacterial consortium (DSC_AB) was constructed with two key members (Sphingobium sp. AM and Burkholderia sp. Bk), of a natural PAH degrading consortium. DSC_AB showed significantly enhanced degradation of PAHs and toxic intermediary metabolites relative to the axenic cultures, indicating the existence of synergistic relationships. Metaproteomic and gene-expression analyses were applied to obtain a view of bacterial performance during phenanthrene removal. Overexpression of the Bk genes, naph, biph, tol and sal and the AM gene, ahdB, in DSC_AB relative to axenic cultures, demonstrated that both strains are actively participating in degradation, which gave evidence of cross-feeding. Several proteins related to stress response were under-expressed in DSC_AB relative to axenic cultures, indicating that the division of labour reduces cellular stress, increasing the efficiency of degradation. This is the one of the first works revealing bacterial relationships during PAH removal in a synthetic consortium applying an omics approach. Our findings could be used to develop criteria for evaluating the potential effectiveness of synthetic bacterial consortia in bioremediation.

Evidence for wastewaters as environments where mobile antibiotic resistance genes emerge

The emergence and spread of mobile antibiotic resistance genes (ARGs) in pathogens have become a serious threat to global health. Still little is known about where ARGs gain mobility in the first place. Here, we aimed to collect evidence indicating where such initial mobilization events of clinically relevant ARGs may have occurred. We found that the majority of previously identified origin species did not carry the mobilizing elements that likely enabled intracellular mobility of the ARGs, suggesting a necessary interplay between different bacteria. Analyses of a broad range of metagenomes revealed that wastewaters and wastewater-impacted environments had by far the highest abundance of both origin species and corresponding mobilizing elements. Most origin species were only occasionally detected in other environments. Co-occurrence of origin species and corresponding mobilizing elements were rare in human microbiota. Our results identify wastewaters and wastewater-impacted environments as plausible arenas for the initial mobilization of resistance genes.

The role of oxidative stress in genome destabilization and adaptive evolution of bacteria

The review is devoted to bacterial genome destabilization by oxidative stress. The article discusses the main groups of substances causing such stress. Stress regulons involved in destabilization of genetic material and mechanisms enhancing mutagenesis, bacterial genome rearrangements, and horizontal gene transfer, induced by oxidative damage to cell components are also considered. Based on the analysis of publications, it can be claimed that rapid development of new food substrates and ecological niches by microorganisms occurs due to acceleration of genetic changes induced by oxidative stress, mediated by several stress regulons (SOS, RpoS and RpoE) and under selective pressure. The authors conclude that non-lethal oxidative stress is probably-one of the fundamental processes that guide evolution of prokaryotes and a powerful universal trigger for adaptive destabilization of bacterial genome under changing environmental conditions.

Rewiring the respiratory pathway of Lactococcus lactis to enhance extracellular electron transfer

Lactococcus lactis, a lactic acid bacterium with a typical fermentative metabolism, can also use oxygen as an extracellular electron acceptor. Here we demonstrate, for the first time, that L. lactis blocked in NAD⁺ regeneration can use the alternative electron acceptor ferricyanide to support growth. By electrochemical analysis and characterization of strains carrying mutations in the respiratory chain, we pinpoint the essential role of the NADH dehydrogenase and 2-amino-3-carboxy-1,4-naphtoquinone in extracellular electron transfer (EET) and uncover the underlying pathway systematically. Ferricyanide respiration has unexpected effects on L. lactis, e.g., we find that morphology is altered from the normal coccoid to a more rod shaped appearance, and that acid resistance is increased. Using adaptive laboratory evolution (ALE), we successfully enhance the capacity for EET. Whole-genome sequencing reveals the underlying reason for the observed enhanced EET capacity to be a late-stage blocking of menaquinone biosynthesis. The perspectives of the study are numerous, especially within food fermentation and microbiome engineering, where EET can help relieve oxidative stress, promote growth of oxygen sensitive microorganisms and play critical roles in shaping microbial communities.

Natural variation in the hrpL promoter renders the phytopathogen Pseudomonas syringae pv. actinidiae nonpathogenic

The genetic basis underlying loss-of-virulence mutations that arise among natural phytopathogen populations is not well documented. In this study, we examined the virulence of 377 isolates of Pseudomonas syringae pv. actinidiae biovar 3 (Psa3) that were isolated from 76 kiwifruit orchards suffering from bacterial canker disease. Eighty-four nonpathogenic isolates were identified in 40 orchards. A nonpathogenic isolate G166 was found to be defective in hrpL transcription and the downstream type III secretion system (T3SS)-dependent phenotypes. Comparative genomics and complementary expression assay revealed that a single-base "G" insertion in the hrpL promoter blocks gene transcription by reducing promoter activity. The electrophoretic mobility shift assay showed that the genetic variation impairs σ⁵⁴ /promoter binding during gene transcription under hrp-inducing conditions, resulting in lower expression of hrpL. A PCR-restriction fragment length polymorphism assay was performed to trace the evolutionary history of this mutation, which revealed the independent onset of genetic variations in natural Psa3 populations. We also found that nonpathogenic variants outperformed virulent Psa3 bacteria for both epiphytic and apoplast colonization of kiwifruit leaves in mixed inoculations. Our study highlights a novel mechanism for loss of virulence in Psa3 and provides insight into bacterial adaptive evolution under natural settings.

Novel configurations of type I-E CRISPR-Cas system in Corynebacterium striatum clinical isolates

Clustered regularly interspaced short palindromic repeats (CRISPR) are a prokaryotic adaptive immune system that, through Cas proteins, promote the degradation of foreign nucleic acids such as phages and plasmids. We analyzed 10 genomes of Corynebacterium striatum clinical isolates from a public hospital in Rio de Janeiro, Brazil, the most emergent multidrug-resistant Corynebacterium species. All isolates were submitted to antimicrobial susceptibility testing. The occurrence and diversity of the CRISPR system were investigated by bioinformatics tools. Our analysis revealed that the isolates exhibited type I-E gene arrangements, and 3 more multidrug-resistant isolates, alternative type I-E gene arrangements, showing a divergent gene arrangement within the cas operon. Phylogenetic analysis of the cas1 gene of this type I-E CRISPR-Cas system alternative arrangement, termed here type I-E', showed a cluster in a distinct clade of the type I-E CRISPR-Cas system. The systems' guanine-cytosine (GC) content is lower than the genomic DNA's GC content, and mobile genetic elements were found in some isolates near the CRISPR-Cas system. Most CRISPR spacers are unknown indicating that there is a reservoir of unexplored corynebacteriophages and plasmids. Some spacers showed perfect homologies with phage and plasmid sequences. Intact phage regions were found in 3 of our isolates, ranging from 9.1 to 43.8 kb, with regions showing similarity to Rhodococcus and Corynebacterium phages. Our results may contribute to research about the CRISPR-Cas system diversity in C. striatum, where there are no published data to date.

Phylogenetic Relationships among TnpB-Containing Mobile Elements in Six Bacterial Species

Some families of mobile elements in bacterial genomes encode not only a transposase but also an accessory TnpB gene. This gene has been shown to encode an RNA-guided DNA endonuclease, co-evolving with Y1 transposase and serine recombinase in mobile elements IS605 and IS607. In this paper, we reveal the evolutionary relationships among TnpB-containing mobile elements (TCMEs) in well-assembled genomes of six bacterial species: Bacillus cereus, Clostridioides difficile, Deinococcus radiodurans, Escherichia coli, Helicobacter pylori and Salmonella enterica. In total, 9996 TCMEs were identified in 4594 genomes. They belonged to 39 different insertion sequences (ISs). Based on their genetic structures and sequence identities, the 39 TCMEs were classified into three main groups and six subgroups. According to our phylogenetic analysis, TnpBs include two main branches (TnpB-A and TnpB-B) and two minor branches (TnpB-C and TnpB-D). The key TnpB motifs and the associated Y1 and serine recombinases were highly conserved across species, even though their overall sequence identities were low. Substantial variation was observed for the rate of invasion across bacterial species and strains. Over 80% of the genomes of B. cereus, C. difficile, D. radiodurans and E. coli contained TCMEs; however, only 64% of the genomes of H. pylori and 44% of S. enterica genomes contained TCMEs. IS605 showed the largest rate of invasion in these species, while IS607 and IS1341 had a relatively narrow distribution. Co-invasions of IS605, IS607 and IS1341 elements were observed in various genomes. The largest average copy number was observed for IS605b elements in C. difficile. The average copy numbers of most other TCMEs were smaller than four. Our findings have important implications for understanding the co-evolution of TnpB-containing mobile elements and their biological roles in host genome evolution.

Factors associated with plasmid antibiotic resistance gene carriage revealed using large-scale multivariable analysis

Plasmids are major vectors of bacterial antibiotic resistance, but understanding of factors associated with plasmid antibiotic resistance gene (ARG) carriage is limited. We curated > 14,000 publicly available plasmid genomes and associated metadata. Duplicate and replicate plasmids were excluded; where possible, sample metadata was validated externally (BacDive database). Using Generalised Additive Models (GAMs) we assessed the influence of 12 biotic/abiotic factors (e.g. plasmid genetic factors, isolation source, collection date) on ARG carriage, modelled as a binary outcome. Separate GAMs were built for 10 major ARG types. Multivariable analysis indicated that plasmid ARG carriage patterns across time (collection years), isolation sources (human/livestock) and host bacterial taxa were consistent with antibiotic selection pressure as a driver of plasmid-mediated antibiotic resistance. Only 0.42% livestock plasmids carried carbapenem resistance (compared with 12% human plasmids); conversely, tetracycline resistance was enriched in livestock vs human plasmids, reflecting known prescribing practices. Interpreting results using a timeline of ARG type acquisition (determined by literature review) yielded additional novel insights. More recently acquired ARG types (e.g. colistin and carbapenem) showed increases in plasmid carriage during the date range analysed (1994-2019), potentially reflecting recent onset of selection pressure; they also co-occurred less commonly with ARGs of other types, and virulence genes. Overall, this suggests that following acquisition, plasmid ARGs tend to accumulate under antibiotic selection pressure and co-associate with other adaptive genes (other ARG types, virulence genes), potentially re-enforcing plasmid ARG carriage through co-selection.

The enterohemorrhagic Escherichia coli insertion sequence-excision enhancer protein is a DNA polymerase with microhomology-mediated end-joining activity

Bacterial genomes contain an abundance of transposable insertion sequence (IS) elements that are essential for genome evolution and fitness. Among them, IS629 is present in most strains of enterohemorrhagic Escherichia coli O157 and accounts for many polymorphisms associated with gene inactivation and/or genomic deletions. The excision of IS629 from the genome is promoted by IS-excision enhancer (IEE) protein. Despite IEE has been identified in the most pathogenic serotypes of E. coli, its biochemical features that could explain its role in IS excision are not yet understood. We show that IEE is present in >30% of all available E. coli genome assemblies, and is highly conserved and very abundant within enterohemorrhagic, enteropathogenic and enterotoxigenic genomes. In vitro analysis of the recombinant protein from E. coli O157:H7 revealed the presence of a Mn2+-dependent error-prone DNA polymerase activity in its N-terminal archaeo-eukaryotic primase (AEP) domain able to promote dislocations of the primer and template strands. Importantly, IEE could efficiently perform in vitro an end-joining reaction of 3'-single-strand DNA overhangs with ≥4 bp of homology requiring both the N-terminal AEP and C-terminal helicase domains. The proposed role for IEE in the novel IS excision mechanism is discussed.

An IS-mediated, RecA-dependent, bet-hedging strategy in Burkholderia thailandensis

Adaptation to fluctuating environmental conditions is difficult to achieve. Phase variation mechanisms can overcome this difficulty by altering genomic architecture in a subset of individuals, creating a phenotypically heterogeneous population with subpopulations optimized to persist when conditions change, or are encountered, suddenly. We have identified a phase variation system in Burkholderia thailandensis that generates a genotypically and phenotypically heterogeneous population. Genetic analyses revealed that RecA-mediated homologous recombination between a pair of insertion sequence (IS) 2-like elements duplicates a 208.6 kb region of DNA that contains 157 coding sequences. RecA-mediated homologous recombination also resolves merodiploids, and hence copy number of the region is varied and dynamic within populations. We showed that the presence of two or more copies of the region is advantageous for growth in a biofilm, and a single copy is advantageous during planktonic growth. While IS elements are well known to contribute to evolution through gene inactivation, polar effects on downstream genes, and altering genomic architecture, we believe that this system represents a rare example of IS element-mediated evolution in which the IS elements provide homologous sequences for amplification of a chromosomal region that provides a selective advantage under specific growth conditions, thereby expanding the lifestyle repertoire of the species.