RNA atlas of human bacterial pathogens uncovers stress dynamics linked to infection

Regulation, structure, and activity of the Pseudomonas aeruginosa MexXY efflux system

The current crisis in bacterial antibiotic resistance can be attributed to the overuse (or misuse) of these essential medicines in healthcare and agriculture, coupled with the slowed progression of new drug development. In the versatile, opportunistic pathogen Pseudomonas aeruginosa, the Resistance-Nodulation-Division (RND) efflux pump MexXY plays critical roles in both cell physiology and the acquisition of multidrug resistance. The mexXY operon is not constitutively expressed, but this process is instead controlled by a complex network of multiple interconnected regulatory mechanisms. These include induction by several of the pump's ribosome-targeting antibiotic substrates and transcriptional repression and anti-repression processes that are themselves influenced by various cellular factors, processes, or stresses. Although extensive studies of the MexXY complex are currently lacking as compared to other RND efflux pumps such as Escherichia coli AcrAB-TolC, recent studies have provided valuable insights into the MexXY architecture and substrate profiles, including its contribution to clinical resistance. Furthermore, while MexXY primarily associates with the outer membrane protein OprM, emerging evidence suggests that this transporter-periplasmic adaptor pair may also partner with other outer membrane proteins, potentially to alter the efflux substrate profile and activity under specific environmental conditions. In this minireview, we summarize current understanding of MexXY regulation, structure, and substrate selectivity within the context of clinical resistance and as a framework for future efflux pump inhibitor development.

The hidden bacterial microproteome

Microproteins encoded by small open reading frames comprise the "dark matter" of proteomes. Although microproteins have been detected in diverse organisms from all three domains of life, many more remain to be identified, and only a few have been functionally characterized. In this comprehensive study of intergenic small open reading frames (ismORFs, 15-70 codons) in 5,668 bacterial genomes of the family Enterobacteriaceae, we identify 67,297 clusters of ismORFs subject to purifying selection. Expression of tagged Escherichia coli microproteins is detected for 11 of the 16 tested, validating the predictions. Although the ismORFs mainly code for hydrophobic, potentially transmembrane, unstructured, or minimally structured microproteins, some globular folds, oligomeric structures, and possible interactions with proteins encoded by neighboring genes are predicted. Complete information on the predicted microprotein families, including evidence of transcription and translation, and structure predictions are available as an easily searchable resource for investigation of microprotein functions.

Deciphering the Molecular Mechanism of Peracetic Acid Response in Listeria monocytogenes

Peracetic acid (PAA), a strong oxidizing agent, has been widely used as a disinfectant in food processing settings as it does not produce harmful chlorinated by-products. In the present study, the transcriptional response of Listeria monocytogenes to a sub-lethal concentration of PAA (2.5 ppm) was assessed using RNA-sequencing (RNA-seq). Our analysis revealed 12 differentially expressed protein-coding genes, of which nine were upregulated (ohrR, ohrA, rpsN, lmo0637, lmo1973, fur, lmo2492, zurM, and lmo1007), and three were down-regulated (argG, lmo0604 and lmo2156) in PAA-treated samples compared to the control samples. A non-coding small RNA gene (rli32) was also found to be down-regulated. In detail, the organic peroxide toxicity protection (OhrA-OhrR) system, the metal homeostasis genes fur and zurM, the SbrE-regulated lmo0636-lmo0637 operon and a carbohydrate phosphotransferase system (PTS) operon component were induced under exposure of L. monocytogenes to PAA. Hence, this study identified key elements involved in the primary response of L. monocytogenes to oxidative stress caused by PAA, including the expression of the peroxide detoxification system and fine-tuning the levels of redox-active metals in the cell. The investigation of the molecular mechanism of PAA response in L. monocytogenes is of utmost importance for the food industry, as residual PAA can lead to stress tolerance in pathogens.

Phenotypic and transcriptomic changes in Acinetobacter baumannii in rich and minimal growth mediums

Acinetobacter baumannii is an opportunistic pathogen that is often studied in commonly used rich media in laboratories worldwide. Due to the metabolic versatility of A. baumannii, it can be cultured in different growth mediums; however, this can lead to genotypic and phenotypic variations. In this study, we compared phenotypic and transcriptomic changes in A. baumannii ATCC17978-VU cultured in M9 minimal media supplemented with 20 mmol/L sodium succinate and rich lysogeny broth media. Phenotypically, growth was significantly slowed, virulence in Galleria mellonella was attenuated, and susceptibility to a variety of antibiotic classes was reduced when A. baumannii ATCC17978-VU was grown in minimal media versus rich media. Transcriptomic analysis showed differential regulation of >700 genes-including those associated with energy production and ribosomal function-when the two growth conditions were compared, with the majority of the upregulated genes seen in minimal media of unknown function. This study showed that culture media has a profound effect on the phenotype and cellular workings of a bacteria, highlighting the need for more studies of pathogens like A. baumannii ATCC17978-VU in minimal media.

A Universal Method for Fingerprinting Multiplexed Bacteria: Evolving Pruned Sensor Arrays via Machine Learning-Driven Combinatorial Group-Specificity Strategy

Array-based sensing technology holds immense potential for discerning the intricacies of biological systems. Nevertheless, developing a universal strategy for simultaneous identification of diverse types of multianalytes and meeting the diagnostic needs of a range of multiclassified clinical diseases poses substantial challenges. Herein, we introduce a combination method for constructing sensor arrays by assembling two types of group-specific elements. Such a method enables the rapid generation of a library of 100 sensing units, each with dual bacterial targeting capabilities. By employing a three-step screening strategy optimized by machine learning algorithms, various optimal five-element arrays were rapidly obtained for diverse clinical infectious models. Moreover, the pruned arrays successfully identified disparate mixing ratios and quantitative detection of clinically prevalent bacterial strains. Optimized through nine multiclassification algorithms, the top-performing multilayer perceptron (MLP) model demonstrated impressive recognition capabilities, achieving 100% accuracy for diagnosing clinical urinary tract infection (UTI) and 99.4% accuracy for clinical sepsis detection in the test models we collected. Such a combinatorial library construction and screening process should be standard and provides insights into successfully generating powerful high-recognition sensor elements and configuring them into highly discriminative mini-sensor arrays.

Mycobacterial β-carbonic anhydrases: Molecular biology, role in the pathogenesis of tuberculosis and inhibition studies

