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

Exploring the Genome-wide Expression Level of the Bacterial Strain Belonging to Bacillus safensis (MM19) Against Phomopsis viticola

Rhizobacteria have the ability to compete with pathogenic microorganisms and contribute to plant immunity and defense mechanisms. Their growth and survival in the rhizosphere ensure a biological balance in favor of the host plant. The differential gene expression profiles of B. safensis (MM19) revealed significantly increased expression of prominent genes related to thiamine biosynthesis involving various metabolites and enzymes that participate in the suppression of mycelium growth and pathogen inhibition. Correspondingly, the expression of three major genes (HOG1, FUS3, SGI) involved in the virulence of P. viticola was assessed using qPCR analysis. HOG1 was the highest expressed gene in the pathogen when it was co-cultivated with MM19. Based on these findings, we performed molecular docking and dynamics analysis to explore the interaction between HOG1 and thiamine, as well as expression network analysis constructed using Cytoscape. The functional genomic data related to thiamine biosynthesis and the corresponding pathways ensure a priming role in the antagonistic behavior of B. safensis (MM19) against P. viticola as a support for plant immunity.

Legume-rhizobia symbiosis: Translatome analysis

Leguminous plants can establish endosymbiotic relationships with nitrogen-fixing soil rhizobacteria. Bacterial infection and nodule organogenesis are two independent but highly coordinated genetic programs that are active during this interaction. These genetic programs can be regulated along all the stages of gene expression. Most of the studies, for both eukaryotes and prokaryotes, focused on the transcriptional regulation level determining the abundance of mRNAs. However, it has been demonstrated that mRNA levels only sometimes correlate with the abundance or activity of the coded proteins. For this reason, in the past two decades, interest in the role of translational control of gene expression has increased, since the subset of mRNA being actively translated outperforms the information gained only by the transcriptome. In the case of legume-rhizobia interactions, the study of the translatome still needs to be explored further. Therefore, this review aims to discuss the methodologies for analyzing polysome-associated mRNAs at the genome-scale and their contribution to studying translational control to understand the complexity of this symbiotic interaction. Moreover, the Dual RNA-seq approach is discussed for its relevance in the context of a symbiotic nodule, where intricate multi-species gene expression networks occur.

Identification and Characterization of Negative Regulators of Rgg1518 Quorum Sensing in Streptococcus pneumoniae

Streptococcus pneumoniae is an agent of otitis media, septicemia, and meningitis and remains the leading cause of community-acquired pneumonia regardless of vaccine use. Of the various strategies that S. pneumoniae takes to enhance its potential to colonize the human host, quorum sensing (QS) is an intercellular communication process that provides coordination of gene expression at a community level. Numerous putative QS systems are identifiable in the S. pneumoniae genome, but their gene-regulatory activities and contributions to fitness have yet to be fully evaluated. To contribute to assessing regulatory activities of rgg paralogs present in the D39 genome, we conducted transcriptomic analysis of mutants of six QS regulators. Our results find evidence that at least four QS regulators impact the expression of a polycistronic operon (encompassing genes spd_1517 to spd_1513) that is directly controlled by the Rgg/SHP1518 QS system. As an approach to unravel the convergent regulation placed on the spd_1513-1517 operon, we deployed transposon mutagenesis screening in search of upstream regulators of the Rgg/SHP1518 QS system. The screen identified two types of insertion mutants that result in increased activity of Rgg1518-dependent transcription, one type being where the transposon inserted into pepO, an annotated endopeptidase, and the other type being insertions in spxB, a pyruvate oxidase. We demonstrate that pneumococcal PepO degrades SHP1518 to prevent activation of Rgg/SHP1518 QS. Moreover, the glutamic acid residue in the conserved "HExxH" domain is indispensable for the catalytic function of PepO. Finally, we confirmed the metalloendopeptidase property of PepO, which requires zinc ions, but not other ions, to facilitate peptidyl hydrolysis. IMPORTANCE Streptococcus pneumoniae uses quorum sensing to communicate and regulate virulence. In our study, we focused on one Rgg quorum sensing system (Rgg/SHP1518) and found that multiple other Rgg regulators also control it. We further identified two enzymes that inhibit Rgg/SHP1518 signaling and revealed and validated one enzyme's mechanisms for breaking down quorum sensing signaling molecules. Our findings shed light on the complex regulatory network of quorum sensing in Streptococcus pneumoniae.

Multiplexed target enrichment of coding and non-coding transcriptomes enables studying Candida spp. infections from human derived samples

The study of transcriptomic interactions between host and pathogens in in vivo conditions is challenged by the low relative amounts of the pathogen RNA. Yeast opportunistic pathogens of the genus Candida can cause life-threatening systemic infections in immunocompromised patients, and are of growing medical concern. Four phylogenetically diverse species account for over 90% of Candida infections, and their specific interactions with various human tissues are still poorly understood. To enable in vivo transcriptomic analysis in these species, we designed and validated pan-Candida target capture probes to enrich protein-coding and non-coding transcriptomes. The probe-based enrichment approach outperformed enrichment based on differential lysis of host cells, and showed similar enrichment performance as an existing capture design, yet achieving better fidelity of expression levels, enabling species multiplexing and capturing of lncRNAs. In addition, we show that our probe-based enrichment strategy allows robust genotype-based identification of the infecting strain present in the sample.

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

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

