Dual RNA-seq of Nontypeable Haemophilus influenzae and Host Cell Transcriptomes Reveals Novel Insights into Host-Pathogen Cross Talk

Tolerance to Haemophilus influenzae infection in human epithelial cells: Insights from a primary cell-based model

Haemophilus influenzae is a human respiratory pathogen and inhabits the human respiratory tract as its only niche. Despite this, the molecular mechanisms that allow H. influenzae to establish persistent infections of human epithelia are not well understood. Here, we have investigated how H. influenzae adapts to the host environment and triggers the host immune response using a human primary cell-based infection model that closely resembles human nasal epithelia (NHNE). Physiological assays combined with dualRNAseq revealed that NHNE from five healthy donors all responded to H. influenzae infection with an initial, 'unproductive' inflammatory response that included a strong hypoxia signature but did not produce pro-inflammatory cytokines. Subsequently, an apparent tolerance to large extracellular and intraepithelial burdens of H. influenzae developed, with NHNE transcriptional profiles resembling the pre-infection state. This occurred in parallel with the development of intraepithelial bacterial populations, and appears to involve interruption of NFκB signalling. This is the first time that large-scale, persistence-promoting immunomodulatory effects of H. influenzae during infection have been observed, and we were able to demonstrate that only infections with live, but not heat-killed H. influenzae led to immunomodulation and reduced expression of NFκB-controlled cytokines such as IL-1β, IL-36γ and TNFα. Interestingly, NHNE were able to re-activate pro-inflammatory responses towards the end of the 14-day infection, resulting in release of IL-8 and TNFα. In addition to providing first molecular insights into mechanisms enabling persistence of H. influenzae in the host, our data further indicate the presence of infection stage-specific gene expression modules, highlighting fundamental similarities between immune responses in NHNE and canonical immune cells, which merit further investigation.

Geogenic particles induce bronchial susceptibility to non-typeable Haemophilus influenzae

Exposure to geogenic (earth-derived) particulate matter (PM) is linked to an increased prevalence of bronchiectasis and other respiratory infections in Australian Indigenous communities. Experimental studies have shown that the concentration of iron in geogenic PM is associated with the magnitude of respiratory health effects, however, the mechanism is unclear. We investigated the effect of geogenic PM and iron oxide on the invasiveness of non-typeable Haemophilus influenzae (NTHi). Peripheral blood mononuclear cell-derived macrophages or epithelial cell lines (A549 & BEAS-2B) were exposed to whole geogenic PM, their primary constituents (haematite, magnetite or silica) or diesel exhaust particles (DEP). The uptake of bacteria was quantified by flow cytometry and whole genome sequencing (WGS) was performed on NTHi strains. Geogenic PM increased the invasiveness of NTHi in bronchial epithelial cells. Of the primary constituents, haematite also increased NTHi invasion with magnetite and silica having significantly less impact. Furthermore, we observed varying levels of invasiveness amongst NTHi isolates. WGS analysis suggested isolates with more genes associated with heme acquisition were more virulent in BEAS-2B cells. The present study suggests that geogenic particles can increase the susceptibility of bronchial epithelial cells to select bacterial pathogens in vitro, a response primarily driven by haematite content in the dust. This demonstrates a potential mechanism linking exposure to iron-laden geogenic PM and high rates of chronic respiratory infections in remote communities in arid environments.

Dual RNA-seq identifies genes and pathways modulated during Clostridioides difficile colonization

IMPORTANCE

The initial interactions between the colonic epithelium and the bacterium are likely critical in the establishment of Clostridioides difficile infection, one of the major causes of hospital-acquired diarrhea worldwide. Molecular interactions between C. difficile and human gut cells have not been well defined mainly due to the technical challenges of studying cellular host-pathogen interactions with this anaerobe. Here we have examined transcriptional changes occurring in the pathogen and host cells during the initial 24 hours of infection. Our data indicate several changes in metabolic pathways and virulence-associated factors during the initial bacterium-host cell contact and early stages of infection. We describe canonical pathways enriched based on the expression profiles of a dual RNA sequencing in the host and bacterium, and functions of bacterial factors that are modulated during infection. This study thus provides fresh insight into the early C. difficile infection process.

In vivo gene expression profile of Haemophilus influenzae during human pneumonia

Haemophilus influenzae is a major cause of community-acquired pneumonia. While studied extensively in various laboratory models, less is known about the cell function while inside the human lung. We present the first analysis of the global gene expression of H. influenzae while the bacteria are in the lung during pneumonia (in vivo conditions) and contrast it with bacterial isolates that have been cultured under standard laboratory conditions (in vitro conditions). Patients with pneumonia were recruited from emergency departments and intensive care units during 2018-2020 (n = 102). Lower respiratory samples were collected for bacterial culture and RNA extraction. Patient samples with H. influenzae (n = 8) and colonies from bacterial cultures (n = 6) underwent RNA sequencing. The reads were then pseudo-aligned to core and pan genomes created from 15 reference strains. While bacteria cultured in vitro clustered tightly by principal component analysis of core genome (n = 1067) gene expression, bacteria in the patient samples had more diverse transcriptomic signatures and did not group with their lab-cultured counterparts. In total, 328 core genes were significantly differentially expressed between in vitro and in vivo conditions. The most highly upregulated genes in vivo included tbpA and fbpA, which are involved in the acquisition of iron from transferrin, and the stress response gene msrAB. The biosynthesis of nucleotides/purines and molybdopterin-scavenging processes were also significantly enriched in vivo. In contrast, major metabolic pathways and iron-sequestering genes were downregulated under this condition. In conclusion, extensive transcriptomic differences were found between bacteria while in the human lung and bacteria that were cultured in vitro. IMPORTANCE The human-specific pathogen Haemophilus influenzae is generally not well suited for studying in animal models, and most laboratory models are unlikely to approximate the diverse environments encountered by bacteria in the human airways accurately. Thus, we have examined the global gene expression of H. influenzae during pneumonia. Extensive differences in the global gene expression profiles were found in H. influenzae while in the human lung compared to bacteria that were grown in the laboratory. In contrast, the gene expression profiles of isolates collected from different patients were found to cluster together when grown under the same laboratory conditions. Interesting observations were made of how H. influenzae acquires and uses iron and molybdate, endures oxidative stress, and regulates central metabolism while in the lung. Our results indicate important processes during infection and can guide future research on genes and pathways that are relevant in the pathogenesis of H. influenzae pneumonia.

Multiplatform Metabolomics Characterization Reveals Novel Metabolites and Phospholipid Compositional Rules of Haemophilus influenzae Rd KW20

Haemophilus influenzae is a gram-negative bacterium of relevant clinical interest. H. influenzae Rd KW20 was the first organism to be sequenced and for which a genome-scale metabolic model (GEM) was developed. However, current H. influenzae GEMs are unable to capture several aspects of metabolome nature related to metabolite pools. To directly and comprehensively characterize the endometabolome of H. influenzae Rd KW20, we performed a multiplatform MS-based metabolomics approach combining LC-MS, GC-MS and CE-MS. We obtained direct evidence of 15-20% of the endometabolome present in current H. influenzae GEMs and showed that polar metabolite pools are interconnected through correlating metabolite islands. Notably, we obtained high-quality evidence of 18 metabolites not previously included in H. influenzae GEMs, including the antimicrobial metabolite cyclo(Leu-Pro). Additionally, we comprehensively characterized and evaluated the quantitative composition of the phospholipidome of H. influenzae, revealing that the fatty acyl chain composition is largely independent of the lipid class, as well as that the probability distribution of phospholipids is mostly related to the conditional probability distribution of individual acyl chains. This finding enabled us to provide a rationale for the observed phospholipid profiles and estimate the abundance of low-level species, permitting the expansion of the phospholipidome characterization through predictive probabilistic modelling.

In Vivo Genome-Wide Gene Expression Profiling Reveals That Haemophilus influenzae Purine Synthesis Pathway Benefits Its Infectivity within the Airways