Mycobacterium tuberculosis (Mtb), which causes tuberculosis (TB), is still a major global health problem. According to the World Health Organization (WHO), TB still causes more deaths worldwide than any other infectious agent. Drug-sensitive TB is treatable using first-line drugs; treatment of multidrug-resistant (MDR) and extensively drug-resistant (XDR) TB requires second- and third-line drugs. However, due to the long duration of treatment, the noncompliance of patients with different levels of resistance of Mtb to these drugs has worsened the situation. Previously developed anti-TB drugs targeted the replication machinery, protein synthesis, and cell wall biosynthesis pathways of Mtb. Therefore, novel drugs targeting alternate pathways crucial for the survival and pathogenesis of Mtb in the human host are needed. The genome of Mtb encodes three β-carbonic anhydrases (CAs) that are fundamental for pH homeostasis, hypoxia, survival, and pathogenesis. Recently, several studies have shown that the β-CAs of Mtb could be inhibited both in vitro and in vivo using small chemical molecules, suggesting that these enzymes could be novel targets for developing anti-TB compounds that are devoid of resistance by Mtb. In addition, homologs of β-CAs are absent in humans; therefore, drugs developed to target these enzymes might have minimal off-target effects. In this work, we describe the roles of β-CAs in Mtb and discuss bioinformatics and cheminformatics tools used in development and discovery of novel inhibitors of these enzymes. In addition, we summarize the in vitro and in vivo studies demonstrating that the β-CAs of Mtb are indeed druggable targets.

Unraveling bacterial stress responses: implications for next-generation antimicrobial solutions

The accelerated spread of antimicrobial-resistant bacteria has caused a serious health problem and rendered antimicrobial treatments ineffective. Innovative approaches are crucial to overcome the health threat posed by resistant pathogens and prevent the emergence of untreatable infections. Triggering stress responses in bacteria can diminish susceptibility to various antimicrobials by inducing resistance mechanisms. Therefore, a thorough understanding of stress response control, especially in relation to antimicrobial resistance, offers valuable perspectives for innovative and efficient therapeutic approaches to combat antimicrobial resistance. The aim of this study was to evaluate the stress responses of 8 different bacteria by analyzing reporter metabolites, around which significant alterations were observed, using a pathway-driven computational approach. For this purpose, the transcriptomic data that the bacterial pathogens were grown under 11 different stress conditions mimicking the human host environments were integrated with the genome-scale metabolic models of 8 pathogenic species (Enterococcus faecalis OG1R, Escherichia coli EPEC O127:H6 E2348/69, Escherichia coli ETEC H10407, Escherichia coli UPEC 536, Klebsiella pneumoniae MGH 78578, Pseudomonas aeruginosa PAO1, Staphylococcus aureus MRSA252, and Staphylococcus aureus MSSA476). The resulting reporter metabolites were enriched in multiple metabolic pathways, with cofactor biosynthesis being the most important. The results of this study will serve as a guide for the development of antimicrobial agents as they provide a first insight into potential drug targets.

Non-interchangeable functions of efflux transporters of Pseudomonas aeruginosa in survival under infection-associated stress

Pseudomonas aeruginosa is a challenging opportunistic pathogen due to its intrinsic and acquired mechanisms of antibiotic resistance. A large repertoire of efflux transporters actively expels antibiotics, toxins, and metabolites from cells and enables growth of P. aeruginosa in diverse environments. In this study, we analyzed the roles of representative efflux pumps from the Resistance-Nodulation-Division (RND), Major Facilitator Superfamily (MFS), and Small Multidrug Resistance (SMR) families of proteins in the susceptibility of P. aeruginosa to antibiotics and bacterial growth under stresses imposed by human hosts during bacterial infections: an elevated temperature, osmotic stress, low iron, bile salts, and acidic pH. We selected five RND pumps MexAB-OprM, MexEF-OprN, MexCD-OprJ, MuxABC-OpmB, and TriABC-OpmH that differ in their substrate specificities and expression profiles, two MFS efflux pumps PA3136-3137 and PA5158-5160 renamed here into MfsAB and MfsCD-OpmG, respectively, and an SMR efflux transporter PA1540-1541 (MdtJI). We found that the most promiscuous RND pumps such as MexEF-OprN and MexAB-OprM are integrated into diverse survival mechanisms and enable P. aeruginosa growth under various stresses. MuxABC-OpmB and TriABC-OpmH pumps with narrower substrate spectra are beneficial only in the presence of the iron chelator 2,2'-dipyridyl and bile salts, respectively. MFS pumps do not contribute to antibiotic efflux but play orthogonal roles in acidic pH, low iron, and in the presence of bile salts. In contrast, MdtJI protects against polycationic antibiotics but does not contribute to survival under stress. Thus, efflux pumps play specific, non-interchangeable functions in P. aeruginosa cell physiology and bacterial survival under stresses.

IMPORTANCE

The role of multidrug efflux pumps in the intrinsic and clinical levels of antibiotic resistance in Pseudomonas aeruginosa and other gram-negative bacteria is well-established. Their functions in bacterial physiology, however, remain unclear. The P. aeruginosa genome comprises an arsenal of efflux pumps from different protein families, the substrate specificities of which are typically assessed by measuring their impact on susceptibility to antibiotics. In this study, we analyzed how deletions and overproductions of efflux pumps affect P. aeruginosa growth under human-infection-induced stresses. Our results show that the physiological functions of multidrug efflux pumps are non-redundant and essential for the survival of this important human pathogen under stress.

Staphylococcus aureus NAD kinase is required for envelop and antibiotic stress responses

Global burden of infectious diseases and antimicrobial resistance are major public health issues calling for innovative control measures. Bacterial NAD kinase (NADK) is a crucial enzyme for production of NADP(H) and growth. In Staphylococcus aureus, NADK promotes pathogenesis by supporting production of key virulence determinants. Here, we find that knockdown of NADK by CRISPR interference sensitizes S. aureus to osmotic stress and to stresses induced by antibiotics targeting the envelop as well as replication, transcription and translation. Thus, NADK represents a promising target for the development of inhibitors which could be used in combination with current antibiotics.

Group A Streptococcus adaptation to diverse niches: lessons from transcriptomic studies

