How to Predict Molecular Interactions between Species?

Bioinformatic and experimental methods to identify and validate bacterial RNA-human RNA interactions

Ample evidence supports the importance of the microbiota on human health and disease. Recent studies suggest that extracellular vesicles are an important means of bacterial-host communication, in part via the transport of small RNAs (sRNAs). Bacterial sRNAs have been shown to co-precipitate with human and mouse RNA-induced silencing complex, hinting that some may regulate gene expression as eukaryotic microRNAs do. Bioinformatic tools, including those that can incorporate an sRNA's secondary structure, can be used to predict interactions between bacterial sRNAs and human messenger RNAs (mRNAs). Validation of these potential interactions using reproducible experimental methods is essential to move the field forward. This review will cover the evidence of interspecies communication via sRNAs, bioinformatic tools currently available to identify potential bacterial sRNA-host (specifically, human) mRNA interactions, and experimental methods to identify and validate those interactions.

Mathematical models to study the biology of pathogens and the infectious diseases they cause

Mathematical models have many applications in infectious diseases: epidemiologists use them to forecast outbreaks and design containment strategies; systems biologists use them to study complex processes sustaining pathogens, from the metabolic networks empowering microbial cells to ecological networks in the microbiome that protects its host. Here, we (1) review important models relevant to infectious diseases, (2) draw parallels among models ranging widely in scale. We end by discussing a minimal set of information for a model to promote its use by others and to enable predictions that help us better fight pathogens and the diseases they cause.

Unconventional routes to developing insect-resistant crops

Concerns over widespread use of insecticides and heightened insect pest virulence under climate change continue to fuel the need for environmentally safe and sustainable control strategies. However, to develop such strategies, a better understanding of the molecular basis of plant-pest interactions is still needed. Despite decades of research investigating plant-insect interactions, few examples exist where underlying molecular mechanisms are well characterized, and even rarer are cases where this knowledge has been successfully applied to manage harmful agricultural pests. Consequently, the field appears to be static, urgently needing shifts in approaches to identify novel mechanisms by which insects colonize plants and plants avoid insect pressure. In this perspective, we outline necessary steps for advancing holistic methodologies that capture complex plant-insect molecular interactions. We highlight novel and underexploited approaches in plant-insect interaction research as essential routes to translate knowledge of underlying molecular mechanisms into durable pest control strategies, including embracing microbial partnerships, identifying what makes a plant an unsuitable host, capitalizing on tolerance of insect damage, and learning from cases where crop domestication and agronomic practices enhance pest virulence.

The Zinc Nutritional Immunity of Epinephelus coioides Contributes to the Importance of znuC During Pseudomonas plecoglossicida Infection

Previously, the dual RNA-seq was carried out in a Pseudomonas plecoglossicida- Epinephelus coioides infection model to investigate the dynamics of pathogen-host interplay in vivo. ZnuC, a member of ZnuCBA Zn importer, was found transcriptionally up-regulated during infection. Thus, this study aimed to assess its role during the trade-off for Zn between host and P. plecoglossicida. ICP-MS analysis and fluorescent staining showed that Zn was withheld from serum and accumulated in the spleen, with increased Zn uptake in the Golgi apparatus of macrophages after infection. Additionally, growth assay, macrophage infection and animal infection after gene knockout / silencing revealed that znuC was necessary for growth in Zn-limiting conditions, colonization, intracellular viability, immune escape and virulence of P. plecoglossicida. Further analysis with dual RNA-seq revealed associations of host's Zn nutritional immunity genes with bacterial Zn assimilation genes. IL6 and ZIP4 played key roles in this network, and markedly affected znuB expression, intracellular viability and immune escape, as revealed by gene silencing. Moreover, EMSA and GFP reporter gene analysis showed that Fur sensed changes in Fe concentration to regulate znuCBA in P. plecoglossicida. Jointly, these findings suggest a trade-off for Zn between host and P. plecoglossicida, while ZnuC is important for P. plecoglossicida Zn acquisition.

Know your enemy - transcriptome of myxozoan Tetracapsuloides bryosalmonae reveals potential drug targets against proliferative kidney disease in salmonids

