Dual RNA-seq identifies human mucosal immunity protein Mucin-13 as a hallmark of Plasmodium exoerythrocytic infection

Transcriptome profiling of intrahepatocytic Plasmodium and their host hepatocytes based on the infection phase and the zonation of the liver

Intrahepatocytic development is a key stage in human Plasmodium infection, in which sporozoites replicate and transform into merozoites. Due to technological limitations, however, previous gene expression studies on malaria parasite liver infection were mostly conducted in vitro. In order to bridge these gaps, our current study compared the gene expression of in vitro-infected parasites at different time points with that of in vivo-infected parasites and revealed distinct patterns between parasite subpopulations in vitro and in vivo. A joint investigation of the Plasmodium transcriptome and their host transcriptome was carried out to confer a comprehensive analysis of gene expression in the liver stage of Plasmodium infection in vivo, which is similar to the disease setting, and therefore deepen our understanding of parasite and host transcriptional dynamics during intrahepatocytic infection.

Albumin-Based MUC13 Peptide Nanomedicine Suppresses Liver Cancer Stem Cells via JNK-ERK Signaling Pathway-Mediated Autophagy Inhibition

Targeting liver cancer stem cells (LCSCs) is a promising strategy for hepatocellular carcinoma (HCC) therapy. Target selection and corresponding inhibitor screening are of vital importance for eliminating the stemness of LCSCs. Peptide-based agents are hopeful but have long been hindered for in vivo application. Herein, we selected a clinically significant target MUC13 and screened out a suitable peptide for preparation of an albumin-based MUC13 peptide nanomedicine, P3@HSA, which suppressed liver cancer stem cells via JNK-ERK signaling pathway-mediated autophagy inhibition. The selected target MUC13 was highly expressed in LCSCs and associated with the prognosis of liver cancer patients. Encouraged by this observation, we screened the corresponding peptide-based inhibitor P3 for further evaluation. P3 could interact with albumin through the intrinsic hydrophobic force and formed the nanomedicine P3@HSA. The prepared nanomedicine could inhibit LCSCs through JNK-ERK signaling pathway-mediated autophagy inhibition and exert potent antitumor effect both in vitro and in vivo. Together, this study provides a promising peptide-based nanomedicine for high-performance HCC treatment.

Exploiting integrative metabolomics to study host-parasite interactions in Plasmodium infections

Despite years of research, malaria remains a significant global health burden, with poor diagnostic tests and increasing antimalarial drug resistance challenging diagnosis and treatment. While 'single-omics'-based approaches have been instrumental in gaining insight into the biology and pathogenicity of the Plasmodium parasite and its interaction with the human host, a more comprehensive understanding of malaria pathogenesis can be achieved through 'multi-omics' approaches. Integrative methods, which combine metabolomics, lipidomics, transcriptomics, and genomics datasets, offer a holistic systems biology approach to studying malaria. This review highlights recent advances, future directions, and challenges involved in using integrative metabolomics approaches to interrogate the interactions between Plasmodium and the human host, paving the way towards targeted antimalaria therapeutics and control intervention methods.

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

IMPORTANCE

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

Epigenetic Regulation and Chromatin Remodeling in Malaria Parasites

Plasmodium falciparum, the human malaria parasite, infects two hosts and various cell types, inducing distinct morphological and physiological changes in the parasite in response to different environmental conditions. These variations required the parasite to adapt and develop elaborate molecular mechanisms to ensure its spread and transmission. Recent findings have significantly improved our understanding of the regulation of gene expression in P. falciparum. Here, we provide an up-to-date overview of technologies used to highlight the transcriptomic adjustments occurring in the parasite throughout its life cycle. We also emphasize the complementary and complex epigenetic mechanisms regulating gene expression in malaria parasites. This review concludes with an outlook on the chromatin architecture, the remodeling systems, and how this 3D genome organization is critical in various biological processes.

MUC13 Cell Surface Mucin Limits Salmonella Typhimurium Infection by Protecting the Mucosal Epithelial Barrier

BACKGROUND & AIMS

MUC13 cell surface mucin is highly expressed on the mucosal surface throughout the intestine, yet its role against bacterial infection is unknown. We investigated how MUC13 impacts Salmonella typhimurium (S Tm) infection and elucidated its mechanisms of action.

METHODS

Muc13-/- and wild-type littermate mice were gavaged with 2 isogenic strains of S Tm after pre-conditioning with streptomycin. We assessed clinical parameters, cecal histology, local and systemic bacterial load, and proinflammatory cytokines after infection. Cecal enteroids and epithelial cell lines were used to evaluate the mechanism of MUC13 activity after infection. The interaction between bacterial SiiE and MUC13 was assessed by using siiE-deficient Salmonella.

