Cross-species RNA-seq for deciphering host-microbe interactions.

The human body is constantly exposed to microorganisms, which entails manifold interactions between human cells and diverse commensal or pathogenic bacteria. The cellular states of the interacting cells are decisive for the outcome of these encounters such as whether bacterial virulence programmes and host defence or tolerance mechanisms are induced. This Review summarizes how next-generation RNA sequencing (RNA-seq) has become a primary technology to study host-microbe interactions with high resolution, improving our understanding of the physiological consequences and the mechanisms at play. We illustrate how the discriminatory power and sensitivity of RNA-seq helps to dissect increasingly complex cellular interactions in time and space down to the single-cell level. We also outline how future transcriptomics may answer currently open questions in host-microbe interactions and inform treatment schemes for microbial disorders.

Improved RNA stability estimation through Bayesian modeling reveals most Salmonella transcripts have subminute half-lives.

RNA decay is a crucial mechanism for regulating gene expression in response to environmental stresses. In bacteria, RNA-binding proteins (RBPs) are known to be involved in posttranscriptional regulation, but their global impact on RNA half-lives has not been extensively studied. To shed light on the role of the major RBPs ProQ and CspC/E in maintaining RNA stability, we performed RNA sequencing of Salmonella enterica over a time course following treatment with the transcription initiation inhibitor rifampicin (RIF-seq) in the presence and absence of these RBPs. We developed a hierarchical Bayesian model that corrects for confounding factors in rifampicin RNA stability assays and enables us to identify differentially decaying transcripts transcriptome-wide. Our analysis revealed that the median RNA half-life in Salmonella in early stationary phase is less than 1 min, a third of previous estimates. We found that over half of the 500 most long-lived transcripts are bound by at least one major RBP, suggesting a general role for RBPs in shaping the transcriptome. Integrating differential stability estimates with cross-linking and immunoprecipitation followed by RNA sequencing (CLIP-seq) revealed that approximately 30% of transcripts with ProQ binding sites and more than 40% with CspC/E binding sites in coding or 3' untranslated regions decay differentially in the absence of the respective RBP. Analysis of differentially destabilized transcripts identified a role for ProQ in the oxidative stress response. Our findings provide insights into posttranscriptional regulation by ProQ and CspC/E, and the importance of RBPs in regulating gene expression.

CRISPR-based screening of small RNA modulators of bile susceptibility in Bacteroides thetaiotaomicron.

Microbiota-centric interventions are limited by our incomplete understanding of the gene functions of many of its constituent species. This applies in particular to small RNAs (sRNAs), which are emerging as important regulators in microbiota species yet tend to be missed by traditional functional genomics approaches. Here, we establish CRISPR interference (CRISPRi) in the abundant microbiota member Bacteroides thetaiotaomicron for genome-wide sRNA screens. By assessing the abundance of different protospacer-adjacent motifs, we identify the Prevotella bryantii B14 Cas12a as a suitable nuclease for CRISPR screens in these bacteria and generate an inducible Cas12a expression system. Using a luciferase reporter strain, we infer guide design rules and use this knowledge to assemble a computational pipeline for automated gRNA design. By subjecting the resulting guide library to a phenotypic screen, we uncover the sRNA BatR to increase susceptibility to bile salts through the regulation of genes involved in Bacteroides cell surface structure. Our study lays the groundwork for unlocking the genetic potential of these major human gut mutualists and, more generally, for identifying hidden functions of bacterial sRNAs.

Group B Streptococcus transcriptome when interacting with brain endothelial cells.

