Search

Abstract Background High-throughput sequencing of cDNA libraries (RNA-Seq) has proven to be a highly effective approach for studying bacterial transcriptomes. A central challenge in designing RNA-Seq-based experiments is estimating a priori the number of reads per sample needed to detect and quantify thousands of individual transcripts with a large dynamic range of abundance. Results We have conducted a systematic examination of how changes in the number of RNA-Seq reads per sample influences both profiling of a single bacterial transcriptome and the comparison of gene expression among samples. Our findings suggest that the number of reads typically produced in a single lane of the Illumina HiSeq sequencer far exceeds the number needed to saturate the annotated transcriptomes of diverse bacteria growing in monoculture. Moreover, as sequencing depth increases, so too does the detection of cDNAs that likely correspond to spurious transcripts or genomic DNA contamination. Finally, even when dozens of barcoded individual cDNA libraries are sequenced in a single lane, the vast majority of transcripts in each sample can be detected and numerous genes differentially expressed between samples can be identified. Conclusions Our analysis provides a guide for the many researchers seeking to determine the appropriate sequencing depth for RNA-Seq-based studies of diverse bacterial species.

Abstract Background Regulation of bacterial gene expression by small RNAs (sRNAs) have proved to be important for many biological processes. Francisella tularensis is a highly pathogenic Gram-negative bacterium that causes the disease tularaemia in humans and animals. Relatively little is known about the regulatory networks existing in this organism that allows it to survive in a wide array of environments and no sRNA regulators have been identified so far. Results We have used a combination of experimental assays and in silico prediction to identify sRNAs in F. tularensis strain LVS. Using a cDNA cloning and sequencing approach we have shown that F. tularensis expresses homologues of several sRNAs that are well-conserved among diverse bacteria. We have also discovered two abundant putative sRNAs that share no sequence similarity or conserved genomic context with any previously annotated regulatory transcripts. Deletion of either of these two loci led to significant changes in the expression of several mRNAs that likely include the cognate target(s) of these sRNAs. Deletion of these sRNAs did not, however, significantly alter F. tularensis growth under various stress conditions in vitro, its replication in murine cells, or its ability to induce disease in a mouse model of F. tularensis infection. We also conducted a genome-wide in silico search for intergenic loci that suggests F. tularensis encodes several other sRNAs in addition to the sRNAs found in our experimental screen. Conclusion Our findings suggest that F. tularensis encodes a significant number of non-coding regulatory RNAs, including members of well conserved families of structural and housekeeping RNAs and other poorly conserved transcripts that may have evolved more recently to help F. tularensis deal with the unique and diverse set of environments with which it must contend.

Abstract Background Small non-coding RNAs (sRNAs) have emerged as ubiquitous regulatory elements in bacteria and other life domains. However, few sRNAs have been identified outside several well-studied species of gamma-proteobacteria and thus relatively little is known about the role of RNA-mediated regulation in most other bacterial genera. Here we have conducted a computational prediction of putative sRNA genes in intergenic regions (IgRs) of the symbiotic α-proteobacterium S. meliloti 1021 and experimentally confirmed the expression of dozens of these candidate loci in the closely related strain S. meliloti 2011. Results Our first sRNA candidate compilation was based mainly on the output of the sRNAPredictHT algorithm. A thorough manual sequence analysis of the curated list rendered an initial set of 18 IgRs of interest, from which 14 candidates were detected in strain 2011 by Northern blot and/or microarray analysis. Interestingly, the intracellular transcript levels varied in response to various stress conditions. We developed an alternative computational method to more sensitively predict sRNA-encoding genes and score these predicted genes based on several features to allow identification of the strongest candidates. With this novel strategy, we predicted 60 chromosomal independent transcriptional units that, according to our annotation, represent strong candidates for sRNA-encoding genes, including most of the sRNAs experimentally verified in this work and in two other contemporary studies. Additionally, we predicted numerous candidate sRNA genes encoded in megaplasmids pSymA and pSymB. A significant proportion of the chromosomal- and megaplasmid-borne putative sRNA genes were validated by microarray analysis in strain 2011. Conclusion Our data extend the number of experimentally detected S. meliloti sRNAs and significantly expand the list of putative sRNA-encoding IgRs in this and closely related α-proteobacteria. In addition, we have developed a computational method that proved useful to predict sRNA-encoding genes in S. meliloti. We anticipate that this predictive approach can be flexibly implemented in many other bacterial species.

The rising prevalence of antibiotic resistance threatens human health. While more sophisticated strategies for antibiotic discovery are being developed, target elucidation of new chemical entities remains challenging. In the post-genomic era, expression profiling can play an important role in mechanism-of-action (MOA) prediction by reporting on the cellular response to perturbation. However, the broad application of transcriptomics has yet to fulfill its promise of transforming target elucidation due to challenges in identifying the most relevant, direct responses to target inhibition. We developed an unbiased strategy for MOA prediction, called Perturbation-Specific Transcriptional Mapping (PerSpecTM), in which large-throughput expression profiling of wildtype or hypomorphic mutants, depleted for essential targets, enables a computational strategy to address this challenge. We applied PerSpecTM to perform reference-based MOA prediction based on the principle that similar perturbations, whether chemical or genetic, will elicit similar transcriptional responses. Using this approach, we elucidated the MOAs of three new molecules with activity against Pseudomonas aeruginosa by comparing their expression profiles to those of a reference set of antimicrobial compounds with known MOAs. We also show that transcriptional responses to small molecule inhibition resemble those resulting from genetic depletion of essential targets by CRISPRi by PerSpecTM, demonstrating proof-of-concept that correlations between expression profiles of small molecule and genetic perturbations can facilitate MOA prediction when no chemical entities exist to serve as a reference. Empowered by PerSpecTM, this work lays the foundation for an unbiased, readily scalable, systematic reference-based strategy for MOA elucidation that could transform antibiotic discovery efforts.