RESULTS

S Tm-infected Muc13-/- mice had increased disease activity, histologic damage, and higher local and systemic bacterial loads. Mechanistically, we found that S Tm binds to MUC13 through its giant SiiE adhesin and that MUC13 acts as a pathogen-binding decoy shed from the epithelial cell surface after pathogen engagement, limiting bacterial invasion. In addition, MUC13 reduces epithelial cell death and intestinal barrier breakdown by enhancing nuclear factor kappa B signaling during infection, independent of its decoy function.

CONCLUSIONS

We show for the first time that MUC13 plays a critical role in antimicrobial defense against pathogenic S Tm at the intestinal mucosal surface by both acting as a releasable decoy limiting bacterial invasion and reducing pathogen-induced cell death. This further implicates the cell surface mucin family in mucosal defense from bacterial infection.

Human nuclear hormone receptor activity contributes to malaria parasite liver stage development

Chemical genetic approaches have had a transformative impact on discovery of drug targets for malaria but have primarily been used for parasite targets. To identify human pathways required for intrahepatic development of parasite, we implemented multiplex cytological profiling of malaria infected hepatocytes treated with liver stage active compounds. Some compounds, including MMV1088447 and MMV1346624, exhibited profiles similar to cells treated with nuclear hormone receptor (NHR) agonist/antagonists. siRNAs targeting human NHRs, or their signaling partners identified eight genes that were critical for Plasmodium berghei infection. Knockdown of NR1D2, a host NHR, significantly impaired parasite growth by downregulation of host lipid metabolism. Importantly, treatment with MMV1088447 and MMV1346624 but not other antimalarials, phenocopied the lipid metabolism defect of NR1D2 knockdown. Our data underlines the use of high-content imaging for host-cellular pathway deconvolution, highlights host lipid metabolism as a drug-able human pathway and provides new chemical biology tools for studying host-parasite interactions.

Mucin 13 (MUC13) as a candidate biomarker for ovarian cancer detection: potential to complement CA125 in detecting non-serous subtypes

OBJECTIVES

Ovarian cancer is the most lethal gynecological malignancy in developed countries. One of the key associations with the high mortality rate is diagnosis at late stages. This clinical limitation is primarily due to a lack of distinct symptoms and detection at the early stages. The ovarian cancer biomarker, CA125, is mainly effective for identifying serous ovarian carcinomas, leaving a gap in non-serous ovarian cancer detection. Mucin 13 (MUC13) is a transmembrane, glycosylated protein with aberrant expression in malignancies, including ovarian cancer. We explored the potential of MUC13 to complement CA125 as an ovarian cancer biomarker, by evaluating its ability to discriminate serous and non-serous subtypes of ovarian cancer at FIGO stages I-IV from benign conditions.

METHODS

We used our newly developed, high sensitivity ELISA to measure MUC13 protein in a large, well-defined cohort of 389 serum samples from patients with ovarian cancer and benign conditions.

RESULTS

MUC13 and CA125 serum levels were elevated in malignant compared to benign cases (p<0.0001). Receiver-operating characteristic (ROC) curve analysis showed similar area under the curve (AUC) of 0.74 (MUC13) and 0.76 (CA125). MUC13 concentrations were significantly higher in mucinous adenocarcinomas compared to benign controls (p=0.0005), with AUC of 0.80. MUC13 and CA125 showed significant elevation in early-stage cases (stage I-II) in relation to benign controls (p=0.0012 and p=0.014, respectively).

CONCLUSIONS

We report the novel role of MUC13 as a serum ovarian cancer biomarker, where it could complement CA125 for detecting some subtypes of non-serous ovarian carcinoma and early-stage disease.

Transcriptional profiling of hepatocytes infected with the replicative form of the malaria parasite Plasmodium cynomolgi

BACKGROUND

The zoonotic simian parasite Plasmodium cynomolgi develops into replicating schizonts and dormant hypnozoites during the infection of hepatocytes and is used as a model organism to study relapsing malaria. The transcriptional profiling of P. cynomolgi liver stages was previously reported and revealed many important biological features of the parasite but left out the host response to malaria infection.

METHODS

Previously published RNA sequencing data were used to quantify the expression of host genes in rhesus macaque hepatocytes infected with P. cynomolgi in comparison to either cells from uninfected samples or uninfected bystander cells.

