Gene expression: Host-pathogen duels revealed by dual RNA-seq in vivo

Dual RNA-seq study of the dynamics of coding and non-coding RNA expression during Clostridioides difficile infection in a mouse model

Clostridioides difficile is the leading cause of healthcare-associated diarrhea in industrialized countries. Many questions remain to be answered about the mechanisms governing its interaction with the host during infection. Non-coding RNAs (ncRNAs) contribute to shape virulence in many pathogens and modulate host responses; however, their role in C. difficile infection (CDI) has not been explored. To better understand the dynamics of ncRNA expression contributing to C. difficile infectious cycle and host response, we used a dual RNA-seq approach in a conventional murine model. From the pathogen side, this transcriptomic analysis revealed the upregulation of virulence factors, metabolism, and sporulation genes, as well as the identification of 61 ncRNAs differentially expressed during infection that correlated with the analysis of available raw RNA-seq data sets from two independent studies. From these data, we identified 118 potential new transcripts in C. difficile, including 106 new ncRNA genes. From the host side, we observed the induction of several pro-inflammatory pathways, and among the 185 differentially expressed ncRNAs, the overexpression of microRNAs (miRNAs) previously associated to inflammatory responses or unknown long ncRNAs and miRNAs. A particular host gene expression profile could be associated to the symptomatic infection. In accordance, the metatranscriptomic analysis revealed specific microbiota changes accompanying CDI and specific species associated with symptomatic infection in mice. This first adaptation of in vivo dual RNA-seq to C. difficile contributes to unravelling the regulatory networks involved in C. difficile infectious cycle and host response and provides valuable resources for further studies of RNA-based mechanisms during CDI.IMPORTANCEClostridioides difficile is a major cause of nosocomial infections associated with antibiotic therapy classified as an urgent antibiotic resistance threat. This pathogen interacts with host and gut microbial communities during infection, but the mechanisms of these interactions remain largely to be uncovered. Noncoding RNAs contribute to bacterial virulence and host responses, but their expression has not been explored during C. difficile infection. We took advantage of the conventional mouse model of C. difficile infection to look simultaneously to the dynamics of gene expression in pathogen, its host, and gut microbiota composition, providing valuable resources for future studies. We identified a number of ncRNAs that could mediate the adaptation of C. difficile inside the host and the crosstalk with the host immune response. Promising inflammation markers and potential therapeutic targets emerged from this work open new directions for RNA-based and microbiota-modulatory strategies to improve the efficiency of C. difficile infection treatments.

Transcriptomic changes in the microsporidia proliferation and host responses in congenitally infected embryos and larvae

Congenital infection caused by vertical transmission of microsporidia N. bombycis can result in severe economic losses in the silkworm-rearing industry. Whole-transcriptome analyses have revealed non-coding RNAs and their regulatory networks in N. bombycis infected embryos and larvae. However, transcriptomic changes in the microsporidia proliferation and host responses in congenitally infected embryos and larvae remains unclear. Here, we simultaneously compared the transcriptomes of N. bombycis and its host B. mori embryos of 5-day and larvae of 1-, 5- and 10-day during congenital infection. For the transcriptome of N. bombycis, a comparison of parasite expression patterns between congenital-infected embryos and larva showed most genes related to parasite central carbon metabolism were down-regulated in larvae during infection, whereas the majority of genes involved in parasite proliferation and growth were up-regulated. Interestingly, a large number of distinct or shared differentially expressed genes (DEGs) were revealed by the Venn diagram and heat map, many of them were connected to infection related factors such as Ricin B lectin, spore wall protein, polar tube protein, and polysaccharide deacetylase. For the transcriptome of B. mori infected with N. bombycis, beyond numerous DEGs related to DNA replication and repair, mRNA surveillance pathway, RNA transport, protein biosynthesis, and proteolysis, with the progression of infection, a large number of DEGs related to immune and infection pathways, including phagocytosis, apoptosis, TNF, Toll-like receptor, NF-kappa B, Fc epsilon RI, and some diseases, were successively identified. In contrast, most genes associated with the insulin signaling pathway, 2-oxacarboxylic acid metabolism, amino acid biosynthesis, and lipid metabolisms were up-regulated in larvae compared to those in embryos. Furthermore, dozens of distinct and three shared DEGs that were involved in the epigenetic regulations, such as polycomb, histone-lysine-specific demethylases, and histone-lysine-N-methyltransferases, were identified via the Venn diagram and heat maps. Notably, many DEGs of host and parasite associated with lipid-related metabolisms were verified by RT-qPCR. Taken together, simultaneous transcriptomic analyses of both host and parasite genes lead to a better understanding of changes in the microsporidia proliferation and host responses in embryos and larvae in N. bombycis congenital infection.

