A Laboratory Methodology for Dual RNA-Sequencing of Bacteria and their Host Cells In Vitro

Non-Targeted RNA Sequencing: Towards the Development of Universal Clinical Diagnosis Methods for Human and Veterinary Infectious Diseases

Metagenomics offers the potential to replace and simplify classical methods used in the clinical diagnosis of human and veterinary infectious diseases. Metagenomics boasts a high pathogen discovery rate and high specificity, advantages absent in most classical approaches. However, its widespread adoption in clinical settings is still pending, with a slow transition from research to routine use. While longer turnaround times and higher costs were once concerns, these issues are currently being addressed by automation, better chemistries, improved sequencing platforms, better databases, and automated bioinformatics analysis. However, many technical options and steps, each producing highly variable outcomes, have reduced the technology's operational value, discouraging its implementation in diagnostic labs. We present a case for utilizing non-targeted RNA sequencing (NT-RNA-seq) as an ideal metagenomics method for the detection of infectious disease-causing agents in humans and animals. Additionally, to create operational value, we propose to identify best practices for the "core" of steps that are invariably shared among many human and veterinary protocols. Reference materials, sequencing procedures, and bioinformatics standards should accelerate the validation processes necessary for the widespread adoption of this technology. Best practices could be determined through "implementation research" by a consortium of interested institutions working on common samples.

Paired single-cell host profiling with multiplex-tagged bacterial mutants reveals intracellular virulence-immune networks

Encounters between host cells and intracellular bacterial pathogens lead to complex phenotypes that determine the outcome of infection. Single-cell RNA sequencing (scRNA-seq) is increasingly used to study the host factors underlying diverse cellular phenotypes but has limited capacity to analyze the role of bacterial factors. Here, we developed scPAIR-seq, a single-cell approach to analyze infection with a pooled library of multiplex-tagged, barcoded bacterial mutants. Infected host cells and barcodes of intracellular bacterial mutants are both captured by scRNA-seq to functionally analyze mutant-dependent changes in host transcriptomes. We applied scPAIR-seq to macrophages infected with a library of Salmonella Typhimurium secretion system effector mutants. We analyzed redundancy between effectors and mutant-specific unique fingerprints and mapped the global virulence network of each individual effector by its impact on host immune pathways. ScPAIR-seq is a powerful tool to untangle bacterial virulence strategies and their complex interplay with host defense strategies that drive infection outcome.

Identification of Burkholderia cenocepacia non-coding RNAs expressed during Caenorhabditis elegans infection

Small non-coding RNAs (sRNAs) are key regulators of post-transcriptional gene expression in bacteria. Despite the identification of hundreds of bacterial sRNAs, their roles on bacterial physiology and virulence remain largely unknown, as is the case of bacteria of the Burkholderia cepacia complex (Bcc). Bcc is a group of opportunistic pathogens with relatively large genomes that can cause lethal lung infections amongst cystic fibrosis (CF) patients. To characterise sRNAs expressed by Bcc bacteria when infecting a host, the nematode Caenorhabditis elegans was used as an infection model by the epidemic CF strain B. cenocepacia J2315. A total of 108 new and 31 previously described sRNAs with a predicted Rho independent terminator were identified, most of them located on chromosome 1. RIT11b, a sRNA downregulated under C. elegans infection conditions, was shown to directly affect B. cenocepacia virulence, biofilm formation, and swimming motility. RIT11b overexpression reduced the expression of the direct targets dusA and pyrC, involved in biofilm formation, epithelial cell adherence, and chronic infections in other organisms. The in vitro direct interaction of RIT11b with the dusA and pyrC messengers was demonstrated by electrophoretic mobility shift assays. To the best of our knowledge this is the first report on the functional characterization of a sRNA directly involved in B. cenocepacia virulence. KEY POINTS: • 139 sRNAs expressed by B. cenocepacia during C. elegans infection were identified • The sRNA RIT11b affects B. cenocepacia virulence, biofilm formation, and motility • RIT11b directly binds to and regulates dusA and pyrC mRNAs.

