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

Bacterial Persisters and Infection: Past, Present, and Progressing

Persisters are nongrowing, transiently antibiotic-tolerant bacteria within a clonal population of otherwise susceptible cells. Their formation is triggered by environmental cues and involves the main bacterial stress response pathways that allow persisters to survive many harsh conditions, including antibiotic exposure. During infection, bacterial pathogens are exposed to a vast array of stresses in the host and form nongrowing persisters that survive both antibiotics and host immune responses, thereby most likely contributing to the relapse of many infections. While antibiotic persisters have been extensively studied over the last decade, the bulk of the work has focused on how these bacteria survive exposure to drugs in vitro. The ability of persisters to survive their interaction with a host is important yet underinvestigated. In order to tackle the problem of persistence of infections that contribute to the worldwide antibiotic resistance crisis, efforts should be made by scientific communities to understand and merge these two fields of research: antibiotic persisters and host-pathogen interactions. Here we give an overview of the history of the field of antibiotic persistence, report evidence for the importance of persisters in infection, and highlight studies that bridge the two areas.

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.

Transcriptional profiling of human macrophages during infection with Bordetella pertussis

Bordetella pertussis, a strictly human re-emerging pathogen and the causative agent of whooping cough, exploits a broad variety of virulence factors to establish efficient infection. Here, we used RNA sequencing to analyse the changes in gene expression profiles of human THP-1 macrophages resulting from B. pertussis infection. In parallel, we attempted to determine the changes in intracellular B. pertussis-specific transcriptomic profiles resulting from interaction with macrophages. Our analysis revealed that global gene expression profiles in THP-1 macrophages are extensively rewired 6 h post-infection. Among the highly expressed genes, we identified those encoding cytokines, chemokines, and transcription regulators involved in the induction of the M1 and M2 macrophage polarization programmes. Notably, several host genes involved in the control of apoptosis and inflammation which are known to be hijacked by intracellular bacterial pathogens were overexpressed upon infection. Furthermore, in silico analyses identified large temporal changes in expression of specific gene subsets involved in signalling and metabolic pathways. Despite limited numbers of the bacterial reads, we observed reduced expression of majority of virulence factors and upregulation of several transcriptional regulators during infection suggesting that intracellular B. pertussis cells switch from virulent to avirulent phase and actively adapt to intracellular environment, respectively.

Burkholderia pseudomallei pathogenesis and survival in different niches

Burkholderia pseudomallei (Bp) is the causative agent of melioidosis, a disease of the tropics with high clinical mortality rates. To date, no vaccines are approved for melioidosis and current treatment relies on antibiotics. Conversely, common misdiagnosis and high pathogenicity of Bp hamper efforts to fight melioidosis. This bacterium can be isolated from a wide range of niches such as waterlogged fields, stagnant water bodies, salt water bodies and from human and animal clinical specimens. Although extensive studies have been undertaken to elucidate pathogenesis mechanisms of Bp, little is known about how a harmless soil bacterium adapts to different environmental conditions, in particular, the shift to a human host to become a highly virulent pathogen. The bacterium has a large genome encoding an armory of factors that assist the pathogen in surviving under stressful conditions and assuming its role as a deadly intracellular pathogen. This review presents an overview of what is currently known about how the pathogen adapts to different environments. With in-depth understanding of Bp adaptation and survival, more effective therapies for melioidosis can be developed by targeting related genes or proteins that play a major role in the bacteria's survival.

Next-generation technologies for studying host-pathogen interactions: a focus on dual transcriptomics, CRISPR/Cas9 screening and organs-on-chips

Pathogens constantly interact with their hosts and the environment, and therefore have evolved unique virulence mechanisms to target and breach host defense barriers and manipulate host immune response to establish an infection. Advances in technologies that allow genome mining, gene editing such as CRISPR/Cas9, genomic, epigenomic and transcriptomic studies such as dual RNA-seq, coupled with bioinformatics, have accelerated the field of host-pathogen interactions within a broad range of infection models. Underpinning of the molecular changes that accompany invasion of eukaryotic cells with pathogenic microorganisms at the intersection of host, pathogen and their local environment has provided a better understanding of infectious disease mechanisms and antimicrobial strategies. The recent evolution of physiologically relevant three-dimensional (3-D) tissue/organ models and microfluidic organ-on-chip devices also provided a window to a more predictive framework of infectious disease processes. These approaches combined hold the potential to highly impact discovery of novel drug targets and vaccine candidates of the future. Here, we review three of the available and emerging technologies-dual RNA-seq, CRISPR/Cas9 screening and organs-on-chips, applicable to the high throughput study and deciphering of interaction networks between pathogens and their hosts that are critical for the development of novel therapeutics.

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.

Long non-coding RNA molecules in tuberculosis

Tuberculosis (TB), a chronic disease caused by Mycobacterium tuberculosis, is one of the deadliest infectious diseases in the world. Despite significant advances in detection techniques and therapeutic approaches for tuberculosis, there is still no suitable solution for early screening and reducing the number of individuals affected and their effective treatment. Various cellular events can disrupt the development of TB. The basis of these events is dysregulating of genes expression patterns related with specific molecules. Long non-coding RNAs (lncRNAs) are molecules discovered to regulate the expression of protein-coding genes and participate in gene silencing, cell cycle regulation and cellular differentiation processes. Dysregulation of lncRNAs has been found to be associated with many diseases, including cancers and infectious diseases. Thus, the recognition of lncRNAs as novel molecular biomarkers and therapeutic targets for tuberculosis is promising. In the present review, we try to summarize the current findings of lncRNA expression patterns and its role in tuberculosis infection process.

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.

Disentangling sRNA-Seq data to study RNA communication between species

Many organisms exchange small RNAs (sRNAs) during their interactions, that can target or bolster defense strategies in host-pathogen systems. Current sRNA-Seq technology can determine the sRNAs present in any symbiotic system, but there are very few bioinformatic tools available to interpret the results. We show that one of the biggest challenges comes from sequences that map equally well to the genomes of both interacting organisms. This arises due to the small size of the sRNAs compared to large genomes, and because a large portion of sequenced sRNAs come from genomic regions that encode highly conserved miRNAs, rRNAs or tRNAs. Here, we present strategies to disentangle sRNA-Seq data from samples of communicating organisms, developed using diverse plant and animal species that are known to receive or exchange RNA with their symbionts. We show that sequence assembly, both de novo and genome-guided, can be used for these sRNA-Seq data, greatly reducing the ambiguity of mapping reads. Even confidently mapped sequences can be misleading, so we further demonstrate the use of differential expression strategies to determine true parasite-derived sRNAs within host cells. We validate our methods on new experiments designed to probe the nature of the extracellular vesicle sRNAs from the parasitic nematode Heligmosomoides bakeri that get into mouse intestinal epithelial cells.