Dual Transcriptome Analysis Reveals That ChATG8 Is Required for Fungal Development, Melanization and Pathogenicity during the Interaction between Colletotrichum higginsianum and Arabidopsis thaliana

Anthracnose disease of cruciferous plants caused by Colletotrichum higginsianum is a serious fungal disease that affects cruciferous crops such as Chinese cabbage, Chinese flowering cabbage, broccoli, mustard plant, as well as the model plant Arabidopsis thaliana. Dual transcriptome analysis is commonly used to identify the potential mechanisms of interaction between host and pathogen. In order to identify differentially expressed genes (DEGs) in both the pathogen and host, the conidia of wild-type (ChWT) and Chatg8 mutant (Chatg8Δ) strains were inoculated onto leaves of A. thaliana, and the infected leaves of A. thaliana at 8, 22, 40, and 60 h post-inoculation (hpi) were subjected to dual RNA-seq analysis. The results showed that comparison of gene expression between the 'ChWT' and 'Chatg8Δ' samples detected 900 DEGs (306 upregulated and 594 down-regulated) at 8 hpi, 692 DEGs (283 upregulated and 409 down-regulated) at 22 hpi, 496 DEGs (220 upregulated and 276 down-regulated) at 40 hpi, and 3159 DEGs (1544 upregulated and 1615 down-regulated) at 60 hpi. GO and KEGG analyses found that the DEGs were mainly involved in fungal development, biosynthesis of secondary metabolites, plant-fungal interactions, and phytohormone signaling. The regulatory network of key genes annotated in the Pathogen-Host Interactions database (PHI-base) and Plant Resistance Genes database (PRGdb), as well as a number of key genes highly correlated with the 8, 22, 40, and 60 hpi, were identified during the infection. Among the key genes, the most significant enrichment was in the gene encoding the trihydroxynaphthalene reductase (THR1) in the melanin biosynthesis pathway. Both Chatg8Δ and Chthr1Δ strains showed varying degrees of reduction of melanin in appressoria and colonies. The pathogenicity of the Chthr1Δ strain was lost. In addition, six DEGs from C. higginsianum and six DEGs from A. thaliana were selected for real-time quantitative PCR (RT-qPCR) to confirm the RNA-seq results. The information gathered from this study enriches the resources available for research into the role of the gene ChATG8 during the infection of A. thaliana by C. higginsianum, such as potential links between melanin biosynthesis and autophagy, and the response of A. thaliana to different fungal strains, thereby providing a theoretical basis for the breeding of cruciferous green leaf vegetable cultivars with resistance to anthracnose disease.

Dual species dynamic transcripts reveal the interaction mechanisms between Chrysanthemum morifolium and Alternaria alternata

BACKGROUND

Chrysanthemum (Chrysanthemum morifolium) black spot disease caused by Alternaria alternata is one of the plant's most destructive diseases. Dual RNA-seq was performed to simultaneously assess their transcriptomes to analyze the potential interaction mechanism between the two species, i.e., host and pathogen.

RESULTS

C. morifolium and A. alternata were subjected to dual RNA-seq at 1, 12, and 24 h after inoculation, and differential expression genes (DEGs) in both species were identified. This analysis confirmed 153,532 DEGs in chrysanthemum and 14,932 DEGs in A. alternata, which were involved in plant-fungal interactions and phytohormone signaling. Fungal DEGs such as toxin synthesis related enzyme and cell wall degrading enzyme genes played important roles during chrysanthemum infection. Moreover, a series of key genes highly correlated with the early, middle, or late infection stage were identified, together with the regulatory network of key genes annotated in the Plant Resistance Genes database (PRGdb) or Pathogen-Host Interactions database (PHI-base). Highly correlated genes were identified at the late infection stage, expanding our understanding of the interplay between C. morifolium and A. alternata. Additionally, six DEGs each from chrysanthemum and A. alternata were selected for quantitative real-time PCR (qRT-PCR) assays to validate the RNA-seq output.

CONCLUSIONS