Transitioning from Soil to Host: Comparative Transcriptome Analysis Reveals the Burkholderia pseudomallei Response to Different Niches

Burkholderia pseudomallei, a soil and water saprophyte, is responsible for the tropical human disease melioidosis. A hundred years since its discovery, there is still much to learn about B. pseudomallei proteins that are essential for the bacterium's survival in and interaction with the infected host, as well as their roles within the bacterium's natural soil habitat. To address this gap, bacteria grown under conditions mimicking the soil environment were subjected to transcriptome sequencing (RNA-seq) analysis. A dual RNA-seq approach was used on total RNA from spleens isolated from a B. pseudomallei mouse infection model at 5 days postinfection. Under these conditions, a total of 1,434 bacterial genes were induced, with 959 induced in the soil environment and 475 induced in bacteria residing within the host. Genes encoding metabolism and transporter proteins were induced when the bacteria were present in soil, while virulence factors, metabolism, and bacterial defense mechanisms were upregulated during active infection of mice. On the other hand, capsular polysaccharide and quorum-sensing pathways were inhibited during infection. In addition to virulence factors, reactive oxygen species, heat shock proteins, siderophores, and secondary metabolites were also induced to assist bacterial adaptation and survival in the host. Overall, this study provides crucial insights into the transcriptome-level adaptations which facilitate infection by soil-dwelling B. pseudomallei. Targeting novel therapeutics toward B. pseudomallei proteins required for adaptation provides an alternative treatment strategy given its intrinsic antimicrobial resistance and the absence of a vaccine. IMPORTANCE Burkholderia pseudomallei, a soil-dwelling bacterium, is the causative agent of melioidosis, a fatal infectious disease of humans and animals. The bacterium has a large genome consisting of two chromosomes carrying genes that encode proteins with important roles for survival in diverse environments as well as in the infected host. While a general mechanism of pathogenesis has been proposed, it is not clear which proteins have major roles when the bacteria are in the soil and whether the same proteins are key to successful infection and spread. To address this question, we grew the bacteria in soil medium and then in infected mice. At 5 days postinfection, bacteria were recovered from infected mouse organs and their gene expression was compared against that of bacteria grown in soil medium. The analysis revealed a list of genes expressed under soil growth conditions and a different set of genes encoding proteins which may be important for survival, replication, and dissemination in an infected host. These proteins are a potential resource for understanding the full adaptation mechanism of this pathogen. In the absence of a vaccine for melioidosis and with treatment being reliant on combinatorial antibiotic therapy, these proteins may be ideal targets for designing antimicrobials to treat melioidosis.

dRNASb: a systems biology approach to decipher dynamics of host-pathogen interactions using temporal dual RNA-seq data

Infection triggers a dynamic cascade of reciprocal events between host and pathogen wherein the host activates complex mechanisms to recognise and kill pathogens while the pathogen often adjusts its virulence and fitness to avoid eradication by the host. The interaction between the pathogen and the host results in large-scale changes in gene expression in both organisms. Dual RNA-seq, the simultaneous detection of host and pathogen transcripts, has become a leading approach to unravelling complex molecular interactions between the host and the pathogen and is particularly informative for intracellular organisms. The amount of in vitro and in vivo dual RNA-seq data is rapidly growing, which demands computational pipelines to effectively analyse such data. In particular, holistic, systems-level, and temporal analyses of dual RNA-seq data are essential to enable further insights into the host-pathogen transcriptional dynamics and potential interactions. Here, we developed an integrative network-driven bioinformatics pipeline, dRNASb, a systems biology-based computational pipeline to analyse temporal transcriptional clusters, incorporate molecular interaction networks (e.g. protein-protein interactions), identify topologically and functionally key transcripts in host and pathogen, and associate host and pathogen temporal transcriptome to decipher potential between-species interactions. The pipeline is applicable to various dual RNA-seq data from different species and experimental conditions. As a case study, we applied dRNASb to analyse temporal dual RNA-seq data of Salmonella-infected human cells, which enabled us to uncover genes contributing to the infection process and their potential functions and to identify putative associations between host and pathogen genes during infection. Overall, dRNASb has the potential to identify key genes involved in bacterial growth or host defence mechanisms for future uses as therapeutic targets.