A benchmark of hemoglobin blocking during library preparation for mRNA-Sequencing of human blood samples

RNA-Sequencing (RNA-Seq) of peripheral blood can be a valuable source of information for investigating the status and mechanism of diseases. However, blood contains 50–80% unwanted hemoglobin (Hb) transcripts. Lexogen’s QuantSeq mRNA-Seq-Kit for Illumina RNA-Seq features a ‘Globin Block’ (GB) module that depletes Hb cDNAs during library preparation. Here, we aimed to assess GB’s effectiveness and checked for technical biases attributable to GB. Using whole blood total RNA samples of 91 healthy individuals, we sequenced 91 pairs of GB and non-blocked samples (noGB) on Illumina HiSeq2500 and 8 pairs of GB/noGB technical replicates on HiSeq4000. GB reduced the fraction of Hb transcripts from 43% (s.d. 14%) to 8.0% (s.d. 4.3%). From GB samples we detected 1,397 more expressed genes at approximately 11 million reads per RNA-isolate. Enrichment and differential expression analyses did not reveal significant differences for GB and noGB samples with respect to molecular function. In contrast to results from studies that have examined the performance of GB during RNA isolation, we were able to assign GB to corresponding noGB samples (from multiple sequencing runs on HiSeq2500) with at least 89.8% accuracy from the complete correlation matrix of all GB/GB, noGB/noGB and GB/noGB pairs. However, the use of different sequencers (HiSeq2500 vs HiSeq4000) impaired assignment of technical replicates, whereas assignment of GB to corresponding noGB samples worked perfectly when sequencing on one lane on HiSeq4000. Lexogen’s GB RNA-Seq module is a valuable addition during mRNA-Seq library preparation which works even with low amounts of input total RNA (50 ng per sample). GB facilitated the detection of low abundant transcripts and yielded more non-hemoglobin reads, while preserving biological information. We observed that differences in sequencing run and platform have a far greater effect on technical variation than the use of GB.

Atypical Salmonella enterica Serovars in Murine and Human Macrophage Infection Models

Nontyphoidal Salmonella species are globally disseminated pathogens and are the predominant cause of gastroenteritis. The pathogenesis of salmonellosis has been extensively studied using in vivo murine models and cell lines, typically challenged with Salmonella enterica serovar Typhimurium. Although S. enterica serovars Enteritidis and Typhimurium are responsible for most of the human infections reported to the Centers for Disease Control and Prevention (CDC), several other serovars also contribute to clinical cases of salmonellosis. Despite their epidemiological importance, little is known about their infection phenotypes. Here, we report the virulence characteristics and genomes of 10 atypical S. enterica serovars linked to multistate foodborne outbreaks in the United States. We show that the murine RAW 264.7 macrophage model of infection is unsuitable for inferring human-relevant differences in nontyphoidal Salmonella infections, whereas differentiated human THP-1 macrophages allowed these isolates to be further characterized in a more human-relevant context.

Bacterial RNA in extracellular vesicles: A new regulator of host-pathogen interactions?

Extracellular vesicles (EVs) are released by cells from all kingdoms and represent one form of cell-cell interaction. This universal system of communication blurs cells type boundaries, offering an new avenue for pathogens to infect their hosts. EVs carry with them an arsenal of virulence factors that have been the focus of numerous studies. During the last years, the RNA content of EVs has also gained increasing attention, particularly in the context of infection. Secreted RNA in EVs (evRNA) from several bacterial pathogens have been characterised but the exact mechanisms promoting pathogenicity remain elusive. In this review, we evaluate the potential of such transcripts to directly interact with targets in infected cells and, by this, represent a novel angle of host cell manipulation during bacterial infection. This article is part of a Special Issue entitled: RNA and gene control in bacteria edited by Dr. M. Guillier and F. Repoila.

Trans-Acting Small RNAs and Their Effects on Gene Expression in Escherichia coli and Salmonella enterica

The last few decades have led to an explosion in our understanding of the major roles that small regulatory RNAs (sRNAs) play in regulatory circuits and the responses to stress in many bacterial species. Much of the foundational work was carried out with Escherichia coli and Salmonella enterica serovar Typhimurium. The studies of these organisms provided an overview of how the sRNAs function and their impact on bacterial physiology, serving as a blueprint for sRNA biology in many other prokaryotes. They also led to the development of new technologies. In this chapter, we first summarize how these sRNAs were identified, defining them in the process. We discuss how they are regulated and how they act and provide selected examples of their roles in regulatory circuits and the consequences of this regulation. Throughout, we summarize the methodologies that were developed to identify and study the regulatory RNAs, most of which are applicable to other bacteria. Newly updated databases of the known sRNAs in E. coli K-12 and S. enterica Typhimurium SL1344 serve as a reference point for much of the discussion and, hopefully, as a resource for readers and for future experiments to address open questions raised in this review.

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

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

A three-dimensional intestinal tissue model reveals factors and small regulatory RNAs important for colonization with Campylobacter jejuni

The Gram-negative Epsilonproteobacterium Campylobacter jejuni is currently the most prevalent bacterial foodborne pathogen. Like for many other human pathogens, infection studies with C. jejuni mainly employ artificial animal or cell culture models that can be limited in their ability to reflect the in-vivo environment within the human host. Here, we report the development and application of a human three-dimensional (3D) infection model based on tissue engineering to study host-pathogen interactions. Our intestinal 3D tissue model is built on a decellularized extracellular matrix scaffold, which is reseeded with human Caco-2 cells. Dynamic culture conditions enable the formation of a polarized mucosal epithelial barrier reminiscent of the 3D microarchitecture of the human small intestine. Infection with C. jejuni demonstrates that the 3D tissue model can reveal isolate-dependent colonization and barrier disruption phenotypes accompanied by perturbed localization of cell-cell junctions. Pathogenesis-related phenotypes of C. jejuni mutant strains in the 3D model deviated from those obtained with 2D-monolayers, but recapitulated phenotypes previously observed in animal models. Moreover, we demonstrate the involvement of a small regulatory RNA pair, CJnc180/190, during infections and observe different phenotypes of CJnc180/190 mutant strains in 2D vs. 3D infection models. Hereby, the CJnc190 sRNA exerts its pathogenic influence, at least in part, via repression of PtmG, which is involved in flagellin modification. Our results suggest that the Caco-2 cell-based 3D tissue model is a valuable and biologically relevant tool between in-vitro and in-vivo infection models to study virulence of C. jejuni and other gastrointestinal pathogens.