The myxozoan Tetracapsuloides bryosalmonae is a widely spread endoparasite that causes proliferative kidney disease (PKD) in salmonid fish. We developed an in silico pipeline to separate transcripts of T. bryosalmonae from the kidney tissue of its natural vertebrate host, brown trout (Salmo trutta). After stringent filtering, we constructed a partial transcriptome assembly T. bryosalmonae, comprising 3427 transcripts. Based on homology-restricted searches of the assembled parasite transcriptome and Atlantic salmon (Salmo salar) proteome, we identified four protein targets (Endoglycoceramidase, Legumain-like protease, Carbonic anhydrase 2, Pancreatic lipase-related protein 2) for the development of anti-parasitic drugs against T. bryosalmonae. Earlier work of these proteins on parasitic protists and helminths suggests that the identified anti-parasitic drug targets represent promising chemotherapeutic candidates also against T. bryosalmonae, and strengthen the view that the known inhibitors can be effective in evolutionarily distant organisms. In addition, we identified differentially expressed T. bryosalmonae genes between moderately and severely infected fish, indicating an increased abundance of T. bryosalmonae sporogonic stages in fish with low parasite load. In conclusion, this study paves the way for future genomic research in T. bryosalmonae and represents an important step towards the development of effective drugs against PKD.

Humic-acid-driven escape from eye parasites revealed by RNA-seq and target-specific metabarcoding

Background

Next generation sequencing (NGS) technologies are extensively used to dissect the molecular mechanisms of host-parasite interactions in human pathogens. However, ecological studies have yet to fully exploit the power of NGS as a rich source for formulating and testing new hypotheses.

Methods

We studied Eurasian perch (Perca fluviatilis) and its eye parasite (Trematoda, Diplostomidae) communities in 14 lakes that differed in humic content in order to explore host-parasite-environment interactions. We hypothesised that high humic content along with low pH would decrease the abundance of the intermediate hosts (gastropods), thus limiting the occurrence of diplostomid parasites in humic lakes. This hypothesis was initially invoked by whole eye RNA-seq data analysis and subsequently tested using PCR-based detection and a novel targeted metabarcoding approach.

Results

Whole eye transcriptome results revealed overexpression of immune-related genes and the presence of eye parasite sequences in RNA-seq data obtained from perch living in clear-water lakes. Both PCR-based and targeted-metabarcoding approach showed that perch from humic lakes were completely free from diplostomid parasites, while the prevalence of eye flukes in clear-water lakes that contain low amounts of humic substances was close to 100%, with the majority of NGS reads assigned to Tylodelphys clavata.

Conclusions

High intraspecific diversity of T. clavata indicates that massively parallel sequencing of naturally pooled samples represents an efficient and powerful strategy for shedding light on cryptic diversity of eye parasites. Our results demonstrate that perch populations in clear-water lakes experience contrasting eye parasite pressure compared to those from humic lakes, which is reflected by prevalent differences in the expression of immune-related genes in the eye. This study highlights the utility of NGS to discover novel host-parasite-environment interactions and provide unprecedented power to characterize the molecular diversity of cryptic parasites.

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.

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.

Coculturing Bacteria Leads to Reduced Phenotypic Heterogeneities