Reducing costs for DNA and RNA sequencing by sample pooling using a metagenomic approach

DNA and RNA sequencing are widely used techniques to investigate genomic modifications and gene expression. The costs for sequencing dropped dramatically in the last decade. However, due to material and labor intense steps, the sample preparation costs could not keep up with that pace. About 80% of the total costs occur prior to sequencing during DNA/RNA extraction, enrichment steps and subsequent library preparation. In this study, we investigate the potential of pooling different organisms samples prior to DNA/RNA extraction to significantly reduce costs in preparative steps. Similar to the common procedure of ligated DNA tags to pool (c)DNA samples, sequence diversity of different organisms intrinsically provide unique sequences that allow separation of reads after sequencing. With this approach, sample pooling can occur before DNA/RNA isolation and library preparation. We show that pooled sequencing of three related bacterial organisms is possible without loss of data quality at a cost reduction of approx. 50% in DNA- and RNA-seq approaches. Furthermore, we show that this approach is highly efficient down to the level of a shared genus and is, therefore, widely applicable in sequencing facilities and companies with diverse sample pools.

Steroid-Refractory Gut Graft-Versus-Host Disease: What We Have Learned From Basic Immunology and Experimental Mouse Model

Intestinal graft-versus-host disease (Gut-GVHD) is one of the major causes of mortality after allogeneic hematopoietic stem cell transplantation (allo-HSCT). While systemic glucocorticoids (GCs) comprise the first-line treatment option, the response rate for GCs varies from 30% to 50%. The prognosis for patients with steroid-refractory acute Gut-GVHD (SR-Gut-aGVHD) remains dismal. The mechanisms underlying steroid resistance are unclear, and apart from ruxolitinib, there are no approved treatments for SR-Gut-aGVHD. In this review, we provide an overview of the current biological understanding of experimental SR-Gut-aGVHD pathogenesis, the advanced technology that can be applied to the human SR-Gut-aGVHD studies, and the potential novel therapeutic options for patients with SR-Gut-aGVHD.

Comparison of DNA/RNA yield and integrity between PMAP36-mediated and other bacterial lysis methods

Currently, several methods are available for the isolation of bacterial DNA and RNA. However, the diversity and complexity of cell envelope structures limit their efficiency depending on the target bacterial species. In this study, we compared the differences in yield and integrity of RNA prepared from four gram-negative and six gram-positive bacterial species using bead-beating, bacteriolytic protein, and PMAP36-vortexing methods. Similarly, we also compared the efficiency of DNA extraction from Staphylococcus aureus. Physical disruption of bacterial cells showed versatility in breaking cells against all tested species; however, a decrease in the integrity of isolated DNA and RNA was observed. Among membranolytic proteins, PMAP36 showed the most promising results, in terms of both the yield and integrity of the prepared nucleic acids. Our results show that each method has inherent advantages and disadvantages depending on its application. Therefore, the characteristics of each method and target species should be considered before the extraction of bacterial DNA and RNA.

Intestinal organoid cocultures with microbes

Adult-stem-cell-derived organoids model human epithelial tissues ex vivo, which enables the study of host-microbe interactions with great experimental control. This protocol comprises methods to coculture organoids with microbes, particularly focusing on human small intestinal and colon organoids exposed to individual bacterial species. Microinjection into the lumen and periphery of 3D organoids is discussed, as well as exposure of organoids to microbes in a 2D layer. We provide detailed protocols for characterizing the coculture with regard to bacterial and organoid cell viability and growth kinetics. Spatial relationships can be studied by fluorescence live microscopy, as well as scanning electron microscopy. Finally, we discuss considerations for assessing the impact of bacteria on gene expression and mutations through RNA and DNA sequencing. This protocol requires equipment for standard mammalian tissue culture, or bacterial or viral culture, as well as a microinjection device.