Haemophilus influenzae is a human-adapted bacterial pathogen that causes airway infections. Bacterial and host elements associated with the fitness of H. influenzae within the host lung are not well understood. Here, we exploited the strength of in vivo-omic analyses to study host-microbe interactions during infection. We used in vivo transcriptome sequencing (RNA-seq) for genome-wide profiling of both host and bacterial gene expression during mouse lung infection. Profiling of murine lung gene expression upon infection showed upregulation of lung inflammatory response and ribosomal organization genes, and downregulation of cell adhesion and cytoskeleton genes. Transcriptomic analysis of bacteria recovered from bronchoalveolar lavage fluid samples from infected mice showed a significant metabolic rewiring during infection, which was highly different from that obtained upon bacterial in vitro growth in an artificial sputum medium suitable for H. influenzae. In vivo RNA-seq revealed upregulation of bacterial de novo purine biosynthesis, genes involved in non-aromatic amino acid biosynthesis, and part of the natural competence machinery. In contrast, the expression of genes involved in fatty acid and cell wall synthesis and lipooligosaccharide decoration was downregulated. Correlations between upregulated gene expression and mutant attenuation in vivo were established, as observed upon purH gene inactivation leading to purine auxotrophy. Likewise, the purine analogs 6-thioguanine and 6-mercaptopurine reduced H. influenzae viability in a dose-dependent manner. These data expand our understanding of H. influenzae requirements during infection. In particular, H. influenzae exploits purine nucleotide synthesis as a fitness determinant, raising the possibility of purine synthesis as an anti-H. influenzae target. IMPORTANCE In vivo-omic strategies offer great opportunities for increased understanding of host-pathogen interplay and for identification of therapeutic targets. Here, using transcriptome sequencing, we profiled host and pathogen gene expression during H. influenzae infection within the murine airways. Lung pro-inflammatory gene expression reprogramming was observed. Moreover, we uncovered bacterial metabolic requirements during infection. In particular, we identified purine synthesis as a key player, highlighting that H. influenzae may face restrictions in purine nucleotide availability within the host airways. Therefore, blocking this biosynthetic process may have therapeutic potential, as supported by the observed inhibitory effect of 6-thioguanine and 6-mercaptopurine on H. influenzae growth. Together, we present key outcomes and challenges for implementing in vivo-omics in bacterial airway pathogenesis. Our findings provide metabolic insights into H. influenzae infection biology, raising the possibility of purine synthesis as an anti-H. influenzae target and of purine analog repurposing as an antimicrobial strategy against this pathogen.

Identification and initial characterization of Hfq-associated sRNAs in Histophilus somni strain 2336

Small RNAs (sRNA), in association with the global chaperone regulator Hfq, positively or negatively regulate gene expression in bacteria. For this study, Histophilus somni sRNAs that bind to Hfq were identified and then partially characterized. The Hfq-associated sRNAs in H. somni were isolated and identified by co-immunoprecipitation using anti-Hfq antibody, followed by sRNA sequencing. Sequence analysis of the sRNA samples identified 100 putative sRNAs, out of which 16 were present in pathogenic strain 2336, but not in non-pathogenic strain 129Pt. Bioinformatic analyses suggested that the sRNAs HS9, HS79, and HS97 could bind to many genes putatively involved in virulence/biofilm formation. Furthermore, multi-sequence alignment of the sRNA regions in the genome revealed that HS9 and HS97 could interact with sigma 54, which is a transcription factor linked to important bacterial traits, including motility, virulence, and biofilm formation. Northern blotting was used to determine the approximate size, abundance and any processing events attributed to the sRNAs. Selected sRNA candidates were confirmed to bind Hfq, as determined by electrophoretic mobility shift assays using sRNAs synthesized by in vitro transcription and recombinant Hfq. The exact transcriptional start site of the sRNA candidates was determined by RNA ligase-mediated rapid amplification of cDNA ends, followed by cloning and sequencing. This is the first investigation of H. somni sRNAs that show they may have important regulatory roles in virulence and biofilm formation.

Pathogenesis of nontypeable Haemophilus influenzae infections in chronic suppurative lung disease

The respiratory tract antimicrobial defense system is a multilayered defense mechanism that relies upon mucociliary clearance and components of both the innate and adaptive immune systems to protect the lungs from inhaled or aspirated microorganisms. One of these potential pathogens, nontypeable Haemophilus influenzae (NTHi), adopts several, multifaceted redundant strategies to successfully colonize the lower airways and establish a persistent infection. NTHi can impair mucociliary clearance, express multiple multifunctional adhesins for various cell types within the respiratory tract and evade host defenses by surviving within and between cells, forming biofilms, increasing antigenic drift, secreting proteases and antioxidants, and by host-pathogen cross-talk, impair macrophage and neutrophil function. NTHi is recognized as an important pathogen in several chronic lower respiratory disorders, such as protracted bacterial bronchitis, bronchiectasis, cystic fibrosis, and primary ciliary dyskinesia. The persistence of NTHi in human airways, including its capacity to form biofilms, results in chronic infection and inflammation, which can ultimately injure airway wall structures. The complex nature of the molecular pathogenetic mechanisms employed by NTHi is incompletely understood but improved understanding of its pathobiology will be important for developing effective therapies and vaccines, especially given the marked genetic heterogeneity of NTHi and its possession of phase-variable genes. Currently, no vaccine candidates are ready for large phase III clinical trials.

Epithelial immune activation and intracellular invasion by non-typeable Haemophilus influenzae

Type-2 low asthma affects 30-50% of people with severe asthma and includes a phenotype characterized by sputum neutrophilia and resistance to corticosteroids. Airways inflammation in type-2 low asthma or COPD is potentially driven by persistent bacterial colonization of the lower airways by bacteria such as non-encapsulated Haemophilus influenzae (NTHi). Although pathogenic in the lower airways, NTHi is a commensal of the upper airways. It is not known to what extent these strains can invade airway epithelial cells, persist intracellularly and activate epithelial cell production of proinflammatory cytokines, and how this differs between the upper and lower airways. We studied NTHi infection of primary human bronchial epithelial cells (PBECs), primary nasal epithelial cells (NECs) and epithelial cell lines from upper and lower airways. NTHi strains differed in propensity for intracellular and paracellular invasion. We found NTHi was internalized within PBECs at 6 h, but live intracellular infection did not persist at 24 h. Confocal microscopy and flow cytometry showed NTHi infected secretory, ciliated and basal PBECs. Infection of PBECs led to induction of CXCL8, interleukin (IL)-1β, IL-6 and TNF. The magnitude of cytokine induction was independent of the degree of intracellular invasion, either by differing strains or by cytochalasin D inhibition of endocytosis, with the exception of the inflammasome-induced mediator IL-1β. NTHi-induced activation of TLR2/4, NOD1/2 and NLR inflammasome pathways was significantly stronger in NECs than in PBECs. These data suggest that NTHi is internalized transiently by airway epithelial cells and has capacity to drive inflammation in airway epithelial cells.

Salmonella Typhimurium induces genome-wide expression and phosphorylation changes that modulate immune response, intracellular survival and vesicle transport in infected neutrophils

Salmonella Typhimurium is a food-borne pathogen that causes salmonellosis. When in contact with the host, neutrophils are rapidly recruited to act as first line of defense. To better understand the pathogenesis of this infection, we used an in vitro model of neutrophil infection to perform dual RNA-sequencing (both host and pathogen). In addition, and given that many pathogens interfere with kinase-mediated phosphorylation in host signaling, we performed a phosphoproteomic analysis. The immune response was overall diminished in infected neutrophils, mainly JAK/STAT and toll-like receptor signaling pathways. We found decreased expression of proinflammatory cytokine receptor genes and predicted downregulation of the mitogen-activated protein (MAPK) signaling pathway. Also, Salmonella infection inhibited interferons I and II signaling pathways by upregulation of SOCS3 and subsequent downregulation of STAT1 and STAT2. Additionally, phosphorylation of PSMC2 and PSMC4, proteasome regulatory proteins, was decreased in infected neutrophils. Cell viability and survival was increased by p53 signaling, cell cycle arrest and NFkB-proteasome pathways activation. Combined analysis of RNA-seq and phosphoproteomics also revealed inhibited vesicle transport mechanisms mediated by dynein/dynactin and exocyst complexes, involved in ER-to-Golgi transport and centripetal movement of lysosomes and endosomes. Among the overexpressed virulence genes from Salmonella we found potential effectors responsible of these dysregulations, such as spiC, sopD2, sifA or pipB2, all of them involved in intracellular replication. Our results suggest that Salmonella induces (through overexpression of virulence factors) transcriptional and phosphorylation changes that increases neutrophil survival and shuts down immune response to minimize host response, and impairing intracellular vesicle transport likely to keep nutrients for replication and Salmonella-containing vacuole formation and maintenance.