Several bacterial pathogens, including Salmonella enterica, can cause persistent infections in humans by mechanisms that are poorly understood. By comparing genomes of isolates longitudinally collected from 256 prolonged salmonellosis patients, we identified repeated mutations in global regulators, including the barA/sirA two-component regulatory system, across multiple patients and Salmonella serovars. Comparative RNA-seq analysis revealed that distinct mutations in barA/sirA led to diminished expression of Salmonella pathogenicity islands 1 and 4 genes, which are required for Salmonella invasion and enteritis. Moreover, barA/sirA mutants were attenuated in an acute salmonellosis mouse model and induced weaker transcription of host immune responses. In contrast, in a persistent infection mouse model, these mutants exhibited long-term colonization and prolonged shedding. Taken together, these findings suggest that selection of mutations in global virulence regulators facilitates persistent Salmonella infection in humans, by attenuating Salmonella virulence and inducing a weaker host inflammatory response., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Lichtheimia ornata is an emerging opportunistic Mucorales pathogen that is associated with fatal infections in immunocompromised individuals. While these environmentally acquired infections have rarely been reported to date, cases were noted in a recent analysis of coronavirus disease 2019 (COVID-19)-associated mucormycosis in India. Here, we report the annotated genome sequence of the environmental isolate CBS 291.66., Competing Interests: The authors declare no conflict of interest.

Previous urinary tract infections (UTIs) can predispose one to future infections; however, the underlying mechanisms affecting recurrence are poorly understood. We previously found that UTIs in mice cause differential bladder epithelial (urothelial) remodelling, depending on disease outcome, that impacts susceptibility to recurrent UTI. Here we compared urothelial stem cell (USC) lines isolated from mice with a history of either resolved or chronic uropathogenic Escherichia coli (UPEC) infection, elucidating evidence of molecular imprinting that involved epigenetic changes, including differences in chromatin accessibility, DNA methylation and histone modification. Epigenetic marks in USCs from chronically infected mice enhanced caspase-1-mediated cell death upon UPEC infection, promoting bacterial clearance. Increased Ptgs2os2 expression also occurred, potentially contributing to sustained cyclooxygenase-2 expression, bladder inflammation and mucosal wounding-responses associated with severe recurrent cystitis. Thus, UPEC infection acts as an epi-mutagen reprogramming the urothelial epigenome, leading to urothelial-intrinsic remodelling and training of the innate response to subsequent infection., (© 2023. The Author(s).)

The ability to measure neutralizing antibodies on large scale can be important for understanding features of the natural history and epidemiology of infection, as well as an aid in determining the efficacy of interventions, particularly in outbreaks such as the current severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic. Because of the assay's rapid scalability and high efficiency, serology measurements that quantify the presence rather than function of serum antibodies often serve as proxies of immune protection. Here, we report the development of a high-throughput, automated fluorescence-based neutralization assay using SARS-CoV-2 virus to quantify neutralizing antibody activity in patient specimens. We performed large-scale testing of over 19,000 COVID-19 convalescent plasma (CCP) samples from patients who had been infected with SARS-CoV-2 between March and August 2020 across the United States. The neutralization capacity of the samples was moderately correlated with serological measurements of anti-receptor-binding domain (RBD) IgG levels. The neutralizing antibody levels within these convalescent-phase serum samples were highly variable against the original USA-WA1/2020 strain with almost 10% of individuals who had had PCR-confirmed SARS-CoV-2 infection having no detectable antibodies either by serology or neutralization, and ~1/3 having no or low neutralizing activity. Discordance between neutralization and serology measurements was mainly due to the presence of non-IgG RBD isotypes. Meanwhile, natural infection with the earliest SARS-CoV-2 strain USA-WA1/2020 resulted in weaker neutralization of subsequent B.1.1.7 (alpha) and the B.1.351 (beta) variants, with 88% of samples having no activity against the BA.1 (omicron) variant. IMPORTANCE The ability to directly measure neutralizing antibodies on live SARS-CoV-2 virus in individuals can play an important role in understanding the efficacy of therapeutic interventions or vaccines. In contrast to functional neutralization assays, serological assays only quantify the presence of antibodies as a proxy of immune protection. Here, we have developed a high-throughput, automated neutralization assay for SARS-CoV-2 and measured the neutralizing activity of ~19,000 COVID-19 convalescent plasma (CCP) samples collected across the United States between March and August of 2020. These data were used to support the FDA's interpretation of CCP efficacy in patients with SARS-CoV-2 infection and their issuance of emergency use authorization of CCP in 2020.

Background: Unbiased assessment of the risks associated with acquisition of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is critical to informing mitigation efforts during pandemics. The objective of our study was to understand the risk factors for acquiring coronavirus disease 2019 (COVID-19) in a large prospective cohort of adult residents in a large US metropolitan area., Methods: We designed a fully remote longitudinal cohort study involving monthly at-home SARS-CoV-2 polymerase chain reaction (PCR) and serology self-testing and monthly surveys., Results: Between October 2020 and January 2021, we enrolled 10 289 adults reflective of the Boston metropolitan area census data. At study entry, 567 (5.5%) participants had evidence of current or prior SARS-CoV-2 infection. This increased to 13.4% by June 15, 2021. Compared with Whites, Black non-Hispanic participants had a 2.2-fold greater risk of acquiring COVID-19 (hazard ratio [HR], 2.19; 95% CI, 1.91-2.50; P < .001), and Hispanics had a 1.5-fold greater risk (HR, 1.52; 95% CI, 1.32-1.71; P < .016). Individuals aged 18-29, those who worked outside the home, and those living with other adults and children were at an increased risk. Individuals in the second and third lowest disadvantaged neighborhood communities were associated with an increased risk of acquiring COVID-19. Individuals with medical risk factors for severe disease were at a decreased risk of SARS-CoV-2 acquisition., Conclusions: These results demonstrate that race/ethnicity and socioeconomic status are the biggest determinants of acquisition of infection. This disparity is significantly underestimated if based on PCR data alone, as noted by the discrepancy in serology vs PCR detection for non-White participants, and points to persistent disparity in access to testing. Medical conditions and advanced age, which increase the risk for severity of SARS-CoV-2 disease, were associated with a lower risk of COVID-19 acquisition, suggesting the importance of behavior modifications. These findings highlight the need for mitigation programs that overcome challenges of structural racism in current and future pandemics., Competing Interests: Potential conflicts of interest. D.H. is a founder, consultant to, equity holder in, and inventor of technology licensed to Sherlock Biosciences and serves on the Scientific Advisory Board for Proof Diagnostics; both are infectious disease diagnostic companies, but neither company’s technologies was used in this work. All other authors report no potential conflicts. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed., (© The Author(s) 2022. Published by Oxford University Press on behalf of Infectious Diseases Society of America.)