Dual RNASeq Reveals NTHi-Macrophage Transcriptomic Changes During Intracellular Persistence

Nontypeable Haemophilus influenzae (NTHi) is a pathobiont which chronically colonises the airway of individuals with chronic respiratory disease and is associated with poor clinical outcomes. It is unclear how NTHi persists in the airway, however accumulating evidence suggests that NTHi can invade and persist within macrophages. To better understand the mechanisms of NTHi persistence within macrophages, we developed an in vitro model of NTHi intracellular persistence using human monocyte-derived macrophages (MDM). Dual RNA Sequencing was used to assess MDM and NTHi transcriptomic regulation occurring simultaneously during NTHi persistence. Analysis of the macrophage response to NTHi identified temporally regulated transcriptomic profiles, with a specific 'core' profile displaying conserved expression of genes across time points. Gene list enrichment analysis identified enrichment of immune responses in the core gene set, with KEGG pathway analysis revealing specific enrichment of intracellular immune response pathways. NTHi persistence was facilitated by modulation of bacterial metabolic, stress response and ribosome pathways. Levels of NTHi genes bioC, mepM and dps were differentially expressed by intracellular NTHi compared to planktonic NTHi, indicating that the transcriptomic adaption was distinct between the two different NTHi lifestyles. Overall, this study provides crucial insights into the transcriptomic adaptations facilitating NTHi persistence within macrophages. Targeting these reported pathways with novel therapeutics to reduce NTHi burden in the airway could be an effective treatment strategy given the current antimicrobial resistance crisis and lack of NTHi vaccines.

Bacteria-host transcriptional response during endothelial invasion by Staphylococcus aureus

Staphylococcus aureus is the cause of serious vascular infections such as sepsis and endocarditis. These infections are notoriously difficult to treat, and it is believed that the ability of S. aureus to invade endothelial cells and persist intracellularly is a key mechanism for persistence despite ongoing antibiotic treatment. Here, we used dual RNA sequencing to study the simultaneous transcriptional response of S. aureus and human endothelial cells during in vitro infections. We revealed discrete and shared differentially expressed genes for both host and pathogen at the different stages of infection. While the endothelial cells upregulated genes involved in interferon signalling and antigen presentation during late infection, S. aureus downregulated toxin expression while upregulating genes related to iron scavenging. In conclusion, the presented data provide an important resource to facilitate functional investigations into host–pathogen interaction during S. aureus invasive infection and a basis for identifying novel drug target sites.

Sigma 54-Regulated Transcription Is Associated with Membrane Reorganization and Type III Secretion Effectors during Conversion to Infectious Forms of Chlamydia trachomatis

The factors that control the growth and infectious processes for Chlamydia are still poorly understood. This study used recently developed genetic tools to determine the regulon for one of the key transcription factors encoded by Chlamydia, sigma 54. Surrogate and computational analyses provide additional support for the hypothesis that sigma 54 plays a key role in controlling the expression of many components critical to converting and enabling the infectious capability of Chlamydia. These components include those that remodel the membrane for the extracellular environment and incorporation of an arsenal of type III secretion effectors in preparation for infecting new cells.