RESULTS

Although the dataset could not be used to resolve the transcriptional profile of hypnozoite-infected hepatocytes, it provided a snapshot of the host response to liver stage schizonts at 9-10 day post-infection and identified specific host pathways that are modulated during the exo-erythrocytic stage of P. cynomolgi.

CONCLUSIONS

This study constitutes a valuable resource characterizing the hepatocyte response to P. cynomolgi infection and provides a framework to build on future research that aims at understanding hepatocyte-parasite interactions during relapsing malaria infection.

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.

Integrated analysis revealing the role of TET3-mediated MUC13 promoter hypomethylation in hepatocellular carcinogenesis

Aim: To explore the function and underlying mechanism of MUC13 in hepatocellular carcinoma (HCC) oncogenesis. Materials & Methods: Online databases and software were used to perform analyses of expression, methylation and enrichment pathway. Experiments were performed to confirm the results using HCC cells in vitro. Results: MUC13 was upregulated in HCC and liver cancer stem cells (CSCs) and had a positive influence on CSC generation. Further analyses revealed that MUC13 with promoter hypomethylated was regulated by DNA demethylase TET3, which was overexpressed in HCC and liver CSCs. Conclusion: These results strongly suggested that high TET3 expression in liver CSCs may mediate MUC13 upregulation via promoter hypomethylation and thereby contribute to hepatocellular carcinogenesis.

Malaria Resilience in South America: Epidemiology, Vector Biology, and Immunology Insights from the Amazonian International Center of Excellence in Malaria Research Network in Peru and Brazil

The 1990s saw the rapid reemergence of malaria in Amazonia, where it remains an important public health priority in South America. The Amazonian International Center of Excellence in Malaria Research (ICEMR) was designed to take a multidisciplinary approach toward identifying novel malaria control and elimination strategies. Based on geographically and epidemiologically distinct sites in the Northeastern Peruvian and Western Brazilian Amazon regions, synergistic projects integrate malaria epidemiology, vector biology, and immunology. The Amazonian ICEMR's overarching goal is to understand how human behavior and other sociodemographic features of human reservoirs of transmission-predominantly asymptomatically parasitemic people-interact with the major Amazonian malaria vector, Nyssorhynchus (formerly Anopheles) darlingi, and with human immune responses to maintain malaria resilience and continued endemicity in a hypoendemic setting. Here, we will review Amazonian ICEMR's achievements on the synergies among malaria epidemiology, Plasmodium-vector interactions, and immune response, and how those provide a roadmap for further research, and, most importantly, point toward how to achieve malaria control and elimination in the Americas.

In vitro models for human malaria: targeting the liver stage

The Plasmodium liver stage represents a vulnerable therapeutic target to prevent disease progression as the parasite resides in the liver before clinical representation caused by intraerythrocytic development. However, most antimalarial drugs target the blood stage of the parasite's life cycle, and the few drugs that target the liver stage are lethal to patients with a glucose-6-phosphate dehydrogenase deficiency. Furthermore, implementation of in vitro liver models to study and develop novel therapeutics against the liver stage of human Plasmodium species remains challenging. In this review, we focus on the progression of in vitro liver models developed for human Plasmodium spp. parasites, provide a brief review on important assay requirements, and lastly present recommendations to improve models to enhance the discovery process of novel preclinical therapeutics.

Single-cell RNA profiling of Plasmodium vivax-infected hepatocytes reveals parasite- and host- specific transcriptomic signatures and therapeutic targets

The resilience of Plasmodium vivax, the most widely-distributed malaria-causing parasite in humans, is attributed to its ability to produce dormant liver forms known as hypnozoites, which can activate weeks, months, or even years after an initial mosquito bite. The factors underlying hypnozoite formation and activation are poorly understood, as is the parasite's influence on the host hepatocyte. Here, we shed light on transcriptome-wide signatures of both the parasite and the infected host cell by sequencing over 1,000 P. vivax-infected hepatocytes at single-cell resolution. We distinguish between replicating schizonts and hypnozoites at the transcriptional level, identifying key differences in transcripts encoding for RNA-binding proteins associated with cell fate. In infected hepatocytes, we show that genes associated with energy metabolism and antioxidant stress response are upregulated, and those involved in the host immune response downregulated, suggesting both schizonts and hypnozoites alter the host intracellular environment. The transcriptional markers in schizonts, hypnozoites, and infected hepatocytes revealed here pinpoint potential factors underlying dormancy and can inform therapeutic targets against P. vivax liver-stage infection.