Developing New Tools to Fight Human Pathogens: A Journey through the Advances in RNA Technologies

A long scientific journey has led to prominent technological advances in the RNA field, and several new types of molecules have been discovered, from non-coding RNAs (ncRNAs) to riboswitches, small interfering RNAs (siRNAs) and CRISPR systems. Such findings, together with the recognition of the advantages of RNA in terms of its functional performance, have attracted the attention of synthetic biologists to create potent RNA-based tools for biotechnological and medical applications. In this review, we have gathered the knowledge on the connection between RNA metabolism and pathogenesis in Gram-positive and Gram-negative bacteria. We further discuss how RNA techniques have contributed to the building of this knowledge and the development of new tools in synthetic biology for the diagnosis and treatment of diseases caused by pathogenic microorganisms. Infectious diseases are still a world-leading cause of death and morbidity, and RNA-based therapeutics have arisen as an alternative way to achieve success. There are still obstacles to overcome in its application, but much progress has been made in a fast and effective manner, paving the way for the solid establishment of RNA-based therapies in the future.

Differing Responses to Phytophthora cinnamomi Infection in Susceptible and Partially Resistant Persea americana (Mill.) Rootstocks: A Case for the Role of Receptor-Like Kinases and Apoplastic Proteases

The hemibiotrophic plant pathogen Phytophthora cinnamomi Rands is the most devastating pathogen of avocado (Persea americana Mill.) and, as such, causes significant annual losses in the industry. Although the molecular basis of P. cinnamomi resistance in avocado and P. cinnamomi virulence determinants have been the subject of recent research, none have yet attempted to compare the transcriptomic responses of both pathogen and host during their interaction. In the current study, the transcriptomes of both avocado and P. cinnamomi were explored by dual RNA sequencing. The basis for partial resistance was sought by the inclusion of both susceptible (R0.12) and partially resistant (Dusa®) rootstocks sampled at early (6, 12 and 24 hours post-inoculation, hpi) and late time-points (120 hpi). Substantial differences were noted in the number of differentially expressed genes found in Dusa® and R0.12, specifically at 12 and 24 hpi. Here, the partially resistant rootstock perpetuated defense responses initiated at 6 hpi, while the susceptible rootstock abruptly reversed course. Instead, gene ontology enrichment confirmed that R0.12 activated pathways related to growth and development, essentially rendering its response at 12 and 24 hpi no different from that of the mock-inoculated controls. As expected, several classes of P. cinnamomi effector genes were differentially expressed in both Dusa® and R0.12. However, their expression differed between rootstocks, indicating that P. cinnamomi might alter the expression of its effector arsenal based on the rootstock. Based on some of the observed differences, several P. cinnamomi effectors were highlighted as potential candidates for further research. Similarly, the receptor-like kinase (RLK) and apoplastic protease coding genes in avocado were investigated, focusing on their potential role in differing rootstock responses. This study suggests that the basis of partial resistance in Dusa® is predicated on its ability to respond appropriately during the early stages following P. cinnamomi inoculation, and that important components of the first line of inducible defense, apoplastic proteases and RLKs, are likely to be important to the observed outcome.

The hidden world within plants: metatranscriptomics unveils the complexity of wood microbiomes

The importance of plants as complex entities influenced by genomes of the associated microorganisms is now seen as a new source of variability for a more sustainable agriculture, also in the light of ongoing climate change. For this reason, we investigated through metatranscriptomics whether the taxa profile and behaviour of microbial communities associated with the wood of 20-year-old grapevine plants are influenced by the health status of the host. We report for the first time a metatranscriptome from a complex tissue in a real environment, highlighting that this approach is able to define the microbial community better than referenced transcriptomic approaches. In parallel, the use of total RNA enabled the identification of bacterial taxa in healthy samples that, once isolated from the original wood tissue, displayed potential biocontrol activities against a wood-degrading fungal taxon. Furthermore, we revealed an unprecedented high number of new viral entities (~120 new viral species among 180 identified) associated with a single and limited environment and with potential impact on the whole holobiont. Taken together, our results suggest a complex multitrophic interaction in which the viral community also plays a crucial role in raising new ecological questions for the exploitation of microbial-assisted sustainable agriculture.

A perspective on insect-microbe holobionts facing thermal fluctuations in a climate-change context

