Impact of suppression of the SOS response on protein expression in clinical isolates of Escherichia coli under antimicrobial pressure of ciprofloxacin

Introduction/objective

Suppression of the SOS response in combination with drugs damaging DNA has been proposed as a potential target to tackle antimicrobial resistance. The SOS response is the pathway used to repair bacterial DNA damage induced by antimicrobials such as quinolones. The extent of lexA-regulated protein expression and other associated systems under pressure of agents that damage bacterial DNA in clinical isolates remains unclear. The aim of this study was to assess the impact of this strategy consisting on suppression of the SOS response in combination with quinolones on the proteome profile of Escherichia coli clinical strains.

Materials and methods

Five clinical isolates of E. coli carrying different chromosomally- and/or plasmid-mediated quinolone resistance mechanisms with different phenotypes were selected, with E. coli ATCC 25922 as control strain. In addition, from each clinical isolate and control, a second strain was created, in which the SOS response was suppressed by deletion of the recA gene. Bacterial inocula from all 12 strains were then exposed to 1xMIC ciprofloxacin treatment (relative to the wild-type phenotype for each isogenic pair) for 1 h. Cell pellets were collected, and proteins were digested into peptides using trypsin. Protein identification and label-free quantification were done by liquid chromatography-mass spectrometry (LC-MS) in order to identify proteins that were differentially expressed upon deletion of recA in each strain. Data analysis and statistical analysis were performed using the MaxQuant and Perseus software.

Results

The proteins with the lowest expression levels were: RecA (as control), AphA, CysP, DinG, DinI, GarL, PriS, PsuG, PsuK, RpsQ, UgpB and YebG; those with the highest expression levels were: Hpf, IbpB, TufB and RpmH. Most of these expression alterations were strain-dependent and involved DNA repair processes and nucleotide, protein and carbohydrate metabolism, and transport. In isolates with suppressed SOS response, the number of underexpressed proteins was higher than overexpressed proteins.

Conclusion

High genomic and proteomic variability was observed among clinical isolates and was not associated with a specific resistant phenotype. This study provides an interesting approach to identify new potential targets to combat antimicrobial resistance.

Understanding the Impact of Temperate Bacteriophages on Their Lysogens Through Transcriptomics

Temperate phages are found integrated as prophages in the majority of bacterial genomes. Some prophages are cryptic and fixed in the bacterial chromosome, but others are active and can be triggered into a replicative form either spontaneously or by exposure to inducing factors. Prophages are commonly associated with the ability to confer toxin production or other virulence-associated traits on their host cell. More recent studies have shown they can play a much bigger role in altering the physiology of their hosts. The technique described here has enabled us to investigate how prophages affect gene expression in the opportunistic bacterium Pseudomonas aeruginosa. In this work, the growth of the wild-type P. aeruginosa strain PAO1 was compared with that of isogenic lysogens carrying different combinations of prophages from the Liverpool Epidemic Strain (LES) LESB58. In a lysogen culture, a proportion of bacterial cells will be supporting lytic bacteriophage replication (spontaneous induction) with a high level of expression per cell of late phage genes, such as those associated with the assembly of phage particles, thus masking the low-level gene expression associated with lysogen-restricted gene expression. The impact of spontaneous induction can thus obscure prophage gene expression across a lysogen population. Growth profiling experiments were used to identify spontaneous induction, which was minimal during the early exponential growth phase. This study reports how to prepare sample cultures during the early exponential growth phase and how to set up adequate controls despite low cell numbers. These protocols ensure the reliable and reproducible comparison of wild-type and lysogenic bacteria under various conditions, thus improving the transcriptomic profiling of prophage genomes and aiding in the identification of previously unrecognized prophage functions.

Mechanisms of low susceptibility to the disinfectant benzalkonium chloride in a multidrug-resistant environmental isolate of Aeromonas hydrophila

Excessive discharge of quaternary ammonium disinfectants such as benzalkonium chloride (BAC) into aquatic systems can trigger several physiological responses in environmental microorganisms. In this study, we isolated a less-susceptible strain of Aeromonas hydrophila to BAC, designated as INISA09, from a wastewater treatment plant in Costa Rica. We characterized its phenotypic response upon exposure to three different concentrations of BAC and characterized mechanisms related to its resistance using genomic and proteomic approaches. The genome of the strain, mapped against 52 different sequenced A. hydrophila strains, consists of approximately 4.6 Mb with 4,273 genes. We found a massive genome rearrangement and thousands of missense mutations compared to the reference strain A. hydrophila ATCC 7966. We identified 15,762 missense mutations mainly associated with transport, antimicrobial resistance, and outer membrane proteins. In addition, a quantitative proteomic analysis revealed a significant upregulation of several efflux pumps and the downregulation of porins when the strain was exposed to three BAC concentrations. Other genes related to membrane fatty acid metabolism and redox metabolic reactions also showed an altered expression. Our findings indicate that the response of A. hydrophila INISA09 to BAC primarily occurs at the envelop level, which is the primary target of BAC. Our study elucidates the mechanisms of antimicrobial susceptibility in aquatic environments against a widely used disinfectant and will help better understand how bacteria can adapt to biocide pollution. To our knowledge, this is the first study addressing the resistance to BAC in an environmental A. hydrophila isolate. We propose that this bacterial species could also serve as a new model to study antimicrobial pollution in aquatic environments.