Enteric pathogens have evolved numerous strategies to successfully colonize and persist in the human gastrointestinal tract. However, especially for the research of virulence mechanisms of human pathogens, often only limited infection models are available. Here, we have applied and further advanced a tissue-engineered human intestinal tissue model based on an extracellular matrix scaffold reseeded with human cells that can faithfully mimic pathogenesis-determining processes of the zoonotic pathogen Campylobacter jejuni. Our three-dimensional (3D) intestinal infection model allows for the assessment of epithelial barrier function during infection as well as for the quantification of bacterial adherence, internalization, and transmigration. Investigation of C. jejuni mutant strains in our 3D tissue model revealed isolate-specific infection phenotypes, in-vivo relevant infection outcomes, and uncovered the involvement of a small RNA pair during C. jejuni pathogenesis. Overall, our results demonstrate the power of tissue-engineered models for studying host-pathogen interactions, and our model will also be helpful to investigate other gastrointestinal pathogens.

An Advanced Human Intestinal Coculture Model Reveals Compartmentalized Host and Pathogen Strategies during Salmonella Infection

A major obstacle in infection biology is the limited ability to recapitulate human disease trajectories in traditional cell culture and animal models, which impedes the translation of basic research into clinics. Here, we introduce a three-dimensional (3D) intestinal tissue model to study human enteric infections at a level of detail that is not achieved by conventional two-dimensional monocultures. Our model comprises epithelial and endothelial layers, a primary intestinal collagen scaffold, and immune cells. Upon Salmonella infection, the model mimics human gastroenteritis, in that it restricts the pathogen to the epithelial compartment, an advantage over existing mouse models. Application of dual transcriptome sequencing to the Salmonella-infected model revealed the communication of epithelial, endothelial, monocytic, and natural killer cells among each other and with the pathogen. Our results suggest that Salmonella uses its type III secretion systems to manipulate STAT3-dependent inflammatory responses locally in the epithelium without accompanying alterations in the endothelial compartment. Our approach promises to reveal further human-specific infection strategies employed by Salmonella and other pathogens.IMPORTANCE Infection research routinely employs in vitro cell cultures or in vivo mouse models as surrogates of human hosts. Differences between murine and human immunity and the low level of complexity of traditional cell cultures, however, highlight the demand for alternative models that combine the in vivo-like properties of the human system with straightforward experimental perturbation. Here, we introduce a 3D tissue model comprising multiple cell types of the human intestinal barrier, a primary site of pathogen attack. During infection with the foodborne pathogen Salmonella enterica serovar Typhimurium, our model recapitulates human disease aspects, including pathogen restriction to the epithelial compartment, thereby deviating from the systemic infection in mice. Combination of our model with state-of-the-art genetics revealed Salmonella-mediated local manipulations of human immune responses, likely contributing to the establishment of the pathogen's infection niche. We propose the adoption of similar 3D tissue models to infection biology, to advance our understanding of molecular infection strategies employed by bacterial pathogens in their human host.

Whole-genome and time-course dual RNA-Seq analyses reveal chronic pathogenicity-related gene dynamics in the ginseng rusty root rot pathogen Ilyonectria robusta

Ilyonectria robusta causes rusty root rot, the most devastating chronic disease of ginseng. Here, we for the first time report the high-quality genome of the I. robusta strain CD-56. Time-course (36 h, 72 h, and 144 h) dual RNA-Seq analysis of the infection process was performed, and many genes, including candidate effectors, were found to be associated with the progression and success of infection. The gene expression profile of CD-56 showed a trend of initial inhibition and then gradually returned to a profile similar to that of the control. Analyses of the gene expression patterns and functions of pathogenicity-related genes, especially candidate effector genes, indicated that the stress response changed to an adaptive response during the infection process. For ginseng, gene expression patterns were highly related to physiological conditions. Specifically, the results showed that ginseng defenses were activated by CD-56 infection and persisted for at least 144 h thereafter but that the mechanisms invoked were not effective in preventing CD-56 growth. Moreover, CD-56 did not appear to fully suppress plant defenses, even in late stages after infection. Our results provide new insight into the chronic pathogenesis of CD-56 and the comprehensive and complex inducible defense responses of ginseng root to I. robusta infection.

ShinyOmics: collaborative exploration of omics-data

BACKGROUND

Omics-profiling is a collection of increasingly prominent approaches that result in large-scale biological datasets, for instance capturing an organism's behavior and response in an environment. It can be daunting to manually analyze and interpret such large datasets without some programming experience. Additionally, with increasing amounts of data; management, storage and sharing challenges arise.

RESULTS

Here, we present ShinyOmics, a web-based application that allows rapid collaborative exploration of omics-data. By using Tn-Seq, RNA-Seq, microarray and proteomics datasets from two human pathogens, we exemplify several conclusions that can be drawn from a rich dataset. We identify a protease and several chaperone proteins upregulated under aminoglycoside stress, show that antibiotics with the same mechanism of action trigger similar transcriptomic responses, point out the dissimilarity in different omics-profiles, and overlay the transcriptional response on a metabolic network.

CONCLUSIONS

ShinyOmics is easy to set up and customize, and can utilize user supplied metadata. It offers several visualization and comparison options that are designed to assist in novel hypothesis generation, as well as data management, online sharing and exploration. Moreover, ShinyOmics can be used as an interactive supplement accompanying research articles or presentations.

Clinical metagenomics for infectious corneal ulcers: Rags to riches?

The emergence of clinical metagenomics as an unbiased, hypothesis-free approach to diagnostic testing is set to fundamentally alter the way infectious diseases are detected. Long envisioned as the solution to the limitations of culture-based conventional microbiology, next generation sequencing methods will soon mature, and our attention will inevitably turn to how they can be applied to areas of medicine which need it most urgently. In ophthalmology, the demand for this technology is particularly pressing for the care of infectious corneal ulcers, where current diagnostic tests may fail to identify a causative organism in over half of cases. However, the optimism found in the budding discourse surrounding clinical metagenomics belies the reality that clinicians and scientists will soon be inundated by oppressive volumes of sequencing data, much of which will be foreign and unfamiliar. Therefore, our success in translating clinical metagenomics is likely to hinge on how we make sense of these data, and understanding its implications for the interpretation and implementation of sequencing into routine clinical care. In this consortium-led review, we provide an outline of these data-related issues and how they may be used to inform technical workflows, with the hope that we may edge closer to realizing the potential of clinical metagenomics for this important unmet need.

Hybridization-based capture of pathogen mRNA enables paired host-pathogen transcriptional analysis

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

Time-resolved dual RNA-seq of tissue uncovers Pseudomonas plecoglossicida key virulence genes in host-pathogen interaction with Epinephelus coioides