Bacterial meningitis is a life-threatening infection of the central nervous system (CNS) that occurs when bacteria are able to cross the blood-brain barrier (BBB) or the meningeal-cerebrospinal fluid barrier (mBCSFB). The BBB and mBCSFB comprise highly specialized brain endothelial cells (BECs) that typically restrict pathogen entry. Group B Streptococcus (GBS or Streptococcus agalactiae) is the leading cause of neonatal meningitis. Until recently, identification of GBS virulence factors has relied on genetic screening approaches. Instead, we here conducted RNA-seq analysis on GBS when interacting with induced pluripotent stem cell-derived BECs (iBECs) to pinpoint virulence-associated genes. Of the 2,068 annotated protein-coding genes of GBS, 430 transcripts displayed significant changes in expression after interacting with BECs. Notably, we found that the majority of differentially expressed GBS transcripts were downregulated (360 genes) during infection of iBECs. Interestingly, codY, encoding a pleiotropic transcriptional repressor in low-G + C Gram-positive bacteria, was identified as being highly downregulated. We conducted qPCR to confirm the codY downregulation observed via RNA-seq during the GBS-iBEC interaction and obtained codY mutants in three different GBS background parental strains. As anticipated from the RNA-seq results, the [Formula: see text]codY strains were more adherent and invasive in two in vitro BEC models. Together, this demonstrates the utility of RNA-seq during the BEC interaction to identify GBS virulence modulators. IMPORTANCE: Group B Streptococcus (GBS) meningitis remains the leading cause of neonatal meningitis. Research work has identified surface factors and two-component systems that contribute to GBS disruption of the blood-brain barrier (BBB). These discoveries often relied on genetic screening approaches. Here, we provide transcriptomic data describing how GBS changes its transcriptome when interacting with brain endothelial cells. Additionally, we have phenotypically validated these data by obtaining mutants of a select regulator that is highly down-regulated during infection and testing on our BBB model. This work provides the research field with a validated data set that can provide an insight into potential pathways that GBS requires to interact with the BBB and open the door to new discoveries.

Comparative genomics provides structural and functional insights into Bacteroides RNA biology.

Bacteria employ noncoding RNA molecules for a wide range of biological processes, including scaffolding large molecular complexes, catalyzing chemical reactions, defending against phages, and controlling gene expression. Secondary structures, binding partners, and molecular mechanisms have been determined for numerous small noncoding RNAs (sRNAs) in model aerobic bacteria. However, technical hurdles have largely prevented analogous analyses in the anaerobic gut microbiota. While experimental techniques are being developed to investigate the sRNAs of gut commensals, computational tools and comparative genomics can provide immediate functional insight. Here, using Bacteroides thetaiotaomicron as a representative microbiota member, we illustrate how comparative genomics improves our understanding of RNA biology in an understudied gut bacterium. We investigate putative RNA-binding proteins and predict a Bacteroides cold-shock protein homolog to have an RNA-related function. We apply an in silico protocol incorporating both sequence and structural analysis to determine the consensus structures and conservation of nine Bacteroides noncoding RNA families. Using structure probing, we validate and refine these predictions and deposit them in the Rfam database. Through synteny analyses, we illustrate how genomic coconservation can serve as a predictor of sRNA function. Altogether, this work showcases the power of RNA informatics for investigating the RNA biology of anaerobic microbiota members.

Improved bacterial RNA-seq by Cas9-based depletion of ribosomal RNA reads.

A major challenge for RNA-seq analysis of gene expression is to achieve sufficient coverage of informative nonribosomal transcripts. In eukaryotic samples, this is typically achieved by selective oligo(dT)-priming of messenger RNAs to exclude ribosomal RNA (rRNA) during cDNA synthesis. However, this strategy is not compatible with prokaryotes in which functional transcripts are generally not polyadenylated. To overcome this, we adopted DASH (depletion of abundant sequences by hybridization), initially developed for eukaryotic cells, to improve both the sensitivity and depth of bacterial RNA-seq. DASH uses the Cas9 nuclease to remove unwanted cDNA sequences prior to library amplification. We report the design, evaluation, and optimization of DASH experiments for standard bacterial short-read sequencing approaches, including software for automated guide RNA (gRNA) design for Cas9-mediated cleavage in bacterial rDNA sequences. Using these gRNA pools, we effectively removed rRNA reads (56%-86%) in RNA-seq libraries from two different model bacteria, the Gram-negative pathogen Salmonella enterica and the anaerobic gut commensal Bacteroides thetaiotaomicron DASH works robustly, even with subnanogram amounts of input RNA. Its efficiency, high sensitivity, ease of implementation, and low cost (∼$5 per sample) render DASH an attractive alternative to rRNA removal protocols, in particular for material-constrained studies where conventional ribodepletion techniques fail.