CRISPRi screen for enhancing heterologous α-amylase yield in Bacillus subtilis

Yield improvements in cell factories can potentially be obtained by fine-tuning the regulatory mechanisms for gene candidates. In pursuit of such candidates, we performed RNA-sequencing of two α-amylase producing Bacillus strains and predict hundreds of putative novel non-coding transcribed regions. Surprisingly, we found among hundreds of non-coding and structured RNA candidates that non-coding genomic regions are proportionally undergoing the highest changes in expression during fermentation. Since these classes of RNA are also understudied, we targeted the corresponding genomic regions with CRIPSRi knockdown to test for any potential impact on the yield. From differentially expression analysis, we selected 53 non-coding candidates. Although CRISPRi knockdowns target both the sense and the antisense strand, the CRISPRi experiment cannot link causes for yield changes to the sense or antisense disruption. Nevertheless, we observed on several instances with strong changes in enzyme yield. The knockdown targeting the genomic region for a putative antisense RNA of the 3' UTR of the skfA-skfH operon led to a 21% increase in yield. In contrast, the knockdown targeting the genomic regions of putative antisense RNAs of the cytochrome c oxidase subunit 1 (ctaD), the sigma factor sigH, and the uncharacterized gene yhfT decreased yields by 31 to 43%.

The impact of PrsA over-expression on the Bacillus subtilis transcriptome during fed-batch fermentation of alpha-amylase production

The production of the alpha-amylase (AMY) enzyme in Bacillus subtilis at a high rate leads to the accumulation of unfolded AMY, which causes secretion stress. The over-expression of the PrsA chaperone aids enzyme folding and reduces stress. To identify affected pathways and potential mechanisms involved in the reduced growth, we analyzed the transcriptomic differences during fed-batch fermentation between a PrsA over-expressing strain and control in a time-series RNA-seq experiment. We observe transcription in 542 unannotated regions, of which 234 had significant changes in expression levels between the samples. Moreover, 1,791 protein-coding sequences, 80 non-coding genes, and 20 riboswitches overlapping UTR regions of coding genes had significant changes in expression. We identified putatively regulated biological processes via gene-set over-representation analysis of the differentially expressed genes; overall, the analysis suggests that the PrsA over-expression affects ATP biosynthesis activity, amino acid metabolism, and cell wall stability. The investigation of the protein interaction network points to a potential impact on cell motility signaling. We discuss the impact of these highlighted mechanisms for reducing secretion stress or detrimental aspects of PrsA over-expression during AMY production.

The Bacillaceae-1 RNA motif comprises two distinct classes

Non-coding RNAs are key regulatory players in bacteria. Many computationally predicted non-coding RNAs, however, lack functional associations. An example is the Bacillaceae-1 RNA motif, whose Rfam model consists of two hairpin loops. We find the motif conserved in nine of 13 non-pathogenic strains of the genus Bacillus but only in one pathogenic strain. To elucidate functional characteristics, we studied 118 hits of the Rfam model in 11 Bacillus spp. and found two distinct classes based on the ensemble diversity of their RNA secondary structure and the genomic context concerning the ribosomal RNA (rRNA) cluster. Forty hits are associated with the rRNA cluster, of which all 19 hits upstream flanking of 16S rRNA have a reverse complementary structure of low structural diversity. Fifty-two hits have large ensemble diversity, of which 38 are located between two coding genes. For eight hits in Bacillus subtilis, we investigated public expression data under various conditions and observed either the forward or the reverse complementary motif expressed. Five hits are associated with the rRNA cluster. Four of them are located upstream of the 16S rRNA and are not transcriptionally active, but instead, their reverse complements with low structural diversity are expressed together with the rRNA cluster. The three other hits are located between two coding genes in non-conserved genomic loci. Two of them are independently expressed from their surrounding genes and are structurally diverse. In summary, we found that Bacillaceae-1 RNA motifs upstream flanking of ribosomal RNA clusters tend to have one stable structure with the reverse complementary motif expressed in B. subtilis. In contrast, a subgroup of intergenic motifs has the thermodynamic potential for structural switches.