Bacterial pathogen-host interactions are highly dynamic, regulated processes that have been primarily investigated using in vitro assays. The dynamics of bacterial pathogen-host interplay in vivo are poorly understood. Using time-resolved dual RNA-seq in a Pseudomonas plecoglossicida-Epinephelus coioides infection model, we observed that bacterial genes encoding classical virulence factors and host genes involved in immune regulation were dynamically expressed during infection. Using network inferencing, we were able to predict interspecies regulatory networks linking bacterial virulence genes to host immune genes. Together with gene co-expression network analysis of the pathogen, secY was predicted to be a key virulence gene for P. plecoglossicida pathogenicity in the host, fliN was predicted to be a less important virulence gene. The results of bioinformatics prediction were confirmed by animal infection experiments. Our work provides the first paradigm to study dynamic alterations of bacterial pathogen and host interactions based on the elucidation of time-resolved interactive transcriptomes in vivo, and may be developed into a novel and universal method for revealing the true complexity of the bacterial infection process.

Use of gene expression studies to investigate the human immunological response to malaria infection

Background

Transcriptional profiling of the human immune response to malaria has been used to identify diagnostic markers, understand the pathogenicity of severe disease and dissect the mechanisms of naturally acquired immunity (NAI). However, interpreting this body of work is difficult given considerable variation in study design, definition of disease, patient selection and methodology employed. This work details a comprehensive review of gene expression profiling (GEP) of the human immune response to malaria to determine how this technology has been applied to date, instances where this has advanced understanding of NAI and the extent of variability in methodology between studies to allow informed comparison of data and interpretation of results.

Methods

Datasets from the gene expression omnibus (GEO) including the search terms; ‘plasmodium’ or ‘malaria’ or ‘sporozoite’ or ‘merozoite’ or ‘gametocyte’ and ‘Homo sapiens’ were identified and publications analysed. Datasets of gene expression changes in relation to malaria vaccines were excluded.

Results

Twenty-three GEO datasets and 25 related publications were included in the final review. All datasets related to Plasmodium falciparum infection, except two that related to Plasmodium vivax infection. The majority of datasets included samples from individuals infected with malaria ‘naturally’ in the field (n = 13, 57%), however some related to controlled human malaria infection (CHMI) studies (n = 6, 26%), or cells stimulated with Plasmodium in vitro (n = 6, 26%). The majority of studies examined gene expression changes relating to the blood stage of the parasite. Significant heterogeneity between datasets was identified in terms of study design, sample type, platform used and method of analysis. Seven datasets specifically investigated transcriptional changes associated with NAI to malaria, with evidence supporting suppression of the innate pro-inflammatory response as an important mechanism for this in the majority of these studies. However, further interpretation of this body of work was limited by heterogeneity between studies and small sample sizes.

Conclusions

GEP in malaria is a potentially powerful tool, but to date studies have been hypothesis generating with small sample sizes and widely varying methodology. As CHMI studies are increasingly performed in endemic settings, there will be growing opportunity to use GEP to understand detailed time-course changes in host response and understand in greater detail the mechanisms of NAI.

Multiomics analyses reveal that NOD-like signaling pathway plays an important role against Streptococcus agalactiae in the spleen of tilapia

Streptococcus aglactiae(GBS) infection in tilapia is a serious global disease that causes significant production loss. Here, we studied the role of GBS in the spleen and the spleen's response against the pathogen through dual RNA-seq and proteome technology. Animals were divided into three groups: control, virulent treated (HN016), and attenuated treated (YM001). Spleen samples were collected and analysis when a disease outbreak. Dual RNA-seq result showed the virulence factor genes of GBS, included CAMP factor, PGK, OCT, enolase, scpB, Sip, bca, were upregulation. downregulation of GapA, cylE, OCT, scpB, C5AP, rlmB, hly, FBP, in HN016 and YM001. But for proteomic, OCT and bca were downregulation, the others were upregulation. For host transcriptome KEGG analysis showed, the NOD-like receptor signaling pathway (NLRs) and TOLL-like receptor signaling pathway (TLRs) were upreguoation in HN016 infected fish than the control fish; But for proteome KEGG, only the NLRS was up, the TLRS was not change. Compared with YM001 infected fishes, for transcriptome, NLRs and TLRs in infected HN016 fishes were significance rise (p < 0.01); for proteome, the NLRs was up (p < 0.05), but TLRs was no change.Analysis of pathogen-host interaction showed that the peptidoglycan (PNG), CD2, LCK, and host's Zap70 were involved in the regulation of NLRs; PNG, LCK, and ZAP70 were involved in the regulation of TRLs. Conclusion: the virulent strain HN016 and attenuated strainYM001 differed in the quantity of virulence factors. In tilapia's innate immune system, NLRs was the main defense factors, but bacteria avoided the host defense through TLRs.

Arabidopsis halleri: a perennial model system for studying population differentiation and local adaptation

Local adaptation is assumed to occur when populations differ in a phenotypic trait or a set of traits, and such variation has a genetic basis. Here, we introduce Arabidopsis halleri and its life history as a perennial model system to study population differentiation and local adaptation. Studies on altitudinal adaptation have been conducted in two regions: Mt. Ibuki in Japan and the European Alps. Several studies have demonstrated altitudinal adaptation in ultraviolet-B (UV-B) tolerance, leaf water repellency against spring frost and anti-herbivore defences. Studies on population differentiation in A. halleri have also focused on metal hyperaccumulation and tolerance to heavy metal contamination. In these study systems, genome scans to identify candidate genes under selection have been applied. Lastly, we briefly discuss how RNA-Seq can broaden phenotypic space and serve as a link to underlying mechanisms. In conclusion, A. halleri provides us with opportunities to study population differentiation and local adaptation, and relate these to the genetic systems underlying target functional traits.

Two novel transcriptional reporter systems for monitoring Helicobacter pylori stress responses

Helicobacter pylori, a human pathogen linked to many stomach diseases, is well adapted to colonize aggressive gastric environments, and its virulence factors contribute this adaptation. Here, we report the construction of two novel H. pylori vectors, pSv2 and pSv4, carrying a reporter gene fused to the promoters of virulence factor genes for monitoring the response of single H. pylori cells to various stresses. H. pylori cryptic plasmids were modified by the introduction of the Escherichia coli origin of replication, chloramphenicol resistance cassette, and promoterless gfp gene to produce E. coli/H. pylori shuttle vectors. The promoter regions of vacA and ureA genes encoding well-characterized H. pylori virulence factors were fused to the promoterless gfp gene. Recording the GFP fluorescence signal from the genetically modified H. pylori cells immobilized in specifically designed microfluidic devices revealed the response of transcriptional reporter systems to osmotic stress, acidic stress, elevated Ni²⁺ concentration or iron chelation. Our observations validate the utility of the pSv2 and pSv4 vectors to monitor the regulation of virulence factor genes in diverse strains and clinical isolates of H. pylori.