Elucidating Spatially-Resolved Changes in Host Signaling During Plasmodium Liver-Stage Infection

Upon transmission to the human host, Plasmodium sporozoites exit the skin, are taken up by the blood stream, and then travel to the liver where they infect and significantly modify a single hepatocyte. Low infection rates within the liver have made proteomic studies of infected hepatocytes challenging, particularly in vivo, and existing studies have been largely unable to consider how protein and phosphoprotein differences are altered at different spatial locations within the heterogeneous liver. Using digital spatial profiling, we characterized changes in host signaling during Plasmodium yoelii infection in vivo without disrupting the liver tissue. Moreover, we measured alterations in protein expression around infected hepatocytes and identified a subset of CD163⁺ Kupffer cells that migrate towards infected cells during infection. These data offer the first insight into the heterogeneous microenvironment that surrounds the infected hepatocyte and provide insights into how the parasite may alter its milieu to influence its survival and modulate immunity.

Contribution of Transcriptome to Elucidate the Biology of Plasmodium spp

In the present review, we discuss some of the new technologies that have been applied to elucidate how Plasmodium spp escape from the immune system and subvert the host physiology to orchestrate the regulation of its biological pathways. Our manuscript describes how techniques such as microarray approaches, RNA-Seq and single-cell RNA sequencing have contributed to the discovery of transcripts and changed the concept of gene expression regulation in closely related malaria parasite species. Moreover, the text highlights the contributions of high-throughput RNA sequencing for the current knowledge of malaria parasite biology, physiology, vaccine target and the revelation of new players in parasite signaling.

Decrypting the complexity of the human malaria parasite biology through systems biology approaches

The human malaria parasite, Plasmodium falciparum, is a unicellular protozoan responsible for over half a million deaths annually. With a complex life cycle alternating between human and invertebrate hosts, this apicomplexan is notoriously adept at evading host immune responses and developing resistance to all clinically administered treatments. Advances in omics-based technologies, increased sensitivity of sequencing platforms and enhanced CRISPR based gene editing tools, have given researchers access to more in-depth and untapped information about this enigmatic micro-organism, a feat thought to be infeasible in the past decade. Here we discuss some of the most important scientific achievements made over the past few years with a focus on novel technologies and platforms that set the stage for subsequent discoveries. We also describe some of the systems-based methods applied to uncover gaps of knowledge left through single-omics applications with the hope that we will soon be able to overcome the spread of this life-threatening disease.

MUC13 promotes lung cancer development and progression by activating ERK signaling

Mucin 13 (MUC13) is a glycoprotein that is expressed on the cell surface and participates in the tumorigenesis of multiple malignancies, including pancreatic cancer, colorectal cancer and renal cancer. However, to the best of our knowledge, the expression levels and function of MUC13 in lung cancer progression have not yet been demonstrated. Therefore, the present study examined the expression pattern and regulatory role of MUC13 in lung cancer tumorigenesis. The results demonstrated that MUC13 was highly expressed in lung cancer tissues and cell lines compared with that in normal tissues and cell lines. Functionally, knockdown of MUC13 inhibited cell proliferation and enhanced the apoptosis of A549 and NCI-H1650 lung cancer cells. Furthermore, silencing of MUC13 suppressed the migration and invasion of lung cancer cells. Additionally, a xenograft tumor model demonstrated that knockdown of MUC13 delayed the development of the lung cancer xenograft and suppressed the expression of proliferation marker Ki-67 in tumor tissues. Mechanistically, MUC13 activated the ERK signaling pathway by enhancing the phosphorylation of ERK, JNK and p38 in lung cancer tissues compared with that in normal tissues. Knockdown of MUC13 inhibited the phosphorylation of ERK/JNK/p38 in A549 and NCI-H1650 cells. Overall, these findings suggested that MUC13 could act as an oncogenic glycoprotein to accelerate the progression of lung cancer via abnormal activation of the ERK/JNK/p38 signaling pathway and might serve as a therapeutic target for lung cancer treatment.

Analysis of pir gene expression across the Plasmodium life cycle

Background

Plasmodium interspersed repeat (pir) is the largest multigene family in the genomes of most Plasmodium species. A variety of functions for the PIR proteins which they encode have been proposed, including antigenic variation, immune evasion, sequestration and rosetting. However, direct evidence for these is lacking. The repetitive nature of the family has made it difficult to determine function experimentally. However, there has been some success in using gene expression studies to suggest roles for some members in virulence and chronic infection.

Methods