Temperature influences the ecology and evolution of insects and their symbionts by impacting each partner independently and their interactions, considering the holobiont as a primary unit of selection. There are sound data about the responses of these partnerships to constant temperatures and sporadic thermal stress (mostly heat shock). However, the current understanding of the thermal ecology of insect-microbe holobionts remains patchy because the complex thermal fluctuations (at different spatial and temporal scales) experienced by these organisms in nature have often been overlooked experimentally. This may drastically constrain our ability to predict the fate of mutualistic interactions under climate change, which will alter both mean temperatures and thermal variability. Here, we tackle down these issues by focusing on the effects of temperature fluctuations on the evolutionary ecology of insect-microbe holobionts. We propose potentially worth-investigating research avenues to (i) evaluate the relevance of theoretical concepts used to predict the biological impacts of temperature fluctuations when applied to holobionts; (ii) acknowledge the plastic (behavioural thermoregulation, physiological acclimation) and genetic responses (evolution) expressed by holobionts in fluctuating thermal environments; and (iii) explore the potential impacts of previously unconsidered patterns of temperature fluctuations on the outcomes and the dynamic of these insect-microbe associations.

Transcriptome Analyses of Barley Roots Inoculated with Novel Paenibacillus sp. and Erwinia gerundensis Strains Reveal Beneficial Early-Stage Plant-Bacteria Interactions

Plant growth-promoting bacteria can improve host plant traits including nutrient uptake and metabolism and tolerance to biotic and abiotic stresses. Understanding the molecular basis of plant–bacteria interactions using dual RNA-seq analyses provides key knowledge of both host and bacteria simultaneously, leading to future enhancements of beneficial interactions. In this study, dual RNA-seq analyses were performed to provide insights into the early-stage interactions between barley seedlings and three novel bacterial strains (two Paenibacillus sp. strains and one Erwinia gerundensis strain) isolated from the perennial ryegrass seed microbiome. Differentially expressed bacterial and barley genes/transcripts involved in plant–bacteria interactions were identified, with varying species- and strain-specific responses. Overall, transcriptome profiles suggested that all three strains improved stress response, signal transduction, and nutrient uptake and metabolism of barley seedlings. Results also suggested potential improvements in seedling root growth via repressing ethylene biosynthesis in roots. Bacterial secondary metabolite gene clusters producing compounds that are potentially associated with interactions with the barley endophytic microbiome and associated with stress tolerance of plants under nutrient limiting conditions were also identified. The results of this study provided the molecular basis of plant growth-promoting activities of three novel bacterial strains in barley, laid a solid foundation for the future development of these three bacterial strains as biofertilisers, and identified key differences between bacterial strains of the same species in their responses to plants.

Lineage-Specific Rewiring of Core Pathways Predating Innovation of Legume Nodules Shapes Symbiotic Efficiency

The interkingdom coevolution innovated the rhizobium-legume symbiosis. The application of this nitrogen-fixing system in sustainable agriculture is usually impeded by incompatible interactions between partners. However, the progressive evolution of rhizobium-legume compatibility remains elusive. In this work, deletions of rhcV encoding a structural component of the type three secretion system allow related Sinorhizobium strains to nodulate a previously incompatible soybean cultivar (Glycine max). These rhcV mutants show low to medium to high symbiotic efficiency on the same cultivated soybean while being indistinguishable on wild soybean plants (Glycine soja). The dual pantranscriptomics reveals nodule-specific activation of core symbiosis genes of Sinorhizobium and Glycine genes associated with genome duplication events along the chronogram. Unexpectedly, symbiotic efficiency is in line with lineage-dependent transcriptional profiles of core pathways which predate the diversification of Fabaceae and Sinorhizobium. This is supported by further physiological and biochemical experiments. Particularly, low-efficiency nodules show disordered antioxidant activity and low-energy status, which restrict nitrogen fixation activity. Collectively, the ancient core pathways play a crucial role in optimizing the function of later-evolved mutualistic arsenals in the rhizobium-legume coevolution.

IMPORTANCE Significant roles of complex extracellular microbiota in environmental adaptation of eukaryotes in ever-changing circumstances have been revealed. Given the intracellular infection ability, facultative endosymbionts can be considered pioneers within complex extracellular microbiota and are ideal organisms for understanding the early stage of interkingdom adaptation. This work reveals that the later innovation of key symbiotic arsenals and the lineage-specific network rewiring in ancient core pathways, predating the divergence of legumes and rhizobia, underline the progressive evolution of rhizobium-legume compatibility. This insight not only is significant for improving the application benefits of rhizobial inoculants in sustainable agriculture but also advances our general understanding of the interkingdom coevolution which is theoretically explored by all host-microbiota interactions.

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

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

Escalation in the host-pathogen arms race: A host resistance response corresponds to a heightened bacterial virulence response