Flagella disruption in Bacillus subtilis increases amylase production yield

BACKGROUND

Bacillus subtilis is a Gram-positive bacterium used as a cell factory for protein production. Over the last decades, the continued optimization of production strains has increased yields of enzymes, such as amylases, and made commercial applications feasible. However, current yields are still significantly lower than the theoretically possible yield based on the available carbon sources. In its natural environment, B. subtilis can respond to unfavorable growth conditions by differentiating into motile cells that use flagella to swim towards available nutrients.

RESULTS

In this study, we analyze existing transcriptome data from a B. subtilis α-amylase production strain at different time points during a 5-day fermentation. We observe that genes of the fla/che operon, essential for flagella assembly and motility, are differentially expressed over time. To investigate whether expression of the flagella operon affects yield, we performed CRISPR-dCas9 based knockdown of the fla/che operon with sgRNA target against the genes flgE, fliR, and flhG, respectively. The knockdown resulted in inhibition of mobility and a striking 2-threefold increase in α-amylase production yield. Moreover, replacing flgE (required for flagella hook assembly) with an erythromycin resistance gene followed by a transcription terminator increased α-amylase yield by about 30%. Transcript levels of the α-amylase were unaltered in the CRISPR-dCas9 knockdowns as well as the flgE deletion strain, but all manipulations disrupted the ability of cells to swim on agar.

CONCLUSIONS

We demonstrate that the disruption of flagella in a B. subtilis α-amylase production strain, either by CRISPR-dCas9-based knockdown of the operon or by replacing flgE with an erythromycin resistance gene followed by a transcription terminator, increases the production of α-amylase in small-scale fermentation.

The International Virus Bioinformatics Meeting 2022

The International Virus Bioinformatics Meeting 2022 took place online, on 23-25 March 2022, and has attracted about 380 participants from all over the world. The goal of the meeting was to provide a meaningful and interactive scientific environment to promote discussion and collaboration and to inspire and suggest new research directions and questions. The participants created a highly interactive scientific environment even without physical face-to-face interactions. This meeting is a focal point to gain an insight into the state-of-the-art of the virus bioinformatics research landscape and to interact with researchers in the forefront as well as aspiring young scientists. The meeting featured eight invited and 18 contributed talks in eight sessions on three days, as well as 52 posters, which were presented during three virtual poster sessions. The main topics were: SARS-CoV-2, viral emergence and surveillance, virus-host interactions, viral sequence analysis, virus identification and annotation, phages, and viral diversity. This report summarizes the main research findings and highlights presented at the meeting.

BSGatlas: a unified Bacillus subtilis genome and transcriptome annotation atlas with enhanced information access