Here pir gene expression was examined across the life cycle of Plasmodium berghei using publicly available RNAseq data-sets, and at high resolution in the intraerythrocytic development cycle using new data from Plasmodium chabaudi.

Results

Expression of pir genes is greatest in stages of the parasite which invade and reside in red blood cells. The marked exception is that liver merozoites and male gametocytes produce a very large number of pir gene transcripts, notably compared to female gametocytes, which produce relatively few. Within the asexual blood stages different subfamilies peak at different times, suggesting further functional distinctions. Representing a subfamily of its own, the highly conserved ancestral pir gene warrants further investigation due to its potential tractability for functional investigation. It is highly transcribed in multiple life cycle stages and across most studied Plasmodium species and thus is likely to play an important role in parasite biology.

Conclusions

The identification of distinct expression patterns for different pir genes and subfamilies is likely to provide a basis for the design of future experiments to uncover their function.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03979-6.

Investigation of the transcriptomic response in Atlantic salmon (Salmo salar) gill exposed to Paramoeba perurans during early onset of disease

Amoebic Gill Disease (AGD), caused by the protozoan extracellular parasite Paramoeba perurans (P. perurans) is a disease affecting Atlantic salmon (Salmo salar). This study investigated the gill transcriptomic profile of pre-clinical AGD using RNA-sequencing (RNA-seq) technology. RNA-seq libraries generated at 0, 4, 7, 14 and 16 days post infection (dpi) identified 19,251 differentially expressed genes (DEGs) of which 56.2% were up-regulated. DEGs mapped to 224 Gene Ontology (GO) terms including 140 biological processes (BP), 45 cellular components (CC), and 39 molecular functions (MF). A total of 27 reference pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) and 15 Reactome gene sets were identified. The RNA-seq data was validated using real-time, quantitative PCR (qPCR). A host immune response though the activation of complement and the acute phase genes was evident at 7 dpi, with a concurrent immune suppression involving cytokine signalling, notably in interleukins, interferon regulatory factors and tumour necrosis factor-alpha (tnf-α) genes. Down-regulated gene expression with involvement in receptor signalling pathways (NOD-like, Toll-like and RIG-1) were also identified. The results of this study support the theory that P. perurans can evade immune surveillance during the initial stages of gill colonisation through interference of signal transduction pathways.

Immunoprofiles associated with controlled human malaria infection and naturally acquired immunity identify a shared IgA pre-erythrocytic immunoproteome

Knowledge of the Plasmodium falciparum antigens that comprise the human liver stage immunoproteome is important for pre-erythrocytic vaccine development, but, compared with the erythrocytic stage immunoproteome, more challenging to classify. Previous studies of P. falciparum antibody responses report IgG and rarely IgA responses. We assessed IgG and IgA antibody responses in adult sera collected during two controlled human malaria infection (CHMI) studies in malaria-naïve volunteers and in 1- to 6-year-old malaria-exposed Malian children on a 251 P. falciparum antigen protein microarray. IgG profiles in the two CHMI groups were equivalent and differed from Malian children. IgA profiles were robust in the CHMI groups and a subset of Malian children. We describe immunoproteome differences in naïve vs. exposed individuals and report pre-erythrocytic proteins recognized by the immune system. IgA responses detected in this study expand the list of pre-erythrocytic antigens for further characterization as potential vaccine candidates.

An Open-Source Toolkit To Expand Bioinformatics Training in Infectious Diseases

As access to high-throughput sequencing technology has increased, the bottleneck in biomedical research has shifted from data generation to data analysis. Here, we describe a modular and extensible framework for didactic instruction in bioinformatics using publicly available RNA sequencing data sets from infectious disease studies, with a focus on host-parasite interactions. We highlight lessons learned from adapting this course for virtual learners during the coronavirus disease 2019 (COVID-19) pandemic.

Characterization of the Tubovesicular Network in Plasmodium vivax Liver Stage Hypnozoites and Schizonts

Plasmodium is a genus of apicomplexan parasites which replicate in the liver before causing malaria. Plasmodium vivax can also persist in the liver as dormant hypnozoites and cause clinical relapse upon activation, but the molecular mechanisms leading to activation have yet to be discovered. In this study, we use high-resolution microscopy to characterize temporal changes of the P. vivax liver stage tubovesicular network (TVN), a parasitophorous vacuole membrane (PVM)-derived network within the host cytosol. We observe extended membrane clusters, tubules, and TVN-derived vesicles present throughout P. vivax liver stage development. Additionally, we demonstrate an unexpected presence of the TVN in hypnozoites and observe some association of this network to host nuclei. We also reveal that the host water and solute channel aquaporin-3 (AQP3) associates with TVN-derived vesicles and extended membrane clusters. AQP3 has been previously shown to localize to the PVM of P. vivax hypnozoites and liver schizonts but has not yet been shown in association to the TVN. Our results highlight host-parasite interactions occur in both dormant and replicating liver stage P. vivax forms and implicate AQP3 function during this time. Together, these findings enhance our understanding of P. vivax liver stage biology through characterization of the TVN with an emphasis on the presence of this network in dormant hypnozoites.