Group A Streptococcus (GAS) is a major human pathogen, causing diseases ranging from mild superficial infections of the skin and pharyngeal epithelium to severe systemic and invasive diseases. Moreover, post infection auto-immune sequelae arise by a yet not fully understood mechanism. The ability of GAS to cause a wide variety of infections is linked to the expression of a large set of virulence factors and their transcriptional regulation in response to various physiological environments. The use of transcriptomics, among others -omics technologies, in addition to traditional molecular methods, has led to a better understanding of GAS pathogenesis and host adaptation mechanisms. This review focusing on bacterial transcriptomic provides new insight into gene-expression patterns in vitro, ex vivo and in vivo with an emphasis on metabolic shifts, virulence genes expression and transcriptional regulators role.

Prediction of Essential Proteins of Klebsiella pneumoniae using Integrative Bioinformatics and Systems Biology Approach: Unveiling New Avenues for Drug Discovery

Klebsiella pneumoniae is an opportunistic multidrug-resistant bacterial pathogen responsible for various health care-associated infections. The prediction of proteins that are essential for the survival of bacterial pathogens can greatly facilitate the drug development and discovery pipeline toward target identification. To this end, the present study reports a comprehensive computational approach integrating bioinformatics and systems biology-based methods to identify essential proteins of K. pneumoniae involved in vital processes. From the proteome of this pathogen, we predicted a total of 854 essential proteins based on sequence, protein-protein interaction (PPI) and genome-scale metabolic model methods. These predicted essential proteins are involved in vital processes for cellular regulation such as translation, metabolism, and biosynthesis of essential factors, among others. Cluster analysis of the PPI network revealed the highly connected modules involved in the basic functionality of the organism. Further, the predicted consensus set of essential proteins of K. pneumoniae was evaluated by comparing them with existing resources (NetGenes and PATHOgenex) and literature. The findings of this study offer guidance toward understanding cell functionality, thereby facilitating the understanding of pathogen systems and providing a way forward to shortlist potential therapeutic candidates for developing novel antimicrobial agents against K. pneumoniae. In addition, the research strategy presented herein is a fusion of sequence and systems biology-based approaches that offers prospects as a model to predict essential proteins for other pathogens.

RNA cis-regulators are important for Streptococcus pneumoniae in vivo success

Bacteria have evolved complex transcriptional regulatory networks, as well as many diverse regulatory strategies at the RNA level, to enable more efficient use of metabolic resources and a rapid response to changing conditions. However, most RNA-based regulatory mechanisms are not well conserved across different bacterial species despite controlling genes important for virulence or essential biosynthetic processes. Here, we characterize the activity of, and assess the fitness benefit conferred by, twelve cis-acting regulatory RNAs (including several riboswitches and a T-box), in the opportunistic pathogen Streptococcus pneumoniae TIGR4. By evaluating native locus mutants of each regulator that result in constitutively active or repressed expression, we establish that growth defects in planktonic culture are associated with constitutive repression of gene expression, while constitutive activation of gene expression is rarely deleterious. In contrast, in mouse nasal carriage and pneumonia models, strains with either constitutively active and repressed gene expression are significantly less fit than matched control strains. Furthermore, two RNA-regulated pathways, FMN synthesis/transport and pyrimidine synthesis/transport display exceptional sensitivity to mis-regulation or constitutive gene repression in both planktonic culture and in vivo environments. Thus, despite lack of obvious phenotypes associated with constitutive gene expression in vitro, the fitness benefit conferred on bacteria via fine-tuned metabolic regulation through cis-acting regulatory RNAs is substantial in vivo, and therefore easily sufficient to drive the evolution and maintenance of diverse RNA regulatory mechanisms.

The Cryptic Bacterial Microproteome

Microproteins encoded by small open reading frames (smORFs) comprise the "dark matter" of proteomes. Although functional microproteins were identified in diverse organisms from all three domains of life, bacterial smORFs remain poorly characterized. In this comprehensive study of intergenic smORFs (ismORFs, 15-70 codons) in 5,668 bacterial genomes of the family Enterobacteriaceae, we identified 67,297 clusters of ismORFs subject to purifying selection. The ismORFs mainly code for hydrophobic, potentially transmembrane, unstructured, or minimally structured microproteins. Using AlphaFold Multimer, we predicted interactions of some of the predicted microproteins encoded by transcribed ismORFs with proteins encoded by neighboring genes, revealing the potential of microproteins to regulate the activity of various proteins, particularly, under stress. We compiled a catalog of predicted microprotein families with different levels of evidence from synteny analysis, structure prediction, and transcription and translation data. This study offers a resource for investigation of biological functions of microproteins.

HiTaxon: a hierarchical ensemble framework for taxonomic classification of short reads

Motivation

Whole microbiome DNA and RNA sequencing (metagenomics and metatranscriptomics) are pivotal to determining the functional roles of microbial communities. A key challenge in analyzing these complex datasets, typically composed of tens of millions of short reads, is accurately classifying reads to their taxa of origin. While still performing worse relative to reference-based short-read tools in species classification, ML algorithms have shown promising results in taxonomic classification at higher ranks. A recent approach exploited to enhance the performance of ML tools, which can be translated to reference-dependent classifiers, has been to integrate the hierarchical structure of taxonomy within the tool's predictive algorithm.

Results

Here, we introduce HiTaxon, an end-to-end hierarchical ensemble framework for taxonomic classification. HiTaxon facilitates data collection and processing, reference database construction and optional training of ML models to streamline ensemble creation. We show that databases created by HiTaxon improve the species-level performance of reference-dependent classifiers, while reducing their computational overhead. In addition, through exploring hierarchical methods for HiTaxon, we highlight that our custom approach to hierarchical ensembling improves species-level classification relative to traditional strategies. Finally, we demonstrate the improved performance of our hierarchical ensembles over current state-of-the-art classifiers in species classification using datasets comprised of either simulated or experimentally derived reads.

Availability and implementation

HiTaxon is available at: https://github.com/ParkinsonLab/HiTaxon.

Regulation of bacterial gene expression by non-coding RNA: It is all about time!