Chlamydia bacteria are obligate intracellular organisms with a phylum-defining biphasic developmental cycle that is intrinsically linked to its ability to cause disease. The progression of the chlamydial developmental cycle is regulated by the temporal expression of genes predominantly controlled by RNA polymerase sigma (σ) factors. Sigma 54 (σ⁵⁴) is one of three sigma factors encoded by Chlamydia for which the role and regulon are unknown. CtcC is part of a two-component signal transduction system that is requisite for σ⁵⁴ transcriptional activation. CtcC activation of σ⁵⁴ requires phosphorylation, which relieves inhibition by the CtcC regulatory domain and enables ATP hydrolysis by the ATPase domain. Prior studies with CtcC homologs in other organisms have shown that expression of the ATPase domain alone can activate σ⁵⁴ transcription. Biochemical analysis of CtcC ATPase domain supported the idea of ATP hydrolysis occurring in the absence of the regulatory domain, as well as the presence of an active-site residue essential for ATPase activity (E242). Using recently developed genetic approaches in Chlamydia to induce expression of the CtcC ATPase domain, a transcriptional profile was determined that is expected to reflect the σ⁵⁴ regulon. Computational evaluation revealed that the majority of the differentially expressed genes were preceded by highly conserved σ⁵⁴ promoter elements. Reporter gene analyses using these putative σ⁵⁴ promoters reinforced the accuracy of the model of the proposed regulon. Investigation of the gene products included in this regulon supports the idea that σ⁵⁴ controls expression of genes that are critical for conversion of Chlamydia from replicative reticulate bodies into infectious elementary bodies.

Spatially distinct physiology of Bacteroides fragilis within the proximal colon of gnotobiotic mice

A complex microbiota inhabits various microenvironments of the gut, with some symbiotic bacteria having evolved traits to invade the epithelial mucus layer and reside deep within intestinal tissue of animals. Whether these distinct bacterial communities across gut biogeographies exhibit divergent behaviors remains largely unknown. Global transcriptomic analysis to investigate microbial physiology in specific mucosal niches has been hampered technically by overabundance of host RNA. Herein, we employed hybrid selection RNA sequencing (hsRNA-Seq) to enable detailed spatial transcriptomic profiling of a prominent human commensal as it colonizes the colonic lumen, mucus or epithelial tissue of mice. Compared to conventional RNA-Seq, hsRNA-Seq increased reads mapping to the Bacteroides fragilis genome by 48- and 154-fold in mucus and tissue, respectively, allowing for high fidelity comparisons across biogeographic sites. Near the epithelium, B. fragilis up-regulated numerous genes involved in protein synthesis, indicating that bacteria inhabiting the mucosal niche are metabolically active. Further, a specific sulfatase (BF3086) and glycosyl hydrolase (BF3134) were highly induced in mucus and tissue compared to bacteria in the lumen. In-frame deletion of these genes impaired in vitro growth on mucus as a carbon source, as well as mucosal colonization of mice. Mutants in either B. fragilis gene displayed a fitness defect in competing for colonization against bacterial challenge, revealing the importance of site-specific gene expression for robust host-microbial symbiosis. As a versatile tool, hsRNA-Seq can be deployed to explore the in vivo spatial physiology of numerous bacterial pathogens or commensals.

New Pathogenesis Mechanisms and Translational Leads Identified by Multidimensional Analysis of Necrotizing Myositis in Primates