A large part of our current understanding of gene regulation in Gram-positive bacteria is based on Bacillus subtilis, as it is one of the most well studied bacterial model systems. The rapid growth in data concerning its molecular and genomic biology is distributed across multiple annotation resources. Consequently, the interpretation of data from further B. subtilis experiments becomes increasingly challenging in both low- and large-scale analyses. Additionally, B. subtilis annotation of structured RNA and non-coding RNA (ncRNA), as well as the operon structure, is still lagging behind the annotation of the coding sequences. To address these challenges, we created the B. subtilis genome atlas, BSGatlas, which integrates and unifies multiple existing annotation resources. Compared to any of the individual resources, the BSGatlas contains twice as many ncRNAs, while improving the positional annotation for 70 % of the ncRNAs. Furthermore, we combined known transcription start and termination sites with lists of known co-transcribed gene sets to create a comprehensive transcript map. The combination with transcription start/termination site annotations resulted in 717 new sets of co-transcribed genes and 5335 untranslated regions (UTRs). In comparison to existing resources, the number of 5' and 3' UTRs increased nearly fivefold, and the number of internal UTRs doubled. The transcript map is organized in 2266 operons, which provides transcriptional annotation for 92 % of all genes in the genome compared to the at most 82 % by previous resources. We predicted an off-target-aware genome-wide library of CRISPR-Cas9 guide RNAs, which we also linked to polycistronic operons. We provide the BSGatlas in multiple forms: as a website (https://rth.dk/resources/bsgatlas/), an annotation hub for display in the UCSC genome browser, supplementary tables and standardized GFF3 format, which can be used in large scale -omics studies. By complementing existing resources, the BSGatlas supports analyses of the B. subtilis genome and its molecular biology with respect to not only non-coding genes but also genome-wide transcriptional relationships of all genes.

Seawater salt-trapped Pseudomonas aeruginosa survives for years and gets primed for salinity tolerance

Background

In nature, microorganisms have to adapt to long-term stressful conditions often with growth limitations. However, little is known about the evolution of the adaptability of new bacteria to such environments. Pseudomonas aeruginosa, an opportunistic pathogen, after natural evaporation of seawater, was shown to be trapped in laboratory-grown halite crystals and to remain viable after entrapment for years. However, how this bacterium persists and survives in such hypersaline conditions is not understood.

Results

In this study, we aimed to understand the basis of survival, and to characterise the physiological changes required to develop salt tolerance using P. aeruginosa as a model. Several clones of P. aeruginosa were rescued after 14 years in naturally evaporated marine salt crystals. Incubation of samples in nutrient-rich broth allowed re-growth and subsequent plating yielded observable colonies. Whole genome sequencing of the P. aeruginosa isolates confirmed the recovery of the original strain. The re-grown strains, however, showed a new phenotype consisting of an enhanced growth in growing salt concentration compared to the ancestor strain. The intracellular accumulation of K⁺ was elicited by high concentration of Na⁺ in the external medium to maintain the homeostasis. Whole transcriptomic analysis by microarray indicated that 78 genes had differential expression between the parental strain and its derivative clones. Sixty-one transcripts were up-regulated, while 17 were down-regulated. Based on a collection of single-gene knockout mutants and gene ontology analysis, we suggest that the adaptive response in P. aeruginosa to hyper-salinity relies on multiple gene product interactions.

Conclusions

The individual gene contributions build up the observed phenotype, but do not ease the identification of salinity-related metabolic pathways. The long-term inclusion of P. aeruginosa in salt crystals primes the bacteria, mediating a readjustment of the bacterial physiology to growth in higher salt concentrations. Our findings provide a starting point to understand how P. aeruginosa, a relevant environmental and pathogenic bacterium, survives to long-term salt stress.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1499-2) contains supplementary material, which is available to authorized users.

Clay-induced DNA breaks as a path for genetic diversity, antibiotic resistance, and asbestos carcinogenesis

Natural clays and synthetic nanofibres can have a severe impact on human health. After several decades of research, the molecular mechanism of how asbestos induces cancer is not well understood. Different fibres, including asbestos, can penetrate cell membranes and introduce foreign DNA in bacterial and eukaryotic cells. Incubating Escherichia coli under friction forces with sepiolite, a clayey material, or with asbestos, causes double-strand DNA breaks. Antibiotics and clays are used together in animal husbandry, the mutagenic effect of these fibres could be a pathway to antibiotic resistance due to the friction provided by peristalsis of the gut from farm animals in addition to horizontal gene transfer. Moreover, we raise the possibility that the same mechanism could generate bacteria diversity in natural scenarios, playing a role in the evolution of species. Finally, we provide a new model on how asbestos may promote mutagenesis and cancer based on the observed mechanical genotoxicity.

Genetic diversity of two Daphnia-infecting microsporidian parasites, based on sequence variation in the internal transcribed spacer region

Background

Microsporidia are spore-forming obligate intracellular parasites that include both emerging pathogens and economically important disease agents. However, little is known about the genetic diversity of microsporidia. Here, we investigated patterns of geographic population structure, intraspecific genetic variation, and recombination in two microsporidian taxa that commonly infect cladocerans of the Daphnia longispina complex in central Europe. Taken together, this information helps elucidate the reproductive mode and life-cycles of these parasite species.

Methods

Microsporidia-infected Daphnia were sampled from seven drinking water reservoirs in the Czech Republic. Two microsporidia species (Berwaldia schaefernai and microsporidium lineage MIC1) were sequenced at the internal transcribed spacer (ITS) region, using the 454 pyrosequencing platform. Geographical structure analyses were performed applying Fisher’s exact tests, analyses of molecular variance, and permutational MANOVA. To evaluate the genetic diversity of the ITS region, the number of polymorphic sites and Tajima’s and Watterson’s estimators of theta were calculated. Tajima’s D was also used to determine if the ITS in these taxa evolved neutrally. Finally, neighbour similarity score and pairwise homology index tests were performed to detect recombination events.

Results

While there was little variation among Berwaldia parasite strains infecting different host populations, the among-population genetic variation of MIC1 was significant. Likewise, ITS genetic diversity was lower in Berwaldia than in MIC1. Recombination signals were detected only in Berwaldia.

Conclusion