Collectively, data obtained in this study enriches the resources available for research into the interactions that exist between chrysanthemum and A. alternata, thereby providing a theoretical basis for the development of new chrysanthemum cultivars with resistance to pathogen.

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.

Unveil the transcriptional landscape at the Cryptococcus-host axis in mice and nonhuman primates

Pathogens and hosts require rapid modulation of virulence and defense mechanisms at the infection axis, but monitoring such modulations is challenging. In studying the human fungal pathogen Cryptococcus neoformans, mouse and rabbit infection models are often employed to shed light on the disease mechanisms but that may not be clinically relevant. In this study, we developed an animal infection model using the non-human primate cynomolgus monkey Macaca fascicularis. In addition, we systematically profiled and compared transcriptional responses between the infected mice and the cynomolgus monkey, using simultaneous or dual RNA next-generation sequencing. We demonstrated that there are shared but distinct transcriptional profiles between the two models following C. neoformans infection. Specifically, genes involved in immune and inflammatory responses are all upregulated. Osteoclastogenesis and insulin signaling are also significantly co-regulated in both models and disrupting an osteoclastogenesis-associated gene (OC-STAMP) or the insulin-signaling process significantly altered the host tolerance to C. neoformans. Moreover, C. neoformans was shown to activate metal sequestration, dampen the sugar metabolism, and control cell morphology during infection. Taking together, we described the development of a non-human primate model of cryptococcosis that allowed us to perform an in-depth analysis and comparison of transcriptome profiles during infections of two animal models and conceptually identify host genes important in disease responses. This study provides new insights in understanding fungal pathogenesis mechanisms that potentially facilitate the identification of novel drug targets for the treatment of cryptococcal infection.

The host-pathogen interaction is highly dynamic and tightly regulated, and yet is difficult to monitor. Traditional investigations provide valuable information for the understanding of pathogenic microbial biology but are time-consuming and often neglect the host immune responses. In addition, current animal models for studying pathogenic fungi are limited in mimicking the responses from humans. The development of a new Cryptococcus neoformans infection model using nonhuman primates and the utilization of simultaneous RNA sequencing analysis provide fast and clinically relevant research data allowing the identification of novel critical players from both the invading fungus and the host. The data from the current study would not only help to decipher disease mechanisms but also promote the discovery of novel antifungal drug targets.

Dual species transcript profiling during the interaction between banana (Musa acuminata) and the fungal pathogen Fusarium oxysporum f. sp. cubense

Background

Banana wilt disease, caused by Fusarium oxysporum f. sp. cubense Tropical Race 4 (Foc TR4), is one of the most devastating diseases in banana (Musa spp.). Foc is a soil borne pathogen that causes rot of the roots or wilt of leaves by colonizing the xylem vessels. The dual RNA sequencing is used to simultaneously assess the transcriptomes of pathogen and host. This method greatly helps to understand the responses of pathogen and host to each other and discover the potential pathogenic mechanism.

Results

Plantlets of two economically important banana cultivars, Foc TR4 less susceptible cultivar NK and susceptible cultivar BX, were used to research the Foc-banana interaction mechanism. Notably, the infected NK had more significantly up-regulated genes on the respiration machinery including TCA cycle, glyoxylate, glycerol, and glycolysis compared to BX at 27 h post inoculation (hpi). In addition, genes involved in plant-pathogen interaction, starch, sucrose, linolenic acid and sphingolipid metabolisms were uniquely more greatly induced in BX than those in NK during the whole infection. Genes related to the biosynthesis and metabolism of SA and JA were greatly induced in the infected NK; while auxin and abscisic acid metabolisms related genes were strongly stimulated in the infected BX at 27 hpi. Furthermore, most of fungal genes were more highly expressed in the roots of BX than in those of NK. The fungal genes related to pathogenicity, pectin and chitin metabolism, reactive oxygen scavenging played the important roles during the infection of Foc. CCP1 (cytochrome c peroxidase 1) was verified to involve in cellulose utilization, oxidative stress response and pathogenicity of fungus.

Conclusion

The transcriptome indicated that NK had much faster defense response against Foc TR4 than BX and the expression levels of fungal genes were higher in BX than those in NK. The metabolisms of carbon, nitrogen, and signal transduction molecular were differentially involved in pathogen infection in BX and NK. Additionally, the putative virulence associated fungal genes involved in colonization, nutrition acquirement and transport provided more insights into the infection process of Foc TR4 in banana roots.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5902-z) contains supplementary material, which is available to authorized users.