Dual RNA-seq unveils noncoding RNA functions in host-pathogen interactions.

Bacteria express many small RNAs for which the regulatory roles in pathogenesis have remained poorly understood due to a paucity of robust phenotypes in standard virulence assays. Here we use a generic 'dual RNA-seq' approach to profile RNA expression simultaneously in pathogen and host during Salmonella enterica serovar Typhimurium infection and reveal the molecular impact of bacterial riboregulators. We identify a PhoP-activated small RNA, PinT, which upon bacterial internalization temporally controls the expression of both invasion-associated effectors and virulence genes required for intracellular survival. This riboregulatory activity causes pervasive changes in coding and noncoding transcripts of the host. Interspecies correlation analysis links PinT to host cell JAK-STAT signalling, and we identify infection-specific alterations in multiple long noncoding RNAs. Our study provides a paradigm for a sensitive RNA-based analysis of intracellular bacterial pathogens and their hosts without physical separation, as well as a new discovery route for hidden functions of pathogen genes.

CRP-cAMP mediates silencing of Salmonella virulence at the post-transcriptional level.

Invasion of epithelial cells by Salmonella enterica requires expression of genes located in the pathogenicity island I (SPI-1). The expression of SPI-1 genes is very tightly regulated and activated only under specific conditions. Most studies have focused on the regulatory pathways that induce SPI-1 expression. Here, we describe a new regulatory circuit involving CRP-cAMP, a widely established metabolic regulator, in silencing of SPI-1 genes under non-permissive conditions. In CRP-cAMP-deficient strains we detected a strong upregulation of SPI-1 genes in the mid-logarithmic growth phase. Genetic analyses revealed that CRP-cAMP modulates the level of HilD, the master regulator of Salmonella invasion. This regulation occurs at the post-transcriptional level and requires the presence of a newly identified regulatory motif within the hilD 3'UTR. We further demonstrate that in Salmonella the Hfq-dependent sRNA Spot 42 is under the transcriptional repression of CRP-cAMP and, when this transcriptional repression is relieved, Spot 42 exerts a positive effect on hilD expression. In vivo and in vitro assays indicate that Spot 42 targets, through its unstructured region III, the 3'UTR of the hilD transcript. Together, our results highlight the biological relevance of the hilD 3'UTR as a hub for post-transcriptional control of Salmonella invasion gene expression.

An RNA-centric view on gut Bacteroidetes.

Bacteria employ noncoding RNAs to maintain cellular physiology, adapt global gene expression to fluctuating environments, sense nutrients, coordinate their interaction with companion microbes and host cells, and protect themselves against bacteriophages. While bacterial RNA research has made fundamental contributions to biomedicine and biotechnology, the bulk of our knowledge of RNA biology stems from the study of a handful of aerobic model species. In comparison, RNA research is lagging in many medically relevant obligate anaerobic species, in particular the numerous commensal bacteria comprising our gut microbiota. This review presents a guide to RNA-based regulatory mechanisms in the phylum Bacteroidetes, focusing on the most abundant bacterial genus in the human gut, Bacteroides spp. This includes recent case reports on riboswitches, an mRNA leader, cis- and trans-encoded small RNAs (sRNAs) in Bacteroides spp., and a survey of CRISPR-Cas systems across Bacteroidetes. Recent work from our laboratory now suggests the existence of hundreds of noncoding RNA candidates in Bacteroides thetaiotaomicron, the emerging model organism for functional microbiota research. Based on these collective observations, we predict mechanistic and functional commonalities and differences between Bacteroides sRNAs and those of other model bacteria, and outline open questions and tools needed to boost Bacteroidetes RNA research.