Dual RNA sequencing of group B Streptococcus-infected human monocytes reveals new insights into host-pathogen interactions and bacterial evasion of phagocytosis

Streptococcus agalactiae, also known as Group B Streptococcus (GBS) is a frequent cause of infections, including bacteraemia and other acute diseases in adults and immunocompromised individuals. We developed a novel system to study GBS within human monocytes to define the co-transcriptome of intracellular GBS (iGBS) and host cells simultaneously using dual RNA-sequencing (RNA-seq) to better define how this pathogen responds to host cells. Using human U937 monocytes and genome-sequenced GBS reference strain 874,391 in antibiotic protection assays we validated a system for dual-RNA seq based on measures of GBS and monocyte viability to ensure that the bacterial and host cell co-transcriptome reflected mainly intracellular (iGBS) rather than extracellular GBS. Elucidation of the co-transcriptome revealed 1119 dysregulated transcripts in iGBS with most genes, including several that encode virulence factors (e.g., scpB, hvgA, ribD, pil2b) exhibiting activation by upregulated expression. Infection with iGBS resulted in significant remodelling of the monocyte transcriptome, with 7587 transcripts differentially expressed including 7040 up-regulated and 547 down-regulated. qPCR confirmed that the most strongly activated genes included sht, encoding Streptococcal Histidine Triad Protein. An isogenic GBS mutant strain deficient in sht revealed a significant effect of this gene on phagocytosis of GBS and survival of the bacteria during systemic infection in mice. Identification of a novel contribution of sht to GBS virulence shows the co-transcriptome responses elucidated in GBS-infected monocytes help to shape the host-pathogen interaction and establish a role for sht in the response of the bacteria to phagocytic uptake. This study provides comprehension of concurrent transcriptional responses that occur in GBS and human monocytes that shape the host-pathogen interaction.

Pneumolysin promotes host cell necroptosis and bacterial competence during pneumococcal meningitis as shown by whole-animal dual RNA-seq

Pneumolysin is a major virulence factor of Streptococcus pneumoniae that plays a key role in interaction with the host during invasive disease. How pneumolysin influences these dynamics between host and pathogen interaction during early phase of central nervous system infection in pneumococcal meningitis remains unclear. Using a whole-animal in vivo dual RNA sequencing (RNA-seq) approach, we identify pneumolysin-specific transcriptional responses in both S. pneumoniae and zebrafish (Danio rerio) during early pneumococcal meningitis. By functional enrichment analysis, we identify host pathways known to be activated by pneumolysin and discover the importance of necroptosis for host survival. Inhibition of this pathway using the drug GSK'872 increases host mortality during pneumococcal meningitis. On the pathogen's side, we show that pneumolysin-dependent competence activation is crucial for intra-host replication and virulence. Altogether, this study provides new insights into pneumolysin-specific transcriptional responses and identifies key pathways involved in pneumococcal meningitis.

Molecular, cellular and neurological consequences of infection by the neglected human pathogen Nocardia

Background

Nocardia is a facultative intracellular pathogen that infects the lungs and brains of immunocompromised patients with consequences that can be fatal. The incidence of such infections is rising, immunocompetent individuals are also being infected, and there is a need to learn more about this neglected bacterial pathogen and the interaction with its human host.

Results

We have applied dual RNA-seq to assess the global transcriptome changes that occur simultaneously in Nocardia farcinica (N. farcinica) and infected human epithelial alveolar host cells, and have tested a series of mutants in this in vitro system to identify candidate determinants of virulence. Using a mouse model, we revealed the profiles of inflammation-related factors in the lung after intranasal infection and confirmed that nbtB and nbtS are key virulence genes for Nocardia infection in vivo. Regarding the host response to infection, we found that the expression of many histones was dysregulated during the infection of lung cells, indicating that epigenetic modification might play a crucial role in the host during Nocardia infection. In our mouse model, Nocardia infection led to neurological symptoms and we found that 15 of 22 Nocardia clinical strains tested could cause obvious PD-like symptoms. Further experiments indicated that Nocardia infection could activate microglia and drive M1 microglial polarization, promote iNOS and CXCL-10 production, and cause neuroinflammation in the substantia nigra, all of which may be involved in causing PD-like symptoms. Importantly, the deletion of nbtS in N. farcinica completely attenuated the neurological symptoms.

Conclusions

Our data contribute to an in-depth understanding of the characteristics of both the host and Nocardia during infection and provide valuable clues for future studies of this neglected human pathogen, especially those addressing the underlying causes of infection-related neurological symptoms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01452-7.

Nontypeable Haemophilus influenzae infection of pulmonary macrophages drives neutrophilic inflammation in severe asthma

BACKGROUND

Nontypeable Haemophilus influenzae (NTHi) is a respiratory tract pathobiont that chronically colonizes the airways of asthma patients and is associated with severe, neutrophilic disease phenotypes. The mechanism of NTHi airway persistence is not well understood, but accumulating evidence suggests NTHi can persist within host airway immune cells such as macrophages. We hypothesized that NTHi infection of pulmonary macrophages drives neutrophilic inflammation in severe asthma.

METHODS

Bronchoalveolar lavage (BAL) samples from 25 severe asthma patients were assessed by fluorescence in situ hybridisation to quantify NTHi presence. Weighted gene correlation network analysis (WGCNA) was performed on RNASeq data from NTHi-infected monocyte-derived macrophages to identify transcriptomic networks associated with NTHi infection.

RESULTS

NTHi was detected in 56% of BAL samples (NTHi+) and was associated with longer asthma duration (34 vs 22.5 years, p = .0436) and higher sputum neutrophil proportion (67% vs 25%, p = .0462). WGCNA identified a transcriptomic network of immune-related macrophage genes significantly associated with NTHi infection, including upregulation of T17 inflammatory mediators and neutrophil chemoattractants IL1B, IL8, IL23 and CCL20 (all p < .05). Macrophage network genes SGPP2 (p = .0221), IL1B (p = .0014) and GBP1 (p = .0477) were more highly expressed in NTHi+ BAL and moderately correlated with asthma duration (IL1B; rho = 0.41, p = .041) and lower prebronchodilator FEV1/FVC% (GBP1; rho = -0.43, p = .046 and IL1B; rho = -0.42, p = .055).

CONCLUSIONS

NTHi persistence with pulmonary macrophages may contribute to chronic airway inflammation and T17 responses in severe asthma, which can lead to decreased lung function and reduced steroid responsiveness. Identifying therapeutic strategies to reduce the burden of NTHi in asthma could improve patient outcomes.

Ribosomal RNA-Depletion Provides an Efficient Method for Successful Dual RNA-Seq Expression Profiling of a Marine Sponge Holobiont