Rapid and accurate diagnosis of infections is fundamental to individual patient care and public health management. Nucleic acid detection methods are critical to this effort, but are limited either in the breadth of pathogens targeted or by the expertise and infrastructure required. We present here a high-throughput system that enables rapid identification of bacterial pathogens, bCARMEN, which utilizes: (1) modular CRISPR-Cas13-based nucleic acid detection with enhanced sensitivity and specificity; and (2) a droplet microfluidic system that enables thousands of simultaneous, spatially multiplexed detection reactions at nanoliter volumes; and (3) a novel preamplification strategy that further enhances sensitivity and specificity. We demonstrate bCARMEN is capable of detecting and discriminating 52 clinically relevant bacterial species and several key antibiotic resistance genes. We further develop a simple proof of principle workflow using stabilized reagents and cell phone camera optical readout, opening up the possibility of a rapid point-of-care multiplexed bacterial pathogen identification and antibiotic susceptibility testing., (© The Author(s) 2022. Published by Oxford University Press on behalf of the National Academy of Sciences.)

Interspecies interactions are known to activate specialized metabolism in diverse actinomycetes. However, how interspecies cues are sensed and ultimately lead to induction of specialized metabolite biosynthetic gene clusters remains largely unexplored. Using transcriptome sequencing (RNA-seq), we analyzed genes that were transcriptionally induced in the model actinomycete Streptomyces coelicolor during interactions with four different actinomycetes, including genes that encode unusual regulatory systems known as conservons. Deletions in one such system, encoded by the cvn8 genes, led to altered patterns of pigmented antibiotic production by S. coelicolor during interactions. Further transcriptomic analysis of mutants lacking each of the five genes in the cvn8 locus demonstrated that this system is a global regulator of at least four different specialized metabolite biosynthetic pathways. How conservon systems work at the mechanistic level to regulate gene expression is not well understood, although it has been hypothesized that they may function in a way similar to eukaryotic G-protein-coupled receptors. The data presented here indicate that the gene products of the cvnA8 and cvnF8 (SCO6939) genes likely function together in one part of the Cvn8 signaling cascade, while the cvnC8 and cvnD8 gene products likely function together in another part. Importantly, because cvnD8 likely encodes a Ras-like GTPase, these results connect G-protein-mediated signaling to gene regulation in a bacterium. Additionally, deletion of any of the cvn8 genes led to abnormally high expression of an adjacent cryptic lanthipeptide biosynthetic gene cluster, indicating that conservon systems may be fruitful targets for manipulation to activate silent specialized metabolite biosynthetic pathways. IMPORTANCE Interactions between different species of actinomycete bacteria often trigger one of the strains to produce specialized metabolites, such as antibiotics. However, how this induction occurs at the genetic level is poorly understood. Using transcriptomic methods, we show that an unusual regulatory system, known as a conservon system, is responsible for regulating expression of multiple specialized metabolite biosynthetic gene clusters in the organism Streptomyces coelicolor during interactions. Conservon systems are unusual because they appear to employ small GTPases as an important component of their signaling cascades. Small GTPases are common in eukaryotic signaling pathways, but the results presented here are notable since they implicate a system that includes a small GTPase in global gene regulation in a bacterium. Mutants lacking this conservon system also showed abnormally high expression of a gene cluster involved in making an unknown specialized metabolite, suggesting that conservon mutants might be useful for driving natural product discovery.

The omnigenic/polygenic theory, which states that complex traits are not shaped by single/few genes, but by situation-specific large networks, offers an explanation for a major enigma in microbiology: deletion of specific small RNAs (sRNAs) playing key roles in various aspects of bacterial physiology, including virulence and antibiotic resistance, results in surprisingly subtle phenotypes. A recent study uncovered polar accumulation of most sRNAs upon osmotic stress, the majority not known to be involved in the applied stress. Here we show that cells deleted for a handful of pole-enriched sRNAs exhibit fitness defect in several stress conditions, as opposed to single, double, or triple sRNA-knockouts, implying that regulation by sRNA relies on sets of genes. Moreover, analysis of RNA-seq data of Escherichia coli and Salmonella typhimurium exposed to antibiotics and/or infection-relevant conditions reveals the involvement of multiple sRNAs in all cases, in line with the existence of a polygenic plan for sRNA-mediated regulation., Competing Interests: The authors declare no competing interests., (© 2021 The Authors.)

Background: Antibody response duration following severe acute respiratory syndrome coronavirus 2 infection tends to be variable and depends on severity of disease and method of detection., Study Design and Methods: COVID-19 convalescent plasma from 18 donors was collected longitudinally for a maximum of 63-129 days following resolution of symptoms. All the samples were initially screened by the Ortho total Ig test to confirm positivity and subsequently tested with seven additional direct sandwich or indirect binding assays (Ortho, Roche, Abbott, Broad Institute) directed against a variety of antigen targets (S1, receptor binding domain, and nucleocapsid [NC]), along with two neutralization assays (Broad Institute live virus PRNT and Vitalant Research Institute [VRI] Pseudovirus reporter viral particle neutralization [RVPN])., Results: The direct detection assays (Ortho total Ig total and Roche total Ig) showed increasing levels of antibodies over the time period, in contrast to the indirect IgG assays that showed a decline. Neutralization assays also demonstrated declining responses; the VRI RVPN pseudovirus had a greater rate of decline than the Broad PRNT live virus assay., Discussion: These data show that in addition to variable individual responses and associations with disease severity, the detection assay chosen contributes to the heterogeneous results in antibody stability over time. Depending on the scope of the research, one assay may be preferable over another. For serosurveillance studies, direct, double Ag-sandwich assays appear to be the best choice due to their stability; in particular, algorithms that include both S1- and NC-based assays can help reduce the rate of false-positivity and discriminate between natural infection and vaccine-derived seroreactivity., (© 2021 AABB.)