RNA target profiles direct the discovery of virulence functions for the cold-shock proteins CspC and CspE.

The functions of many bacterial RNA-binding proteins remain obscure because of a lack of knowledge of their cellular ligands. Although well-studied cold-shock protein A (CspA) family members are induced and function at low temperature, others are highly expressed in infection-relevant conditions. Here, we have profiled transcripts bound in vivo by the CspA family members of Salmonella enterica serovar Typhimurium to link the constitutively expressed CspC and CspE proteins with virulence pathways. Phenotypic assays in vitro demonstrated a crucial role for these proteins in membrane stress, motility, and biofilm formation. Moreover, double deletion of cspC and cspE fully attenuates Salmonella in systemic mouse infection. In other words, the RNA ligand-centric approach taken here overcomes a problematic molecular redundancy of CspC and CspE that likely explains why these proteins have evaded selection in previous virulence factor screens in animals. Our results highlight RNA-binding proteins as regulators of pathogenicity and potential targets of antimicrobial therapy. They also suggest that globally acting RNA-binding proteins are more common in bacteria than currently appreciated.

Resolving host-pathogen interactions by dual RNA-seq.

The transcriptome is a powerful proxy for the physiological state of a cell, healthy or diseased. As a result, transcriptome analysis has become a key tool in understanding the molecular changes that accompany bacterial infections of eukaryotic cells. Until recently, such transcriptomic studies have been technically limited to analyzing mRNA expression changes in either the bacterial pathogen or the infected eukaryotic host cell. However, the increasing sensitivity of high-throughput RNA sequencing now enables "dual RNA-seq" studies, simultaneously capturing all classes of coding and noncoding transcripts in both the pathogen and the host. In the five years since the concept of dual RNA-seq was introduced, the technique has been applied to a range of infection models. This has not only led to a better understanding of the physiological changes in pathogen and host during the course of an infection but has also revealed hidden molecular phenotypes of virulence-associated small noncoding RNAs that were not visible in standard infection assays. Here, we use the knowledge gained from these recent studies to suggest experimental and computational guidelines for the design of future dual RNA-seq studies. We conclude this review by discussing prospective applications of the technique.

An NK Cell Perforin Response Elicited via IL-18 Controls Mucosal Inflammation Kinetics during Salmonella Gut Infection.

Salmonella Typhimurium (S.Tm) is a common cause of self-limiting diarrhea. The mucosal inflammation is thought to arise from a standoff between the pathogen's virulence factors and the host's mucosal innate immune defenses, particularly the mucosal NAIP/NLRC4 inflammasome. However, it had remained unclear how this switches the gut from homeostasis to inflammation. This was studied using the streptomycin mouse model. S.Tm infections in knockout mice, cytokine inhibition and -injection experiments revealed that caspase-1 (not -11) dependent IL-18 is pivotal for inducing acute inflammation. IL-18 boosted NK cell chemoattractants and enhanced the NK cells' migratory capacity, thus promoting mucosal accumulation of mature, activated NK cells. NK cell depletion and Prf-/- ablation (but not granulocyte-depletion or T-cell deficiency) delayed tissue inflammation. Our data suggest an NK cell perforin response as one limiting factor in mounting gut mucosal inflammation. Thus, IL-18-elicited NK cell perforin responses seem to be critical for coordinating mucosal inflammation during early infection, when S.Tm strongly relies on virulence factors detectable by the inflammasome. This may have broad relevance for mucosal defense against microbial pathogens.

Single-cell RNA-seq: advances and future challenges.

Phenotypically identical cells can dramatically vary with respect to behavior during their lifespan and this variation is reflected in their molecular composition such as the transcriptomic landscape. Single-cell transcriptomics using next-generation transcript sequencing (RNA-seq) is now emerging as a powerful tool to profile cell-to-cell variability on a genomic scale. Its application has already greatly impacted our conceptual understanding of diverse biological processes with broad implications for both basic and clinical research. Different single-cell RNA-seq protocols have been introduced and are reviewed here-each one with its own strengths and current limitations. We further provide an overview of the biological questions single-cell RNA-seq has been used to address, the major findings obtained from such studies, and current challenges and expected future developments in this booming field.