Investigations of host-symbiont interactions can benefit enormously from a complete and reliable holobiont gene expression profiling. The most efficient way to acquire holobiont transcriptomes is to perform RNA-Seq on both host and symbionts simultaneously. However, optimal methods for capturing both host and symbiont mRNAs are still under development, particularly when the host is a eukaryote and the symbionts are bacteria or archaea. Traditionally, poly(A)-enriched libraries have been used to capture eukaryotic mRNA, but the ability of this method to adequately capture bacterial mRNAs is unclear because of the short half-life of the bacterial transcripts. Here, we address this gap in knowledge with the aim of helping others to choose an appropriate RNA-Seq approach for analysis of animal host-bacterial symbiont transcriptomes. Specifically, we compared transcriptome bias, depth and coverage achieved by two different mRNA capture and sequencing strategies applied to the marine demosponge Amphimedon queenslandica holobiont. Annotated genomes of the sponge host and the three most abundant bacterial symbionts, which can comprise up to 95% of the adult microbiome, are available. Importantly, this allows for transcriptomes to be accurately mapped to these genomes, and thus quantitatively assessed and compared. The two strategies that we compare here are (i) poly(A) captured mRNA-Seq (Poly(A)-RNA-Seq) and (ii) ribosomal RNA depleted RNA-Seq (rRNA-depleted-RNA-Seq). For the host sponge, we find no significant difference in transcriptomes generated by the two different mRNA capture methods. However, for the symbiont transcriptomes, we confirm the expectation that the rRNA-depleted-RNA-Seq performs much better than the Poly(A)-RNA-Seq. This comparison demonstrates that RNA-Seq by ribosomal RNA depletion is an effective and reliable method to simultaneously capture gene expression in host and symbionts and thus to analyse holobiont transcriptomes.

Dual RNA-Seq of H5N1 Avian Influenza Virus and Host Cell Transcriptomes Reveals Novel Insights Into Host-Pathogen Cross Talk

H5N1 avian influenza virus (AIV) is a highly pathogenic influenza virus that poses a substantial threat to poultry production and public health. A comprehensive understanding of host-pathogen interactions for AIV requires knowledge of gene expression changes in both the pathogen and the host upon infection. We report the use of dual RNA sequencing technology to uncover trends in gene expression in H5N1 AIV and chickens (DF1 cells) during the course of infection. The expression of all viral genes increased continuously from 0 to 20 h post infection. We also identified 2,762 differentially expressed host genes during infection. Pathway analysis found that genes related to the signaling pathways of DNA replication, T cell activation, NF-kappa B signaling pathway, and RNA degradation were significantly enriched. We demonstrated that the cis-acting lncRNA MSTRG.14019.1 targeted CSE1L and may affect virus replication. This study provides a more comprehensive and detailed understanding of host-virus interactions at the RNA level during the course of H5N1 AIV infection.

Transcriptomic Characterization of Tuberculous Sputum Reveals a Host Warburg Effect and Microbial Cholesterol Catabolism

The crucial transmission phase of tuberculosis (TB) relies on infectious sputum and yet cannot easily be modeled. We applied one-step RNA sequencing (RNA-Seq) to sputum from infectious TB patients to investigate the host and microbial environments underlying transmission of Mycobacterium tuberculosis. In such TB sputa, compared to non-TB controls, transcriptional upregulation of inflammatory responses, including an interferon-driven proinflammatory response and a metabolic shift toward glycolysis, was observed in the host. Among all bacterial sequences in the sputum, approximately 1.5% originated from M. tuberculosis, and its transcript abundance was lower in HIV-1-coinfected patients. Commensal bacterial abundance was reduced in the presence of M. tuberculosis infection. Direct alignment to the genomes of the predominant microbiota species also reveals differential adaptation, whereby firmicutes (e.g., streptococci) displayed a nonreplicating phenotype with reduced transcription of ribosomal proteins and reduced activities of ATP synthases, while Neisseria and Prevotella spp. were less affected. The transcriptome of sputum M. tuberculosis more closely resembled aerobic replication and shared similarity in carbon metabolism to in vitro and in vivo models with significant upregulation of genes associated with cholesterol metabolism and downstream propionate detoxification pathways. In addition, and counter to previous reports on intracellular M. tuberculosis infection in vitro, M. tuberculosis in sputum was zinc, but not iron, deprived, and the phoP loci were also significantly downregulated, suggesting that the pathogen is likely extracellular in location.

The contributions of fliG gene to the pathogenicity of Pseudomonas plecoglossicida and pathogen-host interactions with Epinephelus coioides

Pseudomonas plecoglossicida is a Gram-negative aerobic rod-shaped bacterium with polar flagella. It is the causative agent of visceral white spot disease in cultured fish, resulting in serious economic losses. In our previous study, RNA sequencing showed that the expression of the fliG gene in P. plecoglossicida is significantly up-regulated during infection of orange-spotted grouper (Epinephelus coioides). In this study, four P. plecoglossicida RNA interference (RNAi) mutants were successfully constructed by linking four short hairpin RNAs (shRNAs), which target different sites of the fliG gene, to pCM130/tac, respectively. The mRNA expression levels of the fliG gene in P. plecoglossicida were significantly decreased in four mutants. The shRNA-335 mutant (fliG-RNAi strain) showed the best silencing efficiency (88.2%) and was thus chosen for further analysis. Electron microscopy indicated that the flagella of the fliG-RNAi strain of P. plecoglossicida were shorter and finer than those of the wild type strain. The fliG-RNAi strain also showed significantly decreased mobility, chemotaxis, adhesion, and biofilm formation. Furthermore, compared with wild type strain infection, E. coioides infected with the fliG-RNAi strain exhibited a 0.5-d delay in the time of first death and 55% reduction in accumulated mortality, as well as milder splenic symptoms. RNAi of the fliG gene significantly affected the transcriptomes of both pathogen and host in the infected spleens of E. coioides. KEGG analysis revealed that the flagellar assembly pathway, bacterial chemotaxis pathway, and starch and sucrose metabolism pathway were significantly enriched in the pathogen at 3 days post infection (dpi). In contrast, the complement and coagulation cascade pathway and antigen processing and presentation pathway were significantly enriched in the host at 3 dpi. More immune-related pathways were enriched at 5 dpi and more differentially expressed genes were found in the complement and coagulation cascade and antigen processing and presentation pathways. Cytokine-cytokine receptor interaction, hematopoietic cell lineage, and IgA-producing intestinal immune network pathways were significantly enriched in the host at 5 dpi. These results indicate that fliG is an important virulence gene of P. plecoglossicida and contributes to the pathogenicity of P. plecoglossicida as well as pathogen-host interactions with E. coioides.

Reimagining the future of African brain health: Perspectives for basic research on the pathogenesis of cryptococcal meningitis

Cryptococcal meningitis is a fatal opportunistic infection of the brain and a leading cause of neurological damage and death in immunocompromised individuals. This neglected fungal disease of the brain is a huge burden on the health systems of developing countries, especially in Sub-Saharan Africa, where up to 25% of people living with HIV/AIDS succumb to it. Cryptococcal fungal cells have a predilection for the brain and they are capable of traversing the blood brain barrier and invade the brain where they cause infection, inflammation and a disruption of normal brain function. A robust host neuroimmune response is critical for pathogen clearance and survival, and a good understanding of the mechanisms underlying its development in the host is critical for the development of effective treatments. However, past basic research studies have been focussed on the characteristics of the fungus and its effect on the peripheral immune system; with little attention paid to how it interacts with brain immune cells. This mini review briefly discusses the paucity of basic research data on the neuroimmune response to cryptococcal infection, raises pertinent questions on how the brain cells respond to the fungal infection, and thereafter discusses models, techniques and advanced technologies that could be useful for carrying out high-throughput research on the pathogenesis of cryptococcal meningitis.

Mycoplasma genitalium and M. pneumoniae Regulate a Distinct Set of Protein-Coding Genes in Epithelial Cells

Mycoplasma genitalium and M. pneumoniae are two significant mycoplasmas that infect the urogenital and respiratory tracts of humans. Despite distinct tissue tropisms, they both have similar pathogenic mechanisms and infect/invade epithelial cells in the respective regions and persist within these cells. However, the pathogenic mechanisms of these species in terms of bacterium-host interactions are poorly understood. To gain insights on this, we infected HeLa cells independently with M. genitalium and M. pneumoniae and assessed gene expression by whole transcriptome sequencing (RNA-seq) approach. The results revealed that HeLa cells respond to M. genitalium and M. pneumoniae differently by regulating various protein-coding genes. Though there is a significant overlap between the genes regulated by these species, many of the differentially expressed genes were specific to each species. KEGG pathway and signaling network analyses revealed that the genes specific to M. genitalium are more related to cellular processes. In contrast, the genes specific to M. pneumoniae infection are correlated with immune response and inflammation, possibly suggesting that M. pneumoniae has some inherent ability to modulate host immune pathways.

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.