Isogenic bacterial populations are known to exhibit phenotypic heterogeneity at the single-cell level. Because of difficulties in assessing the phenotypic heterogeneity of a single taxon in a mixed community, the importance of this deeper level of organization remains relatively unknown for natural communities. In this study, we have used membrane-based microcosms that allow the probing of the phenotypic heterogeneity of a single taxon while interacting with a synthetic or natural community. Individual taxa were studied under axenic conditions, as members of a coculture with physical separation, and as a mixed culture. Phenotypic heterogeneity was assessed through both flow cytometry and Raman spectroscopy. Using this setup, we investigated the effect of microbial interactions on the individual phenotypic heterogeneities of two interacting drinking water isolates. Through flow cytometry we have demonstrated that interactions between these bacteria lead to a reduction of their individual phenotypic diversities and that this adjustment is conditional on the bacterial taxon. Single-cell Raman spectroscopy confirmed a taxon-dependent phenotypic shift due to the interaction. In conclusion, our data suggest that bacterial interactions may be a general driver of phenotypic heterogeneity in mixed microbial populations.IMPORTANCE Laboratory studies have shown the impact of phenotypic heterogeneity on the survival and functionality of isogenic populations. Because phenotypic heterogeneity plays an important role in pathogenicity and virulence, antibiotic resistance, biotechnological applications, and ecosystem properties, it is crucial to understand its influencing factors. An unanswered question is whether bacteria in mixed communities influence the phenotypic heterogeneity of their community partners. We found that coculturing bacteria leads to a reduction in their individual phenotypic heterogeneities, which led us to the hypothesis that the individual phenotypic diversity of a taxon is dependent on the community composition.

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.

Dual RNA-Seq reveals the role of a transcriptional regulator gene in pathogen-host interactions between Pseudomonas plecoglossicida and Epinephelus coioides

Pseudomonas plecoglossicida is a highly pathogenic bacterium for maricultured fish and causes serious losses. A transcriptional regulator gene RK21_RS10315 was found up-regulated during the whole infection process, which was confirmed by qRT-PCR. Five shRNA were designed to silence RK21_RS10315 gene, and the gene expression was reduced up to 96.1%. Compared with the counterpart infected with wild type strain, the infection of RK21_RS10315-RNAi strain resulted in the death time delay, and 90% reduction in mortality of Epinephelus coioides, as well as the alleviation in the symptoms of E. coioides spleen. Moreover, compared with the fish infected with wild type strain, the infection of RK21_RS10315-RNAi strain of P. plecoglossicida resulted in a significant change both in transcriptome of spleen of infected E. coioides and P. plecoglossicida. The KEGG analysis showed that genes of 16 immune pathways in E. coioides were affected by the silence of RK21_RS10315 of P. plecoglossicida. Among them, intestinal immune network for IgA production pathway and leukocyte transendothelial migration pathway were more prominent than other pathways. 19 euk-DEMs in these immune pathways had varying degrees of correlation with 19 pro-DEMs, and the expression of ipxA, grpE, yhbJ, truD and suhB from 19 pro-DEMs were predicted more related to RK21_RS10315 in P. plecoglossicida.

Dual Analysis of Virus-Host Interactions: The Case of Ostreid herpesvirus 1 and the Cupped Oyster Crassostrea gigas

Dual analyses of the interactions between Ostreid herpesvirus 1 (OsHV-1) and the bivalve Crassostrea gigas during infection can unveil events critical to the onset and progression of this viral disease and can provide novel strategies for mitigating and preventing oyster mortality. Among the currently used “omics” technologies, dual transcriptomics (dual RNA-seq) coupled with the analysis of viral DNA in the host tissues has greatly advanced the knowledge of genes and pathways mostly contributing to host defense responses, expression profiles of annotated and unknown OsHV-1 open reading frames (ORFs), and viral genome variability. In addition to dual RNA-seq, proteomics and metabolomics analyses have the potential to add complementary information, needed to understand how a malacoherpesvirus can redirect and exploit the vital processes of its host. This review explores our current knowledge of “omics” technologies in the study of host-pathogen interactions and highlights relevant applications of these fields of expertise to the complex case of C gigas infections by OsHV-1, which currently threaten the mollusk production sector worldwide.

Two's company: studying interspecies relationships with dual RNA-seq

Organisms do not exist isolated from each other, but constantly interact. Cells can sense the presence of interaction partners by a range of receptors and, via complex regulatory networks, specifically react by changing the expression of many of their genes. Technological advances in next-generation sequencing over the recent years now allow us to apply RNA sequencing to two species at the same time (dual RNA-seq), and thus to directly study the gene expression of two interacting species without the need to physically separate cells or RNA. In this review, we give an overview over the latest studies in interspecies interactions made possible by dual RNA-seq, ranging from pathogenic to symbiotic relationships. We summarize state-of-the-art experimental techniques, bioinformatic data analysis and data interpretation, while also highlighting potential problems and pitfalls starting from the selection of meaningful time points and number of reads to matters of rRNA depletion. A short outlook on new trends in the field of dual RNA-seq concludes this review, looking at sequencing of non-coding RNAs during host-pathogen interactions and the prediction of molecular interspecies interactions networks.