The zig-zag model of host-pathogen interaction describes the relative strength of defense response across a spectrum of pathogen-induced plant phenotypes. A stronger defense response results in increased resistance. Here, we investigate the strength of pathogen virulence during disease and place these findings in the context of the zig-zag model. Xanthomonas vasicola pv. holcicola (Xvh) causes sorghum bacterial leaf streak. Despite being widespread, this disease has not been described in detail at the molecular level. We divided diverse sorghum genotypes into three groups based on disease symptoms: water-soaked lesions, red lesions, and resistance. Bacterial growth assays confirmed that these three phenotypes represent a range of resistance and susceptibility. To simultaneously reveal defense and virulence responses across the spectrum of disease phenotypes, we performed dual RNA-seq on Xvh-infected sorghum. Consistent with the zig-zag model, the expression of plant defense-related genes was strongest in the resistance interaction. Surprisingly, bacterial virulence genes related to the type III secretion system (T3SS) and type III effectors (T3Es) were also most highly expressed in the resistance interaction. This expression pattern was observed at multiple time points within the sorghum-Xvh pathosystem. Further, a similar expression pattern was observed in Arabidopsis infected with Pseudomonas syringae for effector-triggered immunity via AvrRps4 but not AvrRpt2. Specific metabolites were able to repress the Xvh virulence response in vitro and in planta suggesting a possible signaling mechanism. Taken together, these findings reveal multiple permutations of the continually evolving host-pathogen arms race from the perspective of host defense and pathogen virulence responses.

The arms race between plants and pathogens is a complex process. To dissect the plant defense and pathogen virulence responses simultaneously, we used sorghum and Xanthomonas vasicola pv. holcicola, as a model pathosystem. We performed dual RNA-seq on infected sorghum with a range of disease phenotypes. Our characterization of this pathosystem demonstrates that genes related to the plant defense and pathogen virulence responses are most highly induced during a resistance interaction. We observed a similar pattern of escalation in Arabidopsis infected with Pseudomonas syringae. These observations support a conceptual model of fluidity between the different stages of plant immunity and pathogen virulence.

Omics and interspecies interaction

Interspecies interactions are key determinants in biofilm behavior, ecology, and architecture. The cellular responses of microorganisms to each other at transcriptional, proteomic, and metabolomic levels ultimately determine the characteristics of biofilm and the corresponding implications for health and disease. Advances in omics technologies have revolutionized our understanding of microbial community composition and their activities as a whole. Large-scale analyses of the complex interaction between the many microbial species residing within a biofilm, however, are currently still hampered by technical and bioinformatics challenges. Thus, studies of interspecies interactions have largely focused on the transcriptional and proteomic changes that occur during the contact of a few prominent species, such as Porphyromonas gingivalis, Streptococcus mutans, Candida albicans, and a few others, with selected partner species. Expansion of available tools is necessary to grow the revealing, albeit limited, insight these studies have provided into a profound understanding of the nature of individual microbial responses to the presence of others. This will allow us to answer important questions including: Which intermicrobial interactions orchestrate the myriad of cooperative, synergistic, antagonistic, manipulative, and other types of relationships and activities in the complex biofilm environment, and what are the implications for oral health and disease?

Triple RNA-Seq Reveals Synergy in a Human Virus-Fungus Co-infection Model

High-throughput RNA sequencing (RNA-seq) is routinely applied to study diverse biological processes; however, when performed separately on interacting organisms, systemic noise intrinsic to RNA extraction, library preparation, and sequencing hampers the identification of cross-species interaction nodes. Here, we develop triple RNA-seq to simultaneously detect transcriptomes of monocyte-derived dendritic cells (moDCs) infected with the frequently co-occurring pulmonary pathogens Aspergillus fumigatus and human cytomegalovirus (CMV). Comparing expression patterns after co-infection with those after single infections, our data reveal synergistic effects and mutual interferences between host responses to the two pathogens. For example, CMV attenuates the fungus-mediated activation of pro-inflammatory cytokines through NF-κB (nuclear factor κB) and NFAT (nuclear factor of activated T cells) cascades, while A. fumigatus impairs viral clearance by counteracting viral nucleic acid-induced activation of type I interferon signaling. Together, the analytical power of triple RNA-seq proposes molecular hubs in the differential moDC response to fungal/viral single infection or co-infection that contribute to our understanding of the etiology and, potentially, clearance of post-transplant infections.

Dual RNA-sequencing analyses of a coral and its native symbiont during the establishment of symbiosis

Despite the ecological significance of the mutualistic relationship between Symbiodiniaceae and reef-building corals, the molecular interactions during establishment of this relationship are not well understood. This is particularly true of the transcriptional changes that occur in the symbiont. In the current study, a dual RNA-sequencing approach was used to better understand transcriptional changes on both sides of the coral-symbiont interaction during the colonization of Acropora tenuis by a compatible Symbiodiniaceae strain (Cladocopium goreaui; ITS2 type C1). Comparison of transcript levels of the in hospite symbiont 3, 12, 48 and 72 hr after exposure to those of the same strain in culture revealed that extensive and generalized down-regulation of symbiont gene expression occurred during the infection process. Included in this "symbiosis-derived transcriptional repression" were a range of stress response and immune-related genes. In contrast, a suite of symbiont genes implicated in metabolism was upregulated in the symbiotic state. The coral data support the hypothesis that immune-suppression and arrest of phagosome maturation play important roles during the establishment of compatible symbioses, and additionally imply the involvement of some SCRiP family members in the colonization process. Consistent with previous ecological studies, the transcriptomic data suggest that active translocation of metabolites to the host may begin early in the colonization process, and thus that the mutualistic relationship can be established at the larval stage. This dual RNA-sequencing study provides insights into the transcriptomic remodelling that occurs in C. goreaui during transition to a symbiotic lifestyle and the novel coral genes implicated in symbiosis.