Dual RNA-seq analysis of in vitro infection multiplicity and RNA depletion methods in Chlamydia-infected epithelial cells

Dual RNA-seq experiments examining viral and bacterial pathogens are increasing, but vary considerably in their experimental designs, such as infection rates and RNA depletion methods. Here, we have applied dual RNA-seq to Chlamydia trachomatis infected epithelial cells to examine transcriptomic responses from both organisms. We compared two time points post infection (1 and 24 h), three multiplicity of infection (MOI) ratios (0.1, 1 and 10) and two RNA depletion methods (rRNA and polyA). Capture of bacterial-specific RNA were greatest when combining rRNA and polyA depletion, and when using a higher MOI. However, under these conditions, host RNA capture was negatively impacted. Although it is tempting to use high infection rates, the implications on host cell survival, the potential reduced length of infection cycles and real world applicability should be considered. This data highlights the delicate nature of balancing host-pathogen RNA capture and will assist future transcriptomic-based studies to achieve more specific and relevant infection-related biological insights.

Learning from -omics strategies applied to uncover Haemophilus influenzae host-pathogen interactions: Current status and perspectives

Haemophilus influenzae has contributed to key bacterial genome sequencing hallmarks, as being not only the first bacterium to be genome-sequenced, but also starring the first genome-wide analysis of chromosomes directly transformed with DNA from a divergent genotype, and pioneering Tn-seq methodologies. Over the years, the phenomenal and constantly evolving development of -omic technologies applied to a whole range of biological questions of clinical relevance in the H. influenzae-host interplay, has greatly moved forward our understanding of this human-adapted pathogen, responsible for multiple acute and chronic infections of the respiratory tract. In this way, essential genes, virulence factors, pathoadaptive traits, and multi-layer gene expression regulatory networks with both genomic and epigenomic complexity levels are being elucidated. Likewise, the unstoppable increasing whole genome sequencing information underpinning H. influenzae great genomic plasticity, mainly when referring to non-capsulated strains, poses major challenges to understand the genomic basis of clinically relevant phenotypes and even more, to clearly highlight potential targets of clinical interest for diagnostic, therapeutic or vaccine development. We review here how genomic, transcriptomic, proteomic and metabolomic-based approaches are great contributors to our current understanding of the interactions between H. influenzae and the human airways, and point possible strategies to maximize their usefulness in the context of biomedical research and clinical needs on this human-adapted bacterial pathogen.

Genetic interaction analysis in microbial pathogens: unravelling networks of pathogenesis, antimicrobial susceptibility and host interactions

Genetic interaction (GI) analysis is a powerful genetic strategy that analyzes the fitness and phenotypes of single- and double-gene mutant cells in order to dissect the epistatic interactions between genes, categorize genes into biological pathways, and characterize genes of unknown function. GI analysis has been extensively employed in model organisms for foundational, systems-level assessment of the epistatic interactions between genes. More recently, GI analysis has been applied to microbial pathogens and has been instrumental for the study of clinically important infectious organisms. Here, we review recent advances in systems-level GI analysis of diverse microbial pathogens, including bacterial and fungal species. We focus on important applications of GI analysis across pathogens, including GI analysis as a means to decipher complex genetic networks regulating microbial virulence, antimicrobial drug resistance and host-pathogen dynamics, and GI analysis as an approach to uncover novel targets for combination antimicrobial therapeutics. Together, this review bridges our understanding of GI analysis and complex genetic networks, with applications to diverse microbial pathogens, to further our understanding of virulence, the use of antimicrobial therapeutics and host-pathogen interactions. .

Burkholderia cenocepacia transcriptome during the early contacts with giant plasma membrane vesicles derived from live bronchial epithelial cells

Burkholderia cenocepacia is known for its capacity of adherence and interaction with the host, causing severe opportunistic lung infections in cystic fibrosis patients. In this work we produced Giant Plasma Membrane Vesicles (GPMVs) from a bronchial epithelial cell line and validated their use as a cell-like alternative to investigate the steps involved in the adhesion process of B. cenocepacia. RNA-sequencing was performed and the analysis of the B. cenocepacia K56-2 transcriptome after the first contacts with the surface of host cells allowed the recognition of genes implicated in bacterial adaptation and virulence-associated functions. The sensing of host membranes led to a transcriptional shift that caused a cascade of metabolic and physiological adaptations to the host specific environment. Many of the differentially expressed genes encode proteins related with central metabolic pathways, transport systems, cellular processes, and virulence traits. The understanding of the changes in gene expression that occur in the early steps of infection can uncover new proteins implicated in B. cenocepacia-host cell adhesion, against which new blocking agents could be designed to control the progression of the infectious process.

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.

DualSeqDB: the host-pathogen dual RNA sequencing database for infection processes

Despite antibiotic resistance being a matter of growing concern worldwide, the bacterial mechanisms of pathogenesis remain underexplored, restraining our ability to develop new antimicrobials. The rise of high-throughput sequencing technology has made available a massive amount of transcriptomic data that could help elucidate the mechanisms underlying bacterial infection. Here, we introduce the DualSeqDB database, a resource that helps the identification of gene transcriptional changes in both pathogenic bacteria and their natural hosts upon infection. DualSeqDB comprises nearly 300 000 entries from eight different studies, with information on bacterial and host differential gene expression under in vivo and in vitro conditions. Expression data values were calculated entirely from raw data and analyzed through a standardized pipeline to ensure consistency between different studies. It includes information on seven different strains of pathogenic bacteria and a variety of cell types and tissues in Homo sapiens, Mus musculus and Macaca fascicularis at different time points. We envisage that DualSeqDB can help the research community in the systematic characterization of genes involved in host infection and help the development and tailoring of new molecules against infectious diseases. DualSeqDB is freely available at http://www.tartaglialab.com/dualseq.

Non-typeable Haemophilus influenzae chronic colonization in chronic obstructive pulmonary disease (COPD)

Haemophilus influenzae is the most common cause of bacterial infection in the lungs of chronic obstructive pulmonary disease (COPD) patients and contributes to episodes of acute exacerbation which are associated with increased hospitalization and mortality. Due to the ability of H. influenzae to adhere to host epithelial cells, initial colonization of the lower airways can progress to a persistent infection and biofilm formation. This is characterized by changes in bacterial behaviour such as reduced cellular metabolism and the production of an obstructive extracellular matrix (ECM). Herein we discuss the multiple mechanisms by which H. influenzae contributes to the pathogenesis of COPD. In particular, mechanisms that facilitate bacterial adherence to host airway epithelial cells, biofilm formation, and microbial persistence through immune system evasion and antibiotic tolerance will be discussed.

The combined analysis as the best strategy for Dual RNA-Seq mapping

In Dual RNA-Seq experiments the simultaneous extraction of RNA and analysis of gene expression data from both interacting organisms could be a challenge. One alternative is separating the reads during in silico data analysis. There are two main mapping methods used: sequential and combined. Here we present a combined approach in which the libraries were aligned to a concatenated genome to sort the reads before mapping them to the respective annotated genomes. A comparison of this method with the sequential analysis was performed. Two RNA-Seq libraries available in public databases consisting of a eukaryotic (Zea mays) and a prokaryotic (Herbaspirillum seropediceae) organisms were mixed to simulate a Dual RNA-Seq experiment. Libraries from real Dual RNA-Seq experiments were also used. The sequential analysis consistently attributed more reads to the first reference genome used in the analysis (due to cross-mapping) than the combined approach. More importantly, the combined analysis resulted in lower numbers of cross-mapped reads. Our results highlight the necessity of combining the reference genomes to sort reads previously to the counting step to avoid losing information in Dual RNA-Seq experiments. Since most studies first map the RNA-Seq libraries to the eukaryotic genome, much prokaryotic information has probably been lost.

