Bioinformatic analysis of bacteria and host cell dual RNA-sequencing experiments

The Fur-like regulatory protein MAP3773c modulates key metabolic pathways in Mycobacterium avium subsp. paratuberculosis under in-vitro iron starvation

Johne's disease (JD) is a chronic enteric infection of dairy cattle worldwide. Mycobacterium avium subsp. paratuberculosis (MAP), the causative agent of JD, is fastidious often requiring eight to sixteen weeks to produce colonies in culture-a major hurdle in the diagnosis and therefore in implementation of optimal JD control measures. A significant gap in knowledge is the comprehensive understanding of the metabolic networks deployed by MAP to regulate iron both in-vitro and in-vivo. The genome of MAP carries MAP3773c, a putative metal regulator, which is absent in all other mycobacteria. The role of MAP3773c in intracellular iron regulation is poorly understood. In the current study, a field isolate (K-10) and an in-frame MAP3773c deletion mutant (ΔMAP3773c) derived from K-10, were exposed to iron starvation for 5, 30, 60, and 90 min and RNA-Seq was performed. A comparison of transcriptional profiles between K-10 and ΔMAP3773c showed 425 differentially expressed genes (DEGs) at 30 min time post-iron restriction. Functional analysis of DEGs in ΔMAP3773c revealed that pantothenate (Pan) biosynthesis, polysaccharide biosynthesis and sugar metabolism genes were downregulated at 30 min post-iron starvation whereas ATP-binding cassette (ABC) type metal transporters, putative siderophore biosynthesis, PPE and PE family genes were upregulated. Pathway analysis revealed that the MAP3773c knockout has an impairment in Pan and Coenzyme A (CoA) biosynthesis pathways suggesting that the absence of those pathways likely affect overall metabolic processes and cellular functions, which have consequences on MAP survival and pathogenesis.

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.

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.

Dual transcriptomic analysis reveals metabolic changes associated with differential persistence of human pathogenic bacteria in leaves of Arabidopsis and lettuce

Understanding the molecular determinants underlying the interaction between the leaf and human pathogenic bacteria is key to provide the foundation to develop science-based strategies to prevent or decrease the pathogen contamination of leafy greens. In this study, we conducted a dual RNA-sequencing analysis to simultaneously define changes in the transcriptomic profiles of the plant and the bacterium when they come in contact. We used an economically relevant vegetable crop, lettuce (Lactuca sativa L. cultivar Salinas), and a model plant, Arabidopsis thaliana Col-0, as well as two pathogenic bacterial strains that cause disease outbreaks associated with fresh produce, Escherichia coli O157:H7 and Salmonella enterica serovar Typhimurium 14028s (STm 14028s). We observed commonalities and specificities in the modulation of biological processes between Arabidopsis and lettuce and between O157:H7 and STm 14028s during early stages of the interaction. We detected a larger alteration of gene expression at the whole transcriptome level in lettuce and Arabidopsis at 24 h post inoculation with STm 14028s compared to that with O157:H7. In addition, bacterial transcriptomic adjustments were substantially larger in Arabidopsis than in lettuce. Bacterial transcriptome was affected at a larger extent in the first 4 h compared to the subsequent 20 h after inoculation. Overall, we gained valuable knowledge about the responses and counter-responses of both bacterial pathogen and plant host when these bacteria are residing in the leaf intercellular space. These findings and the public genomic resources generated in this study are valuable for additional data mining.

Molecular insights into the mechanisms of susceptibility of Labeo rohita against oomycete Aphanomyces invadans

Aphanomyces invadans, the causative agent of epizootic ulcerative syndrome, is one of the most destructive pathogens of freshwater fishes. To date, the disease has been reported from over 160 fish species in 20 countries and notably, this is the first non-salmonid disease that has resulted in major impacts globally. In particular, Indian major carps (IMCs) are highly susceptible to this disease. To increase our knowledge particularly with regards to host immune response against A. invadans infection in a susceptible host, the gene expression profile in head kidney of A. invadans-infected and control rohu, Labeo rohita was investigated using RNA sequencing. Time course analysis of RNA-Seq data revealed 5608 differentially expressed genes, involved among others in Antigen processing and presentation, Leukocyte transendothelial migration, IL-17 signaling, Chemokine signaling, C-type lectin receptor signaling and Toll-like receptor signaling pathways. In the affected pathways, a number of immune genes were found to be downregulated, suggesting an immune evasion strategy of A. invadans in establishing the infection. The information generated in this study offers first systematic mechanistic understanding of the host–pathogen interaction that might underpin the development of new management strategies for this economically devastating fish-pathogenic oomycete A. invadans.

Early Transcriptional Landscapes of Chlamydia trachomatis-Infected Epithelial Cells at Single Cell Resolution