CROSSMAPPER: estimating cross-mapping rates and optimizing experimental design in multi-species sequencing studies

Motivation

Numerous sequencing studies, including transcriptomics of host-pathogen systems, sequencing of hybrid genomes, xenografts, mixed species systems, metagenomics and meta-transcriptomics, involve samples containing genetic material from divergent organisms. A crucial step in these studies is identifying from which organism each sequencing read originated, and the experimental design should be directed to minimize biases caused by cross-mapping of reads to incorrect source genomes. Additionally, pooling of sufficiently different genetic material into a single sequencing library could significantly reduce experimental costs but requires careful planning and assessment of the impact of cross-mapping. Having these applications in mind we designed Crossmapper, the first to our knowledge tool able to assess cross-mapping prior to sequencing, therefore allowing optimization of experimental design.

Results

Using any combination of reference genomes, Crossmapper performs read simulation and back-mapping of those reads to the pool of references, quantifies and reports the cross-mapping rates for each organism. Crossmapper performs these analyses with numerous user-specified parameters, including, among others, read length, read layout, coverage, mapping parameters, genomic or transcriptomic data. Additionally, it outputs the results in highly interactive and publication-ready reports. This allows the user to perform multiple comparisons at once and choose the experimental setup minimizing cross-mapping rates. Moreover, Crossmapper can be used for resource optimization in sequencing facilities by pooling different samples into one sequencing library.

Availability and implementation

Crossmapper is a command line tool implemented in Python 3.6 and available as a conda package, allowing effortless installation. The source code, detailed information and a step-by-step tutorial is available at our GitHub page https://github.com/Gabaldonlab/crossmapper.

Supplementary information

Supplementary data are available at Bioinformatics online.

Soil microbiota influences clubroot disease by modulating Plasmodiophora brassicae and Brassica napus transcriptomes

The contribution of surrounding plant microbiota to disease development has led to the 'pathobiome' concept, which represents the interaction between the pathogen, the host plant and the associated biotic microbial community, resulting or not in plant disease. The aim herein is to understand how the soil microbial environment may influence the functions of a pathogen and its pathogenesis, and the molecular response of the plant to the infection, with a dual-RNAseq transcriptomics approach. We address this question using Brassica napus and Plasmodiophora brassicae, the pathogen responsible for clubroot. A time-course experiment was conducted to study interactions between P. brassicae, two B. napus genotypes and three soils harbouring high, medium or low microbiota diversities and levels of richness. The soil microbial diversity levels had an impact on disease development (symptom levels and pathogen quantity). The P. brassicae and B. napus transcriptional patterns were modulated by these microbial diversities, these modulations being dependent on the host genotype plant and the kinetic time. The functional analysis of gene expressions allowed the identification of pathogen and plant host functions potentially involved in the change of plant disease level, such as pathogenicity-related genes (NUDIX effector) in P. brassicae and plant defence-related genes (glucosinolate metabolism) in B. napus.

In vivo dual RNA-seq reveals that neutrophil recruitment underlies differential tissue tropism of Streptococcus pneumoniae

Streptococcus pneumoniae is a genetically diverse human-adapted pathogen commonly carried asymptomatically in the nasopharynx. We have recently shown that a single nucleotide polymorphism (SNP) in the raffinose pathway regulatory gene rafR accounts for a difference in the capacity of clonally-related strains to cause localised versus systemic infection. Using dual RNA-seq, we show that this SNP affects expression of bacterial genes encoding multiple sugar transporters, and fine-tunes carbohydrate metabolism, along with extensive rewiring of host transcriptional responses to infection, particularly expression of genes encoding cytokine and chemokine ligands and receptors. The data predict a crucial role for differential neutrophil recruitment (confirmed by in vivo neutrophil depletion and IL-17 neutralization) indicating that early detection of bacteria by the host in the lung environment is crucial for effective clearance. Thus, dual RNA-seq provides a powerful tool for understanding complex host-pathogen interactions and reveals how a single bacterial SNP can drive differential disease outcomes.

Minhas, Aprianto et al. apply dual RNA seq to a set of related Streptococcus pneumoniae strains to find that differential neutrophil recruitment explains different tissue tropism of these strains. This study highlights the power of dual RNA-seq in investigating how a single bacterial SNP determines the host’s disease outcomes.

Investigations of the mechanisms of interactions between four non-conventional species with Saccharomyces cerevisiae in oenological conditions