Host-Parasite Interaction of Atlantic salmon (Salmo salar) and the Ectoparasite Neoparamoeba perurans in Amoebic Gill Disease

Marine farmed Atlantic salmon (Salmo salar) are susceptible to recurrent amoebic gill disease (AGD) caused by the ectoparasite Neoparamoeba perurans over the growout production cycle. The parasite elicits a highly localized response within the gill epithelium resulting in multifocal mucoid patches at the site of parasite attachment. This host-parasite response drives a complex immune reaction, which remains poorly understood. To generate a model for host-parasite interaction during pathogenesis of AGD in Atlantic salmon the local (gill) and systemic transcriptomic response in the host, and the parasite during AGD pathogenesis was explored. A dual RNA-seq approach together with differential gene expression and system-wide statistical analyses of gene and transcription factor networks was employed. A multi-tissue transcriptomic data set was generated from the gill (including both lesioned and non-lesioned tissue), head kidney and spleen tissues naïve and AGD-affected Atlantic salmon sourced from an in vivo AGD challenge trial. Differential gene expression of the salmon host indicates local and systemic upregulation of defense and immune responses. Two transcription factors, znfOZF-like and znf70-like, and their associated gene networks significantly altered with disease state. The majority of genes in these networks are candidates for mediators of the immune response, cellular proliferation and invasion. These include Aurora kinase B-like, rho guanine nucleotide exchange factor 25-like and protein NDNF-like inhibited. Analysis of the N. perurans transcriptome during AGD pathology compared to in vitro cultured N. perurans trophozoites, as a proxy for wild type trophozoites, identified multiple gene candidates for virulence and indicates a potential master regulatory gene system analogous to the two-component PhoP/Q system. Candidate genes identified are associated with invasion of host tissue, evasion of host defense mechanisms and formation of the mucoid lesion. We generated a novel model for host-parasite interaction during AGD pathogenesis through integration of host and parasite functional profiles. Collectively, this dual transcriptomic study provides novel molecular insights into the pathology of AGD and provides alternative theories for future research in a step towards improved management of AGD.

Dual RNA Sequencing Meta-analysis in Plasmodium Infection Identifies Host-Parasite Interactions

Dual RNA sequencing (RNA-Seq) is the simultaneous transcriptomic analysis of interacting symbionts, for example, in malaria. Potential cross-species interactions identified by correlated gene expression might highlight interlinked signaling, metabolic, or gene regulatory pathways in addition to physically interacting proteins. Often, malaria studies address one of the interacting organisms-host or parasite-rendering the other "contamination." Here we perform a meta-analysis using such studies for cross-species expression analysis. We screened experiments for gene expression from host and Plasmodium. Out of 171 studies in Homo sapiens, Macaca mulatta, and Mus musculus, we identified 63 potential studies containing host and parasite data. While 16 studies (1,950 samples) explicitly performed dual RNA-Seq, 47 (1,398 samples) originally focused on one organism. We found 915 experimental replicates from 20 blood studies to be suitable for coexpression analysis and used orthologs for meta-analysis across different host-parasite systems. Centrality metrics from the derived gene expression networks correlated with gene essentiality in the parasites. We found indications of host immune response to elements of the Plasmodium protein degradation system, an antimalarial drug target. We identified well-studied immune responses in the host with our coexpression networks, as our approach recovers known broad processes interlinked between hosts and parasites in addition to individual host and parasite protein associations. The set of core interactions represents commonalities between human malaria and its model systems for prioritization in laboratory experiments. Our approach might also allow insights into the transferability of model systems for different pathways in malaria studies.IMPORTANCE Malaria still causes about 400,000 deaths a year and is one of the most studied infectious diseases. The disease is studied in mice and monkeys as lab models to derive potential therapeutic intervention in human malaria. Interactions between Plasmodium spp. and its hosts are either conserved across different host-parasite systems or idiosyncratic to those systems. Here we use correlation of gene expression from different RNA-Seq studies to infer common host-parasite interactions across human, mouse, and monkey studies. First, we find a set of very conserved interactors, worth further scrutiny in focused laboratory experiments. Second, this work might help assess to which extent experiments and knowledge on different pathways can be transferred from models to humans for potential therapy.

