Genome-wide liver transcriptomic profiling of a malaria mouse model reveals disturbed immune and metabolic responses
Parasites & Vectors volume 16, Article number: 40 (2023)
The liver is responsible for a range of functions in vertebrates, such as metabolism and immunity. In malaria, the liver plays a crucial role in the interaction between the parasite and host. Although malarial hepatitis is a common clinical complication of severe malaria, other malaria-related liver changes have been overlooked during the blood stage of the parasite life-cycle, in contrast to the many studies that have focused on parasite invasion of and replication in the liver during the hepatic stage of the parasite.
A rodent model of malaria was established using Plasmodium yoelii strain 17XL, a lethal strain of rodent malaria, for liver transcriptomic profiling.
Differentially expressed messenger RNAs were associated with innate and adaptive immune responses, while differentially expressed long noncoding RNAs were enriched in the regulation of metabolism-related pathways, such as lipid metabolism. The coexpression network showed that host genes were related to cellular transport and tissue remodeling. Hub gene analysis of P. yoelii indicated that ubiquitination genes that were coexpressed with the host were evolutionarily conserved.
Our analysis yielded evidence of activated immune responses, aberrant metabolic processes and tissue remodeling changes in the livers of mice with malaria during the blood stage of the parasite, which provided a systematic outline of liver responses during Plasmodium infection.
Malaria is a widespread infectious disease that is prevalent mainly in tropical and subtropical regions of the world where it is a major health burden . Systemic changes occur in organs of the infected organism, such as severe anemia and impaired microvascular perfusion, which are central to the development of severe malaria [2,3,4]. Severe malaria in turn is associated with ischemia/hypoxia, metabolic disturbances, hyperlactatemia and disturbance of consciousness, leading to high mortality [5,6,7]. Investigating the involvement and responses of important organs responsible for metabolism and immunity, such as the liver, may contribute to a systematic understanding of the pathogenesis of severe malaria and provide new perspectives on the treatment and care of critically ill patients [8, 9].
The liver is an important organ for maintaining immune and metabolic homeostasis. As part of the immune response, the liver plays a critical role not only in host defense but also in liver injury and repair . During the liver stage of plasmodial infection, bone marrow-derived cells in the liver interact with the parasite, modulating the replication boost and the differentiation of sporozoites into merozoites [11,12,13]. During the blood stage of the infection, the liver is a main target organ for Plasmodium, as demonstrated by the activated immune responses and various liver injuries associated with high parasite loads observed in several malaria mouse models [14,15,16,17]. On the other hand, as a central metabolic organ, the liver integrates systemic metabolism to maintain glucose and fatty acid homeostasis. Impaired hepatic gluconeogenesis and hepatic lactate clearance have been reported to directly affect blood lactate concentrations, resulting in death-related acidosis in patients with malaria . A positive association between high Plasmodium-infected erythrocyte (pRBC) load of the liver and jaundice, hepatomegaly and elevated liver enzymes has been observed in different malaria patient cohorts [18, 19]. Liver dysfunction in these patients is related to other organ dysfunctions and poor outcomes [19, 20], suggesting a potential role of the liver in malaria responses.
Compared with the liver stage of the parasite life-cycle, involvement of the liver in the blood stage of the parasite has been overlooked in most studies, and only a few animal and human studies on hepatic pathologies in malaria have been reported [14,15,16, 21, 22]. Histopathological evidence of the involvement of reactive Kupffer cells, retention of heme pigment, minimal pRBC sequestration and inflammatory responses in the pathogenesis of malaria-induced liver injury has been accumulating [15, 18, 23, 24]. However, systemic identification of molecular events is required to comprehensively demonstrate liver changes and provide a basis for exploring the supportive features of the liver response during infection.
In the study reported here, we explored the global transcriptomic changes in the host liver in Plasmodium yoelii-infected mice. We found that differentially expressed (DE) messenger RNAs (mRNAs) in the liver were mainly responsible for innate and adaptive immune responses. In contrast, DE long noncoding RNAs (lncRNAs) mainly focused on the regulation of fatty acid metabolism and oxidation‒reduction processes. Host genes highly coexpressed with P. yoelii genes were enriched in the processes of cellular transport and tissue remodeling. Additionally, hub genes of P. yoelii in the coexpression network suggested that conserved ubiquitination genes were highly coexpressed with host genes. Collectively, our work reveals that immune responses, metabolic changes and tissue repair are major events that occur in the liver during the blood stage of plasmodial infection and further provides potentially distinct regulatory roles of mRNAs and lncRNAs in liver responses.
The International Guiding Principles for Biomedical Research Involving Animals were strictly followed during animal experiments. Procedures involving vertebrate animals were reviewed and approved by the Animal Care and Use Committee of Peking University Health Science Center with permit number PUIRB-LA2022677.
Establishment of the P. yoelii-mouse model
Four-week-old male ICR/JCL mice (20–25 g) were obtained from the Department of Laboratory Animal Science of Peking University (Beijing, China; permit: SCXK-2012–0015) and housed under specific pathogen-free conditions. Mouse blood containing P. yoelii lethal strain 17XL (70% parasitemia) frozen in liquid nitrogen was thawed quickly for injection. Five mice were infected with 2 × 106 infected red blood cells (RBCs) by the intraperitoneal route, similar to earlier studies [25, 26]. The parasitemia of infected mice was assessed by microscopic examination of Giemsa-stained blood smears post-infection. Five fields were randomly selected from each blood smear for cell counting, and the parasitemia percentage was measured as parasitized RBCs/total RBCs × 100.
RNA extraction and transcriptomic analysis
As oxidative stress and hepatic apoptosis peak at 50%–60% parasitemia , mouse liver tissues were sampled at approximately 50% parasitemia. Total RNA was extracted from liver tissues using TRIzol reagent. The RNA integrity number (RIN) was assessed with an Agilent 2100 Bioanalyzer system (Aglient Technologies, Santa Clara, CA, USA), and extracted RNA with an RIN > 7.0 was retained for sequencing. Ribosomal RNA depletion libraries were prepared and sequenced with IGENECODE (Beijing, China) on the Illumina HiSeq PE150 platform (Illumina, Inc., San Diego, CA, USA).
Read mapping and normalization of RNA sequencing data
RNA sequencing (RNA-seq) data were mapped to the reference genome of Mus musculus (GRCm39) using Hisat2 . Gene expression levels were calculated with Stringtie . Genes of M. musculus with average counts > 50 were kept for differential expression analysis. Circular RNAs (circRNAs) were identified with CIRI2 using the same RNA-seq reference genome for guidance .
Identification and functional enrichment analysis
The R package DEseq2  was used to identify DE RNAs. Genes with adjusted P values < 0.05 and fold changes > 2 were considered to be DE. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of DE genes were carried out with DAVID v6.8 .
Gene coexpression network construction and visualization
An artificial genome was generated by combining the reference genomes of M. musculus (GRCm39) and P. yoelii (GCF_900002385.2) . The same RNA-seq data were mapped to the artificial genome as described above. Mus musculus genes with an average TPM (transcript per million) > 1 were kept, as were Plasmodium yoelii genes that were not expressed in the control group and had an average TPM > 1 in the infected group. The R package WGCNA  was used for gene coexpression network construction. Briefly, filtered and normalized expression data were used for average linkage hierarchical clustering. Network modules were identified using a dynamic tree cut algorithm with a minimum cluster size of 30. The correlation between modules and phenotypes (control and infection) was calculated, and the modules containing P. yoelii genes with a correlation coefficient > 0.8 were selected for downstream GO enrichment analysis. The visualization of special modules was performed with Cytoscape_v3.8.2 .
Homolog identification among Plasmodium spp.
Reference protein sequences of P. yoelii were compared with those of a number of other Plasmodium species, including P. falciparum, P. malariae, P. ovale, P. vivax, P. knowlesi, P. berghei and P. chabaudi, using BLAST 2.10.1 + with a parameter e value = 1e10. Only proteins with > 90% overlapping sequences with other proteins and with overlapping sequences accounting for > 90% of the two protein sequences were considered to be homologous proteins.
Hematoxylin–eosin staining of liver
Liver hematoxylin–eosin (HE) staining of infected and control mice was performed by Wuhan Servicebio Technology Co., Ltd. In brief, paraffin-embedded sections were deparaffinized and rehydrated using ethanol and xylene. HE staining was performed to evaluate inflamed cell infiltration and hepatocellular edema/necrosis/fibrosis.
Real-time quantitative PCR
A total of 1 μg RNA was reverse transcribed into complementary DNA (cDNA) using the HiScript III 1st Strand cDNA Synthesis Kit (+ gDNA wiper) (Vazyme, Nanjing, China) according to the manufacturer's instructions. Then, real-time quantitative PCR (RT‒qPCR) analysis was performed using the Hieff UNICON® Universal Blue qPCR SYBR Green Master Mix (Yeasen Biotechnology Co., Ltd., Shanghai, China) according to the manufacturer's instructions. The cycling program consisted of an initial incubation at 95 °C for 2 min, followed by 40 cycles of 10 s at 95 °C and 30 s at 60 °C. Each gene was analyzed in biological triplicates and technical triplicates with the glyceraldehyde 3-phosphate dehydrogenase gene (GAPDH) as the internal reference. All primers used in RT‒qPCR assays are listed in Additional file 1: Table S1).
Differential expression profiling of whole transcriptomes between infected and control groups
The P. yoelii-rodent model described in section Establishment of the P. yoelii-mouse model was established with noninoculated mice as the control group. Giemsa-stained blood smears of all five infected mice were examined until parasitemia reached 50% (Additional file 2: Figure S1). Liver tissues were resected from five infected mice and three controls for histological examination and RNA extraction. Hepatic sinusoidal expansion, hepatic cell edema, irregular arrangement of hepatic cords, isolated lipid droplets and the presence of macrophages engulfing a large number of malaria pigment particles were observed in the HE-stained liver sections of infected mice (Additional file 3: Figure S2). Total RNA was extracted from liver tissues for ribosomal RNA (rRNA) depletion high-throughput sequencing, and an average of 7.6 Gb of data were generated from each sample (Additional file 4: Table S2). RNA-seq of infected and control samples showed a high correlation within and between groups, indicating favorable data quality (Fig. 1a). Principal component analysis (PCA) indicated that biological replicates clustered together, while strong separation among infected and control samples distinguished the transcriptomes of the two groups (Fig. 1b).
To further explore the differences in the transcriptomes between the infected and control groups, we identified 4875 DE genes among 19,066 stably expressed genes using DEseq2 with thresholds of fold change > 2 and adjusted P value < 0.05. A total of 3975 mRNAs and 493 lncRNAs were found to be DE, accounting for > 90% of DE genes (Fig. 1c). In the infected group compared with the control group, 2537 mRNAs and 233 lncRNAs were upregulated and 1438 mRNAs and 260 lncRNAs were downregulated (Fig. 1d, e). We also identified 477 circRNAs among samples using CIRI2 (Fig. 1f), but no marked change in expression was found due to the low number of overlapping circRNAs of samples without RNase R enrichment. The top 10 DE mRNAs and lncRNAs were validated by RT‒qPCR (Fig. 1g, h).
GO term and KEGG pathway analysis of DE mRNAs underscored the complexity of immune responses
Hierarchical clustering of DE mRNAs (Fig. 2a) and DE lncRNAs (Fig. 2b) indicated a clearcut discrimination between infected and control conditions, which suggested latent biological changes during infection. GO and KEGG analyses were performed to further investigate the crucial molecular events and processes. A series of classical liver injury pathways were significantly detected , including tumor necrosis factor (TNF)-related pathways, natural killer cell-related pathways, natural killer T (NKT) cell-related pathways, NF-kappa B-related pathways and interleukin-6 (IL6)-related pathways. Genes in those pathways were significantly enriched in DE gene sets (P < 0.001, Fisher’s exact test; detailed statistical information on 6 pathways is provided in Additional file 5: Table S3). Both innate and adaptive immune responses were enriched, suggesting comprehensive inflammatory activation due to parasite defense responses.
The significantly enriched biological process (BP) and molecular function (MF) GO terms involved immunological signatures, such as “innate immune response,” “positive regulation of T-cell proliferation,” “positive regulation of inflammatory response,” “chemokine receptor activity” and “cytokine activity” (Fig. 2c, d). Notably, most cellular component (CC) GO terms were related to membrane and extracellular interaction processes (Fig. 2e), implying that those genes were associated with excessive accumulation of extracellular matrix, as was reported in previous studies of liver damage in malaria .
Likewise, the significantly enriched KEGG pathways were related to immune processes (Fig. 2f). The uncovered pathways “cytokine−cytokine receptor interaction,” “cell adhesion molecules,” “NF−kappa B signaling pathway” and “TNF signaling pathway” indicated the links between cytokines and inflammatory responses, which have been reported in other Plasmodium–rodent systems and human malaria patients [36,37,38]. Intriguingly, the pathways “toxo–plasmosis” and “leishmaniasis”, which are involved in similar intracellular parasitic infections affecting internal organs, including the liver, were enriched. The autoimmune disease terms “rheumatoid arthritis” and “inflammatory bowel disease (IBD)” were unexpectedly enriched KEGG pathways, suggesting underlying cross-talk of immune responses between intracellular parasitic infections and autoimmunity due to the imbalance of pro-/anti-inflammatory responses [39, 40].
GO term and KEGG pathway analysis of DE lncRNA-targeted genes highlighted metabolic disturbances
To explore the potential regulatory roles of DE lncRNAs, we predicted the targeted mRNAs of 493 DE lncRNAs by correlation of expression levels. A total of 1,599 mRNAs of 349 DE lncRNAs were identified as lncRNA-targeted mRNAs (correlation > 0.99 and adjusted p value < 0.01). Remarkably, metabolism-associated pathways were highlighted in the functional analysis of targeted mRNAs.
In BP enrichment terms, “fatty acid biosynthetic process,” “fatty acid metabolic process,” “oxidation−reduction process,” “steroid biosynthetic process” and “oxidoreductase activity” were the pivotal metabolic pathways associated with infection (Fig. 3a, b). CC and MF enrichment terms indicated the active status of the intracellular membrane and enzyme catalysis. These results suggested that lncRNA-targeted mRNAs predominantly contributed to metabolic flexibility during infection, such as in cellular lipid metabolism and oxidation reduction reactions (Fig. 3c).
Consistent with the GO enrichment results, KEGG analysis also revealed similar metabolism-related pathways, including “regulation of lipolysis in adipocytes” and “sphingolipid signaling pathway” (Fig. 3d). Additionally, hypoxia/ischemia response-associated pathways were detected. The term “HIF–1 signaling pathway,” which has been reported to be essential in the pathogenesis of severe malaria, highlights the hypoxic response because of pRBC sequestration and vessel occlusion . The “cGMP−PKG signaling pathway” was found to be involved in protection against ischemia‒reperfusion injury , suggesting that similar multifocal damage occurred in the liver (Additional file 3: Fig. S2), which might be a consequence of clogging of microvessels.
The host-parasite gene coexpression network indicated tissue remodeling-like changes
During the infection process, the interaction between P. yoelii and the host results in large-scale changes in gene expression within both organisms. We subjected Plasmodium-containing liver samples to RNA-seq, with the aim to explore the simultaneous expression profiles of host and parasite, thus obtaining a novel perspective on host-parasite interactions. Through an artificial genome based on the reference genomes of both M. musculus and P. yoelii, we filtered 13,702 genes (386 parasite genes and 13,316 host genes) in dual RNA-seq. The PCA showed obvious differences between the infected and control groups. We then constructed a gene coexpression network, and five modules that were most highly related to infection were selected (Fig. 4a).
To minimize the background noise of the host-parasite cotranscriptome, we considered only P. yoelii genes in the top 1000 soft-connectivity of each module to be potential host-influencing genes. The five selected infection-related modules contained a total of 60 eligible genes (Fig. 4a). In the module “pink,” which contained the most parasite genes, the most highly coexpressed host genes included the endosome homeostasis-related gene SNX10 , the apoptosis/inflammation-related gene Casp4 , the liver fibrosis-related gene Jam2 , the ischemia/reperfusion regulating gene Tbk1  and the lymphocyte activation-related gene Cd28  (Fig. 4b), indicating elevated inflammation levels and activated immune responses that were consistent with the previous results of functional enrichment (Fig. 2c, f).
Further functional enrichment analyses of all highly coexpressed host genes revealed 67 overlapping terms that were detected in at least three modules (Fig. 4c; Additional file 6: Table S4). A predominant portion of these terms was associated with dynamic changes in the cellular transport system, including “cytoskeleton,” “cell cortex,” “endocytosis,” “Golgi apparatus,” “endosome,” “early endosome,” “intracellular membrane-bounded organelle,” “mitochondrial outer membrane,” “protein transport” and “extracellular exosome”, reflecting a complex intercellular and extracellular interaction network. Moreover, terms related to cell‒cell connections were also significantly enriched, such as “cadherin binding involved in “cell‒cell adhesion,” “cell‒cell adhesion,”, “focal adhesion” and “cell‒cell adherens junction”, which highlighted the altered tissue structure and histological environment during liver damage. The terms “angiogenesis”, “regulation of cell proliferation” and “apoptotic process” demonstrated that the tissue remodeling processes were similar to those of chronic liver diseases , which might be the adaptive repair of tissue injury.
Immunity-related genes were enriched in “positive regulation of inflammatory response,” “innate immune response,” “response to cytokine,” “positive regulation of NF-kappa B transcription factor activity,” “positive regulation of I-kappa B kinase/NF-kappa B signaling” and “TNF signaling” pathways. Specifically, the metabolic term “biosynthesis of antibiotics” was uncovered in the KEGG pathway, which might be related to microbial translocation occurring in the early stage of liver diseases associated with Plasmodium burden .
Conserved and nonconserved coexpressed genes of P. yoelii
The 60 hub genes of P. yoelii identified in the coexpression network were mainly composed of energy biogenesis-related genes, FAM genes, Ras family genes and transcription/translation-related genes. We next concentrated on the homologous genes of other Plasmodium species to enhance our understanding of these protein-coding genes and thus shed light on the pathogenesis and virulence of Plasmodium spp. The Plasmodium genomes of the human parasites P. falciparum , P. malariae , P. ovale , P. vivax  and P. knowlesi  as well as the rodent parasites P. berghei  and P. chabaudi  were included for homology identification. Taking an e-value < 1 × 10–10 as the threshold, only P. yoelii proteins with > 90% overlapping sequences with other Plasmodium proteins and overlapping sequences accounting for > 90% of the two Plasmodium protein sequences were kept for downstream analysis (Additional file 7: Table S5).
Of these 60 hub genes, three failed homology screening, which included genes encoding an uncharacterized protein, a Plasmodium exported protein of unknown function and an S-antigen. S-antigen is assumed to contribute to the extensive malaria antigenic diversity of the asexual blood stage of the parasite, which is considered useful in the stereotypic analysis of natural parasite populations . As expected, more conserved genes of P. yoelii were found in rodent parasites P. berghei and P. chabaudi than in human parasites. Ten P. yoelii genes were identified as conserved among seven Plasmodium spp. (Fig. 5), including four histone genes, the 40S ribosomal protein S27 gene, the translation initiation factor eIF-1A gene, the tubulin subunit beta gene, the centrin-1 gene and the ubiquitin-conjugating enzyme E2 gene. Tubulin and centrin are cytoskeleton proteins with key roles in cell division, such as merozoite biogenesis [56, 57]. Ubiquitin-conjugating enzyme E2 and the less conserved ubiquitin-conjugating enzyme E2 N are involved in ATP-dependent ubiquitination, which has been acknowledged to be involved in the growth and adaptation of apicomplexan parasites .
Jaundice, hepatocyte dysfunction and hepatic encephalopathy are common severe manifestations in patients with malaria [19, 59,60,61] and are diseases believed to be associated with complex complications. However, the direct roles of the liver in the blood stage of the parasite remain unclear. In our study, using a P. yoelii-rodent model system, we demonstrated the molecular characteristics of the infected mouse liver during the blood stage of P. yoelii infection, when activated immune responses, disrupted lipid metabolism, extensive hypoxia/ischemia stress and accompanying tissue remodeling were significant biological events.
Infection-specific perturbations in lipid metabolism have been demonstrated in malaria patients, and these could contribute to poor clinical outcomes [62,63,64]. The intertwined metabolic pathways of the host and parasite are preferentially studied integrally, but this approach leads to the loss of detailed information on specially involved organs. Our results show that the lipid metabolism processes of hepatocytes were also disturbed. Whether there exists a correlation between aberrant liver metabolism and disordered metabolomics remains to be clarified.
Coinfection with parasites and bacteria is a noteworthy clinical issue in malaria-endemic regions and is associated with concomitant sepsis and poor prognosis [65,66,67]. Approximately 6% of children with severe P. falciparum malaria were observed to suffer from invasive bacterial diseases, such as nontyphoidal Salmonella [65, 68, 69]. In our study, the bacterial immune response was detected in DE mRNA and host-parasite coexpression functional analysis, suggesting that a potential concomitant bacterial invasion occurred in Plasmodium infection. The liver has a unique vascular system within the gastrointestinal tract, as the majority of the liver's blood supply comes from the gut via the portal vein. The gastrointestinal tract microenvironment has been demonstrated to be the target of Plasmodium parasites . A malaria infection can cause disruption of the vascular and epithelial barriers in the gut, leading to dysfunction of gut microbiota homeostasis. Long-term liver damage was found to be associated with altered gut microbes in a malaria mouse model , while increased abundances of some probiotics were found to decrease the parasite burden . The resolution of this three-way crosstalk between host-plasmodium-microbiota may help to elucidate immune responses, anti-parasitic responses and regulatory roles of the gut microbiota, exploiting the gut microbiota as novel targets for antimalarial treatment.
The dual RNA-seq approach allows parallel analysis of the host and pathogen transcriptomes, providing profound insights into a comprehensive understanding of host‒pathogen interactions. In our study, hHost genes highly coexpressed with P. yoelii genes were mainly focused on processes of tissue remodeling, which was similar to the structural reconstruction that is observed in chronic liver diseases . We hypothesized that sinusoid disorder might worsen hepatocyte hypoxia and disrupt portal circulation, as occurs in other peripheral organs , which leads to further liver dysfunction and provides shelter for parasites. This immune escape might be an explanation for the higher incidence of complications in patients with malarial hepatitis, but more experiments are needed to prove this hypothesis.
Conserved genes reflect common responses, while nonconserved genes indicate unique changes in species. In our study, three out of 60 P. yoelii genes failed homology screening, including the S-antigen gene. Despite the unknown function and polymorphism of the S-antigen system among strains, this system is capable of eliciting antibodies with neutralizing properties , which may warrant revisiting in future vaccine studies of P. yoelii. Conserved genes related to cytoskeleton ubiquitination are coexpressed with host genes in the parasite-host network. Several microtubule inhibitors are reported to block the growth and development of malaria parasites at certain concentrations [75, 76]. The presented selective toxicity is partly due to different tubulin affinities, rendering microtubule inhibitors with different species affinities a promising anti-malaria approach. The type and chain length of ubiquitination further determines whether the protein is assigned to play further roles in multicellular processes, such as transcription. Ubiquitination of the E3 ligase NEDD4 has been reported to result in downregulation of autophagy-related genes in exoerythrocytic forms (EEFs) during the liver stage of the parasite in the Plasmodium berghei ANKA model, potentially promoting survival of and completion of hepatocyte infection by EEFs . Due to the significant structural conservation of several ubiquitination genes, drugs targeting the ubiquitin system could potentially be developed into antimalarial agents with promise for overcoming resistance [78,79,80,81].
The P. chabaudi–rodent model is another malaria mouse model that has been used in recent years to study liver pathology, with abnormal metabolism and activated immune responses reported using this model. Hepatic cytochrome 450 (CYP) activity has been reported to be regulated in the P. chabaudi-mouse model, which may further interfere with drug metabolism . We noted that in our study, in addition to identifying both Cyp2a5 upregulation and Cyp1a2 downregulation, we found that the gene expression levels of several genes of the CYP family 2, subgroup family c (Cyp2c23, Cyp2c29, Cyp2c40, Cyp2c50, Cyp2c55, Cyp2c67, Cyp2c78, Cyp2c69, Cyp2c70) were also significantly decreased under infection conditions. Similar to our enriched immune pathways, activated immune responses, such as pro-inflammatory IL-1α, anti-inflammatory factors IL-22, TNF-α and hepatic hemozoin levels, were demonstrated in previous studies to be involved in the development of liver pathology in the P. chabaudi-mouse model [14, 21, 22, 83]. Further transcriptomic studies performed on P. chabaudi-mouse models showed that innate immune responses and adaptive immune responses play distinct roles in the acute and recrudescence stages of the parasite [84, 85]. Considering the differences between P. chabaudi and P. yoelii , detailed liver changes need to be studied and compared to different clinical liver damage manifestations in human malaria.
Taken together, our study presents a systematic overview of biological changes in the P. yoelii-infected mouse liver during the blood stage of the parasite, in which activated immune responses, aberrant metabolism processes and tissue remodeling changes were the dominant molecular characteristics. However, the small amount of P. yoelii in the liver and the restricted sequencing depth limited our research to partial parasite genes in the coexpression network analysis; hence, multi-organ studies with deeper depth are required for further insight into host-parasite interactions. We expect that our findings will advance current understanding of malaria complications and facilitate the design of novel interventions to improve the clinical outcomes of severe malaria.
Activated immune responses, aberrant metabolic processes and tissue remodeling changes were the dominant molecular characteristics in the livers of P. yoelii-infected mice during the blood stage of the parasite. Future work will require the use of experimental approaches to demonstrate these changes as well as the roles of the liver in malaria.
Kyoto Encyclopedia of Genes and Genomes
Long noncoding RNA
Principal component analysis
Plasmodium-infected red blood cells
Rogerson SJ, Desai M, Mayor A, Sicuri E, Taylor SM, van Eijk AM. Burden, pathology, and costs of malaria in pregnancy: new developments for an old problem. Lancet Infect Dis. 2018;18:e107–18.
Bush MA, Anstey NM, Yeo TW, Florence SM, Granger DL, Mwaikambo ED, et al. Vascular dysfunction in malaria: understanding the role of the endothelial glycocalyx. Front Cell Dev Biol. 2021;9:751251.
Cowman AF, Healer J, Marapana D, Marsh K. Malaria: biology and disease. Cell. 2016;167:610–24.
Dondorp AM, Pongponratn E, White NJ. Reduced microcirculatory flow in severe falciparum malaria: pathophysiology and electron-microscopic pathology. Acta Trop. 2004;89:309–17.
Wilairatana P, Mala W, Kotepui M, Kotepui KU. Alteration of blood lactate levels in severe falciparum malaria: a systematic review and meta-analysis. Biology. 2021;10:1085.
Possemiers H, Vandermosten L, Van den Steen PE. Etiology of lactic acidosis in malaria. PLoS Pathog. 2021;17:e1009122.
White NJ. Anaemia and malaria. Malar J. 2018;17:371.
Namazzi R, Opoka R, Datta D, Bangirana P, Batte A, Berrens Z, et al. Acute kidney injury interacts with coma, acidosis, and impaired perfusion to significantly increase risk of death in children with severe malaria. Clin Infect Dis. 2022;75:1511–9.
Shaikh MS, Ali B, Janjua M, Akbar A, Haider SA, Moiz B, et al. Plasmodium in the bone marrow: case series from a hospital in Pakistan, 2007–2015. Malar J. 2021;20:254.
Robinson MW, Harmon C, O’Farrelly C. Liver immunology and its role in inflammation and homeostasis. Cell Mol Immunol. 2016;13:267–76.
Baer K, Roosevelt M, Clarkson AB, van Rooijen N, Schnieder T, Frevert U. Kupffer cells are obligatory for Plasmodium yoelii sporozoite infection of the liver. Cell Microbiol. 2007;9:397–412.
Liehl P, Zuzarte-Luís V, Chan J, Zillinger T, Baptista F, Carapau D, et al. Host-cell sensors for Plasmodium activate innate immunity against liver-stage infection. Nat Med. 2014;20:47–53.
Miller JL, Sack BK, Baldwin M, Vaughan AM, Kappe SHI. Interferon-mediated innate immune responses against malaria parasite liver stages. Cell Rep. 2014;7:436–47.
de Menezes MN, Salles ÉM, Vieira F, Amaral EP, Zuzarte-Luís V, Cassado A, et al. IL-1α promotes liver inflammation and necrosis during blood-stage Plasmodium chabaudi malaria. Sci Rep. 2019;9:7575.
Dey S, Bindu S, Goyal M, Pal C, Alam A, Iqbal MS, et al. Impact of intravascular hemolysis in malaria on liver dysfunction: involvement of hepatic free heme overload, NF-κB activation, and neutrophil infiltration. J Biol Chem. 2012;287:26630–46.
Kim J, Wang S, Lee C, Sung S, Shin Y, Song KS, et al. Blood-stage Plasmodium berghei ANKA infection promotes hepatic fibrosis by enhancing hedgehog signaling in mice. Cell Physiol Biochem. 2018;50:1414–28.
Medeiros MM, da Silva HB, Reis AS, Barboza R, Thompson J, Lima MRDI, et al. Liver accumulation of Plasmodium chabaudi-infected red blood cells and modulation of regulatory T cell and dendritic cell responses. PLoS ONE. 2013;8:e81409.
Prommano O, Chaisri U, Turner GDH, Wilairatana P, Ferguson DJP, Viriyavejakul P, et al. A quantitative ultrastructural study of the liver and the spleen in fatal falciparum malaria. Southeast Asian J Trop Med Public Health. 2005;36:1359–70.
Abro AH, Ustadi AM, Abro HA, Abdou AS, Younis NJ, Akaila SI. Jaundice with hepatic dysfunction in P. falciparum malaria. J Coll Physicians Surg Pak. 2009;19:363–6.
Kaeley N, Ahmad S, Shirazi N, Bhatia R, Bhat NK, Srivastava S, et al. Malarial hepatopathy: a 6-year retrospective observational study from Uttarakhand, North India. Trans R Soc Trop Med Hyg. 2017;111:220–5.
Mastelic B, do Rosario APF, Veldhoen M, Renauld JC, Jarra W, Sponaas A-M, et al. IL-22 protects against liver pathology and lethality of an experimental blood-stage malaria infection. Front Immunol. 2012;3:85.
Deroost K, Lays N, Pham T-T, Baci D, Van den Eynde K, Komuta M, et al. Hemozoin induces hepatic inflammation in mice and is differentially associated with liver pathology depending on the Plasmodium strain. PLoS ONE. 2014;9:e113519.
Rupani AB, Amarapurkar AD. Hepatic changes in fatal malaria: an emerging problem. Ann Trop Med Parasitol. 2009;103:119–27.
Brugat T, Cunningham D, Sodenkamp J, Coomes S, Wilson M, Spence PJ, et al. Sequestration and histopathology in Plasmodium chabaudi malaria are influenced by the immune response in an organ-specific manner. Cell Microbiol. 2014;16:687–700.
Pandey AV, Joshi R, Tekwani BL, Singh RL, Chauhan VS. Synthetic peptides corresponding to a repetitive sequence of malarial histidine rich protein bind haem and inhibit haemozoin formation in vitro. Mol Biochem Parasitol. 1997;90:281–7.
Guha M, Kumar S, Choubey V, Maity P, Bandyopadhyay U. Apoptosis in liver during malaria: role of oxidative stress and implication of mitochondrial pathway. FASEB J. 2006;20:1224–6.
Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019;37:907–15.
Pertea M, Kim D, Pertea GM, Leek JT, Salzberg SL. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat Protoc. 2016;11:1650–67.
Gao Y, Zhang J, Zhao F. Circular RNA identification based on multiple seed matching. Brief Bioinform. 2018;19:803–10.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.
Carlton JM, Angiuoli SV, Suh BB, Kooij TW, Pertea M, Silva JC, et al. Genome sequence and comparative analysis of the model rodent malaria parasite Plasmodium yoelii yoelii. Nature. 2002;419:512–9.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Tacke F, Luedde T, Trautwein C. Inflammatory pathways in liver homeostasis and liver injury. Clin Rev Allergy Immunol. 2009;36:4–12.
Viriyavejakul P, Khachonsaksumet V, Punsawad C. Liver changes in severe Plasmodium falciparum malaria: histopathology, apoptosis and nuclear factor kappa B expression. Malar J. 2014;13:106.
Deroost K, Pham T-T, Opdenakker G, Van den Steen PE. The immunological balance between host and parasite in malaria. FEMS Microbiol Rev. 2016;40:208–57.
Gowda DC, Wu X. Parasite recognition and signaling mechanisms in innate immune responses to malaria. Front Immunol. 2018;9:3006.
Jafari AA, Keikha M, Mirmoeeni S, Rahimi MT, Jafari R. Parasite-based interventions in systemic lupus erythematosus (SLE): a systematic review. Autoimmun Rev. 2021;20:102896.
Kot K, Łanocha-Arendarczyk N, Ptak M, Łanocha A, Kalisińska E, Kosik-Bogacka D. Pathomechanisms in the kidneys in selected protozoan parasitic infections. Int J Mol Sci. 2021;22:4209.
Boldt ABW, van Tong H, Grobusch MP, Kalmbach Y, Dzeing Ella A, Kombila M, et al. The blood transcriptome of childhood malaria. eBioMedicine. 2019;40:614–25.
Burley DS, Ferdinandy P, Baxter GF. Cyclic GMP and protein kinase-G in myocardial ischaemia-reperfusion: opportunities and obstacles for survival signaling. Br J Pharmacol. 2007;152:855–69.
You Y, Li W-Z, Zhang S, Hu B, Li Y-X, Li H-D, et al. SNX10 mediates alcohol-induced liver injury and steatosis by regulating the activation of chaperone-mediated autophagy. J Hepatol. 2018;69:129–41.
Khanova E, Wu R, Wang W, Yan R, Chen Y, French SW, et al. Pyroptosis by Caspase11/4-Gasdermin-D pathway in alcoholic hepatitis in mice and patients. Hepatology. 2018;67:1737–53.
Hintermann E, Bayer M, Conti CB, Fuchs S, Fausther M, Leung PS, et al. Junctional adhesion molecules JAM-B and JAM-C promote autoimmune-mediated liver fibrosis in mice. J Autoimmun. 2018;91:83–96.
Lv P, Li C, Wang M, Ren J, Zhang Y, Fu G. TANK-binding kinase 1 alleviates myocardial ischemia/reperfusion injury through regulating apoptotic pathway. Biochem Biophys Res Commun. 2020;528:574–9.
Spengler U, Leifeld L, Braunschweiger I, Dumoulin FL, Lechmann M, Sauerbruch T. Anomalous expression of costimulatory molecules B7–1, B7–2 and CD28 in primary biliary cirrhosis. J Hepatol. 1997;26:31–6.
Ucuzian AA, Gassman AA, East AT, Greisler HP. Molecular mediators of angiogenesis. J Burn Care Res. 2010;31:158–75.
Walker GT, Yang G, Tsai JY, Rodriguez JL, English BC, Faber F, et al. Malaria parasite infection compromises colonization resistance to an enteric pathogen by reducing gastric acidity. Sci Adv. 2021;7:6232.
Böhme U, Otto TD, Sanders M, Newbold CI, Berriman M. Progression of the canonical reference malaria parasite genome from 2002–2019. Wellcome Open Res. 2019;4:58.
Rutledge GG, Böhme U, Sanders M, Reid AJ, Cotton JA, Maiga-Ascofare O, et al. Plasmodium malariae and P. ovale genomes provide insights into malaria parasite evolution. Nature. 2017;542:101–4.
Carlton JM, Adams JH, Silva JC, Bidwell SL, Lorenzi H, Caler E, et al. Comparative genomics of the neglected human malaria parasite Plasmodium vivax. Nature. 2008;455:757–63.
Pain A, Böhme U, Berry AE, Mungall K, Finn RD, Jackson AP, et al. The genome of the simian and human malaria parasite Plasmodium knowlesi. Nature. 2008;455:799–803.
Otto TD, Böhme U, Jackson AP, Hunt M, Franke-Fayard B, Hoeijmakers WAM, et al. A comprehensive evaluation of rodent malaria parasite genomes and gene expression. BMC Biol. 2014;12:86.
Anders RF, McColl DJ, Coppel RL. Molecular variation in Plasmodium falciparum: polymorphic antigens of asexual erythrocytic stages. Acta Trop. 1993;53:239–53.
Fennell BJ, Al-shatr ZA, Bell A. Isotype expression, post-translational modification and stage-dependent production of tubulins in erythrocytic Plasmodium falciparum. Int J Parasitol. 2008;38:527–39.
Mahajan B, Selvapandiyan A, Gerald NJ, Majam V, Zheng H, Wickramarachchi T, et al. Centrins, cell cycle regulation proteins in human malaria parasite Plasmodium falciparum. J Biol Chem. 2008;283:31871–83.
Ponts N, Yang J, Chung D-WD, Prudhomme J, Girke T, Horrocks P, et al. Deciphering the ubiquitin-mediated pathway in apicomplexan parasites: a potential strategy to interfere with parasite virulence. PLoS ONE. 2008;3:e2386.
Kochar DK, Agarwal P, Kochar SK, Jain R, Rawat N, Pokharna RK, et al. Hepatocyte dysfunction and hepatic encephalopathy in Plasmodium falciparum malaria. QJM. 2003;96:505–12.
Asma U-e, Taufiq F, Khan W. Prevalence and clinical manifestations of malaria in Aligarh India. Korean J Parasitol. 2014;52:621–9.
Chaparro-Narváez PE, Lopez-Perez M, Rengifo LM, Padilla J, Herrera S, Arévalo-Herrera M. Clinical and epidemiological aspects of complicated malaria in Colombia, 2007–2013. Malar J. 2016;15:269.
Olszewski KL, Morrisey JM, Wilinski D, Burns JM, Vaidya AB, Rabinowitz JD, et al. Host-parasite interactions revealed by Plasmodium falciparum metabolomics. Cell Host Microbe. 2009;5:191–9.
Gardinassi LG, Arévalo-Herrera M, Herrera S, Cordy RJ, Tran V, Smith MR, et al. Integrative metabolomics and transcriptomics signatures of clinical tolerance to Plasmodium vivax reveal activation of innate cell immunity and T cell signaling. Redox Biol. 2018;17:158–70.
Abdrabou W, Dieng MM, Diawara A, Sermé SS, Almojil D, Sombié S, et al. Metabolome modulation of the host adaptive immunity in human malaria. Nat Metab. 2021;3:1001–16.
Church J, Maitland K. Invasive bacterial co-infection in African children with Plasmodium falciparum malaria: a systematic review. BMC Med. 2014;12:31.
Nyein PP, Aung NM, Kyi TT, Htet ZW, Anstey NM, Kyi MM, et al. High frequency of clinically significant bacteremia in adults hospitalized with falciparum malaria. Open Forum Infect Dis. 2016;3:ofw028.
Phu NH, Day NPJ, Tuan PQ, Mai NTH, Chau TTH, Van Chuong L, et al. Concomitant bacteremia in adults with severe falciparum malaria. Clin Infect Dis. 2020;71:e465–70.
Bassat Q, Guinovart C, Sigaúque B, Mandomando I, Aide P, Sacarlal J, et al. Severe malaria and concomitant bacteraemia in children admitted to a rural Mozambican hospital. Trop Med Int Health. 2009;14:1011–9.
Were T, Davenport GC, Hittner JB, Ouma C, Vulule JM, Ong’echa JM, et al. Bacteremia in Kenyan children presenting with malaria. J Clin Microbiol. 2011;49:671–6.
Coban C, Lee MSJ, Ishii KJ. Tissue-specific immunopathology during malaria infection. Nat Rev Immunol. 2018;18:266–78.
Denny JE, Powers JB, Castro HF, Zhang J, Joshi-Barve S, Campagna SR, et al. Differential sensitivity to Plasmodium yoelii infection in C57BL/6 mice impacts gut-liver axis homeostasis. Sci Rep. 2019;9:3472.
Villarino NF, LeCleir GR, Denny JE, Dearth SP, Harding CL, Sloan SS, et al. Composition of the gut microbiota modulates the severity of malaria. Proc Natl Acad Sci USA. 2016;113:2235–40.
Paternostro C, David E, Novo E, Parola M. Hypoxia, angiogenesis and liver fibrogenesis in the progression of chronic liver diseases. World J Gastroenterol. 2010;16:281–8.
Fu Y, Ding Y, Zhou T-l, Ou Q-Y, Xu W-Y. Comparative histopathology of mice infected with the 17XL and 17XNL strains of Plasmodium yoelii. J Parasitol. 2012;98:310–5.
Fennell BJ, Carolan S, Pettit GR, Bell A. Effects of the antimitotic natural product dolastatin 10, and related peptides, on the human malarial parasite Plasmodium falciparum. J Antimicrob Chemother. 2003;51:833–41.
Naughton JA, Hughes R, Bray P, Bell A. Accumulation of the antimalarial microtubule inhibitors trifluralin and vinblastine by Plasmodium falciparum. Biochem Pharmacol. 2008;75:1580–7.
Zheng H, Lu X, Li K, Zhu F, Zhao C, Liu T, et al. ATG ubiquitination is required for circumsporozoite protein to subvert host innate immunity against rodent malaria liver stage. Front Immunol. 2022;13:815936.
Ponts N, Saraf A, Chung D-WD, Harris A, Prudhomme J, Washburn MP, et al. Unraveling the ubiquitome of the human malaria parasite. J Biol Chem. 2011;286:40320–30.
Chung D-WD, Ponts N, Prudhomme J, Rodrigues EM, Le Roch KG. Characterization of the ubiquitylating components of the human malaria parasite’s protein degradation pathway. PLoS ONE. 2012;7:e43477.
Green JL, Wu Y, Encheva V, Lasonder E, Prommaban A, Kunzelmann S, et al. Ubiquitin activation is essential for schizont maturation in Plasmodium falciparum blood-stage development. PLoS Pathog. 2020;16:e1008640.
Kumar P, Kumar P, Mandal D, Velayutham R. The emerging role of Deubiquitinases (DUBs) in parasites: a foresight review. Front Cell Infect Microbiol. 2022;12:985178.
De-Oliveira ACAX, Carvalho RS, Paixão FHM, Tavares HS, Gueiros LS, Siqueira CM, et al. Up- and down-modulation of liver cytochrome P450 activities and associated events in two murine malaria models. Malar J. 2010;9:81.
Seixas E, Oliveira P, Moura Nunes JF, Coutinho A. An experimental model for fatal malaria due to TNF-α-dependent hepatic damage. Parasitology. 2008;135:683–90.
Chen S, Gao Y, Fan Y, Guo S, Zhou J, Liu T, et al. The dynamic change of immune responses between acute and recurrence stages of rodent malaria infection. Front Microbiol. 2022;13:844975.
Kumarasingha R, Ioannidis LJ, Abeysekera W, Studniberg S, Wijesurendra D, Mazhari R, et al. Transcriptional memory-like imprints and enhanced functional activity in γδ T cells following resolution of malaria infection. Front Immunol. 2020;11:582358.
Pattaradilokrat S, Wu J, Xu F, Su X-Z. The origins, isolation, and biological characterization of rodent malaria parasites. Parasitol Int. 2022;91:102636.
This work was supported by the Fundamental Research Funds for the Central Universities (BMU2019XK007).
Ethics approval and consent to participate
Procedures involving vertebrate animals were reviewed and approved by the Animal Care and Use Committee of Peking University Health Science Center (permit number PUIRB-LA2022677).
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Primers of the top 10 mRNAs and top 10 lncRNAs used for RT‒qPCR validation.
Blood smears of five infected mice on day 5 post-infection.
Parasitemia and histopathology of mouse liver. A HE staining of control mouse liver. B, C, D HE staining of P. yoelii-infected mice.
Detailed information on RNA-seq.
GO terms related to liver injury and Fisher's exact test of term enrichment.
Overlapping GO and KEGG terms of coexpression modules.
Homologous genes among Plasmodium spp.
About this article
Cite this article
Hu, X., Zhao, J., Zhao, J. et al. Genome-wide liver transcriptomic profiling of a malaria mouse model reveals disturbed immune and metabolic responses. Parasites Vectors 16, 40 (2023). https://doi.org/10.1186/s13071-023-05672-w