A fundamental goal of contemporary biomedical research is to understand the molecular basis of disease pathogenesis and exploit this information to develop targeted and more-effective therapies. Necrotizing myositis caused by the bacterial pathogen Streptococcus pyogenes is a devastating human infection with a high mortality rate and few successful therapeutic options. We used dual transcriptome sequencing (RNA-seq) to analyze the transcriptomes of S. pyogenes and host skeletal muscle recovered contemporaneously from infected nonhuman primates. The in vivo bacterial transcriptome was strikingly remodeled compared to organisms grown in vitro, with significant upregulation of genes contributing to virulence and altered regulation of metabolic genes. The transcriptome of muscle tissue from infected nonhuman primates (NHPs) differed significantly from that of mock-infected animals, due in part to substantial changes in genes contributing to inflammation and host defense processes. We discovered significant positive correlations between group A streptococcus (GAS) virulence factor transcripts and genes involved in the host immune response and inflammation. We also discovered significant correlations between the magnitude of bacterial virulence gene expression in vivo and pathogen fitness, as assessed by previously conducted genome-wide transposon-directed insertion site sequencing (TraDIS). By integrating the bacterial RNA-seq data with the fitness data generated by TraDIS, we discovered five new pathogen genes, namely, S. pyogenes 0281 (Spy0281 [dahA]), ihk-irr, slr, isp, and ciaH, that contribute to necrotizing myositis and confirmed these findings using isogenic deletion-mutant strains. Taken together, our study results provide rich new information about the molecular events occurring in severe invasive infection of primate skeletal muscle that has extensive translational research implications.IMPORTANCE Necrotizing myositis caused by Streptococcus pyogenes has high morbidity and mortality rates and relatively few successful therapeutic options. In addition, there is no licensed human S. pyogenes vaccine. To gain enhanced understanding of the molecular basis of this infection, we employed a multidimensional analysis strategy that included dual RNA-seq and other data derived from experimental infection of nonhuman primates. The data were used to target five streptococcal genes for pathogenesis research, resulting in the unambiguous demonstration that these genes contribute to pathogen-host molecular interactions in necrotizing infections. We exploited fitness data derived from a recently conducted genome-wide transposon mutagenesis study to discover significant correlation between the magnitude of bacterial virulence gene expression in vivo and pathogen fitness. Collectively, our findings have significant implications for translational research, potentially including vaccine efforts.

A survey of Fusobacterium nucleatum genes modulated by host cell infection

Here, we report comprehensive transcriptomic profiles from Fusobacterium nucleatum under conditions that mimic the first stages of bacterial infection in a highly differentiated adenocarcinoma epithelial cell line. Our transcriptomic in vitro adenocarcinoma approach allows us to measure the expression dynamics and regulation of bacterial virulence and response factors in real time, and is a novel strategy for clarifying the role of F. nucleatum infection in colorectal cancer (CRC) progression. Our data show that: (i) infection alters metabolic and functional pathways in F. nucleatum, allowing the bacterium to adapt to the host-imposed milieu; (ii) infection also stimulates the expression of genes required to help induce and promote a hypoxic and inflammatory microenvironment in the host; and (iii) F. nucleatum invasion occurs by a haematogenous route of infection. Our study identifies novel gene targets from F. nucleatum that are activated during invasion and which may aid in determining how this species invades and promotes disease within the human gastrointestinal tract. These invasion-specific genes may be useful as biomarkers for CRC progression in a host and could also assist in the development of new diagnostic tools and treatments (such as vaccines or small molecule drug targets), which will be able to combat infection and inflammation in the host while circumventing the potential problem of F. nucleatum tolerization.

Transcriptional Profiling of Host Cell Responses to Virulent Haemophilus parasuis: New Insights into Pathogenesis

Haemophilus parasuis is the causative agent of Glässer’s disease in pigs. H. parasuis can cause vascular damage, although the mechanism remains unclear. In this study, we investigated the host cell responses involved in the molecular pathway interactions in porcine aortic vascular endothelial cells (PAVECs) induced by H. parasuis using RNA-Seq. The transcriptome results showed that when PAVECs were infected with H. parasuis for 24 h, 281 differentially expressed genes (DEGs) were identified; of which, 236 were upregulated and 45 downregulated. The 281 DEGs were involved in 136 KEGG signaling pathways that were organismal systems, environmental information processing, metabolism, cellular processes, and genetic information processing. The main pathways were the Rap1, FoxO, and PI3K/Akt signaling pathways, and the overexpressed genes were determined and verified by quantitative reverse transcription polymerase chain reaction. In addition, 252 genes were clustered into biological processes, molecular processes, and cellular components. Our study provides new insights for understanding the interaction between bacterial and host cells, and analyzed, in detail, the possible mechanisms that lead to vascular damage induced by H. parasuis. This may lead to development of novel therapeutic targets to control H. parasuis infection.