Pseudogymnoascus destructans transcriptome changes during white-nose syndrome infections

White nose syndrome (WNS) is caused by the psychrophilic fungus Pseudogymnoascus destructans that can grow in the environment saprotrophically or parasitically by infecting hibernating bats. Infections are pathological in many species of North American bats, disrupting hibernation and causing mortality. To determine what fungal pathways are involved in infection of living tissue, we examined fungal gene expression using RNA-Seq. We compared P. destructans gene expression when grown in culture to that during infection of a North American bat species, Myotis lucifugus, that shows high WNS mortality. Cultured P. destructans was grown at 10 to 14 C and P. destructans growing in vivo was presumably exposed to temperatures ranging from 4 to 8 C during torpor and up to 37 C during periodic arousals. We found that when P. destructans is causing WNS, the most significant differentially expressed genes were involved in heat shock responses, cell wall remodeling, and micronutrient acquisition. These results indicate that this fungal pathogen responds to host-pathogen interactions by regulating gene expression in ways that may contribute to evasion of host responses. Alterations in fungal cell wall structures could allow P. destructans to avoid detection by host pattern recognition receptors and antibody responses. This study has also identified several fungal pathways upregulated during WNS infection that may be candidates for mitigating infection pathology. By identifying host-specific pathogen responses, these observations have important implications for host-pathogen evolutionary relationships in WNS and other fungal diseases.

Resolving host-pathogen interactions by dual RNA-seq

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

Dual-species transcriptional profiling during systemic candidiasis reveals organ-specific host-pathogen interactions

Candida albicans is a common cause of life-threatening fungal bloodstream infections. In the murine model of systemic candidiasis, the kidney is the primary target organ while the fungal load declines over time in liver and spleen. To better understand these organ-specific differences in host-pathogen interaction, we performed gene expression profiling of murine kidney, liver and spleen and determined the fungal transcriptome in liver and kidney. We observed a delayed transcriptional immune response accompanied by late induction of fungal stress response genes in the kidneys. In contrast, early upregulation of the proinflammatory response in the liver was associated with a fungal transcriptome resembling response to phagocytosis, suggesting that phagocytes contribute significantly to fungal control in the liver. Notably, C. albicans hypha-associated genes were upregulated in the absence of visible filamentation in the liver, indicating an uncoupling of gene expression and morphology and a morphology-independent effect by hypha-associated genes in this organ. Consistently, integration of host and pathogen transcriptional data in an inter-species gene regulatory network indicated connections of C. albicans cell wall remodelling and metabolism to the organ-specific immune responses.

Integrating Large-Scale Data and RNA Technology to Protect Crops from Fungal Pathogens

With a rapidly growing human population it is expected that plant science researchers and the agricultural community will need to increase food productivity using less arable land. This challenge is complicated by fungal pathogens and diseases, many of which can severely impact crop yield. Current measures to control fungal pathogens are either ineffective or have adverse effects on the agricultural enterprise. Thus, developing new strategies through research innovation to protect plants from pathogenic fungi is necessary to overcome these hurdles. RNA sequencing technologies are increasing our understanding of the underlying genes and gene regulatory networks mediating disease outcomes. The application of invigorating next generation sequencing strategies to study plant-pathogen interactions has and will provide unprecedented insight into the complex patterns of gene activity responsible for crop protection. However, questions remain about how biological processes in both the pathogen and the host are specified in space directly at the site of infection and over the infection period. The integration of cutting edge molecular and computational tools will provide plant scientists with the arsenal required to identify genes and molecules that play a role in plant protection. Large scale RNA sequence data can then be used to protect plants by targeting genes essential for pathogen viability in the production of stably transformed lines expressing RNA interference molecules, or through foliar applications of double stranded RNA.