Modeling Relapsing Malaria: Emerging Technologies to Study Parasite-Host Interactions in the Liver

Recent studies of liver stage malaria parasite-host interactions have provided exciting new insights on the cross-talk between parasite and its mammalian (predominantly rodent) host. We review the latest state of the art and and zoom in on new technologies that will provide the tools necessary to investigate host-parasite interactions of relapsing parasites. Interactions between hypnozoites and hepatocytes are particularly interesting because the parasite can remain in a quiescent state for prolonged periods of time and triggers for reactivation have not been irrefutably identified. If we learn more about the cross-talk between hypnozoite and host we may be able to identify factors that encourage waking up these dormant parasite reservoirs and help to achieve the total eradication of malaria.

Role of Heat Shock Proteins in Immune Modulation in Malaria

Malaria is one of the major parasitic killer diseases worldwide. Severe cases of malaria are mostly in children under the age of 5 years due to their naïve immune system and in pregnant women with weakened immune responses. Inflammatory immune responses against the parasite involve complement activation as well as the antibody and effector cell-mediated immune system. However, after an infection with Plasmodium falciparum (P. falciparum), the most dangerous malaria species, the host-derived immunity is often insufficient to completely inhibit the infection cycles of the parasite in red blood cells for yet unknown reasons. In the present chapter we aim to elucidate the role of the host's and the parasite's heat shock proteins (HSPs) in the development of a novel anti-malaria therapeutic approach.

Evolutionarily Selected Overexpression of the Cytokine BAFF Enhances Mucosal Immune Response Against P. falciparum

We have previously shown that a variant of the TNFSF13B gene that we called BAFF-var increases the production of the cytokine BAFF, upregulating humoral immunity and increasing the risk for certain autoimmune diseases. In addition, genetic population signatures revealed that BAFF-var was evolutionarily advantageous, most likely by increasing resistance to malaria infection, which is a prime candidate for selective pressure. To evaluate whether the increased soluble BAFF (sBAFF) production confers protection, we experimentally assessed the role of BAFF-var in response to malaria antigens. Lysates of erythrocytes infected with Plasmodium falciparum (iRBCs) or left uninfected (uRBCs, control) were used to treat peripheral blood mononuclear cells (PBMCs) with distinct BAFF genotypes. The PBMCs purified from BAFF-var donors and treated with iRBCs showed different levels of specific cells, immunoglobulins, and cytokines as compared with BAFF-WT. In particular, a relevant differential effect on mucosal immunity B subpopulations have been observed. These findings point to specific immune cells and molecules through which the evolutionary selected BAFF-var may have improved fitness during P. falciparum infection.

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.

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.

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.

RNA-Seq Analysis Illuminates the Early Stages of Plasmodium Liver Infection

The LS of Plasmodium infection is an asymptomatic yet necessary stage for producing blood-infective parasites, the causative agents of malaria. Blocking the liver stage of the life cycle can prevent clinical malaria, but relatively less is known about the parasite’s biology at this stage. Using the rodent model P. berghei, we investigated whole-transcriptome changes occurring as early as 2 hpi of hepatocytes. The transcriptional profiles of early time points (2, 4, 12, and 18 hpi) have not been accessible before due to the technical challenges associated with liver-stage infections. Our data now provide insights into these early parasite fluxes that may facilitate establishment of infection, transformation, and replication in the liver.

The apicomplexan parasites Plasmodium spp. are the causative agents of malaria, a disease that poses a significant global health burden. Plasmodium spp. initiate infection of the human host by transforming and replicating within hepatocytes. This liver stage (LS) is poorly understood compared to other Plasmodium life stages, which has hindered our ability to target these parasites for disease prevention. We conducted an extensive transcriptome sequencing (RNA-Seq) analysis throughout the Plasmodium berghei LS, covering as early as 2 h postinfection (hpi) and extending to 48 hpi. Our data revealed that hundreds of genes are differentially expressed at 2 hpi and that multiple genes shown to be important for later infection are upregulated as early as 12 hpi. Using hierarchical clustering along with coexpression analysis, we identified clusters functionally enriched for important liver-stage processes such as interactions with the host cell and redox homeostasis. Furthermore, some of these clusters were highly correlated to the expression of ApiAP2 transcription factors, while showing enrichment of mostly uncharacterized DNA binding motifs. This finding indicates potential LS targets for these transcription factors, while also hinting at alternative uncharacterized DNA binding motifs and transcription factors during this stage. Our work presents a window into the previously undescribed transcriptome of Plasmodium upon host hepatocyte infection to enable a comprehensive view of the parasite’s LS. These findings also provide a blueprint for future studies that extend hypotheses concerning LS gene function in P. berghei to human-infective Plasmodium parasites.