Commensal and pathogenic bacteria continuously evolve to survive in diverse ecological niches by efficiently coordinating gene expression levels in their ever-changing environments. Regulation through the RNA transcript itself offers a faster and more cost-effective way to adapt than protein-based mechanisms and can be leveraged for diagnostic or antimicrobial purposes. However, RNA can fold into numerous intricate, not always functional structures that both expand and obscure the plethora of roles that regulatory RNAs serve within the cell. Here, we review the current knowledge of bacterial non-coding RNAs in relation to their folding pathways and interactions. We posit that co-transcriptional folding of these transcripts ultimately dictates their downstream functions. Elucidating the spatiotemporal folding of non-coding RNAs during transcription therefore provides invaluable insights into bacterial pathogeneses and predictive disease diagnostics. Finally, we discuss the implications of co-transcriptional folding andapplications of RNAs for therapeutics and drug targets.

Co-PATHOgenex web application for assessing complex stress responses in pathogenic bacteria

IMPORTANCE

Unveiling gene co-expression networks in bacterial pathogens has the potential for gaining insights into their adaptive strategies within the host environment. Here, we developed Co-PATHOgenex, an interactive and user-friendly web application that enables users to construct networks from gene co-expressions using custom-defined thresholds (https://avicanlab.shinyapps.io/copathogenex/). The incorporated search functions and visualizations within the tool simplify the usage and facilitate the interpretation of the analysis output. Co-PATHOgenex also includes stress stimulons for various bacterial species, which can help identify gene products not previously associated with a particular stress condition.

Acinetobacter baumannii Survival under Infection-Associated Stresses Depends on the Expression of Resistance-Nodulation-Division and Major Facilitator Superfamily Efflux Pumps

Multidrug efflux transporters are major contributors to the antibiotic resistance of Acinetobacter baumannii in clinical settings. Previous studies showed that these transporters are tightly integrated into the physiology of A. baumannii and have diverse functions. However, for many of the efflux pumps, such functions remain poorly defined. In this study, we characterized two putative drug efflux pumps, AmfAB and AmfCD (Acinetobacter Major Facilitator), that are homologous to EmrAB-like transporters from Escherichia coli and other Gram-negative bacteria. These pumps comprise the Major Facilitator Superfamily (MFS) transporters AmfB and AmfD and the periplasmic membrane fusion proteins AmfA and AmfC, respectively. We inactivated and overproduced these pumps in the wild-type ATCC 17978 strain and its derivative strains lacking the major efflux pumps from the Resistance-Nodulation-Division (RND) superfamily and characterized antibiotic susceptibilities and growth of the strains under stresses typical during human infections. We found that neither AmfAB nor AmfCD contribute to the antibiotic non-susceptibility phenotypes of A. baumannii. The two pumps, however, are critical for the adaptation and growth of the bacterium under acidic stress, whereas AmfCD also contributes to growth under conditions of low iron, high temperature, and in the presence of bile salts. These functions are dependent on the presence of the RND pumps, the inactivation of which further diminishes A. baumannii survival and growth. Our results suggest that MFS transporters contribute to stress survival by affecting the permeability properties of the A. baumannii cell envelope.

Fitness and transcriptomic analysis of pathogenic Vibrio parahaemolyticus in seawater at different shellfish harvesting temperatures

IMPORTANCE

Given the involvement of Vibrio parahaemolyticus (Vp) in a wide range of seafood outbreaks, a systematical characterization of Vp fitness and transcriptomic changes at temperatures of critical importance for seafood production and storage is needed. In this study, one of each virulent Vp strain (tdh+ and trh+) was tested. While no difference in survival behavior of the two virulent strains was observed at 10°C, the tdh+ strain had a faster growth rate than the trh+ strain at 30°C. Transcriptomic analysis showed that a significantly higher number of genes were upregulated at 30°C than at 10°C. The majority of differentially expressed genes of Vp at 30°C were annotated to functional categories supporting cellular growth. At 10°C, the downregulation of the biofilm formation and histidine metabolism indicates that the current practice of storing seafood at low temperatures not only protects seafood quality but also ensures seafood safety.

The genomes of Scedosporium between environmental challenges and opportunism

Emerging fungal pathogens are a global challenge for humankind. Many efforts have been made to understand the mechanisms underlying pathogenicity in bacteria, and OMICs techniques are largely responsible for those advancements. By contrast, our limited understanding of opportunism and antifungal resistance is preventing us from identifying, limiting and interpreting the emergence of fungal pathogens. The genus Scedosporium (Microascaceae) includes fungi with high tolerance to environmental pollution, whilst some species can be considered major human pathogens, such as Scedosporium apiospermum and Scedosporium boydii. However, unlike other fungal pathogens, little is known about the genome evolution of these organisms. We sequenced two novel genomes of Scedosporium aurantiacum and Scedosporium minutisporum isolated from extreme, strongly anthropized environments. We compared all the available Scedosporium and Microascaceae genomes, that we systematically annotated and characterized ex novo in most cases. The genomes in this family were integrated in a Phylum-level comparison to infer the presence of putative, shared genomic traits in filamentous ascomycetes with pathogenic potential. The analysis included the genomes of 100 environmental and clinical fungi, revealing poor evolutionary convergence of putative pathogenicity traits. By contrast, several features in Microascaceae and Scedosporium were detected that might have a dual role in responding to environmental challenges and allowing colonization of the human body, including chitin, melanin and other cell wall related genes, proteases, glutaredoxins and magnesium transporters. We found these gene families to be impacted by expansions, orthologous transposon insertions, and point mutations. With RNA-seq, we demonstrated that most of these anciently impacted genomic features responded to the stress imposed by an antifungal compound (voriconazole) in the two environmental strains S. aurantiacum MUT6114 and S. minutisporum MUT6113. Therefore, the present genomics and transcriptomics investigation stands on the edge between stress resistance and pathogenic potential, to elucidate whether fungi were pre-adapted to infect humans. We highlight the strengths and limitations of genomics applied to opportunistic human pathogens, the multifactoriality of pathogenicity and resistance to drugs, and suggest a scenario where pressures other than anthropic contributed to forge filamentous human pathogens.

Desiccation induces varied responses within a soil bacterial genus

Desiccation impacts a suite of physiological processes in microbes by elevating levels of damaging reactive oxygen species and inducing DNA strand breaks. In response to desiccation-induced stress, microbes have evolved specialized mechanisms to help them survive. Here, we performed a 128-day lab desiccation experiment on nine strains from three clades of an abundant soil bacterium, Curtobacterium. We sequenced RNA from each strain at three time points to investigate their response. Curtobacterium was highly resistant to desiccation, outlasting both Escherichia coli and a famously DNA damage-resistant bacterium, Deinococcus radiodurans. However, within the genus, there were also 10-fold differences in survival rates among strains. Transcriptomic profiling revealed responses shared within the genus including up-regulation of genes involved in DNA damage repair, osmolyte production, and efflux pumps, but also up-regulation of pathways and genes unique to the three clades. For example, trehalose synthesis gene otsB, the chaperone groEL, and the oxygen scavenger katA were all found in either one or two clades but not the third. Here, we provide evidence of considerable variation in closely related strains, and further elucidation of the phylogenetic conservation of desiccation tolerance remains an important goal for microbial ecologists.

TRS: a method for determining transcript termini from RNAtag-seq sequencing data

In bacteria, determination of the 3' termini of transcripts plays an essential role in regulation of gene expression, affecting the functionality and stability of the transcript. Several experimental approaches were developed to identify the 3' termini of transcripts, however, these were applied only to a limited number of bacteria and growth conditions. Here we present a straightforward approach to identify 3' termini from widely available RNA-seq data without the need for additional experiments. Our approach relies on the observation that the RNAtag-seq sequencing protocol results in overabundance of reads mapped to transcript 3' termini. We present TRS (Termini by Read Starts), a computational pipeline exploiting this property to identify 3' termini in RNAtag-seq data, and show that the identified 3' termini are highly reliable. Since RNAtag-seq data are widely available for many bacteria and growth conditions, our approach paves the way for studying bacterial transcription termination in an unprecedented scope.

Models of Hfq interactions with small non-coding RNA in Gram-negative and Gram-positive bacteria

Hfq is required by many Gram-negative bacteria to chaperone the interaction between small non-coding RNA (sRNA) and mRNA to facilitate annealing. Conversely and despite the presence of Hfq in many Gram-positive bacteria, sRNAs in Gram-positive bacteria bind the mRNA target independent of Hfq. Details provided by the Hfq structures from both Gram-negative and Gram-positive bacteria have demonstrated that despite a conserved global structure of the protein, variations of residues on the binding surfaces of Hfq results in the recognition of different RNA sequences as well as the ability of Hfq to facilitate the annealing of the sRNA to the mRNA target. Additionally, a subset of Gram-negative bacteria has an extended C-terminal Domain (CTD) that has been shown to affect the stability of the Hfq hexamer and increase the rate of release of the annealed sRNA-mRNA product. Here we review the structures of Hfq and biochemical data that have defined the interactions of the Gram-negative and Gram-positive homologues to highlight the similarities and differences in the interactions with RNA. These interactions provided a deeper understanding of the how Hfq functions to facilitate the annealing of sRNA-mRNA, the selectivity of the interactions with RNA, and the role of the CTD of Hfq in the interactions with sRNA.

Microgravity and evasion of plant innate immunity by human bacterial pathogens

Spaceflight microgravity and modeled-microgravity analogs (MMA) broadly alter gene expression and physiology in both pathogens and plants. Research elucidating plant and bacterial responses to normal gravity or microgravity has shown the involvement of both physiological and molecular mechanisms. Under true and simulated microgravity, plants display differential expression of pathogen-defense genes while human bacterial pathogens exhibit increased virulence, antibiotic resistance, stress tolerance, and reduced LD50 in animal hosts. Human bacterial pathogens including Salmonella enterica and E. coli act as cross-kingdom foodborne pathogens by evading and suppressing the innate immunity of plants for colonization of intracellular spaces. It is unknown if evasion and colonization of plants by human pathogens occurs under microgravity and if there is increased infection capability as demonstrated using animal hosts. Understanding the relationship between microgravity, plant immunity, and human pathogens could prevent potentially deadly outbreaks of foodborne disease during spaceflight. This review will summarize (1) alterations to the virulency of human pathogens under microgravity and MMA, (2) alterations to plant physiology and gene expression under microgravity and MMA, (3) suppression and evasion of plant immunity by human pathogens under normal gravity, (4) studies of plant-microbe interactions under microgravity and MMA. A conclusion suggests future study of interactions between plants and human pathogens under microgravity is beneficial to human safety, and an investment in humanity's long and short-term space travel goals.

Heterogeneity of the group B streptococcal type VII secretion system and influence on colonization of the female genital tract

Type VIIb secretion systems (T7SSb) in Gram-positive bacteria facilitate physiology, interbacterial competition, and/or virulence via EssC ATPase-driven secretion of small ɑ-helical proteins and toxins. Recently, we characterized T7SSb in group B Streptococcus (GBS), a leading cause of infection in newborns and immunocompromised adults. GBS T7SS comprises four subtypes based on variation in the C-terminus of EssC and the repertoire of downstream effectors; however, the intraspecies diversity of GBS T7SS and impact on GBS-host interactions remains unknown. Bioinformatic analysis indicates that GBS T7SS loci encode subtype-specific putative effectors, which have low interspecies and inter-subtype homology but contain similar domains/motifs and therefore may serve similar functions. We further identify orphaned GBS WXG100 proteins. Functionally, we show that GBS T7SS subtype I and III strains secrete EsxA in vitro and that in subtype I strain CJB111, esxA1 appears to be differentially transcribed from the T7SS operon. Furthermore, we observe subtype-specific effects of GBS T7SS on host colonization, as CJB111 subtype I but not CNCTC 10/84 subtype III T7SS promotes GBS vaginal colonization. Finally, we observe that T7SS subtypes I and II are the predominant subtypes in clinical GBS isolates. This study highlights the potential impact of T7SS heterogeneity on host-GBS interactions.

Mechanobiology as a tool for addressing the genotype-to-phenotype problem in microbiology

The central hypothesis of the genotype-phenotype relationship is that the phenotype of a developing organism (i.e., its set of observable attributes) depends on its genome and the environment. However, as we learn more about the genetics and biochemistry of living systems, our understanding does not fully extend to the complex multiscale nature of how cells move, interact, and organize; this gap in understanding is referred to as the genotype-to-phenotype problem. The physics of soft matter sets the background on which living organisms evolved, and the cell environment is a strong determinant of cell phenotype. This inevitably leads to challenges as the full function of many genes, and the diversity of cellular behaviors cannot be assessed without wide screens of environmental conditions. Cellular mechanobiology is an emerging field that provides methodologies to understand how cells integrate chemical and physical environmental stress and signals, and how they are transduced to control cell function. Biofilm forming bacteria represent an attractive model because they are fast growing, genetically malleable and can display sophisticated self-organizing developmental behaviors similar to those found in higher organisms. Here, we propose mechanobiology as a new area of study in prokaryotic systems and describe its potential for unveiling new links between an organism's genome and phenome.

Single-cell transcriptomics and data analyses for prokaryotes-Past, present and future concepts

Transcriptomics, or more specifically mRNA sequencing, is a powerful tool to study gene expression at the single-cell level (scRNA-seq) which enables new insights into a plethora of biological processes. While methods for single-cell RNA-seq in eukaryotes are well established, application to prokaryotes is still challenging. Reasons for that are rigid and diverse cell wall structures hampering lysis, the lack of polyadenylated transcripts impeding mRNA enrichment, and minute amounts of RNA requiring amplification steps before sequencing. Despite those obstacles, several promising scRNA-seq approaches for bacteria have been published recently, albeit difficulties in the experimental workflow and data processing and analysis remain. In particular, bias is often introduced by amplification which makes it difficult to distinguish between technical noise and biological variation. Future optimization of experimental procedures and data analysis algorithms are needed for the improvement of scRNA-seq but also to aid in the emergence of prokaryotic single-cell multi-omics. to help address 21st century challenges in the biotechnology and health sector.

Phenotyping of Methicillin-Resistant Staphylococcus aureus Using a Ratiometric Sensor Array

Chemical tools capable of classifying multidrug-resistant bacteria (superbugs) can facilitate early-stage disease diagnosis and help guide precision therapy. Here, we report a sensor array that permits the facile phenotyping of methicillin-resistant Staphylococcus aureus (MRSA), a clinically common superbug. The array consists of a panel of eight separate ratiometric fluorescent probes that provide characteristic vibration-induced emission (VIE) profiles. These probes bear a pair of quaternary ammonium salts in different substitution positions around a known VIEgen core. The differences in the substituents result in varying interactions with the negatively charged cell walls of bacteria. This, in turn, dictates the molecular conformation of the probes and affects their blue-to-red fluorescence intensity ratios (ratiometric changes). Within the sensor array, the differences in the ratiometric changes for the probes result in "fingerprints" for MRSA of different genotypes. This allows them to be identified using principal component analysis (PCA) without the need for cell lysis and nucleic acid isolation. The results obtained with the present sensor array agree well with those obtained using polymerase chain reaction (PCR) analysis.

Impact of Temperature and Oxygen Availability on Gene Expression Patterns of Mycobacterium ulcerans

Buruli ulcer disease is a neglected tropical disease caused by the environmental pathogen Mycobacterium ulcerans. The M. ulcerans major virulence factor is mycolactone, a lipid cytotoxic compound whose genes are carried on a plasmid. Although an exact reservoir and mode(s) of transmission are unknown, data provide evidence of both. First, Buruli ulcer incidence and M. ulcerans presence have been linked to slow-moving water with low oxygen. M. ulcerans has also been suggested to be sensitive to UV due to termination in crtI, encoding a phytoene dehydrogenase, required for carotenoid production. Further, M. ulcerans has been shown to cause disease following puncture but not when introduced to open abrasion sites, suggesting that puncture is necessary for transmission and pathology. Despite these findings, the function and modulation of mycolactone and other genes in response to dynamic abiotic conditions such as UV, temperature, and oxygen have not been shown. In this study, we investigated modulation of mycolactone and other genes on exposure to changing UV and oxygen microenvironmental conditions. Mycolactone expression was downregulated on exposure to the single stress high temperature and did not change significantly with exposure to UV; however, it was upregulated when exposed to microaerophilic conditions. Mycolactone expression was downregulated under combined stresses of high temperature and low oxygen, but there was upregulation of several stress response genes. Taken together, results suggest that temperature shapes M. ulcerans metabolic response more so than UV exposure or oxygen requirements. These data help to define the environmental niche of M. ulcerans and metabolic responses during initial human infection. IMPORTANCE Buruli ulcer is a debilitating skin disease caused by the environmental pathogen Mycobacterium ulcerans. M. ulcerans produces a toxic compound, mycolactone, which leads to tissue necrosis and ulceration. Barriers to preventing Buruli ulcer include an incomplete understanding of M. ulcerans reservoirs, how the pathogen is transmitted, and under what circumstances mycolactone and other M. ulcerans genes are expressed and produced in its natural environment and in the host. We conducted a study to investigate M. ulcerans gene expression under several individual or combined abiotic conditions. Our data showed that mycolactone expression was downregulated under combined stresses of high temperature and low oxygen but there was upregulation of several stress response genes. These data are among only a few studies measuring modulation of mycolactone and other M. ulcerans genes that could be involved in pathogen fitness in its natural environment and virulence while within the host.

Yersiniomics, a Multi-Omics Interactive Database for Yersinia Species

The genus Yersinia includes a large variety of nonpathogenic and life-threatening pathogenic bacteria, which cause a broad spectrum of diseases in humans and animals, such as plague, enteritis, Far East scarlet-like fever (FESLF), and enteric redmouth disease. Like most clinically relevant microorganisms, Yersinia spp. are currently subjected to intense multi-omics investigations whose numbers have increased extensively in recent years, generating massive amounts of data useful for diagnostic and therapeutic developments. The lack of a simple and centralized way to exploit these data led us to design Yersiniomics, a web-based platform allowing straightforward analysis of Yersinia omics data. Yersiniomics contains a curated multi-omics database at its core, gathering 200 genomic, 317 transcriptomic, and 62 proteomic data sets for Yersinia species. It integrates genomic, transcriptomic, and proteomic browsers, a genome viewer, and a heatmap viewer to navigate within genomes and experimental conditions. For streamlined access to structural and functional properties, it directly links each gene to GenBank, the Kyoto Encyclopedia of Genes and Genomes (KEGG), UniProt, InterPro, IntAct, and the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and each experiment to Gene Expression Omnibus (GEO), the European Nucleotide Archive (ENA), or the Proteomics Identifications Database (PRIDE). Yersiniomics provides a powerful tool for microbiologists to assist with investigations ranging from specific gene studies to systems biology studies. IMPORTANCE The expanding genus Yersinia is composed of multiple nonpathogenic species and a few pathogenic species, including the deadly etiologic agent of plague, Yersinia pestis. In 2 decades, the number of genomic, transcriptomic, and proteomic studies on Yersinia grew massively, delivering a wealth of data. We developed Yersiniomics, an interactive web-based platform, to centralize and analyze omics data sets on Yersinia species. The platform allows user-friendly navigation between genomic data, expression data, and experimental conditions. Yersiniomics will be a valuable tool to microbiologists.

An integrated transcriptomics-functional genomics approach reveals a small RNA that modulates Bacteroides thetaiotaomicron sensitivity to tetracyclines

Gene expression plasticity allows bacteria to adapt to diverse environments, tie their metabolism to available nutrients, and cope with stress. This is particularly relevant in a niche as dynamic and hostile as the human intestinal tract, yet transcriptional networks remain largely unknown in gut Bacteroides spp. Here, we map transcriptional units and profile their expression levels in Bacteroides thetaiotaomicron over a suite of 15 defined experimental conditions that are relevant in vivo , such as variation of temperature, pH, and oxygen tension, exposure to antibiotic stress, and growth on simple carbohydrates or on host mucin-derived glycans. Thereby, we infer stress- and carbon source-specific transcriptional regulons, including conditional expression of capsular polysaccharides and polysaccharide utilization loci, and expand the annotation of small regulatory RNAs (sRNAs) in this organism. Integrating this comprehensive expression atlas with transposon mutant fitness data, we identify conditionally important sRNAs. One example is MasB, whose inactivation led to increased bacterial tolerance of tetracyclines. Using MS2 affinity purification coupled with RNA sequencing, we predict targets of this sRNA and discuss their potential role in the context of the MasB-associated phenotype. Together, this transcriptomic compendium in combination with functional sRNA genomics-publicly available through a new iteration of the 'Theta-Base' web browser (www.helmholtz-hiri.de/en/datasets/bacteroides-v2)-constitutes a valuable resource for the microbiome and sRNA research communities alike.

Heterogeneity of the group B streptococcal type VII secretion system and influence on colonization of the female genital tract

Type VIIb secretion systems (T7SSb) in Gram-positive bacteria facilitate physiology, interbacterial competition, and/or virulence via EssC ATPase-driven secretion of small ɑ-helical proteins and toxins. Recently, we characterized T7SSb in group B Streptococcus (GBS), a leading cause of infection in newborns and immunocompromised adults. GBS T7SS comprises four subtypes based on variation in the C-terminus of EssC and the repertoire of downstream effectors; however, the intra-species diversity of GBS T7SS and impact on GBS-host interactions remains unknown. Bioinformatic analysis indicates that GBS T7SS loci encode subtype-specific putative effectors, which have low inter-species and inter-subtype homology but contain similar domains/motifs and therefore may serve similar functions. We further identify orphaned GBS WXG100 proteins. Functionally, we show that GBS T7SS subtype I and III strains secrete EsxA in vitro and that in subtype I strain CJB111, esxA1 appears to be differentially transcribed from the T7SS operon. Further, we observe subtype-specific effects of GBS T7SS on host colonization, as subtype I but not subtype III T7SS promotes GBS vaginal persistence. Finally, we observe that T7SS subtypes I and II are the predominant subtypes in clinical GBS isolates. This study highlights the potential impact of T7SS heterogeneity on host-GBS interactions.

Making sense of drug-efflux transporters in the physiological environment

Bacterial drug-efflux transporters act synergistically with diffusion barriers of cellular membranes and other resistance mechanisms to protect cells from antibiotics and toxic metabolites. Their critical roles in clinical antibiotic and multidrug resistance are well established. In addition, a large body of evidence has been accumulated in support of their important contributions to bacterial growth and proliferation during infections. However, how these diverse functions of drug transporters are integrated at the level of bacterial cell physiology remains unclear. This opinion briefly summarizes the current understanding of substrate specificities and physiological roles of drug-efflux pumps from Resistance-Nodulation-Division (RND) superfamily of proteins in two ESKAPE pathogens Pseudomonas aeruginosa and Acinetobacter baumannii. Based on the analysis of phenotypic and transcriptomic studies in vitro and in vivo, we propose that RND pumps of Gram-negative bacteria fall into three categories: constitutively expressed, regulated, and silent. These three categories of efflux pumps participate in different physiological programs, which are not involved in the central metabolism and bacterial growth.

The impacts of viral infection and subsequent antimicrobials on the microbiome-resistome of growing pigs

BACKGROUND

Antimicrobials are used in food-producing animals for purposes of preventing, controlling, and/or treating infections. In swine, a major driver of antimicrobial use is porcine reproductive and respiratory syndrome (PRRS), which is caused by a virus that predisposes infected animals to secondary bacterial infections. Numerous antimicrobial protocols are used to treat PRRS, but we have little insight into how these treatment schemes impact antimicrobial resistance (AMR) dynamics within the fecal microbiome of commercial swine. The aim of this study was to determine whether different PRRS-relevant antimicrobial treatment protocols were associated with differences in the fecal microbiome and resistome of growing pigs. To accomplish this, we used a metagenomics approach to characterize and compare the longitudinal wean-to-market resistome and microbiome of pigs challenged with PRRS virus and then exposed to different antimicrobial treatments, and a group of control pigs not challenged with PRRS virus and having minimal antimicrobial exposure. Genomic DNA was extracted from pen-level composite fecal samples from each treatment group and subjected to metagenomic sequencing and microbiome-resistome bioinformatic and statistical analysis. Microbiome-resistome profiles were compared over time and between treatment groups.

RESULTS

Fecal microbiome and resistome compositions both changed significantly over time, with a dramatic and stereotypic shift between weaning and 9 days post-weaning (dpw). Antimicrobial resistance gene (ARG) richness and diversity were significantly higher at earlier time points, while microbiome richness and diversity were significantly lower. The post-weaning shift was characterized by transition from a Bacteroides-dominated enterotype to Lactobacillus- and Streptococcus-dominated enterotypes. Both the microbiome and resistome stabilized by 44 dpw, at which point the trajectory of microbiome-resistome maturation began to diverge slightly between the treatment groups, potentially due to physical clustering of the pigs. Challenge with PRRS virus seemed to correspond to the re-appearance of many very rare and low-abundance ARGs within the feces of challenged pigs. Despite very different antimicrobial exposures after challenge with PRRS virus, resistome composition remained largely similar between the treatment groups. Differences in ARG abundance between the groups were mostly driven by temporal changes in abundance that occurred prior to antimicrobial exposures, with the exception of ermG, which increased in the feces of treated pigs, and was significantly more abundant in the feces of these pigs compared to the pigs that did not receive post-PRRS antimicrobials.

CONCLUSIONS

The fecal microbiome-resistome of growing pigs exhibited a stereotypic trajectory driven largely by weaning and physiologic aging of the pigs. Events such as viral illness, antimicrobial exposures, and physical grouping of the pigs exerted significant yet relatively minor influence over this trajectory. Therefore, the AMR profile of market-age pigs is the culmination of the life history of the individual pigs and the populations to which they belong. Disease status alone may be a significant driver of AMR in market-age pigs, and understanding the interaction between disease processes and antimicrobial exposures on the swine microbiome-resistome is crucial to developing effective, robust, and reproducible interventions to control AMR. Video Abstract.

NAD kinase promotes Staphylococcus aureus pathogenesis by supporting production of virulence factors and protective enzymes

Nicotinamide adenine dinucleotide phosphate (NADPH) is the primary electron donor for reductive reactions that are essential for the biosynthesis of major cell components in all organisms. Nicotinamide adenine dinucleotide kinase (NADK) is the only enzyme that catalyzes the synthesis of NADP(H) from NAD(H). While the enzymatic properties and physiological functions of NADK have been thoroughly studied, the role of NADK in bacterial pathogenesis remains unknown. Here, we used CRISPR interference to knock down NADK gene expression to address the role of this enzyme in Staphylococcus aureus pathogenic potential. We find that NADK inhibition drastically decreases mortality of zebrafish infected with S. aureus. Furthermore, we show that NADK promotes S. aureus survival in infected macrophages by protecting bacteria from antimicrobial defense mechanisms. Proteome-wide data analysis revealed that production of major virulence-associated factors is sustained by NADK. We demonstrate that NADK is required for expression of the quorum-sensing response regulator AgrA, which controls critical S. aureus virulence determinants. These findings support a key role for NADK in bacteria survival within innate immune cells and the host during infection.

Extensively Drug-Resistant Klebsiella pneumoniae Counteracts Fitness and Virulence Costs That Accompanied Ceftazidime-Avibactam Resistance Acquisition

The ability of extensively drug-resistant (XDR) Klebsiella pneumoniae to rapidly acquire resistance to novel antibiotics is a global concern. Moreover, Klebsiella clonal lineages that successfully combine resistance and hypervirulence have increasingly occurred during the last years. However, the underlying mechanisms of counteracting fitness costs that accompany antibiotic resistance acquisition remain largely unexplored. Here, we investigated whether and how an XDR sequence type (ST)307 K. pneumoniae strain developed resistance against the novel drug combination ceftazidime-avibactam (CAZ-AVI) using experimental evolution. In addition, we performed in vitro and in vivo assays, molecular modeling, and bioinformatics to identify resistance-conferring processes and explore the resulting decrease in fitness and virulence. The subsequent amelioration of the initial costs was also addressed. We demonstrate that distinct mutations of the major nonselective porin OmpK36 caused CAZ-AVI resistance that persists even upon following a second experimental evolution without antibiotic selection pressure and that the Klebsiella strain compensates the resulting fitness and virulence costs. Furthermore, the genomic and transcriptomic analyses suggest the envelope stress response regulator rpoE and associated RpoE-regulated genes as drivers of this compensation. This study verifies the crucial role of OmpK36 in CAZ-AVI resistance and shows the rapid adaptation of a bacterial pathogen to compensate fitness- and virulence-associated resistance costs, which possibly contributes to the emergence of successful clonal lineages.

IMPORTANCE Extensively drug-resistant Klebsiella pneumoniae causing major outbreaks and severe infections has become a significant challenge for health care systems worldwide. Rapid resistance development against last-resort therapeutics like ceftazidime-avibactam is a significant driver for the accelerated emergence of such pathogens. Therefore, it is crucial to understand what exactly mediates rapid resistance acquisition and how bacterial pathogens counteract accompanying fitness and virulence costs. By combining bioinformatics with in vitro and in vivo phenotypic approaches, this study revealed the critical role of mutations in a particular porin channel in ceftazidime-avibactam resistance development and a major metabolic regulator for ameliorating fitness and virulence costs. These results highlight underlying mechanisms and contribute to the understanding of factors important for the emergence of successful bacterial pathogens.

Advances in Non-Coding RNA Sequencing

Non-coding RNAs (ncRNAs) comprise a set of abundant and functionally diverse RNA molecules. Since the discovery of the first ncRNA in the 1960s, ncRNAs have been shown to be involved in nearly all steps of the central dogma of molecular biology. In recent years, the pace of discovery of novel ncRNAs and their cellular roles has been greatly accelerated by high-throughput sequencing. Advances in sequencing technology, library preparation protocols as well as computational biology helped to greatly expand our knowledge of which ncRNAs exist throughout the kingdoms of life. Moreover, RNA sequencing revealed crucial roles of many ncRNAs in human health and disease. In this review, we discuss the most recent methodological advancements in the rapidly evolving field of high-throughput sequencing and how it has greatly expanded our understanding of ncRNA biology across a large number of different organisms.

Nucleoid-associated proteins shape chromatin structure and transcriptional regulation across the bacterial kingdom

Genome architecture has proven to be critical in determining gene regulation across almost all domains of life. While many of the key components and mechanisms of eukaryotic genome organization have been described, the interplay between bacterial DNA organization and gene regulation is only now being fully appreciated. An increasing pool of evidence has demonstrated that the bacterial chromosome can reasonably be thought of as chromatin, and that bacterial chromosomes contain transcriptionally silent and transcriptionally active regions analogous to heterochromatin and euchromatin, respectively. The roles played by histones in eukaryotic systems appear to be shared across a range of nucleoid-associated proteins (NAPs) in bacteria, which function to compact, structure, and regulate large portions of bacterial chromosomes. The broad range of extant NAPs, and the extent to which they differ from species to species, has raised additional challenges in identifying and characterizing their roles in all but a handful of model bacteria. Here we review the regulatory roles played by NAPs in several well-studied bacteria and use the resulting state of knowledge to provide a working definition for NAPs, based on their function, binding pattern, and expression levels. We present a screening procedure which can be applied to any species for which transcriptomic data are available. Finally, we note that NAPs tend to play two major regulatory roles - xenogeneic silencers and developmental regulators - and that many unrecognized potential NAPs exist in each bacterial species examined.