Fermentation by microorganisms is a key step in the production of traditional food products such as bread, cheese, beer and wine. In these fermentative ecosystems, microorganisms interact in various ways, namely competition, predation, commensalism and mutualism. Traditional wine fermentation is a complex microbial process performed by Saccharomyces and non-Saccharomyces (NS) yeast species. To better understand the different interactions occurring within wine fermentation, isolated yeast cultures were compared with mixed co-cultures of one reference strain of S. cerevisiae with one strain of four NS yeast species (Metschnikowia pulcherrima, M. fructicola, Hanseniaspora opuntiae and H. uvarum). In each case, we studied population dynamics, resource consumed and metabolites produced from central carbon metabolism. This phenotyping of competition kinetics allowed us to confirm the main mechanisms of interaction between strains of four NS species. S. cerevisiae competed with H. uvarum and H. opuntiae for resources although both Hanseniaspora species were characterized by a strong mortality either in mono or mixed fermentations. M. pulcherrima and M. fructicola displayed a negative interaction with the S. cerevisiae strain tested, with a decrease in viability in co-culture. Overall, this work highlights the importance of measuring specific cell populations in mixed cultures and their metabolite kinetics to understand yeast-yeast interactions. These results are a first step towards ecological engineering and the rational design of optimal multi-species starter consortia using modeling tools. In particular the originality of this paper is for the first times to highlight the joint-effect of different species population dynamics on glycerol production and also to discuss on the putative role of lipid uptake on the limitation of some non-conventional species growth although interaction processes.

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

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

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.

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.

The Transcriptional Profile of Trichophyton rubrum Co-Cultured with Human Keratinocytes Shows New Insights about Gene Modulation by Terbinafine

The dermatophyte Trichophyton rubrum is the main causative agent of dermatophytoses worldwide. Although a superficial mycosis, its incidence has been increasing especially among diabetic and immunocompromised patients. Terbinafine is commonly used for the treatment of infections caused by dermatophytes. However, cases of resistance of T. rubrum to this allylamine were reported even with the efficacy of this drug. The present study is the first to evaluate the effect of terbinafine using a co-culture model of T. rubrum and human keratinocytes, mimicking a fungus-host interaction, in conjunction with RNA-seq technique. Our data showed the repression of several genes involved in the ergosterol biosynthesis cascade and the induction of genes encoding major facilitator superfamily (MFS)- and ATP-binding cassette superfamily (ABC)-type membrane transporter which may be involved in T. rubrum mechanisms of resistance to this drug. We observed that some genes reported in the scientific literature as candidates of new antifungal targets were also modulated. In addition, we found the modulation of several genes that are hypothetical in T. rubrum but that possess known orthologs in other dermatophytes. Taken together, the results indicate that terbinafine can act on various targets related to the physiology of T. rubrum other than its main target of ergosterol biosynthetic pathway.

Transcription Profiling Reveals Cooperative Metabolic Interactions in a Microbial Cheese-Ripening Community Composed of Debaryomyces hansenii, Brevibacterium aurantiacum, and Hafnia alvei

Ripening cultures containing fungi and bacteria are widely used in smear-ripened cheese production processes, but little is known about the biotic interactions of typical ripening microorganisms at the surface of cheese. We developed a lab-scale mini-cheese model to investigate the biotic interactions of a synthetic community that was composed of Debaryomyces hansenii, Brevibacterium aurantiacum, and Hafnia alvei, three species that are commonly used for smear-ripened cheese production. Transcriptomic analyses of cheese samples produced with different combinations of these three species revealed potential mechanisms of biotic interactions concerning iron acquisition, proteolysis, lipolysis, sulfur metabolism, and D-galactonate catabolism. A strong mutualistic interaction was observed between H. alvei and B. aurantiacum. We propose an explanation of this positive interaction in which B. aurantiacum would benefit from siderophore production by H. alvei, and the latter would be stimulated by the energy compounds liberated from caseins and triglycerides through the action of the proteases and lipases secreted by B. aurantiacum. In the future, it would be interesting to take the iron acquisition systems of cheese-associated strains into account for the purpose of improving the selection of the ripening culture components and their association in mixed cultures.

Clostridioides difficile LuxS mediates inter-bacterial interactions within biofilms

The anaerobic gut pathogen, Clostridioides difficile, forms adherent biofilms that may play an important role in recurrent C. difficile infections. The mechanisms underlying C. difficile community formation and inter-bacterial interactions are nevertheless poorly understood. C. difficile produces AI-2, a quorum sensing molecule that modulates biofilm formation across many bacterial species. We found that a strain defective in LuxS, the enzyme that mediates AI-2 production, is defective in biofilm development in vitro. Transcriptomic analyses of biofilms formed by wild type (WT) and luxS mutant (luxS) strains revealed a downregulation of prophage loci in the luxS mutant biofilms compared to the WT. Detection of phages and eDNA within biofilms may suggest that DNA release by phage-mediated cell lysis contributes to C. difficile biofilm formation. In order to understand if LuxS mediates C. difficile crosstalk with other gut species, C. difficile interactions with a common gut bacterium, Bacteroides fragilis, were studied. We demonstrate that C. difficile growth is significantly reduced when co-cultured with B. fragilis in mixed biofilms. Interestingly, the absence of C. difficile LuxS alleviates the B. fragilis-mediated growth inhibition. Dual species RNA-sequencing analyses from single and mixed biofilms revealed differential modulation of distinct metabolic pathways for C. difficile WT, luxS and B. fragilis upon co-culture, indicating that AI-2 may be involved in induction of selective metabolic responses in B. fragilis. Overall, our data suggest that C. difficile LuxS/AI-2 utilises different mechanisms to mediate formation of single and mixed species communities.