Repertoire and Diversity of Toxin - Antitoxin Systems of Crohn's Disease-Associated Adherent-Invasive Escherichia coli. New Insight of T his Emergent E. coli Pathotype

Adherent-invasive Escherichia coli (AIEC) corresponds to an E. coli pathovar proposed as a possible agent trigger associated to Crohn's disease. It is characterized for its capacity to adhere and to invade epithelial cells, and to survive and replicate inside macrophages. Mechanisms that allow intestinal epithelium colonization, and host factors that favor AIEC persistence have been partly elucidated. However, bacterial factors involved in AIEC persistence are currently unknown. Toxin-antitoxin (TA) systems are recognized elements involved in bacterial persistence, in addition to have a role in stabilization of mobile genetic elements and stress response. The aim of this study was to elucidate the repertoire and diversity of TA systems in the reference AIEC NRG857c strain and to compare it with AIEC strains whose genomes are available at databases. In addition, toxin expression levels under in vitro stress conditions found by AIEC through the intestine and within the macrophage were measured. Our results revealed that NRG857c encodes at least 33 putative TA systems belonging to types I, II, IV, and V, distributed around all the chromosome, and some in close proximity to genomic islands. A TA toxin repertoire marker of the pathotype was not found and the repertoire of 33 TA toxin genes described here was exclusive of the reference strains, NRG857c and LF82. Most toxin genes were upregulated in the presence of bile salts and acidic pH, as well as within the macrophage. However, different transcriptional responses were detected between reference strains (NRG857c and HM605), recalling the high diversity associated to this pathotype. To our knowledge this is the first analysis of TA systems associated to AIEC and it has revealed new insight associated to this emergent E. coli pathotype.

Microbial extracellular RNAs and their roles in human diseases

IMPACT STATEMENT

The number of commensal bacteria in the body surpasses the number of actual human cells. Thus, various interactions between microbes and human cells constitute an inevitable phenomenon. Recent evidence has led to bacterial extracellular RNAs (exRNAs) being proposed as good candidates for microbe-host inter-kingdom communication tools as they can modulate the expression of host genes. However, research findings on the relevance of interactions between extracellular RNA and human diseases are still in their infancy. Nevertheless, substantial data suggest that microbial exRNAs are implicated in various human diseases both at local and distant sites. By exploring various scenarios for the involvement of microbial exRNAs in human diseases, we may better understand the role of exRNAs as "communication signals" for diseases and thereby develop novel therapeutic strategies by using them and their carrier extracellular vesicles.

Organism dual RNA-seq reveals the importance of BarA/UvrY in Vibrio parahaemolyticus virulence

Elucidation of host-pathogen interaction is essential for developing effective strategies to combat bacterial infection. Dual RNA-Seq using cultured cells or tissues/organs as the host of pathogen has emerged as a novel strategy to understand the responses concurrently from both pathogen and host at cellular level. However, bacterial infection mostly causes systematic responses from the host at organism level where the interplay is urgently to be understood but inevitably being neglected by the current practice. Here, we developed an approach that simultaneously monitor the genome-wide infection-linked transcriptional alterations in both pathogenic Vibrio parahaemolyticus and the infection host nematode Caenorhabditis elegans. Besides the dynamic alterations in transcriptomes of both C. elegans and V. parahaemolyticus during infection, we identify a two-component system, BarA/UvrY, that is important for virulence in host. BarA/UvrY not only controls the virulence factors in V. parahaemolyticus including Type III and Type VI secretion systems, but also attenuates innate immune responses in C. elegans, including repression on the MAP kinase-mediated cascades. Thus, our study exemplifies the use of dual RNA-Seq at organism level to uncover previously unrecognized interplay between host and pathogen.

Bacterial genetics and molecular pathogenesis in the age of high throughput DNA sequencing

When Stanley Falkow introduced Molecular Koch's Postulates (Falkow, 1988) as a conceptual framework to identify microbial factors that contributed to disease, he reaffirmed the prominent role that the basic principles of genetic analysis should play in defining genotype-phenotype associations in microbial pathogens. In classical bacterial genetics the nature of mutations is inferred through cis-trans complementation and by indirectly mapping their relative position and physical distance through recombination frequencies - all of which were made possible by the genetic tools of the day: natural transformations, conjugation and transduction. Unfortunately, many of these genetic tools are not always available to study pathogenic bacteria. The recombinant DNA revolution in the 1980s launched the field of molecular pathogenesis as genes could be treated as physical units that could be cut, spliced and transplanted from one microbe to another and thus not only 'prove' that an individual gene complemented a virulence defect in a mutant strain but also could impart pathogenic properties to otherwise benign microbes. The recombinant DNA revolution also enabled the generation of newer versions of genetic tools to generate mutations and engineer microbial genomes. The last decade has ushered in next generation sequencing technologies as a new powerful tool for bacterial genetics. The routine and inexpensive sequencing of microbial genomes has increased the number and phylogenetic scope of microbes that are amenable to functional characterization and experimentation. In this review, we highlight some salient advances in this rapidly evolving area.

Panel 3: Genomics, precision medicine and targeted therapies

OBJECTIVE

To review the most recent advances in human and bacterial genomics as applied to pathogenesis and clinical management of otitis media.

DATA SOURCES

PubMed articles published since the last meeting in June 2015 up to June 2019.

REVIEW METHODS

A panel of experts in human and bacterial genomics of otitis media was formed. Each panel member reviewed the literature in their respective fields and wrote draft reviews. The reviews were shared with all panel members, and a merged draft was created. The panel met at the 20th International Symposium on Recent Advances in Otitis Media in June 2019, discussed the review and refined the content. A final draft was made, circulated, and approved by the panel members.

CONCLUSION

Trans-disciplinary approaches applying pan-omic technologies to identify human susceptibility to otitis media and to understand microbial population dynamics, patho-adaptation and virulence mechanisms are crucial to the development of novel, personalized therapeutics and prevention strategies for otitis media.

IMPLICATIONS FOR PRACTICE

In the future otitis media prevention strategies may be augmented by mucosal immunization, combination vaccines targeting multiple pathogens, and modulation of the middle ear microbiome. Both treatment and vaccination may be tailored to an individual's otitis media phenotype as defined by molecular profiles obtained by using rapidly developing techniques in microbial and host genomics.

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.

Novel Approaches for Systems Biology of Metabolism-Oriented Pathogen-Human Interactions: A Mini-Review

Pathogenic microorganisms exploit host metabolism for sustained survival by rewiring its metabolic interactions. Therefore, several metabolic changes are induced in both pathogen and host cells in the course of infection. A systems-based approach to elucidate those changes includes the integrative use of genome-scale metabolic networks and molecular omics data, with the overall goal of better characterizing infection mechanisms for novel treatment strategies. This review focuses on novel aspects of metabolism-oriented systems-based investigation of pathogen-human interactions. The reviewed approaches are the generation of dual-omics data for the characterization of metabolic signatures of pathogen-host interactions, the reconstruction of pathogen-host integrated genome-scale metabolic networks, which has a high potential to be applied to pathogen-gut microbiota interactions, and the structure-based analysis of enzymes playing role in those interactions. The integrative use of those approaches will pave the way for the identification of novel biomarkers and drug targets for the prediction and prevention of infectious diseases.

Borrelia burgdorferi hijacks cellular metabolism of immune cells: Consequences for host defense

Changes in cellular metabolism have proven to be important factors in driving cell behavior. It has been shown that cellular metabolism of immune cells changes when exposed to or infected by several pathogens: while this is often an adaptation of the host cells to the infection, sometimes it represents a mechanism through which the pathogens evade immune activation. Borrelia burgdorferi sensu lato, the causative agent of Lyme borreliosis, is a pathogen that highly depends on the host to survive, as the bacterium lacks many central metabolic pathways to generate its own nutrients. It is therefore quite likely that the bacterium interacts with host cells to obtain these metabolites and thereby affects metabolism in the host. Previously, several studies have assessed metabolic pathways in B. burgdorferi s.l. and how it adapts to its different host species. However, few studies have looked into how the interaction with the bacterium might affect the host cell metabolism. In this review we present the major metabolic pathways activated during Lyme borreliosis, viewed from both bacterium and host metabolism, and we discuss how these pathways interact with each other, and how they influence pathogenesis of Lyme borreliosis.