Advances in transcriptome sequencing allow for simultaneous interrogation of differentially expressed genes from multiple species originating from a single RNA sample, termed dual or multi-species transcriptomics. Compared to single-species differential expression analysis, the design of multi-species differential expression experiments must account for the relative abundances of each organism of interest within the sample, often requiring enrichment methods and yielding differences in total read counts across samples. The analysis of multi-species transcriptomics datasets requires modifications to the alignment, quantification, and downstream analysis steps compared to the single-species analysis pipelines. We describe best practices for multi-species transcriptomics and differential gene expression.

In this era of rising antibiotic resistance, in contrast to our increasing understanding of mechanisms that cause resistance, our understanding of mechanisms that influence the propensity to evolve resistance remains limited. Here, we identified genetic factors that facilitate the evolution of resistance to carbapenems, the antibiotic of 'last resort', in Klebsiella pneumoniae , the major carbapenem-resistant species. In clinical isolates, we found that high-level transposon insertional mutagenesis plays an important role in contributing to high-level resistance frequencies in several major and emerging carbapenem-resistant lineages. A broader spectrum of resistance-conferring mutations for select carbapenems such as ertapenem also enables higher resistance frequencies and, importantly, creates stepping-stones to achieve high-level resistance to all carbapenems. These mutational mechanisms can contribute to the evolution of resistance, in conjunction with the loss of systems that restrict horizontal resistance gene uptake, such as the CRISPR-Cas system. Given the need for greater antibiotic stewardship, these findings argue that in addition to considering the current efficacy of an antibiotic for a clinical isolate in antibiotic selection, considerations of future efficacy are also important. The genetic background of a clinical isolate and the exact antibiotic identity can and should also be considered as they are determinants of a strain's propensity to become resistant. Together, these findings thus provide a molecular framework for understanding acquisition of carbapenem resistance in K. pneumoniae with important implications for diagnosing and treating this important class of pathogens., Competing Interests: PM, LH, AP, HL, CE, AM, RB, AE, JL, DH No competing interests declared, (© 2021, Ma et al.)

The genomes of gut Bacteroidales contain numerous invertible regions, many of which contain promoters that dictate phase-variable synthesis of surface molecules such as polysaccharides, fimbriae, and outer surface proteins. Here, we characterize a different type of phase-variable system of Bacteroides fragilis, a Type I restriction modification system (R-M). We show that reversible DNA inversions within this R-M locus leads to the generation of eight specificity proteins with distinct recognition sites. In vitro grown bacteria have a different proportion of specificity gene combinations at the expression locus than bacteria isolated from the mammalian gut. By creating mutants, each able to produce only one specificity protein from this region, we identified the R-M recognition sites of four of these S-proteins using SMRT sequencing. Transcriptome analysis revealed that the locked specificity mutants, whether grown in vitro or isolated from the mammalian gut, have distinct transcriptional profiles, likely creating different phenotypes, one of which was confirmed. Genomic analyses of diverse strains of Bacteroidetes from both host-associated and environmental sources reveal the ubiquity of phase-variable R-M systems in this phylum., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Dual transcriptional profiling of host and bacteria during infection is challenging due to the low abundance of bacterial mRNA. We report Pathogen Hybrid Capture (PatH-Cap), a method to enrich for bacterial mRNA and deplete bacterial rRNA simultaneously from dual RNA-seq libraries using transcriptome-specific probes. By addressing both the differential RNA content of the host relative to the infecting bacterium and the overwhelming abundance of uninformative structural RNAs (rRNA, tRNA) of both species in a single step, this approach enables analysis of very low-input RNA samples. By sequencing libraries before (pre-PatH-Cap) and after (post-PatH-Cap) enrichment, we achieve dual transcriptional profiling of host and bacteria, respectively, from the same sample. Importantly, enrichment preserves relative transcript abundance and increases the number of unique bacterial transcripts per gene in post-PatH-Cap libraries compared to pre-PatH-Cap libraries at the same sequencing depth, thereby decreasing the sequencing depth required to fully capture the transcriptional profile of the infecting bacteria. We demonstrate that PatH-Cap enables the study of low-input samples including single eukaryotic cells infected by 1-3 Pseudomonas aeruginosa bacteria and paired host-pathogen temporal gene expression analysis of Mycobacterium tuberculosis infecting macrophages. PatH-Cap can be applied to the study of a range of pathogens and microbial species, and more generally, to lowly-abundant species in mixed populations.

Multidrug resistant organisms are a serious threat to human health 1,2 . Fast, accurate antibiotic susceptibility testing (AST) is a critical need in addressing escalating antibiotic resistance, since delays in identifying multidrug resistant organisms increase mortality 3,4 and use of broad-spectrum antibiotics, further selecting for resistant organisms. Yet current growth-based AST assays, such as broth microdilution 5 , require several days before informing key clinical decisions. Rapid AST would transform the care of patients with infection while ensuring that our antibiotic arsenal is deployed as efficiently as possible. Growth-based assays are fundamentally constrained in speed by doubling time of the pathogen, and genotypic assays are limited by the ever-growing diversity and complexity of bacterial antibiotic resistance mechanisms. Here we describe a rapid assay for combined genotypic and phenotypic AST through RNA detection, GoPhAST-R, that classifies strains with 94-99% accuracy by coupling machine learning analysis of early antibiotic-induced transcriptional changes with simultaneous detection of key genetic resistance determinants to increase accuracy of resistance detection, facilitate molecular epidemiology and enable early detection of emerging resistance mechanisms. This two-pronged approach provides phenotypic AST 24-36 h faster than standard workflows, with <4 h assay time on a pilot instrument for hybridization-based multiplexed RNA detection implemented directly from positive blood cultures.

RNA localization in eukaryotes is a mechanism to regulate transcripts fate. Conversely, bacterial transcripts were not assumed to be specifically localized. We previously demonstrated that E. coli mRNAs may localize to where their products localize in a translation-independent manner, thus challenging the transcription-translation coupling extent. However, the scope of RNA localization in bacteria remained unknown. Here, we report the distribution of the E. coli transcriptome between the membrane, cytoplasm, and poles by combining cell fractionation with deep-sequencing (Rloc-seq). Our results reveal asymmetric RNA distribution on a transcriptome-wide scale, significantly correlating with proteome localization and prevalence of translation-independent RNA localization. The poles are enriched with stress-related mRNAs and small RNAs, the latter becoming further enriched upon stress in an Hfq-dependent manner. Genome organization may play a role in localizing membrane protein-encoding transcripts. Our results show an unexpected level of intricacy in bacterial transcriptome organization and highlight the poles as hubs for regulation., (Copyright © 2019 Elsevier Inc. All rights reserved.)