High-resolution analysis of the pneumococcal transcriptome under a wide range of infection-relevant conditions

Streptococcus pneumoniae is an opportunistic human pathogen that typically colonizes the nasopharyngeal passage and causes lethal disease in other host niches, such as the lung or the meninges. The expression and regulation of pneumococcal genes at different life-cycle stages, such as commensal or pathogenic, are not entirely understood. To chart the transcriptional responses of S. pneumoniae, we used RNA-seq to quantify the relative abundance of the transcriptome under 22 different infection-relevant conditions. The data demonstrated a high level of dynamic expression and, strikingly, all annotated pneumococcal genomic features were expressed in at least one of the studied conditions. By computing the correlation values of every pair of genes across all studied conditions, we created a co-expression matrix that provides valuable information on both operon structure and regulatory processes. The co-expression data are highly consistent with well-characterized operons and regulons, such as the PyrR, ComE and ComX regulons, and have allowed us to identify a new member of the competence regulon. Lastly, we created an interactive data center named PneumoExpress (https://veeninglab.com/pneumoexpress) that enables users to access the expression data as well as the co-expression matrix in an intuitive and efficient manner, providing a valuable resource to the pneumococcal research community.

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.

Candida parapsilosis: from Genes to the Bedside

Patients with suppressed immunity are at the highest risk for hospital-acquired infections. Among these, invasive candidiasis is the most prevalent systemic fungal nosocomial infection. Over recent decades, the combined prevalence of non-albicans Candida species outranked Candida albicans infections in several geographical regions worldwide, highlighting the need to understand their pathobiology in order to develop effective treatment and to prevent future outbreaks. Candida parapsilosis is the second or third most frequently isolated Candida species from patients. Besides being highly prevalent, its biology differs markedly from that of C. albicans, which may be associated with C. parapsilosis' increased incidence. Differences in virulence, regulatory and antifungal drug resistance mechanisms, and the patient groups at risk indicate that conclusions drawn from C. albicans pathobiology cannot be simply extrapolated to C. parapsilosis Such species-specific characteristics may also influence their recognition and elimination by the host and the efficacy of antifungal drugs. Due to the availability of high-throughput, state-of-the-art experimental tools and molecular genetic methods adapted to C. parapsilosis, genome and transcriptome studies are now available that greatly contribute to our understanding of what makes this species a threat. In this review, we summarize 10 years of findings on C. parapsilosis pathogenesis, including the species' genetic properties, transcriptome studies, host responses, and molecular mechanisms of virulence. Antifungal susceptibility studies and clinician perspectives are discussed. We also present regional incidence reports in order to provide an updated worldwide epidemiology summary.

Transcriptome Sequencing Approaches to Elucidate Host-Microbe Interactions in Opportunistic Human Fungal Pathogens

Infections caused by opportunistic human fungal pathogens are a source of increasing medical concern, due to their growing incidence, the emergence of novel pathogenic species, and the lack of effective diagnostics tools. Fungal pathogens are phylogenetically diverse, and their virulence mechanisms can differ widely across species. Despite extensive efforts, the molecular bases of virulence in pathogenic fungi and their interactions with the human host remain poorly understood for most species. In this context, next-generation sequencing approaches hold the promise of helping to close this knowledge gap. In particular, high-throughput transcriptome sequencing (RNA-Seq) enables monitoring the transcriptional profile of both host and microbes to elucidate their interactions and discover molecular mechanisms of virulence and host defense. Here, we provide an overview of transcriptome sequencing techniques and approaches, and survey their application in studying the interplay between humans and fungal pathogens. Finally, we discuss novel RNA-Seq approaches in studying host-pathogen interactions and their potential role in advancing the clinical diagnostics of fungal infections.

Systems Biology and Machine Learning in Plant-Pathogen Interactions

Systems biology is an inclusive approach to study the static and dynamic emergent properties on a global scale by integrating multiomics datasets to establish qualitative and quantitative associations among multiple biological components. With an abundance of improved high throughput -omics datasets, network-based analyses and machine learning technologies are playing a pivotal role in comprehensive understanding of biological systems. Network topological features reveal most important nodes within a network as well as prioritize significant molecular components for diverse biological networks, including coexpression, protein-protein interaction, and gene regulatory networks. Machine learning techniques provide enormous predictive power through specific feature extraction from biological data. Deep learning, a subtype of machine learning, has plausible future applications because a domain expert for feature extraction is not needed in this algorithm. Inspired by diverse domains of biology, we here review classic systems biology techniques applied in plant immunity thus far. We also discuss additional advanced approaches in both graph theory and machine learning, which may provide new insights for understanding plant-microbe interactions. Finally, we propose a hybrid approach in plant immune systems that harnesses the power of both network biology and machine learning, with a potential to be applicable to both model systems and agronomically important crop plants.

Dynamic interactions between RNA viruses and human hosts unravelled by a decade of next generation sequencing

BACKGROUND

Next generation sequencing (NGS) methods have significantly contributed to a paradigm shift in genomic research for nearly a decade now. These methods have been useful in studying the dynamic interactions between RNA viruses and human hosts.

SCOPE OF THE REVIEW

In this review, we summarise and discuss key applications of NGS in studying the host - pathogen interactions in RNA viral infections of humans with examples.

MAJOR CONCLUSIONS

Use of NGS to study globally relevant RNA viral infections have revolutionized our understanding of the within host and between host evolution of these viruses. These methods have also been useful in clinical decision-making and in guiding biomedical research on vaccine design.

GENERAL SIGNIFICANCE

NGS has been instrumental in viral genomic studies in resolving within-host viral genomic variants and the distribution of nucleotide polymorphisms along the full-length of viral genomes in a high throughput, cost effective manner. In the future, novel advances such as long read, single molecule sequencing of viral genomes and simultaneous sequencing of host and pathogens may become the standard of practice in research and clinical settings. This will also bring on new challenges in big data analysis.

Concept of microbial gatekeepers: Positive guys?

Microorganisms contribute diverse and fundamental roles in biogeochemical processes. In a given microbial community, individuals interact with one another to form complex regulatory networks in which gatekeepers contribute disproportional roles in sustaining stability, dynamics, and function. Owing to the ecological and functional importance of microbial gatekeeper, this review provides an overview on its history, identification, roles, application in biological sciences, and clinical diagnostics. The roles of microbial gatekeepers can be beneficial or detrimental, which depends on our purpose. As the field is rather new, some limitations are raised, and further efforts devoted to solving these concerns are proposed. Collectively, gatekeepers provide promising targets for sustaining and re-establishing a desired microbial community. However, substantial obstacles, such as factors governing gatekeepers, must be overcome to manipulate gatekeepers as positive guys.

Toward Cell Type-Specific In Vivo Dual RNA-Seq

Dual RNA-seq has emerged as a genome-wide expression profiling technique, simultaneously measuring RNA transcript levels in a given host and its pathogen during an infection. Recently, the method was transferred from cell culture to in vivo models of bacterial infections; however, specific host cell-type resolution has not yet been achieved. Here we present a detailed protocol that describes the application of Dual RNA-seq to murine colonocytes isolated from mice infected with the enteric pathogen Salmonella Typhimurium. At day 5 after oral infection, the mice were humanely euthanized, their colons extracted, and colonocytes isolated and fixed. Upon antibody staining of cell type-specific surface markers, the fraction of Salmonella-invaded colonocytes was collected by fluorescence-activated cell sorting based on a fluorescent signal emitted by the internalized bacteria. Total RNA was extracted from cells enriched by this method, and ribosomal transcripts from host and microbial cells were removed prior to cDNA synthesis and library generation. We compared different protocols for library preparation and discuss their respective advantages and caveats when applied to minute RNA amounts that constitute an inherent challenge for in vivo transcriptomics. Our results introduce an ultralow input protocol that holds promise for cell type-specific in vivo Dual RNA-seq for charting gene expression of a bacterial pathogen within its respective in vivo niche, along with the consequent host response.

Dual RNA-Seq Analysis of Trichophyton rubrum and HaCat Keratinocyte Co-Culture Highlights Important Genes for Fungal-Host Interaction

The dermatophyte Trichophyton rubrum is the major fungal pathogen of skin, hair, and nails that uses keratinized substrates as the primary nutrients during infection. Few strategies are available that permit a better understanding of the molecular mechanisms involved in the interaction of T. rubrum with the host because of the limitations of models mimicking this interaction. Dual RNA-seq is a powerful tool to unravel this complex interaction since it enables simultaneous evaluation of the transcriptome of two organisms. Using this technology in an in vitro model of co-culture, this study evaluated the transcriptional profile of genes involved in fungus-host interactions in 24 h. Our data demonstrated the induction of glyoxylate cycle genes, ERG6 and TERG_00916, which encodes a carboxylic acid transporter that may improve the assimilation of nutrients and fungal survival in the host. Furthermore, genes encoding keratinolytic proteases were also induced. In human keratinocytes (HaCat) cells, the SLC11A1, RNASE7, and CSF2 genes were induced and the products of these genes are known to have antimicrobial activity. In addition, the FLG and KRT1 genes involved in the epithelial barrier integrity were inhibited. This analysis showed the modulation of important genes involved in T. rubrum–host interaction, which could represent potential antifungal targets for the treatment of dermatophytoses.