Dual and Triple Epithelial Coculture Model Systems with Donor-Derived Microbiota and THP-1 Macrophages To Mimic Host-Microbe Interactions in the Human Sinonasal Cavities

Despite the relevance of the resident microbiota in sinonasal health and disease and the need for cross talk between immune and epithelial cells in the upper respiratory tract, these parameters have not been combined in a single in vitro model system. We have developed a coculture system of differentiated respiratory epithelium and natural nasal microbiota and incorporated an immune component. As indicated by absence of cytotoxicity and stable cytokine profiles and epithelial integrity, nasal microbiota from human origin appeared to be well tolerated by host cells, while microbial community composition remained representative for that of the human (sino)nasal cavity. Importantly, the introduction of macrophage-like cells enabled us to obtain a differential readout from the epithelial cells dependent on the donor microbial background to which the cells were exposed. We conclude that both model systems offer the means to investigate host-microbe interactions in the upper respiratory tract in a more representative way.

The epithelium of the human sinonasal cavities is colonized by a diverse microbial community, modulating epithelial development and immune priming and playing a role in respiratory disease. Here, we present a novel in vitro approach enabling a 3-day coculture of differentiated Calu-3 respiratory epithelial cells with a donor-derived bacterial community, a commensal species (Lactobacillus sakei), or a pathobiont (Staphylococcus aureus). We also assessed how the incorporation of macrophage-like cells could have a steering effect on both epithelial cells and the microbial community. Inoculation of donor-derived microbiota in our experimental setup did not pose cytotoxic stress on the epithelial cell layers, as demonstrated by unaltered cytokine and lactate dehydrogenase release compared to a sterile control. Epithelial integrity of the differentiated Calu-3 cells was maintained as well, with no differences in transepithelial electrical resistance observed between coculture with donor-derived microbiota and a sterile control. Transition of nasal microbiota from in vivo to in vitro conditions maintained phylogenetic richness, and yet a decrease in phylogenetic and phenotypic diversity was noted. Additional inclusion and coculture of THP-1-derived macrophages did not alter phylogenetic diversity, and yet donor-independent shifts toward higher Moraxella and Mycoplasma abundance were observed, while phenotypic diversity was also increased. Our results demonstrate that coculture of differentiated airway epithelial cells with a healthy donor-derived nasal community is a viable strategy to mimic host-microbe interactions in the human upper respiratory tract. Importantly, including an immune component allowed us to study host-microbe interactions in the upper respiratory tract more in depth.

IMPORTANCE Despite the relevance of the resident microbiota in sinonasal health and disease and the need for cross talk between immune and epithelial cells in the upper respiratory tract, these parameters have not been combined in a single in vitro model system. We have developed a coculture system of differentiated respiratory epithelium and natural nasal microbiota and incorporated an immune component. As indicated by absence of cytotoxicity and stable cytokine profiles and epithelial integrity, nasal microbiota from human origin appeared to be well tolerated by host cells, while microbial community composition remained representative for that of the human (sino)nasal cavity. Importantly, the introduction of macrophage-like cells enabled us to obtain a differential readout from the epithelial cells dependent on the donor microbial background to which the cells were exposed. We conclude that both model systems offer the means to investigate host-microbe interactions in the upper respiratory tract in a more representative way.

Roles of aminoacyl-tRNA synthetases in immune regulation and immune diseases

Aminoacyl-tRNA synthetases (ARSs) play a vital role in protein synthesis by linking amino acids to their cognate transfer RNAs (tRNAs). This typical function has been well recognized over the past few decades. However, accumulating evidence reveals that ARSs are involved in a wide range of physiological and pathological processes apart from translation. Strikingly, certain ARSs are closely related to different types of immune responses. In this review, we address the infection and immune responses induced by pathogen ARSs, as well as the potential anti-infective compounds that target pathogen ARSs. Meanwhile, we describe the functional mechanisms of ARSs in the development of immune cells. In addition, we focus on the roles of ARSs in certain immune diseases, such as autoimmune diseases, infectious diseases, and tumor immunity. Although our knowledge of ARSs in the immunological context is still in its infancy, research in this field may provide new ideas for the treatment of immune-related diseases.

Expression of the Nontypeable Haemophilus influenzae Type IV Pilus Is Stimulated by Coculture with Host Respiratory Tract Epithelial Cells

The type IV pilus (Tfp) of nontypeable Haemophilus influenzae (NTHI) mediates adherence, colonization, motility, and biofilm formation, and the major protein subunit, PilA, is a promising vaccine candidate. Thus, it is crucial to understand how Tfp expression is regulated within the microenvironments of the human nasopharynx, which NTHI colonizes asymptomatically, and the more distal regions of the respiratory tract where NTHI-induced diseases occur. Here, we examined the effects of coculture of NTHI with human airway epithelial cells and heme availability on Tfp expression at temperatures typical of the human nasopharynx (34°C) or warmer anatomical sites during infection (37°C). Tfp expression was estimated by pilA promoter activity, pilA gene expression, and relative abundances of PilA and pilin protein. The results revealed that at both temperatures, NTHI cocultured with airway epithelial cells demonstrated significantly greater expression of pilA, PilA/pilin protein, and likely, fully assembled Tfp than NTHI cultured on an abiotic surface. Because NTHI is a heme auxotroph, we hypothesized that availability of heme from host cells might be a signal for Tfp expression. Thereby, we cultured NTHI in iron-limited medium, and we observed that supplementation with heme significantly increased pilA promoter activity. Collectively, our data suggested that NTHI Tfp expression was stimulated by soluble factor(s) released by epithelial cells, which are present in all microenvironments of the respiratory tract. The expression of this target antigen under conditions that mimic the human airway strongly supports the rationale for the use of PilA as a vaccine immunogen to prevent NTHI-induced diseases of the respiratory tract.

Cellular Immune Response against Nontypeable Haemophilus influenzae Infecting the Preinflamed Middle Ear of the Junbo Mouse

Nontypeable Haemophilus influenzae (NTHi) is a major pathogen causing acute otitis media (AOM). The pathology of AOM increases during long-term infection in the middle ear (ME), but the host cellular immune response to bacterial infection in this inflamed environment is poorly understood. Using the Junbo mouse, a characterized NTHi infection model, we analyzed the cellular response to NTHi infection in the Junbo mouse middle ear fluid (MEF). NTHi infection increased the total cell number and significantly decreased the proportion of live cells in the MEF at day 1, and this further decreased gradually on each day up to day 7. Flow cytometry analysis showed that neutrophils were the dominant immune cell population in the MEF and that NTHi infection significantly increased their proportion whereas it decreased the monocyte, macrophage, and dendritic cell proportions. Neutrophil and macrophage numbers increased in blood and spleen after NTHi infection. The T-cell population was dominated by T-helper (Th) cells in noninoculated MEF, and the effector Th (CD44⁺) cell population increased at day 2 of NTHi infection with an increase in IL-12p40 levels. Sustained NTHi infection up to 3 days increased the transforming growth factor β levels, decreasing the effector cell population and increasing the T-regulatory (T-reg) cell population. In the preinflamed ME environment of the Junbo mouse, neutrophils are the first responder to NTHi infection followed by T-reg immune suppressive cells. These data indicate that sustained NTHi infection in the ME induces the immune suppressive response by inducing the T-reg cell population and reducing immune cell infiltration, thus promoting longer-term infection.

Molecular Signatures of Non-typeable Haemophilus influenzae Lung Adaptation in Pediatric Chronic Lung Disease