Dual RNA-Seq Unveils the Role of the Pseudomonas plecoglossicida fliA Gene in Pathogen-Host Interaction with Larimichthys crocea

In the present study, Larimichthys crocea and Pseudomonas plecoglossicida were selected as a host-pathogen interaction model for teleosts and prokaryotic pathogens. Five shRNAs were designed and synthesized to silence the fliA gene, all of which resulted in pronounced reductions in fliA mRNA; the mutant strain with the best silencing efficiency of 92.16% was chosen for subsequent analysis. A significant decrease in motility, intracellular survival and escape was observed for the fliA-RNAi strain of P. plecoglossicida, whereby silencing of the fliA gene led to a 30% decrease in mortality and a four-day delay in the onset of infection in L. crocea. Moreover, silencing of P. plecoglossicida fliA resulted in a significant change in both the pathogen and host transcriptome in the spleens of infected L. crocea. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of pathogen transcriptome data showed that silencing fliA resulted in downregulation of 18 flagellum-related genes; KEGG analysis of host transcriptome data revealed that infection with the fliA-RNAi strain caused upregulation of 47 and downregulation of 106 immune-related genes. These pathogen-host interactions might facilitate clearance of P. plecoglossicida by L. crocea, with a significant decrease in fliA-RNAi P. plecoglossicida strain virulence in L. crocea.

Dual RNA-seq uncovers the immune response of Larimichthys crocea to the secY gene of Pseudomonas plecoglossicida from the perspective of host-pathogen interactions

Pseudomonas plecoglossicida is a Gram-negative aerobic bacterium that causes high mortality and serious economic losses in some commercial marine fish. Expression of secY was found to be significantly upregulated at 18 °C compared to 28 °C by RNA-seq and qRT-PCR. All five tested recombinant vectors (pCM130/tac + shRNA) significantly reduced secY mRNA levels in P. plecoglossicida. The recombinant vector encoding shRNA-1165 exhibited the best gene-silencing efficiency, 82.4% and was used to create an RNAi strain for further studies. Compared with the wildtype strain, infections of Larimichthys crocea with the RNAi strain resulted in a 2-day delay in onset time and a 35% reduction in mortality, as well as the alleviation of spleen symptoms. The spleens of L. crocea infected by the wild type or RNAi strain of P. plecoglossicida were subjected to dual RNA-seq at 2 dpi. Compared with the wildtype strain, infection of P. plecoglossicida with the RNAi strain resulted in significant changes in the transcriptomes of both host and pathogen. KEGG analysis showed that the complement and coagulation cascade and the Toll-like receptor signalling pathway were the most enriched host pathways. In the pathogen, genes of the "Sec secretion system" were significantly downregulated. This downregulation of "Sec secretion system" genes hindered the secretion of bacterial proteins and reduced the virulence of P. plecoglossicida. Thus, it was easier for L. crocea to clear the RNAi strain of P. plecoglossicida, and the immune response was similarly reduced. The results indicated that secY was a virulence gene of P. plecoglossicida and played roles in the host-pathogen interactions of L. crocea and P. plecoglossicida.

Mechanistic insight into the roles of Pseudomonas plecoglossicida clpV gene in host-pathogen interactions with Larimichthys crocea by dual RNA-seq

Large yellow croaker (Larimichthys crocea) is an economical important farmed fish in China. "Visceral White Spot Disease" caused by Pseudomonas plecoglossicida is a disease with a high mortality rate in cage-cultured L. crocea in recent years and resulted in heavy economy lossess. The dual RNA-seq results of previous study showed that the expression of clpV gene in P. plecoglossicida was significantly up-regulated during infection. RNAi significantly reduced the expression of clpV in P. plecoglossicida with maximum silencing efficiency of 96.1%. Compared with the wild type strain, infection of clpV-RNAi strain resulted in a delayed onset time and a 25% reduction in mortality of L. crocea, as well as lessening the symptoms of the spleen. The results of dual RNA-seq of L. crocea infected by clpV-RNAi strain of P. plecoglossicida changed considerably, compared with the counterpart infected with the wild strain. The KEGG enrichment analysis showed that Cytokine-cytokine receptor interaction, Toll-like receptor signaling pathway, C-type lectin receptor signaling pathway and MAPK signaling pathway of L. crocea were most affected by the silence of clpV in P. plecoglossicida. RNAi of clpV resulted in the downregulation of genes in flagella assembly pathway and a weaker immune response of host.

Transcriptomic profiling of interacting nasal staphylococci species reveals global changes in gene and non-coding RNA expression

Interspecies interactions between bacterial pathogens and the commensal microbiota can influence disease outcome. In the nasal cavities, Staphylococcus epidermidis has been shown to be a determining factor for Staphylococcus aureus colonization and biofilm formation. However, the interaction between S. epidermidis and S. aureus has mainly been described by phenotypic analysis, and little is known about how this interaction modulates gene expression. This study aimed to determine the interactome of nasal S. aureus and S. epidermidis isolates to understand the molecular effect of interaction. After whole-genome sequencing of two nasal staphylococcal isolates, an agar-based RNA sequencing setup was utilized to identify interaction-induced transcriptional alterations in surface-associated populations. Our results revealed differential expression of several virulence genes in both species. We also identified putative non-coding RNAs (ncRNAs) and, interestingly, detected a putative ncRNA transcribed antisense to esp, the serine protease of S. epidermidis, that has previously been shown to inhibit nasal colonization of S. aureus. In our study, the gene encoding Esp and the antisense ncRNA are both downregulated during interaction with S. aureus. Our findings contribute to a better understanding of pathogen physiology in the context of interactions with the commensal microbiota, and may provide targets for future therapeutics.

The cellular microbiology of Salmonellae interactions with macrophages

Salmonellae are important enteric pathogens that cause gastroenteritis and systemic illnesses. Macrophages are important components of both the innate and acquired immune system, acting as phagocytes with significant antimicrobial killing activities that present antigen to the adaptive immune system. Macrophages can also be cultured from a variety of sites as primary cells, and the study of the survival and interactions of Salmonellae with these cells is a very early model of infection and cellular microbiology. This review traces the history of discoveries made using Salmonellae infection of macrophages and addresses the possibility of future research in this area, in particular with regards to understanding the complexity of individual bacteria and macrophage cell variability and how such heterogeneity may alter the outcome of infection.

The Small RNA PinT Contributes to PhoP-Mediated Regulation of the Salmonella Pathogenicity Island 1 Type III Secretion System in Salmonella enterica Serovar Typhimurium