A mucosal infectious disease episode can render the host either more or less susceptible to recurrent infection, but the specific mechanisms that tip the balance remain unclear. We investigated this question in a mouse model of recurrent urinary tract infection and found that a prior bladder infection resulted in an earlier onset of tumor necrosis factor-alpha (TNFɑ)-mediated bladder inflammation upon subsequent bacterial challenge, relative to age-matched naive mice. However, the duration of TNFɑ signaling activation differed according to whether the first infection was chronic (Sensitized) or self-limiting (Resolved). TNFɑ depletion studies revealed that transient early-phase TNFɑ signaling in Resolved mice promoted clearance of bladder-colonizing bacteria via rapid recruitment of neutrophils and subsequent exfoliation of infected bladder cells. In contrast, sustained TNFɑ signaling in Sensitized mice prolonged damaging inflammation, worsening infection. This work reveals how TNFɑ signaling dynamics can be rewired by a prior infection to shape diverse susceptibilities to future mucosal infections., Competing Interests: LY, VO, JL, DD, NB, MC, MC, ER, SH, TH No competing interests declared, (© 2019, Yu et al.)

Mycobacterium tuberculosis ( Mtb ) killed more people in 2017 than any other single infectious agent. This dangerous pathogen is able to withstand stresses imposed by the immune system and tolerate exposure to antibiotics, resulting in persistent infection. The global tuberculosis (TB) epidemic has been exacerbated by the emergence of mutant strains of Mtb that are resistant to frontline antibiotics. Thus, both phenotypic drug tolerance and genetic drug resistance are major obstacles to successful TB therapy. Using a chemical approach to identify compounds that block stress and drug tolerance, as opposed to traditional screens for compounds that kill Mtb , we identified a small molecule, C10, that blocks tolerance to oxidative stress, acid stress, and the frontline antibiotic isoniazid (INH). In addition, we found that C10 prevents the selection for INH-resistant mutants and restores INH sensitivity in otherwise INH-resistant Mtb strains harboring mutations in the katG gene, which encodes the enzyme that converts the prodrug INH to its active form. Through mechanistic studies, we discovered that C10 inhibits Mtb respiration, revealing a link between respiration homeostasis and INH sensitivity. Therefore, by using C10 to dissect Mtb persistence, we discovered that INH resistance is not absolute and can be reversed., Competing Interests: Conflict of interest statement: C.L.S., S.J.H., and F.A. have ownership interests in Quretech Bio AB, which licenses C10.

Rapid bacterial identification remains a critical challenge in infectious disease diagnostics. We developed a novel molecular approach to detect and identify a wide diversity of bacterial pathogens in a single, simple assay, exploiting the conservation, abundance, and rich phylogenetic content of ribosomal RNA in a rapid fluorescent hybridization assay that requires no amplification or enzymology. Of 117 isolates from 64 species across 4 phyla, this assay identified bacteria with >89% accuracy at the species level and 100% accuracy at the family level, enabling all critical clinical distinctions. In pilot studies on primary clinical specimens, including sputum, blood cultures, and pus, bacteria from 5 different phyla were identified.

Toxin synthesis and endospore formation are two of the most critical factors that determine the outcome of infection by Clostridioides difficile. The two major toxins, TcdA and TcdB, are the principal factors causing damage to the host. Spores are the infectious form of C. difficile, permit survival of the bacterium during antibiotic treatment and are the predominant cell form that leads to recurrent infection. Toxin production and sporulation have their own specific mechanisms of regulation, but they share negative regulation by the global regulatory protein CodY. Determining the extent of such regulation and its detailed mechanism is important for understanding the linkage between two apparently independent biological phenomena and raises the possibility of creating new ways of limiting infection. The work described here shows that a codY null mutant of a hypervirulent (ribotype 027) strain is even more virulent than its parent in a mouse model of infection and that the mutant expresses most sporulation genes prematurely during exponential growth phase. Moreover, examining the expression patterns of mutants producing CodY proteins with different levels of residual activity revealed that expression of the toxin genes is dependent on total CodY inactivation, whereas most sporulation genes are turned on when CodY activity is only partially diminished. These results suggest that, in wild-type cells undergoing nutrient limitation, sporulation genes can be turned on before the toxin genes., Competing Interests: The authors have declared that no competing interests exist.

One key to the success of Mycobacterium tuberculosis as a pathogen is its ability to reside in the hostile environment of the human macrophage. Bacteria adapt to stress through a variety of mechanisms, including the use of small regulatory RNAs (sRNAs), which posttranscriptionally regulate bacterial gene expression. However, very little is currently known about mycobacterial sRNA-mediated riboregulation. To date, mycobacterial sRNA discovery has been performed primarily in log-phase growth, and no direct interaction between any mycobacterial sRNA and its targets has been validated. Here, we performed large-scale sRNA discovery and expression profiling in M. tuberculosis during exposure to five pathogenically relevant stresses. From these data, we identified a subset of sRNAs that are highly induced in multiple stress conditions. We focused on one of these sRNAs, ncRv11846, here renamed mycobacterial regulatory sRNA in iron (MrsI). We characterized the regulon of MrsI and showed in mycobacteria that it regulates one of its targets, bfrA , through a direct binding interaction. MrsI mediates an iron-sparing response that is required for optimal survival of M. tuberculosis under iron-limiting conditions. However, MrsI is induced by multiple host-like stressors, which appear to trigger MrsI as part of an anticipatory response to impending iron deprivation in the macrophage environment., Competing Interests: The authors declare no conflict of interest., (Copyright © 2018 the Author(s). Published by PNAS.)