Genetic tests showed that parasite populations could have expanded recently after a bottleneck or that the ITS could be under negative selection in both microsporidia species. Recombination analyses might indicate cryptic sex in Berwaldia and pure asexuality in MIC1. The differences observed between the two microsporidian species present an exciting opportunity to study the genetic basis of microsporidia-Daphnia coevolution in natural populations, and to better understand reproduction in these parasites.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1584-4) contains supplementary material, which is available to authorized users.

Daphnia parasite dynamics across multiple Caullerya epidemics indicate selection against common parasite genotypes

Studies of parasite population dynamics in natural systems are crucial for our understanding of host-parasite coevolutionary processes. Some field studies have reported that host genotype frequencies in natural populations change over time according to parasite-driven negative frequency-dependent selection. However, the temporal patterns of parasite genotypes have rarely been investigated. Moreover, parasite-driven negative frequency-dependent selection is contingent on the existence of genetic specificity between hosts and parasites. In the present study, the population dynamics and host-genotype specificity of the ichthyosporean Caullerya mesnili, a common endoparasite of Daphnia water fleas, were analysed based on the observed sequence variation in the first internal transcribed spacer (ITS1) of the ribosomal DNA. The Daphnia population of lake Greifensee (Switzerland) was sampled and subjected to parasite screening and host genotyping during C. mesnili epidemics of four consecutive years. The ITS1 of wild-caught C. mesnili-infected Daphnia was sequenced using the 454 pyrosequencing platform. The relative frequencies of C. mesnili ITS1 sequences differed significantly among years: the most abundant C. mesnili ITS1 sequence decreased and rare sequences increased over the course of the study, a pattern consistent with negative frequency-dependent selection. However, only a weak signal of host-genotype specificity between C. mesnili and Daphnia genotypes was detected. Use of cutting edge genomic techniques will allow further investigation of the underlying micro-evolutionary relationships within the Daphnia-C. mesnili system.

The Quantification of Representative Sequences pipeline for amplicon sequencing: case study on within-population ITS1 sequence variation in a microparasite infecting Daphnia

Next generation sequencing (NGS) platforms are replacing traditional molecular biology protocols like cloning and Sanger sequencing. However, accuracy of NGS platforms has rarely been measured when quantifying relative frequencies of genotypes or taxa within populations. Here we developed a new bioinformatic pipeline (QRS) that pools similar sequence variants and estimates their frequencies in NGS data sets from populations or communities. We tested whether the estimated frequency of representative sequences, generated by 454 amplicon sequencing, differs significantly from that obtained by Sanger sequencing of cloned PCR products. This was performed by analysing sequence variation of the highly variable first internal transcribed spacer (ITS1) of the ichthyosporean Caullerya mesnili, a microparasite of cladocerans of the genus Daphnia. This analysis also serves as a case example of the usage of this pipeline to study within-population variation. Additionally, a public Illumina data set was used to validate the pipeline on community-level data. Overall, there was a good correspondence in absolute frequencies of C. mesnili ITS1 sequences obtained from Sanger and 454 platforms. Furthermore, analyses of molecular variance (amova) revealed that population structure of C. mesnili differs across lakes and years independently of the sequencing platform. Our results support not only the usefulness of amplicon sequencing data for studies of within-population structure but also the successful application of the QRS pipeline on Illumina-generated data. The QRS pipeline is freely available together with its documentation under GNU Public Licence version 3 at http://code.google.com/p/quantification-representative-sequences.

Chromosomal replication dynamics and interaction with the β sliding clamp determine orientation of bacterial transposable elements

Insertion sequences (ISs) are small transposable elements widespread in bacterial genomes, where they play an essential role in chromosome evolution by stimulating recombination and genetic flow. Despite their ubiquity, it is unclear how ISs interact with the host. Here, we report a survey of the orientation patterns of ISs in bacterial chromosomes with the objective of gaining insight into the interplay between ISs and host chromosomal functions. We find that a significant fraction of IS families present a consistent and family-specific orientation bias with respect to chromosomal DNA replication, especially in Firmicutes. Additionally, we find that the transposases of up to nine different IS families with different transposition pathways interact with the β sliding clamp, an essential replication factor, suggesting that this is a widespread mechanism of interaction with the host. Although we find evidence that the interaction with the β sliding clamp is common to all bacterial phyla, it also could explain the observed strong orientation bias found in Firmicutes, because in this group β is asymmetrically distributed during synthesis of the leading or lagging strands. Besides the interaction with the β sliding clamp, other asymmetries also play a role in the biased orientation of some IS families. The utilization of the highly conserved replication sliding clamps suggests a mechanism for host regulation of IS proliferation and also a universal platform for IS dispersal and transmission within bacterial populations and among phylogenetically distant species.