Salmonella enterica serovar Typhimurium induces inflammatory diarrhea and bacterial uptake into intestinal epithelial cells using the Salmonella pathogenicity island 1 (SPI1) type III secretion system (T3SS). HilA activates transcription of the SPI1 structural components and effector proteins. Expression of hilA is activated by HilD, HilC, and RtsA, which act in a complex feed-forward regulatory loop. Many environmental signals and other regulators are integrated into this regulatory loop, primarily via HilD. After the invasion of Salmonella into host intestinal epithelial cells or during systemic replication in macrophages, the SPI T3SS is no longer required or expressed. We have shown that the two-component regulatory system PhoPQ, required for intracellular survival, represses the SPI1 T3SS mostly by controlling the transcription of hilA and hilD Here we show that PinT, one of the PhoPQ-regulated small RNAs (sRNAs), contributes to this regulation by repressing hilA and rtsA translation. PinT base pairs with both the hilA and rtsA mRNAs, resulting in translational inhibition of hilA, but also induces degradation of the rts transcript. PinT also indirectly represses expression of FliZ, a posttranslational regulator of HilD, and directly represses translation of ssrB, encoding the primary regulator of the SPI2 T3SS. Our in vivo mouse competition assays support the concept that PinT controls a series of virulence genes at the posttranscriptional level in order to adapt Salmonella from the invasion stage to intracellular survival.IMPORTANCE Salmonella is one of the most important food-borne pathogens, infecting over one million people in the United States every year. These bacteria use a needle-like device to interact with intestinal epithelial cells, leading to invasion of the cells and induction of inflammatory diarrhea. A complex regulatory network controls expression of the invasion system in response to numerous environmental signals. Here we explore the molecular mechanisms by which the small RNA PinT contributes to this regulation, facilitating inactivation of the system after invasion. PinT controls several important virulence systems in Salmonella, tuning the transition between different stages of infection.

Dual RNA-Seq uncovers the function of an ABC transporter gene in the host-pathogen interaction between Epinephelus coioides and Pseudomonas plecoglossicida

As an important pathogen in aquaculture, Pseudomonas plecoglossicida has caused heavy losses. The expression of an ABC transporter gene-L321_23611 of P. plecoglossicida at 18 °C was found significant higher than those at 28 °C by RNA-seq and qRT-PCR. RNAi significantly reduced the content of L321_23611 mRNA in P. plecoglossicida with a maximal decrease of 89.2%. Compared with the wild type strain, the infection of L321_23611-RNAi strain resulted in the reduction in mortality and the onset time delay of a kind of marine teleosts, Epinephelus coioides. The results of dual RNA-seq showed that the RNAi of L321_23611 resulted in a significant change in both pathogen and host transcriptome in the spleens of infected E. coioides. The result of GO and KEGG analysis from dual RNA-seq data showed both host genes of chemokine signaling pathway, coagulation and complement system, hematopoietic cell lineage pathway as well as hemoglobin complex GO term and pathogenic genes of bacterial-type flagellum-dependent cell mortality GO term and flagellar assembly, biosynthesis of amino acids and lysine biosynthesis systems pathways were mainly affected by L321_23611 gene of P. plecoglossicida. The results indicated that: 1. ABC transporter gene-L321_23611 was a virulent gene of P. plecoglossicida. 2. Both the activation of the host immune pathways and depression of pathogenic virulence-related pathways facilitated E. coioides to remove L321_23611-RNAi strain than the wild type strain of P. plecoglossicida.

Functional expansion of a TCA cycle operon mRNA by a 3' end-derived small RNA

Global RNA profiling studies in bacteria have predicted the existence of many of small noncoding RNAs (sRNAs) that are processed off mRNA 3′ ends to regulate other mRNAs via the RNA chaperones Hfq and ProQ. Here, we present targets of SdhX (RybD), an Hfq-dependent sRNA that is generated by RNase E mediated 3′ processing of the ∼10 000-nt mRNA of the TCA cycle operon sdhCDAB-sucABCD in enteric bacteria. An in silico search predicted ackA mRNA, which encodes acetate kinase, as a conserved primary target of SdhX. Through base pairing, SdhX represses AckA synthesis during growth of Salmonella on acetate. Repression can be achieved by a naturally occurring 38-nucleotide SdhX variant, revealing the shortest functional Hfq-associated sRNA yet. Salmonella SdhX also targets the mRNAs of fumB (anaerobic fumarase) and yfbV, a gene of unknown function adjacent to ackA. Instead, through a slightly different seed sequence, SdhX can repress other targets in Escherichia coli, namely katG (catalase) and fdoG (aerobic formate dehydrogenase). This study illustrates how a key operon from central metabolism is functionally connected to other metabolic pathways through a 3′ appended sRNA, and supports the notion that mRNA 3′UTRs are a playground for the evolution of regulatory RNA networks in bacteria.

PhoP-Mediated Repression of the SPI1 Type 3 Secretion System in Salmonella enterica Serovar Typhimurium

Salmonella must rapidly adapt to various niches in the host during infection. Relevant virulence factors must be appropriately induced, and systems that are detrimental in a particular environment must be turned off. Salmonella infects intestinal epithelial cells using a type 3 secretion system (T3SS) encoded on Salmonella pathogenicity island 1 (SPI1). The system is controlled by three AraC-like regulators, HilD, HilC, and RtsA, which form a complex feed-forward loop to activate expression of hilA, encoding the main transcriptional regulator of T3SS structural genes. This system is tightly regulated, with many of the activating signals acting at the level of hilD translation or HilD protein activity. Once inside the phagosomes of epithelial cells, or in macrophages during systemic stages of disease, the SPI1 T3SS is no longer required or expressed. Here, we show that the PhoPQ two-component system, critical for intracellular survival, appears to be the primary mechanism by which Salmonella shuts down the SPI1 T3SS. PhoP negatively regulates hilA through multiple distinct mechanisms: direct transcriptional repression of the hilA promoter, indirect transcriptional repression of both the hilD and rtsA promoters, and activation of the small RNA (sRNA) PinT. Genetic analyses and electrophoretic mobility shift assays suggest that PhoP specifically binds the hilA promoter to block binding of activators HilD, HilC, and RtsA as a mechanism of repression.IMPORTANCE Salmonella is one of the most common foodborne pathogens, causing an estimated 1.2 million illnesses per year in the United States. A key step in infection is the activation of the bacterial invasion machinery, which induces uptake of the bacterium into epithelial cells and leads to induction of inflammatory diarrhea. Upon entering the vacuolar compartments of host cells, Salmonella senses an environmental transition and represses the invasion machinery with a two-component system relevant for survival within the vacuole. This adaptation to specific host niches is an important example of how signals are integrated for survival of the pathogen.

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.

The effect of a LysR-type transcriptional regulator gene of Pseudomonas plecoglossicida on the immune responses of Epinephelus coioides