Non-typeable Haemophilus influenzae (NTHi), an opportunistic pathogen of the upper airways of healthy children, can infect the lower airways, driving chronic lung disease. However, the molecular basis underpinning NTHi transition from a commensal to a pathogen is not clearly understood. Here, we performed comparative genomic and transcriptomic analyses of 12 paired, isogenic NTHi strains, isolated from the nasopharynx (NP) and bronchoalveolar lavage (BAL) of 11 children with chronic lung disease, to identify convergent molecular signatures associated with lung adaptation. Comparative genomic analyses of the 12 NP-BAL pairs demonstrated that five were genetically identical, with the remaining seven differing by only 1 to 3 mutations. Within-patient transcriptomic analyses identified between 2 and 58 differentially expressed genes in 8 of the 12 NP-BAL pairs, including pairs with no observable genomic changes. Whilst no convergence was observed at the gene level, functional enrichment analysis revealed significant under-representation of differentially expressed genes belonging to Coenzyme metabolism, Function unknown, Translation, ribosomal structure, and biogenesis Cluster of Orthologous Groups categories. In contrast, Carbohydrate transport and metabolism, Cell motility and secretion, Intracellular trafficking and secretion, and Energy production categories were over-represented. This observed trend amongst genetically unrelated NTHi strains provides evidence of convergent transcriptional adaptation of NTHi to pediatric airways that deserves further exploration. Understanding the pathoadaptative mechanisms that NTHi employs to infect and persist in the lower pediatric airways is essential for devising targeted diagnostics and treatments aimed at minimizing disease severity, and ultimately, preventing NTHi lung infections and subsequent chronic lung disease in children.

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

Bacterial pathogens subvert host cells by manipulating cellular pathways for survival and replication; in turn, host cells respond to the invading pathogen through cascading changes in gene expression. Deciphering these complex temporal and spatial dynamics to identify novel bacterial virulence factors or host response pathways is crucial for improved diagnostics and therapeutics. Dual RNA sequencing (dRNA-Seq) has recently been developed to simultaneously capture host and bacterial transcriptomes from an infected cell. This approach builds on the high sensitivity and resolution of RNA sequencing technology and is applicable to any bacteria that interact with eukaryotic cells, encompassing parasitic, commensal or mutualistic lifestyles. Several laboratory protocols have been presented that outline the collection, extraction and sequencing of total RNA for dRNA-Seq experiments, but there is relatively little guidance available for the detailed bioinformatic analyses required. This protocol outlines a typical dRNA-Seq experiment, based on a Chlamydia trachomatis-infected host cell, with a detailed description of the necessary bioinformatic analyses with currently available software tools.

clpV is a key virulence gene during in vivo Pseudomonas plecoglossicida infection

Interaction between bacterial pathogen and aquatic animal host is exceedingly complex, which involves large dynamic changes in gene expression during different stages of the disease. However, research on identifying key virulence genes based on the dynamics of gene expression changes of a one-sided bacterial pathogen in tissue has not been reported so far across different stages of infectious disease. The clpV for the T6SS of Pseudomonas plecoglossicida was identified for a candidate for key virulence gene based on dynamic changes of gene expression. For the Epinephelus coioides infected using clpV-RNAi strain, no deaths were observed up to 20 dpi. The spleens, kidneys and livers of all the E. coioides that received clpV-RNAi strain failed to develop visible nodules at 5-8 dpi, with the swelling gradually disappearing. The burdens of clpV-RNAi strain in the spleen and blood were greatly reduced at most of the time points after injection, and the burdens of clpV-RNAi strain in the head kidneys and trunk kidneys also had a sharp reduction from 72 to 120 hpi. This paper provides a new insight into the discovery of key virulence genes of pathogens in infected tissue systems.

Determination of an Interaction Network between an Extracellular Bacterial Pathogen and the Human Host

Dual RNA sequencing (RNA-seq) offers the promise of determining an interactome at a transcriptional level between a bacterium and the host but has yet to be done on any bacterial infection in human tissue. We performed dual RNA-seq and metabolomics analyses on wounded and infected sites following experimental infection of the arm with H. ducreyi. Our results suggest that H. ducreyi survives in an abscess by utilizing l-ascorbate as an alternative carbon source, possibly taking advantage of host ascorbic acid recycling, and that H. ducreyi also adapts by upregulating genes involved in anaerobic metabolism and inorganic ion and nutrient transport. To our knowledge, this is the first description of an interaction network between a bacterium and the human host at a site of infection.

A major gap in understanding infectious diseases is the lack of information about molecular interaction networks between pathogens and the human host. Haemophilus ducreyi causes the genital ulcer disease chancroid in adults and is a leading cause of cutaneous ulcers in children in the tropics. We developed a model in which human volunteers are infected on the upper arm with H. ducreyi until they develop pustules. To define the H. ducreyi and human interactome, we determined bacterial and host transcriptomic and host metabolomic changes in pustules. We found that in vivoH. ducreyi transcripts were distinct from those in the inocula, as were host transcripts in pustule and wounded control sites. Many of the upregulated H. ducreyi genes were found to be involved in ascorbic acid and anaerobic metabolism and inorganic ion/nutrient transport. The top 20 significantly expressed human pathways showed that all were involved in immune responses. We generated a bipartite network for interactions between host and bacterial gene transcription; multiple positively correlated networks contained H. ducreyi genes involved in anaerobic metabolism and host genes involved with the immune response. Metabolomic studies showed that pustule and wounded samples had different metabolite compositions; the top ion pathway involved ascorbate and aldarate metabolism, which correlated with the H. ducreyi transcriptional response and upregulation of host genes involved in ascorbic acid recycling. These data show that an interactome exists between H. ducreyi and the human host and suggest that H. ducreyi exploits the metabolic niche created by the host immune response.

Dual RNA-Seq Unveils Pseudomonas plecoglossicida htpG Gene Functions During Host-Pathogen Interactions With Epinephelus coioides

Pseudomonas plecoglossicida is a temperature-dependent opportunistic pathogen which is associated with a variety of diseases in fish. During the development of "white nodules" disease, the expression of htpG in P. plecoglossicida was found to be significantly up-regulated at its virulent temperature of 18°C. The infection of htpG-RNAi strain resulted in the onset time delay, reduction in mortality and infection symptoms in spleen of Epinephelus coioides, and affected the bacterial tissue colonization. In order to reveal the effect of htpG silencing of P. plecoglossicida on the virulence regulation in P. plecoglossicida and immune response in E. coioides, dual RNA-seq was performed and a pathogen-host integration network was constructed. Our results showed that infection induced the expression of host genes related to immune response, but attenuated the expression of bacterial virulence genes. Novel integration was found between host immune genes and bacterial virulence genes, while IL6, IL1R2, IL1B, and TLR5 played key roles in the network. Further analysis with GeneMANIA indicated that flgD and rplF might play key roles during the htpG-dependent virulence regulation, which was in accordance with the reduced biofilm production, motility and virulence in htpG-RNAi strain. Meanwhile, IL6 and IL1B were found to play key roles during the defense against P. plecoglossicida, while CELA2, TRY, CPA1, CPA2, and CPB1 were important targets for P. plecoglossicida attacking to the host.

Dual RNA-seq reveals the effect of the flgM gene of Pseudomonas plecoglossicida on the immune response of Epinephelus coioides

Pseudomonas plecoglossicida is an important and highly pathogenic bacterium for aquaculture and causes serious losses. The expression level of flgM was found to be significantly upregulated post-infection compared with in vitro results, which was confirmed by quantitative real-time PCR. RNAi significantly reduced the expression level of flgM mRNA of P. plecoglossicida. Compared with infection with the wild-type strain, infection with the flgM-RNAi strain resulted in a delay in death and a 75% reduction in the mortality of Epinephelus coioides, followed by alleviation of the symptoms in E. coioides spleen. Moreover, compared with infection with the wild-type strain, infection with the flgM-RNAi strain of P. plecoglossicida resulted in a significant change in the transcriptome of the spleens of infected E. coioides and P. plecoglossicida. KEGG analysis for E. coioides showed that genes of 17 immune pathways were most affected by flgM-RNAi of P. plecoglossicida. Among them, the expression of mhc2, zap70, rhoh, tlr2, ca79a, hcst and cd32 in E. coioides spleen was predicted to be negatively related to flgM in P. plecoglossicida but positively related to genes involved in communication, metabolism and motility.