Organoids for Liver Stage Malaria Research

Plasmodium parasites cause malaria and are maintained between Anopheles mosquitoes and mammalian hosts in a complex life cycle. Malaria parasites occupy tissue niches that can be difficult to access, and models to study them can be challenging to recapitulate experimentally, particularly for Plasmodium species that infect humans. 2D culture models provide extremely beneficial tools to investigate Plasmodium biology but they have limitations. More complex 3D structural networks, such as organoids, have unveiled new avenues for developing more physiological tissue models, and their application to malaria research offers great promise. Here, we review current models for studying Plasmodium infection with a key focus on the obligate pre-erythrocytic stage that culminates in blood infection, causing malaria, and discuss how organoids should fulfil an important and unmet need.

Advances in omics-based methods to identify novel targets for malaria and other parasitic protozoan infections

A major advance in antimalarial drug discovery has been the shift towards cell-based phenotypic screening, with notable progress in the screening of compounds against the asexual blood stage, liver stage, and gametocytes. A primary method for drug target deconvolution in Plasmodium falciparum is in vitro evolution of compound-resistant parasites followed by whole-genome scans. Several of the most promising antimalarial drug targets, such as translation elongation factor 2 (eEF2) and phenylalanine tRNA synthetase (PheRS), have been identified or confirmed using this method. One drawback of this method is that if a mutated gene is uncharacterized, a substantial effort may be required to determine whether it is a drug target, a drug resistance gene, or if the mutation is merely a background mutation. Thus, the availability of high-throughput, functional genomic datasets can greatly assist with target deconvolution. Studies mapping genome-wide essentiality in P. falciparum or performing transcriptional profiling of the host and parasite during liver-stage infection with P. berghei have identified potentially druggable pathways. Advances in mapping the epigenomic regulation of the malaria parasite genome have also enabled the identification of key processes involved in parasite development. In addition, the examination of the host genome during infection has identified novel gene candidates associated with susceptibility to severe malaria. Here, we review recent studies that have used omics-based methods to identify novel targets for interventions against protozoan parasites, focusing on malaria, and we highlight the advantages and limitations of the approaches used. These approaches have also been extended to other protozoan pathogens, including Toxoplasma, Trypanosoma, and Leishmania spp., and these studies highlight how drug discovery efforts against these pathogens benefit from the utilization of diverse omics-based methods to identify promising drug targets.

Transcriptome analysis of Plasmodium berghei during exo-erythrocytic development

BACKGROUND

The complex life cycle of malaria parasites requires well-orchestrated stage specific gene expression. In the vertebrate host the parasites grow and multiply by schizogony in two different environments: within erythrocytes and within hepatocytes. Whereas erythrocytic parasites are well-studied in this respect, relatively little is known about the exo-erythrocytic stages.

METHODS

In an attempt to fill this gap, genome wide RNA-seq analyses of various exo-erythrocytic stages of Plasmodium berghei including sporozoites, samples from a time-course of liver stage development and detached cells were performed. These latter contain infectious merozoites and represent the final step in exo-erythrocytic development.

RESULTS

The analysis represents the complete transcriptome of the entire life cycle of P. berghei parasites with temporal detailed analysis of the liver stage allowing comparison of gene expression across the progression of the life cycle. These RNA-seq data from different developmental stages were used to cluster genes with similar expression profiles, in order to infer their functions. A comparison with published data from other parasite stages confirmed stage-specific gene expression and revealed numerous genes that are expressed differentially in blood and exo-erythrocytic stages. One of the most exo-erythrocytic stage-specific genes was PBANKA_1003900, which has previously been annotated as a "gametocyte specific protein". The promoter of this gene drove high GFP expression in exo-erythrocytic stages, confirming its expression profile seen by RNA-seq.

CONCLUSIONS

The comparative analysis of the genome wide mRNA expression profiles of erythrocytic and different exo-erythrocytic stages could be used to improve the understanding of gene regulation in Plasmodium parasites and can be used to model exo-erythrocytic stage metabolic networks toward the identification of differences in metabolic processes during schizogony in erythrocytes and hepatocytes.