Dual 3'Seq using deepSuperSAGE uncovers transcriptomes of interacting Salmonella enterica Typhimurium and human host cells.

BACKGROUND: The interaction of eukaryotic host and prokaryotic pathogen cells is linked to specific changes in the cellular proteome, and consequently to infection-related gene expression patterns of the involved cells. To simultaneously assess the transcriptomes of both organisms during their interaction we developed dual 3'Seq, a tag-based sequencing protocol that allows for exact quantification of differentially expressed transcripts in interacting pro- and eukaryotic cells without prior fixation or physical disruption of the interaction. RESULTS: Human epithelial cells were infected with Salmonella enterica Typhimurium as a model system for invasion of the intestinal epithelium, and the transcriptional response of the infected host cells together with the differential expression of invading and intracellular pathogen cells was determined by dual 3'Seq coupled with the next-generation sequencing-based transcriptome profiling technique deepSuperSAGE (deep Serial Analysis of Gene Expression). Annotation to reference transcriptomes comprising the operon structure of the employed S. enterica Typhimurium strain allowed for in silico separation of the interacting cells including quantification of polycistronic RNAs. Eighty-nine percent of the known loci are found to be transcribed in prokaryotic cells prior or subsequent to infection of the host, while 75% of all protein-coding loci are represented in the polyadenylated transcriptomes of human host cells. CONCLUSIONS: Dual 3'Seq was alternatively coupled to MACE (Massive Analysis of cDNA ends) to assess the advantages and drawbacks of a library preparation procedure that allows for sequencing of longer fragments. Additionally, the identified expression patterns of both organisms were validated by qRT-PCR using three independent biological replicates, which confirmed that RELB along with NFKB1 and NFKB2 are involved in the initial immune response of epithelial cells after infection with S. enterica Typhimurium.

Differential activation and functional specialization of miR-146 and miR-155 in innate immune sensing.

Many microRNAs (miRNAs) are co-regulated during the same physiological process but the underlying cellular logic is often little understood. The conserved, immunomodulatory miRNAs miR-146 and miR-155, for instance, are co-induced in many cell types in response to microbial lipopolysaccharide (LPS) to feedback-repress LPS signalling through Toll-like receptor TLR4. Here, we report that these seemingly co-induced regulatory RNAs dramatically differ in their induction behaviour under various stimuli strengths and act non-redundantly through functional specialization; although miR-146 expression saturates at sub-inflammatory doses of LPS that do not trigger the messengers of inflammation markers, miR-155 remains tightly associated with the pro-inflammatory transcriptional programmes. Consequently, we found that both miRNAs control distinct mRNA target profiles; although miR-146 targets the messengers of LPS signal transduction components and thus downregulates cellular LPS sensitivity, miR-155 targets the mRNAs of genes pervasively involved in pro-inflammatory transcriptional programmes. Thus, miR-155 acts as a broad limiter of pro-inflammatory gene expression once the miR-146 dependent barrier to LPS triggered inflammation has been breached. Importantly, we also report alternative miR-155 activation by the sensing of bacterial peptidoglycan through cytoplasmic NOD-like receptor, NOD2. We predict that dose-dependent responses to environmental stimuli may involve functional specialization of seemingly co-induced miRNAs in other cellular circuitries as well.

Transgene regulation in plants by alternative splicing of a suicide exon.