As an important pathogen in aquaculture, Pseudomonas plecoglossicida has caused heavy losses. It was determined with RNA-seq that the expression of a LysR-type transcriptional regulator gene (L321_20267) of P. plecoglossicida at 18 °C was significantly higher than that at 28 °C, which was verified by quantitative real-time PCR (qRT-PCR). RNAi significantly reduced the content of L321_20267 mRNA in P. plecoglossicida, with a maximal decrease of 90.63%. Compared with the wild-type strain, infection with the L321_20267-RNAi strain resulted in a 50% reduction in mortality and an onset time delay of Epinephelus coioides, as well as alleviation of the symptoms in E. coioides spleens. Compared with the wild-type strain of P. plecoglossicida, the L321_20267-RNAi strain resulted in a significant change in the spleen transcriptome of infected E. coioides. The results of GO and KEGG analysis showed that genes of serine hydrolase activity, the antigen processing and presentation pathway, the B cell receptor signalling pathway and the chemokine signalling pathway were most affected by the L321_20267 gene of P. plecoglossicida. Meanwhile, the immune genes were related to different numbers of miRNAs and lncRNAs, and some miRNAs were related to more than one gene. The results indicated that 1. L321_20267 is a virulence gene of P. plecoglossicida; 2. the upregulation of the immune pathways facilitated E. coioides to remove the L321_20267-RNAi strain compared with the wild-type strain of P. plecoglossicida; and 3. the immune genes were regulated by miRNA and lncRNA in a complex manner.

Induced Pluripotent Stem Cell-Derived Brain Endothelial Cells as a Cellular Model to Study Neisseria meningitidis Infection

Meningococcal meningitis is a severe central nervous system infection that occurs when Neisseria meningitidis (Nm) penetrates brain endothelial cells (BECs) of the meningeal blood-cerebrospinal fluid barrier. As a human-specific pathogen, in vivo models are greatly limited and pose a significant challenge. In vitro cell models have been developed, however, most lack critical BEC phenotypes limiting their usefulness. Human BECs generated from induced pluripotent stem cells (iPSCs) retain BEC properties and offer the prospect of modeling the human-specific Nm interaction with BECs. Here, we exploit iPSC-BECs as a novel cellular model to study Nm host-pathogen interactions, and provide an overview of host responses to Nm infection. Using iPSC-BECs, we first confirmed that multiple Nm strains and mutants follow similar phenotypes to previously described models. The recruitment of the recently published pilus adhesin receptor CD147 underneath meningococcal microcolonies could be verified in iPSC-BECs. Nm was also observed to significantly increase the expression of pro-inflammatory and neutrophil-specific chemokines IL6, CXCL1, CXCL2, CXCL8, and CCL20, and the secretion of IFN-γ and RANTES. For the first time, we directly observe that Nm disrupts the three tight junction proteins ZO-1, Occludin, and Claudin-5, which become frayed and/or discontinuous in BECs upon Nm challenge. In accordance with tight junction loss, a sharp loss in trans-endothelial electrical resistance, and an increase in sodium fluorescein permeability and in bacterial transmigration, was observed. Finally, we established RNA-Seq of sorted, infected iPSC-BECs, providing expression data of Nm-responsive host genes. Altogether, this model provides novel insights into Nm pathogenesis, including an impact of Nm on barrier properties and tight junction complexes, and suggests that the paracellular route may contribute to Nm traversal of BECs.

Ro60 and Y RNAs: structure, functions, and roles in autoimmunity

Ro60, also known as SS-A or TROVE2, is an evolutionarily conserved RNA-binding protein that is found in most animal cells, approximately 5% of sequenced prokaryotic genomes and some archaea. Ro60 is present in cells as both a free protein and as a component of a ribonucleoprotein complex, where its best-known partners are members of a class of noncoding RNAs called Y RNAs. Structural and biochemical analyses have revealed that Ro60 is a ring-shaped protein that binds Y RNAs on its outer surface. In addition to Y RNAs, Ro60 binds misfolded and aberrant noncoding RNAs in some animal cell nuclei. Although the fate of these defective Ro60-bound noncoding RNAs in animal cells is not well-defined, a bacterial Ro60 ortholog functions with 3' to 5' exoribonucleases to assist structured RNA degradation. Studies of Y RNAs have revealed that these RNAs regulate the subcellular localization of Ro60, tether Ro60 to effector proteins and regulate the access of other RNAs to its central cavity. As both mammalian cells and bacteria lacking Ro60 are sensitized to ultraviolet irradiation, Ro60 function may be important during exposure to some environmental stressors. Here we summarize the current knowledge regarding the functions of Ro60 and Y RNAs in animal cells and bacteria. Because the Ro60 RNP is a clinically important target of autoantibodies in patients with rheumatic diseases such as Sjogren's syndrome, systemic lupus erythematosus, and neonatal lupus, we also discuss potential roles for Ro60 RNPs in the initiation and pathogenesis of systemic autoimmune rheumatic disease.

Uncovering complex molecular networks in host-pathogen interactions using systems biology

Interactions between pathogens and their hosts can induce complex changes in both host and pathogen states to privilege pathogen survival or host clearance of the pathogen. To determine the consequences of specific host-pathogen interactions, a variety of techniques in microbiology, cell biology, and immunology are available to researchers. Systems biology that enables unbiased measurements of transcriptomes, proteomes, and other biomolecules has become increasingly common in the study of host-pathogen interactions. These approaches can be used to generate novel hypotheses or to characterize the effects of particular perturbations across an entire biomolecular network. With proper experimental design and complementary data analysis tools, high-throughput omics techniques can provide novel insights into the mechanisms that underlie processes from phagocytosis to pathogen immune evasion. Here, we provide an overview of the suite of biochemical approaches for high-throughput analyses of host-pathogen interactions, analytical frameworks for understanding the resulting datasets, and a vision for the future of this exciting field.

The Interplay between Salmonella enterica Serovar Typhimurium and the Intestinal Mucosa during Oral Infection

Bacterial infection results in a dynamic interplay between the pathogen and its host. The underlying interactions are multilayered, and the cellular responses are modulated by the local environment. The intestine is a particularly interesting tissue regarding host-pathogen interaction. It is densely colonized by commensal microbes and a portal of entry for ingested pathogens. This necessitates constant monitoring of microbial stimuli in order to maintain homeostasis during encounters with benign microbiota and to trigger immune defenses in response to bacterial pathogens. Homeostasis is maintained by physical barriers (the mucus layer and epithelium), chemical defenses (antimicrobial peptides), and innate immune responses (NLRC4 inflammasome), which keep the bacteria from reaching the sterile lamina propria. Intestinal pathogens represent potent experimental tools to probe these barriers and decipher how pathogens can circumvent them. The streptomycin mouse model of oral Salmonella enterica serovar Typhimurium infection provides a well-characterized, robust experimental system for such studies. Strikingly, each stage of the gut tissue infection poses a different set of challenges to the pathogen and requires tight control of virulence factor expression, host response modulation, and cooperation between phenotypic subpopulations. Therefore, successful infection of the intestinal tissue relies on a delicate and dynamic balance between responses of the pathogen and its host. These mechanisms can be deciphered to their full extent only in realistic in vivo infection models.

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.