Chlamydia are Gram-negative obligate intracellular bacterial pathogens responsible for a variety of disease in humans and animals worldwide. Chlamydia trachomatis causes trachoma in disadvantaged populations, and is the most common bacterial sexually transmitted infection in humans, causing reproductive tract disease. Antibiotic therapy successfully treats diagnosed chlamydial infections, however asymptomatic infections are common. High-throughput transcriptomic approaches have explored chlamydial gene expression and infected host cell gene expression. However, these were performed on large cell populations, averaging gene expression profiles across all cells sampled and potentially obscuring biologically relevant subsets of cells. We generated a pilot dataset, applying single cell RNA-Seq (scRNA-Seq) to C. trachomatis infected and mock-infected epithelial cells to assess the utility, pitfalls and challenges of single cell approaches applied to chlamydial biology, and to potentially identify early host cell biomarkers of chlamydial infection. Two hundred sixty-four time-matched C. trachomatis-infected and mock-infected HEp-2 cells were collected and subjected to scRNA-Seq. After quality control, 200 cells were retained for analysis. Two distinct clusters distinguished 3-h cells from 6- and 12-h. Pseudotime analysis identified a possible infection-specific cellular trajectory for Chlamydia-infected cells, while differential expression analyses found temporal expression of metallothioneins and genes involved with cell cycle regulation, innate immune responses, cytoskeletal components, lipid biosynthesis and cellular stress. We find that changes to the host cell transcriptome at early times of C. trachomatis infection are readily discernible by scRNA-Seq, supporting the utility of single cell approaches to identify host cell biomarkers of chlamydial infection, and to further deconvolute the complex host response to infection.

Stress responses, M2 macrophages, and a distinct microbial signature in fatal intestinal acute graft-versus-host disease

Steroid-refractory intestinal acute graft-versus-host disease (aGVHD) is a frequently fatal condition with little known about mechanisms driving failed steroid responses in gut mucosa. To uncover novel molecular insights in steroid-refractory aGVHD, we compared gene expression profiles of rectosigmoid biopsies from patients at diagnosis of clinical stage 3-4 lower intestinal aGVHD (N=22), to repeat biopsies when the patients became steroid refractory (N=22), and normal controls (N=10). We also performed single gene analyses of factors associated with tolerance (programmed death ligand-1 [PDL1], indoleamine 2,3 dioxygenase [IDO1], and T cell immunoreceptor with Ig and ITIM domains [TIGIT]) and found that significantly higher expression levels of these aGVHD inhibitory genes (PDL1, IDO1, TIGIT) at aGVHD onset became decreased in the steroid-refractory state. We examined genes triggered by microbial ligands to stimulate gut repair, amphiregulin (AREG) and the aryl hydrocarbon receptor (AhR), and found that both AREG and AhR gene expression levels were increased at aGVHD onset and remained elevated in steroid-refractory aGVHD. We also identified higher expression levels of metallothioneines, metal-binding enzymes induced in stress responses, and M2 macrophage genes in steroid-refractory aGVHD. We observed no differences in T-cell subsets between onset and steroid-refractory aGVHD. Patients with a rapidly fatal course showed greater DNA damage and a distinct microbial signature at aGVHD onset, whereas patients with more prolonged survival exhibited a gene expression profile consistent with activation of Smoothened. Our results extend the paradigm beyond T cell-centric therapies for steroid-refractory GI aGVHD and highlight new mechanisms for therapeutic exploration.

Dual RNA-Seq of Chlamydia and Host Cells

During the infection of a host cell by a bacterial pathogen, a cascading series of gene expression changes occurs as each organism manipulates or responds to the other via defense or survival strategies. Unraveling this complex interplay is key for our understanding of bacterial virulence and host response pathways for the development of novel therapeutics. Dual RNA sequencing (dual RNA-Seq) has recently been developed to simultaneously capture host and bacterial transcriptomes from an infected cell. Leveraging the sensitivity and resolution allowed by RNA-seq, dual RNA-Seq can be applied to any bacteria-eukaryotic host interaction. We pioneered dual RNA-Seq to simultaneously capture Chlamydia and host expression profiles during an in vitro infection as proof of principle. Here we provide a detailed laboratory protocol and bioinformatics analysis guidelines for dual RNA-seq experiments focusing on Chlamydia as the organism of interest.

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

Dual RNA-Sequencing leverages established next-generation sequencing (NGS)-enabled RNA-Seq approaches to measure genome-wide transcriptional changes of both an infecting bacteria and host cells. By simultaneously investigating both organisms from the same biological sample, dual RNA-Seq can provide unique insight into bacterial infection processes and reciprocal host responses at once. However, the difficulties involved in handling both prokaryotic and eukaryotic material require distinct, optimized procedures. We previously developed and applied dual RNA-Seq to measure prokaryotic and eukaryotic expression profiles of human cells infected with bacteria, using in vitro Chlamydia-infected epithelial cells as proof of principle. Here we provide a detailed laboratory protocol for in vitro dual RNA-Seq that is readily adaptable to any host-bacteria system of interest.