Compared to transcriptional activation, other mechanisms of gene regulation have not been widely exploited for the control of transgenes. One barrier to the general use and application of alternative splicing is that splicing-regulated transgenes have not been shown to be reliably and simply designed. Here, we demonstrate that a cassette bearing a suicide exon can be inserted into a variety of open reading frames (ORFs), generating transgenes whose expression is activated by exon skipping in response to a specific protein inducer. The surprisingly minimal sequence requirements for the maintenance of splicing fidelity and regulation indicate that this splicing cassette can be used to regulate any ORF containing one of the amino acids Glu, Gln or Lys. Furthermore, a single copy of the splicing cassette was optimized by rational design to confer robust gene activation with no background expression in plants. Thus, conditional splicing has the potential to be generally useful for transgene regulation.

CRISPR Interference-Based Functional Small RNA Genomics.

Small RNAs (sRNAs) are versatile regulators universally present in species across the prokaryotic kingdom, yet their functional characterization remains a major bottleneck. Gene inactivation through random transposon insertion has proven extremely valuable in discovering hidden gene functions. However, this approach is biased toward long genes and usually results in the underrepresentation of sRNA mutants. In contrast, CRISPR interference (CRISPRi) harnesses guide RNAs to recruit cleavage-deficient Cas nucleases to specific DNA loci. The ensuing steric hindrance inhibits RNA polymerase assembly at-or migration along-predefined genes, allowing for targeted knockdown screens without major length bias. In this chapter, we provide a detailed protocol for CRISPRi-based functional screening of bacterial sRNAs. Using the abundant microbiota species Bacteroides thetaiotaomicron as a model, we describe the design and generation of a guide library targeting the full intergenic sRNA repertoire of this organism and its application to identify sRNA knockdown-associated fitness effects. Our protocol is generic and thus suitable for the systematic assessment of sRNA-associated phenotypes in a wide range of bacterial species and experimental conditions. We expect CRISPRi-based functional genomics to boost sRNA research in understudied bacterial taxa, for instance, members of the gut microbiota.

Non-professional efferocytosis of Salmonella-infected intestinal epithelial cells in the neonatal host.

The intestinal epithelium is the first line of defense against enteric pathogens. Removal of infected cells by exfoliation prevents mucosal translocation and systemic infection in the adult host, but is less commonly observed in the neonatal intestine. Instead, here, we describe non-professional efferocytosis of Salmonella-infected enterocytes by neighboring epithelial cells in the neonatal intestine. Intestinal epithelial stem cell organoid cocultures of neonatal and adult cell monolayers with damaged enterocytes replicated this observation, confirmed the age-dependent ability of intestinal epithelial cells for efferocytosis, and identified the involvement of the "eat-me" signals and adaptors phosphatidylserine and C1q as well as the "eat-me" receptors integrin-αv (CD51) and CD36 in cellular uptake. Consistent with this, massive epithelial cell membrane protrusions and CD36 accumulation at the contact site with apoptotic cells were observed in the infected neonatal host in vivo. Efferocytosis of infected small intestinal enterocytes by neighboring epithelial cells may represent a previously unrecognized mechanism of neonatal antimicrobial host defense to maintain barrier integrity.

An expanded transcriptome atlas for Bacteroides thetaiotaomicron reveals a small RNA that modulates tetracycline sensitivity.

Plasticity in gene expression allows bacteria to adapt to diverse environments. This is particularly relevant in the dynamic niche of the human intestinal tract; however, transcriptional networks remain largely unknown for gut-resident bacteria. Here we apply differential RNA sequencing (RNA-seq) and conventional RNA-seq to the model gut bacterium Bacteroides thetaiotaomicron to map transcriptional units and profile their expression levels across 15 in vivo-relevant growth conditions. We infer stress- and carbon source-specific transcriptional regulons and expand the annotation of small RNAs (sRNAs). Integrating this expression atlas with published transposon mutant fitness data, we predict conditionally important sRNAs. These include MasB, which downregulates tetracycline tolerance. Using MS2 affinity purification and RNA-seq, we identify a putative MasB target and assess its role in the context of the MasB-associated phenotype. These data-publicly available through the Theta-Base web browser ( http://micromix.helmholtz-hiri.de/bacteroides/ )-constitute a valuable resource for the microbiome community.