Immunometabolism at the interface between macrophages and pathogens

It is generally regarded that the progression of an infection within host macrophages is the consequence of a failed immune response. However, recent appreciation of macrophage heterogeneity, with respect to both development and metabolism, indicates that the reality is more complex. Different lineages of tissue-resident macrophages respond divergently to microbial, environmental and immunological stimuli. The emerging picture that the developmental origin of macrophages determines their responses to immune stimulation and to infection stresses the importance of in vivo infection models. Recent investigations into the metabolism of infecting microorganisms and host macrophages indicate that their metabolic interface can be a major determinant of pathogen growth or containment. This Review focuses on the integration of data from existing studies, the identification of challenges in generating and interpreting data from ongoing studies and a discussion of the technologies and tools that are required to best address future questions in the field.

Transcriptomic basis and evolution of the ant nurse-larval social interactome

Development is often strongly regulated by interactions among close relatives, but the underlying molecular mechanisms are largely unknown. In eusocial insects, interactions between caregiving worker nurses and larvae regulate larval development and resultant adult phenotypes. Here, we begin to characterize the social interactome regulating ant larval development by collecting and sequencing the transcriptomes of interacting nurses and larvae across time. We find that the majority of nurse and larval transcriptomes exhibit parallel expression dynamics across larval development. We leverage this widespread nurse-larva gene co-expression to infer putative social gene regulatory networks acting between nurses and larvae. Genes with the strongest inferred social effects tend to be peripheral elements of within-tissue regulatory networks and are often known to encode secreted proteins. This includes interesting candidates such as the nurse-expressed giant-lens, which may influence larval epidermal growth factor signaling, a pathway known to influence various aspects of insect development. Finally, we find that genes with the strongest signatures of social regulation tend to experience relaxed selective constraint and are evolutionarily young. Overall, our study provides a first glimpse into the molecular and evolutionary features of the social mechanisms that regulate all aspects of social life.

Transcriptomic basis and evolution of the ant nurse-larval social interactome

Development is often strongly regulated by interactions among close relatives, but the underlying molecular mechanisms are largely unknown. In eusocial insects, interactions between caregiving worker nurses and larvae regulate larval development and resultant adult phenotypes. Here, we begin to characterize the social interactome regulating ant larval development by collecting and sequencing the transcriptomes of interacting nurses and larvae across time. We find that the majority of nurse and larval transcriptomes exhibit parallel expression dynamics across larval development. We leverage this widespread nurse-larva gene co-expression to infer putative social gene regulatory networks acting between nurses and larvae. Genes with the strongest inferred social effects tend to be peripheral elements of within-tissue regulatory networks and are often known to encode secreted proteins. This includes interesting candidates such as the nurse-expressed giant-lens, which may influence larval epidermal growth factor signaling, a pathway known to influence various aspects of insect development. Finally, we find that genes with the strongest signatures of social regulation tend to experience relaxed selective constraint and are evolutionarily young. Overall, our study provides a first glimpse into the molecular and evolutionary features of the social mechanisms that regulate all aspects of social life.

Integrative Analyses of Long Non-coding RNA and mRNA Involved in Piglet Ileum Immune Response to Clostridium perfringens Type C Infection

Long non-coding RNAs (lncRNAs) have been shown to play important roles in regulating host immune and inflammatory responses to bacterial infection. Infection with Clostridium perfringens (C. perfringens), a food-borne zoonotic pathogen, can lead to a series of inflammatory diseases in human and piglet, greatly challenging the healthy development of global pig industry. However, the roles of lncRNAs involved in piglet immune response against C. perfringens type C infection remain unknown. In this study, the regulatory functions of ileum lncRNAs and mRNAs were investigated in piglet immune response to C. perfringens type C infection among resistance (IR), susceptibility (IS) and sham-inoculation (control, IC) groups. A total of 480 lncRNAs and 3,669 mRNAs were significantly differentially expressed, the differentially expressed lncRNAs and mRNAs in the IR and IS groups were enriched in various pathways of ABC transporters, olfactory transduction, PPAR signaling pathway, chemokine signaling pathway and Toll-like receptor signaling pathway, involving in regulating piglet immune responses and resistance during infection. There were 212 lncRNAs and 505 target mRNAs found to have important association with C. perfringens infectious diseases, furthermore, 25 dysregulated lncRNAs corresponding to 13 immune-related target mRNAs were identified to play potential roles in defense against bacterial infection. In conclusion, the results improve our understanding on the characteristics of lncRNAs and mRNAs on regulating host immune response against C. perfringens type C infection, which will provide a reference for future research into exploring C. perfringens-related diseases in human.

Clinical significance of lnc-AC145676.2.1-6 and lnc-TGS1-1 and their variants in western Chinese tuberculosis patients

OBJECTIVES

Tuberculosis (TB) remains a global public health problem. Recent studies have implicated long non-coding RNAs and their variants as possibly playing important roles in TB. The aim of this study was to assess the clinical relevance of lnc-AC145676.2.1-6 and lnc-TGS1-1 and their variants in a western Chinese population.

METHODS

This case-control study included 467 TB patients and 473 healthy controls from West China Hospital. The expression levels of lnc-AC145676.2.1-6 and lnc-TGS1-1 were analyzed by reverse transcriptase quantitative real-time PCR. Single-nucleotide polymorphism genotyping was performed using a custom-designed 2 × 48-Plex SNPscan kit.

RESULTS

It was observed that lnc-AC145676.2.1-6 and lnc-TGS1-1 expression levels were both obviously down-regulated in TB patients. In addition, a lower expression level of lnc-TGS1-1 was associated with the presence of thrombocytopenia in TB patients during anti-TB treatment, and the homozygous CC genotype of rs4737420 correlated with a decreased risk of leukopenia, compared with individuals with the T allele (TT/CT genotype), in the dominant mode.

CONCLUSIONS

For the first time, potential TB-associated promoting effects were identified for the decreased expression levels of lnc-AC145676.2.1-6 and lnc-TGS-1, while lnc-TGS1-1 and its variant rs4737420 may be predictive indicators of anti-TB drug-induced adverse drug reactions. Larger validation studies on different populations are warranted to confirm these findings.