The epidemic pathogen Vibrio cholerae senses and responds to different external stresses it encounters in the aquatic environment and in the human host. One stress that V. cholerae encounters in the host is exposure to antimicrobial peptides on mucosal surfaces. We used massively parallel cDNA sequencing (RNA-Seq) to quantitatively identify the transcriptome of V. cholerae grown in the presence and absence of sub-lethal concentrations of the antimicrobial peptide polymyxin B. We evaluated the transcriptome of both wild type V. cholerae and a mutant carrying a deletion of vc1639, a putative sensor kinase of an uncharacterized two-component system, under these conditions. In addition to many previously uncharacterized pathways responding with elevated transcript levels to polymyxin B exposure, we confirmed the predicted elevated transcript levels of a previously described LPS modification system in response to polymyxin B exposure. Additionally, we identified the V. cholerae homologue of visP (ygiW) as a regulatory target of VC1639. VisP is a conserved periplasmic protein implicated in lipid A modification in Salmonellae. This study provides the first systematic analysis of the transcriptional response of Vibrio cholerae to polymyxin B, raising important questions for further study regarding mechanisms used by V. cholerae to sense and respond to envelope stress.

Growing evidence suggests that microbes can influence the efficacy of cancer therapies. By studying colon cancer models, we found that bacteria can metabolize the chemotherapeutic drug gemcitabine (2',2'-difluorodeoxycytidine) into its inactive form, 2',2'-difluorodeoxyuridine. Metabolism was dependent on the expression of a long isoform of the bacterial enzyme cytidine deaminase (CDD L ), seen primarily in Gammaproteobacteria. In a colon cancer mouse model, gemcitabine resistance was induced by intratumor Gammaproteobacteria, dependent on bacterial CDD L expression, and abrogated by cotreatment with the antibiotic ciprofloxacin. Gemcitabine is commonly used to treat pancreatic ductal adenocarcinoma (PDAC), and we hypothesized that intratumor bacteria might contribute to drug resistance of these tumors. Consistent with this possibility, we found that of the 113 human PDACs that were tested, 86 (76%) were positive for bacteria, mainly Gammaproteobacteria., (Copyright © 2017 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

This study was conducted to investigate the role of iron deprivation in the persistence of Mycobacterium tuberculosis We present evidence of iron restriction in human necrotic granulomas and demonstrate that under iron starvation M. tuberculosis persists, refractive to antibiotics and capable of restarting replication when iron is made available. Transcriptomics and metabolomic analyses indicated that the persistence of M. tuberculosis under iron starvation is dependent on strict control of endogenous Fe utilization and is associated with upregulation of pathogenicity and intrinsic antibiotic resistance determinants. M. tuberculosis mutants compromised in their ability to survive Fe starvation were identified. The findings of this study advance the understanding of the physiological settings that may underpin the chronicity of human tuberculosis (TB) and are relevant to the design of effective antitubercular therapies. IMPORTANCE One-third of the world population may harbor persistent M. tuberculosis , causing an asymptomatic infection that is refractory to treatment and can reactivate to become potentially lethal tuberculosis disease. However, little is known about the factors that trigger and maintain M. tuberculosis persistence in infected individuals. Iron is an essential nutrient for M. tuberculosis growth. In this study, we show, first, that in human granulomas the immune defense creates microenvironments in which M. tuberculosis likely experiences drastic Fe deprivation and, second, that Fe-starved M. tuberculosis is capable of long-term persistence without growth. Together, these observations suggest that Fe deprivation in the lung might trigger a state of persistence in M. tuberculosis and promote chronic TB. We also identified vulnerabilities of iron-restricted persistent M. tuberculosis , which can be exploited for the design of new antitubercular therapies., (Copyright © 2017 Kurthkoti et al.)

Rapid, inexpensive, and sensitive nucleic acid detection may aid point-of-care pathogen detection, genotyping, and disease monitoring. The RNA-guided, RNA-targeting clustered regularly interspaced short palindromic repeats (CRISPR) effector Cas13a (previously known as C2c2) exhibits a "collateral effect" of promiscuous ribonuclease activity upon target recognition. We combine the collateral effect of Cas13a with isothermal amplification to establish a CRISPR-based diagnostic (CRISPR-Dx), providing rapid DNA or RNA detection with attomolar sensitivity and single-base mismatch specificity. We use this Cas13a-based molecular detection platform, termed Specific High-Sensitivity Enzymatic Reporter UnLOCKing (SHERLOCK), to detect specific strains of Zika and Dengue virus, distinguish pathogenic bacteria, genotype human DNA, and identify mutations in cell-free tumor DNA. Furthermore, SHERLOCK reaction reagents can be lyophilized for cold-chain independence and long-term storage and be readily reconstituted on paper for field applications., (Copyright © 2017, American Association for the Advancement of Science.)

Rifampicin, which inhibits bacterial RNA polymerase, provides one of the most effective treatments for tuberculosis. Inhibition of the transcription termination factor Rho is used to treat some bacterial infections, but its importance varies across bacteria. Here we show that Rho of Mycobacterium tuberculosis functions to both define the 3' ends of mRNAs and silence substantial fragments of the genome. Brief inactivation of Rho affects over 500 transcripts enriched for genes of foreign DNA elements and bacterial virulence factors. Prolonged inactivation of Rho causes extensive pervasive transcription, a genome-wide increase in antisense transcripts, and a rapid loss of viability of replicating and non-replicating M. tuberculosis in vitro and during acute and chronic infection in mice. Collectively, these data suggest that inhibition of Rho may provide an alternative strategy to treat tuberculosis with an efficacy similar to inhibition of RNA polymerase.

Urinary tract infections (UTIs) are caused by uropathogenic Escherichia coli (UPEC) strains. In contrast to many enteric E. coli pathogroups, no genetic signature has been identified for UPEC strains. We conducted a high-resolution comparative genomic study using E. coli isolates collected from the urine of women suffering from frequent recurrent UTIs. These isolates were genetically diverse and varied in their urovirulence, that is, their ability to infect the bladder in a mouse model of cystitis. We found no set of genes, including previously defined putative urovirulence factors (PUFs), that were predictive of urovirulence. In addition, in some patients, the E. coli strain causing a recurrent UTI had fewer PUFs than the supplanted strain. In competitive experimental infections in mice, the supplanting strain was more efficient at colonizing the mouse bladder than the supplanted strain. Despite the lack of a clear genomic signature for urovirulence, comparative transcriptomic and phenotypic analyses revealed that the expression of key conserved functions during culture, such as motility and metabolism, could be used to predict subsequent colonization of the mouse bladder. Together, our findings suggest that UTI risk and outcome may be determined by complex interactions between host susceptibility and the urovirulence potential of diverse bacterial strains., (Copyright © 2017, American Association for the Advancement of Science.)

Antibiotic resistance arising via chromosomal mutations is typically specific to a particular antibiotic or class of antibiotics. We have identified mutations in genes encoding ribosomal components in Mycobacterium smegmatis that confer resistance to several structurally and mechanistically unrelated classes of antibiotics and enhance survival following heat shock and membrane stress. These mutations affect ribosome assembly and cause large-scale transcriptomic and proteomic changes, including the downregulation of the catalase KatG, an activating enzyme required for isoniazid sensitivity, and upregulation of WhiB7, a transcription factor involved in innate antibiotic resistance. Importantly, while these ribosomal mutations have a fitness cost in antibiotic-free medium, in a multidrug environment they promote the evolution of high-level, target-based resistance. Further, suppressor mutations can then be easily acquired to restore wild-type growth. Thus, ribosomal mutations can serve as stepping-stones in an evolutionary path leading to the emergence of high-level, multidrug resistance.

Recurrent bacterial infections are a significant burden worldwide, and prior history of infection is often a significant risk factor for developing new infections. For urinary tract infection (UTI), a history of two or more episodes is an independent risk factor for acute infection. However, mechanistic knowledge of UTI pathogenesis has come almost exclusively from studies in naive mice. Here we show that, in mice, an initial Escherichia coli UTI, whether chronic or self-limiting, leaves a long-lasting molecular imprint on the bladder tissue that alters the pathophysiology of subsequent infections, affecting host susceptibility and disease outcome. In bladders of previously infected versus non-infected, antibiotic-treated mice, we found (1) an altered transcriptome and defects in cell maturation, (2) a remodelled epithelium that confers resistance to intracellular bacterial colonization, and (3) changes to cyclooxygenase-2-dependent inflammation. Furthermore, in mice with a history of chronic UTI, cyclooxygenase-2-dependent inflammation allowed a variety of clinical E. coli isolates to circumvent intracellular colonization resistance and cause severe recurrent UTI, which could be prevented by cyclooxygenase-2 inhibition or vaccination. This work provides mechanistic insight into how a history of infection can impact the risk for developing recurrent infection and has implications for the development of therapeutics for recurrent UTI.

Borrelia burgdorferi, the agent of Lyme disease, differentially expresses numerous genes and proteins as it cycles between mammalian hosts and tick vectors. Insights on regulatory mechanisms have been provided by earlier studies that examined B. burgdorferi gene expression patterns during cultivation. However, prior studies examined bacteria at only a single time point of cultivation, providing only a snapshot of what is likely a dynamic transcriptional program driving B. burgdorferi adaptations to changes during culture growth phases. To address that concern, we performed RNA sequencing (RNA-Seq) analysis of B. burgdorferi cultures at early-exponential, mid-exponential, and early-stationary phases of growth. We found that expression of nearly 18% of annotated B. burgdorferi genes changed significantly during culture maturation. Moreover, genome-wide mapping of the B. burgdorferi transcriptome in different growth phases enabled insight on transcript boundaries, operon structures, and identified numerous putative non-coding RNAs. These RNA-Seq data are discussed and presented as a resource for the community of researchers seeking to better understand B. burgdorferi biology and pathogenesis., Competing Interests: C. A. Brissette, J. Seshu, and B. Stevenson are Academic Editors of PLOS ONE. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

The ability to simultaneously characterize the bacterial and host expression programs during infection would facilitate a comprehensive understanding of pathogen-host interactions. Although RNA sequencing (RNA-seq) has greatly advanced our ability to study the transcriptomes of prokaryotes and eukaryotes separately, limitations in existing protocols for the generation and analysis of RNA-seq data have hindered simultaneous profiling of host and bacterial pathogen transcripts from the same sample. Here we provide a detailed protocol for simultaneous analysis of host and bacterial transcripts by RNA-seq. Importantly, this protocol details the steps required for efficient host and bacteria lysis, barcoding of samples, technical advances in sample preparation for low-yield sample inputs and a computational pipeline for analysis of both mammalian and microbial reads from mixed host-pathogen RNA-seq data. Sample preparation takes 3 d from cultured cells to pooled libraries. Data analysis takes an additional day. Compared with previous methods, the protocol detailed here provides a sensitive, facile and generalizable approach that is suitable for large-scale studies and will enable the field to obtain in-depth analysis of host-pathogen interactions in infection models.

The global regulator CodY controls the expression of dozens of metabolism and virulence genes in the opportunistic pathogen Staphylococcus aureus in response to the availability of isoleucine, leucine and valine (ILV), and GTP. Using RNA-Seq transcriptional profiling and partial activity variants, we reveal that S. aureus CodY activity grades metabolic and virulence gene expression as a function of ILV availability, mediating metabolic reorganization and controlling virulence factor production in vitro. Strains lacking CodY regulatory activity produce a PIA-dependent biofilm, but development is restricted under conditions that confer partial CodY activity. CodY regulates the expression of thermonuclease (nuc) via the Sae two-component system, revealing cascading virulence regulation and factor production as CodY activity is reduced. Proteins that mediate the host-pathogen interaction and subvert the immune response are shut off at intermediate levels of CodY activity, while genes coding for enzymes and proteins that extract nutrients from tissue, that kill host cells, and that synthesize amino acids are among the last genes to be derepressed. We conclude that S. aureus uses CodY to limit host damage to only the most severe starvation conditions, providing insight into one potential mechanism by which S. aureus transitions from a commensal bacterium to an invasive pathogen., (© 2016 John Wiley & Sons Ltd.)

Although RNA-seq is a powerful tool, the considerable time and cost associated with library construction has limited its utilization for various applications. RNAtag-Seq, an approach to generate multiple RNA-seq libraries in a single reaction, lowers time and cost per sample, and it produces data on prokaryotic and eukaryotic samples that are comparable to those generated by traditional strand-specific RNA-seq approaches.

Vibrio parahaemolyticus is the leading worldwide cause of seafood-associated gastroenteritis, yet little is known regarding its intraintestinal gene expression or physiology. To date, in vivo analyses have focused on identification and characterization of virulence factors--e.g. a crucial Type III secretion system (T3SS2)--rather than genome-wide analyses of in vivo biology. Here, we used RNA-Seq to profile V. parahaemolyticus gene expression in infected infant rabbits, which mimic human infection. Comparative transcriptomic analysis of V. parahaemolyticus isolated from rabbit intestines and from several laboratory conditions enabled identification of mRNAs and sRNAs induced during infection and of regulatory factors that likely control them. More than 12% of annotated V. parahaemolyticus genes are differentially expressed in the intestine, including the genes of T3SS2, which are likely induced by bile-mediated activation of the transcription factor VtrB. Our analyses also suggest that V. parahaemolyticus has access to glucose or other preferred carbon sources in vivo, but that iron is inconsistently available. The V. parahaemolyticus transcriptional response to in vivo growth is far more widespread than and largely distinct from that of V. cholerae, likely due to the distinct ways in which these diarrheal pathogens interact with and modulate the environment in the small intestine., (© The Author(s) 2014. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Unlabelled: Chaperone-usher pathway (CUP) pili are extracellular organelles produced by Gram-negative bacteria that mediate bacterial pathogenesis. Small-molecule inhibitors of CUP pili, termed pilicides, were rationally designed and shown to inhibit type 1 or P piliation. Here, we show that pilicide ec240 decreased the levels of type 1, P, and S piliation. Transcriptomic and proteomic analyses using the cystitis isolate UTI89 revealed that ec240 dysregulated CUP pili and decreased motility. Paradoxically, the transcript levels of P and S pilus genes were increased during growth in ec240, even though the level of P and S piliation decreased. In contrast, the most downregulated transcripts after growth in ec240 were from the type 1 pilus genes. Type 1 pilus expression is controlled by inversion of the fimS promoter element, which can oscillate between phase on and phase off orientations. ec240 induced the fimS phase off orientation, and this effect was necessary for the majority of ec240's inhibition of type 1 piliation. ec240 increased levels of the transcriptional regulators SfaB and PapB, which were shown to induce the fimS promoter phase off orientation. Furthermore, the effect of ec240 on motility was abolished in the absence of the SfaB, PapB, SfaX, and PapX regulators. In contrast to the effects of ec240, deletion of the type 1 pilus operon led to increased S and P piliation and motility. Thus, ec240 dysregulated several uropathogenic Escherichia coli (UPEC) virulence factors through different mechanisms and independent of its effects on type 1 pilus biogenesis and may have potential as an antivirulence compound., Importance: CUP pili and flagella play active roles in the pathogenesis of a variety of Gram-negative bacterial infections, including urinary tract infections mediated by UPEC. These are extremely common infections that are often recurrent and increasingly caused by antibiotic-resistant organisms. Preventing piliation and motility through altered regulation and assembly of these important virulence factors could aid in the development of novel therapeutics. This study increases our understanding of the regulation of these virulence factors, providing new avenues by which to target their expression., (Copyright © 2014 Greene et al.)

Many small, noncoding RNAs (sRNAs) in bacteria act as posttranscriptional regulators of messenger RNAs. TargetRNA2 is a web server that identifies mRNA targets of sRNA regulatory action in bacteria. As input, TargetRNA2 takes the sequence of an sRNA and the name of a sequenced bacterial replicon. When searching for targets of RNA regulation, TargetRNA2 uses a variety of features, including conservation of the sRNA in other bacteria, the secondary structure of the sRNA, the secondary structure of each candidate mRNA target and the hybridization energy between the sRNA and each candidate mRNA target. TargetRNA2 outputs a ranked list of likely regulatory targets for the input sRNA. When evaluated on a comprehensive set of sRNA-target interactions, TargetRNA2 was found to be both accurate and efficient in identifying targets of sRNA regulatory action. Furthermore, TargetRNA2 has the ability to integrate RNA-seq data, if available. If an sRNA is differentially expressed in two or more RNA-seq experiments, TargetRNA2 considers co-differential gene expression when searching for regulatory targets, significantly improving the accuracy of target identifications. The TargetRNA2 web server is freely available for use at http://cs.wellesley.edu/∼btjaden/TargetRNA2., (© The Author(s) 2014. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Global regulators that bind strategic metabolites allow bacteria to adapt rapidly to dynamic environments by coordinating the expression of many genes. We report an approach for determining gene regulation hierarchy using the regulon of the Bacillus subtilis global regulatory protein CodY as proof of principle. In theory, this approach can be used to measure the dynamics of any bacterial transcriptional regulatory network that is affected by interaction with a ligand. In B. subtilis, CodY controls dozens of genes, but the threshold activities of CodY required to regulate each gene are unknown. We hypothesized that targets of CodY are differentially regulated based on varying affinity for the protein's many binding sites. We used RNA sequencing to determine the transcription profiles of B. subtilis strains expressing mutant CodY proteins with different levels of residual activity. In parallel, we quantified intracellular metabolites connected to central metabolism. Strains producing CodY variants F71Y, R61K, and R61H retained varying degrees of partial activity relative to the WT protein, leading to gene-specific, differential alterations in transcript abundance for the 223 identified members of the CodY regulon. Using liquid chromatography coupled to MS, we detected significant increases in branched-chain amino acids and intermediates of arginine, proline, and glutamate metabolism, as well as decreases in pyruvate and glycerate as CodY activity decreased. We conclude that a spectrum of CodY activities leads to programmed regulation of gene expression and an apparent rerouting of carbon and nitrogen metabolism, suggesting that during changes in nutrient availability, CodY prioritizes the expression of specific pathways.

Bacterial small RNAs (sRNAs) are short transcripts that typically do not encode proteins and often act as regulators of gene expression through a variety of mechanisms. Regulatory sRNAs have been identified in many species, including Mycobacterium tuberculosis, the causative agent of tuberculosis. Here, we use a computational algorithm to predict sRNA candidates in the mycobacterial species M. smegmatis and M. bovis BCG and confirmed the expression of many sRNAs using Northern blotting. Thus, we have identified 17 and 23 novel sRNAs in M. smegmatis and M. bovis BCG, respectively. We have also applied a high-throughput technique (Deep-RACE) to map the 5' and 3' ends of many of these sRNAs and identified potential regulators of sRNAs by analysis of existing ChIP-seq datasets. The sRNAs identified in this work likely contribute to the unique biology of mycobacteria.