A comprehensive rhythmicity analysis of host proteins and immune factors involved in malaria pathogenesis to decipher the importance of host circadian clock in malaria

Background Circadian rhythms broadly impact human health by regulating our daily physiological and metabolic processes. The circadian clocks substantially regulate our immune responses and susceptibility to infections. Malaria parasites have intrinsic molecular oscillations and coordinate their infection cycle with host rhythms. Considering the cyclical nature of malaria, a clear understanding of the circadian regulations in malaria pathogenesis and host responses is of immense importance. Methods We have thoroughly investigated the transcript level rhythmic patterns in blood proteins altered in falciparum and vivax malaria and malaria-related immune factors in mice, baboons, and humans by analyzing datasets from published literature and comprehensive databases. Using the Metascape and DAVID platforms, we analyzed Gene Ontology terms and physiological pathways associated with the rhythmic malaria-associated host immune factors. Results We observed that almost 50% of the malaria-associated host immune factors are rhythmic in mice and humans. Overlapping rhythmic genes identified in mice, baboons, and humans, exhibited enrichment (Q < 0.05, fold-enrichment > 5) of multiple physiological pathways essential for host immune and defense response, including cytokine production, leukocyte activation, cellular defense, and response, regulation of kinase activity, B-cell receptor signaling pathway, and cellular response to cytokine stimulus. Conclusions Our analysis indicates a robust circadian regulation on multiple interconnected host response pathways and immunological networks in malaria, evident from numerous rhythmic genes involved in those pathways. Host immune rhythms play a vital role in the temporal regulation of host-parasite interactions and defense machinery in malaria.


Introduction
Circadian rhythm is a biological mechanism that governs the cyclic variations in a wide diversity of physiological, metabolic, immunological, and behavioral processes over approximately 24 hours (1,2). The circadian clock tightly regulates our immune system, including the functioning of natural killer cells, neutrophils, monocytes, eosinophils, and macrophages, most likely to provide a steady defense in accordance with diurnal peaks of pathogenic infections (3,4). Cellular immunity against pathogenic microbes temporally controlled by the core clock circuit leads to fluctuations in our susceptibility to infections at different times of the day (5). In this vein, multiple studies have demonstrated that the time of day of infection has a vital influence on the progression and severity of infectious diseases (6)(7)(8).
The cyclical nature of malaria is well documented, as the infection cycles for Plasmodium species that infect humans last a multiple of approximately 24 hours (9). Importantly, recent studies have demonstrated that malaria parasites have intrinsic molecular oscillations and effectively synchronize their rhythms with the host (10,11). Malaria parasites control the intra-erythrocytic developmental cycle (IDC) progression through coordination with host timekeeping machinery (12). Synchronous maturation of Plasmodium in vivo is coordinated and regulated by the rhythmic changes in melatonin concentration in our body (13). If malaria parasites' cell cycle timing is perturbed relative to the host's circadian rhythm, it causes reduced replication and transmission of the pathogen (14, 15).
Among the five species of Plasmodium that can infect humans and cause malaria, Plasmodium falciparum and Plasmodium vivax are the primary pathogens of this protozoan infection. Severe P. falciparum and P. vivax infections alter many vital physiological pathways, including cytokine signaling, complement systems, inflammation, platelet hemostasis, and coagulation cascades in humans (16)(17)(18)(19). Various pathogens, including the malaria parasites, disrupt host circadian rhythms as a virulence trait to favor their replication and dissemination (20). Robust circadian timekeeping machinery in the host may reduce vulnerability towards any acute infections or severity of their clinical manifestations. This study systematically investigated the rhythmic patterns in serum, plasma, and tissue proteins altered in malaria and malaria-related immune factors by analyzing data from published literature and comprehensive databases ( Figure 1). We also mapped the diurnal rhythmicity in vital physiological pathways and host immunity cascades involved in malaria. Our findings indicate that several components within the host physiological pathways and immunological networks activated in malaria exhibit robust 24-hour rhythms. A detailed understanding of host rhythms may open new opportunities for developing clock-based efficient prophylactic and therapeutic strategies for malaria.

Documentation of host proteins altered in falciparum and vivax malaria
Information on differentially abundant serum or plasma proteins in two primary clinical forms of malaria, P. falciparum and P. vivax infections was retrieved from relevant published literature (16,17,19,(21)(22)(23)(24)(25)(26)(27)(28)(29)(30) (Supplementary Tables 1, 2). We also included a study describing the altered proteome of the frontal lobe in cerebral falciparum malaria patients (18). Human proteins altered in non-severe and severe falciparum and vivax malaria patients were recorded for subsequent analyses (Figure 1). UniProt IDs of the altered proteins from the published literature were converted into official gene symbols using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) gene conversion tool (https://david.ncifcrf.gov/home.jsp) (31).

Cataloging of immunity-related genes in malaria
Genes involved in immune responses in malaria were searched using the keywords "Malaria immunity genes," "Host immune response in malaria," and "Host defense mechanism in malaria" in PubMed and Google Scholar. Moreover, potential malaria-related immunity genes were searched in the Molecular Signatures Database (MSigDB) (https://www.gsea-msigdb.org/gsea/msigdb/index.jsp) using the keyword "Plasmodium." MSigDB is an extensively used, wide-ranging database of genes for the execution of gene set enrichment analysis (32). We removed the duplicates after combining the two data sets (entities from the published literature and MSigDB). A list of 392 malaria-related immunity genes was generated for further validation (Supplementary Table 3). The immunogenic identity of the identified 392 genes was verified further in the InnateDB database (https://www.innatedb.com) ( Figure 1). InnateDB is a comprehensive centralized resource for immune function-related genes (33). Venny 2.1 was used to analyze overlaps among different datasets (https://bioinfogp.cnb.csic.es/tools/ venny/).

Rhythmicity analysis of the malariaassociated host factors
Circadian expression profile database (CircaDB) (http:// circadb.hogeneschlab.org/) was used as a comprehensive resource to investigate the rhythmicity of the malaria-associated immunity genes and host factors in mice and humans. CircaDB provides information on the transcript level oscillations in various tissues of mice and humans (34, 35). Genes with robust oscillations in humans and mice were defined using period length and FDR thresholds (Period 24 ± 3 hours; JTK Q < 0.1). Rhythmicity of the malaria-associated immunity genes and host factors in a primate (baboon, Papio anubis)" was investigated in the diurnal transcriptomics datasets reported by Mure et al. (36). The study analyzed the diurnal transcriptome profiles in 64 tissues collected every 2 hours over a 24-hour cycle, where a p-value < 0.05 (MetaCycle algorithm) statistical threshold was used to define rhythmicity. Subsequently, the rhythmic malaria-associated immunity genes and host factors were compared in these three mammals (humans, mice, and baboons), and their overlaps were Schematic showing the analysis pipeline for investigating rhythmicity of host proteins and immune factors involved in malaria pathogenesis. At the onset, we collected data on altered serum, plasma, and tissue proteins in falciparum and vivax malaria and malaria-related immunity genes from published literature and databases. Consequently, we investigated 24-hour rhythmicity in mice, primates (baboons), and humans using the CircaDB database and published literature. Further, we performed Gene Ontology, pathway, and interaction analysis on the rhythmic host proteins and immune factors involved in malaria pathogenesis to interpret the importance of the host circadian clock in malaria.
visualized in a Circos plot (37). We also used Venny 2.1 to analyze overlaps among three datasets (humans, mice, and baboons) and plot the Venn diagrams.

Gene Ontology and pathway enrichment analysis
The pathway and Gene Ontology (GO)-Biological Process (BP) enrichment analysis of the rhythmic genes was carried out using the Metascape platform (https://metascape.org/gp/ index.html), which combines multiple databases (GO Biological Processes, Kyoto encyclopedia of genes and genomes (KEGG) pathway, Reactome gene sets, Comprehensive resource of mammalian protein complexes (CORUM), WikiPathways, and Protein analysis through evolutionary relationships (PANTHER) pathway (38). The top network clusters were generated by choosing terms based on p-value, enrichment factor, and minimum count thresholds. Grouped terms with a similarity > 0.3 were connected by edges and visualized using Cytoscape (https://cytoscape.org/) (39). Further, using the DAVID database, we examined the KEGG pathways of the overlapping genes (40). We visualized the physiological pathways associated with the rhythmic genes through a bubble plot representation using the SR Plot platform (https:// www.bioinformatics.com.cn/en).
Enrichment of transcription regulator using the Transcriptional regulatory relationships unraveled by sentence-based text mining (TRRUST) database (https://www.grnpedia.org/trrust/), and based on kappa-statistical similarities, significant terms were hierarchically clustered into a tree. TRRUST is a comprehensive reference database of human transcriptional regulatory networks (41). A threshold of 0.3 kappa score was applied to divide the tree into term clusters, and the lowest p-value was selected to visualize the top 20 terms into a heatmap/dendrogram.

Protein-protein interaction network analysis
Protein-protein interaction (PPI) enrichment analysis was performed utilizing the Metascape platform using the following databases (38): Search tool for the retrieval of interacting genes/ proteins (STRING), Biological general repository for interaction datasets (BioGrid), OmniPath, and InWeb-IM. Proteins that form physical interactions following the default analysis parameters with at least one other member of the list were used in the analysis (38). The Molecular complex detection (MCODE) algorithm was applied using best-scoring by p-value threshold to identify densely connected network components (42). The final cluster was visualized by Cytoscape using Y-files radial layout (39).

Rhythmic patterns in serum, plasma, and tissue proteins altered in falciparum malaria
We analyzed possible rhythmic patterns of 603 differentially abundant proteins in falciparum malaria (FM) at the transcription level using the CircaDB platform (Supplementary Table 1). Among these 603 malaria-associated host factors, we identified 198 (32.8%) and 265 (43.9%) rhythmic candidates (Period 24 ± 3 h, JTK Q < 0.1) in humans and mice, respectively. Using a diurnal transcriptome dataset (36), we also investigated the expression patterns of these malaria-associated host factors in the baboon. Non-human primates models are vital to malaria research for understanding its pathophysiology. We identified 322 (53.3%) rhythmic malaria-associated host factors in the primate model ( Figure 2A). Furthermore, we conducted a comparative analysis of rhythmic candidates identified in severe falciparum malaria (SFM) and non-severe falciparum malaria (NSFM); we identified 329 rhythmic candidates in SFM and 191 rhythmic candidates in NSFM across mice, humans, and primates (Supplementary Table 1

Rhythmicity in multiple vital physiological pathways associated with host factors altered in falciparum malaria
We investigated the physiological pathways associated with rhythmic malaria-associated host factors in mice, baboons, and humans, encompassing SFM and NSFM. Employing Metascape analysis, we examined the rhythmic candidates specific to SFM and NSFM to understand the dysregulated pathways with malaria severity. Analysis of rhythmic host proteins identified in SFM revealed enrichment of multiple physiological pathways vital for host-pathogen interaction and immune response, including complement and coagulation cascades (Q = 6.  Table 4). Notably, we observed complement and coagulation cascades, neutrophil degranulation, acute inflammatory response, response to inorganic substance, extracellular matrix organization, and integrin cell surface interactions as overlapping pathways associated with the rhythmic host proteins identified in both SFM and NSFM ( Figures 2F, G).

Analysis of rhythmic patterns in serum and plasma proteins altered in vivax malaria
We also investigated the rhythmic patterns of 345 serum or plasma proteins that exhibited differential abundance in vivax malaria (VM) (Supplementary Table 2). We identified 128 (37%) and 108 (31.3%) rhythmic candidates (Period 24 ± 3 h, JTK Q < 0.1) among these malaria-associated host factors in mice and humans, respectively. Moreover, our analysis revealed 179 (51.8%) rhythmic malaria-associated host factors in primates ( Figure 3A). Subsequently, we compared the rhythmic human proteins identified in FM and VM, which indicated 73 overlapping rhythmic candidates ( Figure 3B).
Alike FM, we precisely investigated the rhythmic host proteins in different severity levels of VM. We identified 126 rhythmic candidates in severe vivax malaria (SVM), while non-severe vivax malaria (NSVM) was associated with 164 rhythmic candidates ( Figures 3C, D). Cross-comparing the rhythmicity data across three species, i.e., mice, humans, and baboons, we identified 20 and 28 overlapping rhythmic candidates in SVM and NSVM, respectively. We have shown the rhythmic expression profiles of a few representative transcripts in SVM and NSVM, robustly oscillating in mice, humans, and baboons ( Figures 3E, F).

Physiological pathways associated with the rhythmic host factors altered in vivax malaria
We intended to map the various physiological pathways associated with the rhythmic host factors identified in SVM and NSVM (Supplementary Table 4 Figure 3H). Four pathways, namely complement and coagulation cascades, neutrophil degranulation, regulation of vesicle-mediated transport, and the acute-phase response, were common between SVM and NSVM.
Inclusive analysis of all the rhythmic host proteins altered in FM and VM indicated enrichment of complement and coagulation cascades, neutrophil degranulation, and vascular endothelial growth factor A (VEGFA) -VEGF receptor 2 (VEGFR2) signaling as the most significant entities ( Table 1). The circadian clock tightly regulates these pathways under healthy conditions, and blood levels of several components of these pathways are altered in P. falciparum and P. vivax infections.

Rhythmic expression of malaria-related immune factors in mammals
Besides the differentially abundant serum or plasma proteins, we investigated the possible rhythmicity in the expression of 307 malaria-associated host immune factors using the CircaDB (Supplementary Table 3). We identified 153 (49.84%) rhythmic immune factors (Period 24 ± 3 h, JTK Q < 0.1) in mice and 131 (42.67%) in humans. Further, analyzing the diurnal transcriptome dataset (36), we identified 218 (71.01%) rhythmic malariaassociated host immune factors in baboons. We observed a significant overlap of rhythmic host immunity genes with known biological functions across all three species represented as a Circos plot ( Figure 4A). Precisely, 71 overlapping rhythmic candidates were identified by comparing the malaria-associated rhythmic host immunity genes of mice, primates, and humans, accounting for 23% of total rhythmic candidates identified in these three species (Figure 1). We show here rhythmic expression patterns for several malaria-associated immune factors (top 15 of the 71 overlapping genes) that show almost 24-hour oscillation at the transcript level (JTK Q < 0.1) ( Figure 4B).

Diurnal rhythmicity in host immunity pathways in malaria
Based on the identified rhythmic malaria-associated host immune factors, we mapped the circadian regulation of diverse immunity pathways that play a crucial role in host response and defense mechanisms in malaria (Supplementary Table 5 Figure 4D). These TFs may regulate the temporal expression of host genes related to malaria immunity, suggesting the involvement of a complex molecular network controlled by the timekeeping machinery (Supplementary Table 6).   PPI analysis using Metascape and MCODE involving the rhythmic host factors related to malaria-associated immunity in mice, primates, and humans identified 11 PPI modules ( Figure 4E, Supplementary Table 7). Further, for the top three PPI modules based on p-values, we applied GO enrichment to the overall MCODE network to get biological relevance from the network component. The top three terms were found to be cell activation (GO:0001775, p-value = 1.15E-17), leukocyte activation (GO:0045321, p-value = 8.53E-17), and positive regulation of immune response (GO:0050778, p-value = 2.83E-16) (Supplementary Table 7). Taken together, our analysis indicates that the circadian clock tightly controls multiple interconnected host responses and immunological networks in malaria, which is evident from the presence of several rhythmic components within those pathways ( Figure 5).

Discussion
An inadequate understanding of the molecular mechanisms driving complex cyclical pathogenesis and multidimensional fatal manifestations in severe malaria impedes the clinical management of this protozoan infection. Investigation of clock-infection biology can substantially impact various aspects of malaria research, including understanding disease pathogenesis, dynamic hostpathogen interactions, and the time-dependent efficacy of therapeutic interventions. Malaria parasites show developmental rhythms during replication in the mammalian hosts and lead to rhythmic fevers due to the synchronous bursting of the parasiteinfected host red blood cells (RBCs) (43). Notably, self-sustained non-canonical redox and metabolic oscillations exist in human RBCs (44,45). Since our immune system shows rhythmic  Rhythmicity of malaria-associated immunity genes in mammalian systems. (A) The Circos plot shows the overlap of malaria-associated rhythmic immunity genes in mice, primates, and humans. The outer arc represents rhythmic genes specific to mice (green), primates (red), and humans (blue). Dark orange indicates overlapping malaria-associated rhythmic genes in mice, primates, and humans on the inner arc. Light orange indicates malaria-associated rhythmic genes unique to each organism. Blue linker lines represent the same enriched ontology terms, while the purple links represent input gene list overlap. (B) Transcript level rhythmic profiles of selected malaria-associated immunity genes in mammalian systems. CircaDB (http://circadb.hogeneschlab.org/) was used for rhythmicity analysis. 24-hour transcript level rhythmic profiles of the immunity genes are shown in different organelles of mice (JTK Q < 0.1, period 24 ± 3 hours). (C) Statistically enriched (Top 15, p-value < 0.01, Enrichment factor > 1.5, Count > 4) ontology terms (GO/KEGG terms and canonical pathways) associated with the overlapping malaria-associated rhythmic immunity genes identified in mice, primates (baboons), and humans. (D) Enrichment analysis of all statistically significant transcription factors (TFs)-target interaction networks associated with malaria-associated rhythmic immunity genes. The heatmap shows the strength of evidence for each TF-target relationship. The intensity of the color represents the strength of the evidence, with brighter colors indicating more robust evidence. (E) The protein-protein interaction (PPI) network of the malaria-associated rhythmic immunity genes created by Metascape shows 11 interacting modules (see Supplementary Table 7 for further details).
variations, investigating circadian control in the host response pathways in malaria is vital for understanding the complex pathogenesis of this parasitic infection. Here we report transcriptlevel rhythmic patterns in malaria-related immune factors and host response pathways in mice, baboons, and humans by analyzing datasets from published literature and wide-ranging databases. Several physiological pathways controlled by the circadian timekeeping machinery, including complement and coagulation cascades, interleukin signaling, acute inflammatory response, neutrophil degranulation, and VEGFA-VEGFR2 signaling, are dysregulated in malaria. Earlier, we have shown that severe manifestations in P. falciparum and P. vivax infections are caused by the modulation of host defense machinery and vital physiological pathways, including cytokine signaling, lipid metabolism, complement pathways, and coagulation cascades, which can be assessed by measuring the alterations in host serum or plasma proteome in malaria patients (16,17,19,46). Other contemporary studies on blood proteome in severe malaria patients also indicated similar alterations in acute phase response, inflammation, and metabolic pathways (21, 24). This study shows that several components of these physiological pathways exhibit robust 24hour rhythms in both nocturnal and diurnal mammals.
It is critical to decipher whether a circadian disruption in the host contributes to dysregulation in these diverse physiological pathways in malaria. For instance, circadian rhythm firmly controls angiogenesis, potentially by regulating VEGF expression and secretion (47). Notably, the core clock transcription factor brain and muscle ARNT-like-1 (Bmal1) regulates the rhythmic expression of VEGF through direct binding to its promoter region (47,48). Elevated levels of VEGF and VEGFR-2 signaling involved in angiogenesis are frequently observed in individuals suffering from malaria and can serve as a potential biomarker of malaria complications (49,50). Such alteration in VEGF in malaria is particularly relevant in cerebral manifestations, where the binding of parasitized erythrocytes to the cerebral endothelium and angiogenic dysregulation contributes to the pathogenesis of malaria (49,51). Thus, circadian disruptions in the host and their subsequent effects on the multidimensional physiological network involved in malaria pathogenesis and complications entail a comprehensive investigation. Rhythmicity information for the elements of host response pathways will help understand possible circadian disruptions in their rhythmic patterns under the febrile or convalescent stages or different severity levels of malaria.
Impairments in liver, brain, and kidney functions in severe P. falciparum infection often cause morbidity and mortality in complicated malaria (52)(53)(54). P. vivax infection also leads to severe malaria involving complicated clinical symptoms, including liver dysfunction, cerebral manifestations, and acute kidney injury Circadian regulation of host response in malaria. The malaria infection-immunity-clock represents the 24-hour oscillations in host proteins and immune factors that play an essential role in host-pathogen interaction. Here we show circadian regulations of the components of four major host response pathways in malaria (FDR < 0.05, fold enrichment > 2.5, Count > 5). The clock symbols indicate the transcript level rhythmicity (JTK Q < 0.1, period 24 ± 3 hours) of that entity involved in that particular physiological pathway [Created with BioRender.com]. (55,56). Our findings indicate that several malaria-induced altered proteins exhibit 24-hour rhythmicity in the liver, brain, and kidney with high statistical significance (JTK Q < 0.05). Numerous studies have highlighted the crucial role of the rhythmic candidates identified in our analysis in malaria progression and severity. For instance, ADAMTS13 deficiency has been observed in severe malaria, particularly in cerebral manifestations, which may contribute to the pathophysiology of complicated malaria by impairing the regulation of von Willebrand factor (VWF) and promoting endothelial dysfunction (57). KLKB1 involves blood coagulation, inflammation, fibrinolysis, and complement activation (58). The cell adhesion molecule VCAM1 attaches to infected RBCs and is involved in their sequestration (59,60). SERPINF2 (a2AP) shows altered levels in falciparum malaria patients, contributing to the development of disseminated intravascular coagulation and increasing the risk of thrombosis (19). Additionally, FKBP5, another highly rhythmic candidate in mammals, influences the host immune system and cerebral manifestations of malaria in mice (61, 62).
Host immune factors, cells, and pathways such as proinflammatory cascade, cytokine levels, IFN-g, and natural killer T (NKT) cells determines the disease state, susceptibility, and fatality in malaria (63,64). Circadian variations in the expression levels of the vital host factors involved in malaria pathophysiology and immunity shown in our study indicate the fluctuations in our defense machinery, which hints at a possible time-dependent vulnerability to the severe manifestations of malaria. Investigating whether disruptions in host circadian rhythms contribute to the severe clinical manifestations frequently observed in malaria patients is an open question. Recently, the "One Health" concept, which emphasizes a trans-disciplinary collaborative approach to strengthen the interrelation between people, animals, and their environment, is gaining popularity in malaria research and elimination programs (65). The circadian aspect of malaria infection biology in the hosts, vector, and pathogen should be emphasized in such a multi-sectoral health strategy.
If validated empirically in preclinical animal models of malaria or malaria patients, our findings can provide novel mechanistic insights into the dynamic host-pathogen interactions in malaria and strengthen the concept of chronotherapeutics for this protozoan infection. In this vein, a systematic investigation of circadian disruptions in non-severe and severe falciparum and vivax malaria patients will be beneficial to decipher if the magnitude of circadian alterations governs the severity of malaria, which could be an imperative continuation of the present study. The rhythmic malaria-related host proteins and immune factors identified in our study could be targeted to understand the overall circadian disruptions in malaria caused by different Plasmodium species. Moreover, recent studies have demonstrated the existence of potential markers for peripheral circadian rhythm in biological fluids, such as blood or saliva. Kramer and colleagues have developed a blood-based assay through multiplex gene expression profiling to determine the internal circadian time in humans (66). In another study, b-arrestin 1 (ARRB1) is identified as a molecular marker of the peripheral circadian rhythm in saliva (67).
Quantitive analysis of the possible alterations in rhythmicity of proteins and metabolites in blood or saliva samples of malaria patients of different severity levels collected in 3-6 hour resolution for at least two environmental day-night cycles can provide an inclusive representation of dysfunctions in the circadian timekeeping system with malaria severity. Rhythmic molecular markers discussed above could be used as an indicator of a healthy circadian clock to evaluate circadian disruptions in malaria patients. In such case-control longitudinal circadian studies, clock parameters in malaria patients need to be compared against age and gender-matched healthy and febrile controls to provide meaningful information regarding the specificity and adversity of circadian aberrations in core clock machinery and clock-regulated pathways in Plasmodium infections.

Conclusion and limitations
The present study indicates a robust circadian control of various host response pathways and immunological networks that can influence the outcome and severity of clinical manifestations in malaria. We observed a significant overlap among the rhythmic host immunity genes across mice, non-human primates, and humans. Host rhythms are essential for the cyclic regulation of host-parasite interactions and defense machinery in malaria. Understanding host rhythms is critical for developing novel therapeutic strategies for malaria based on boosting the body clocks or correcting circadian disruptions in patients suffering from infection with malaria parasites.
There are a couple of limitations of the present study. The rhythmicity of host immunity factors presented here is analyzed only at the transcript levels based on the available datasets in existing literature and databases. Transcript-level observations do not always reflect rhythmicity at the protein levels, and lower correlations between the rhythmic transcriptome and proteome are observed in many genome-scale circadian studies in different model organisms (69)(70)(71). Dysfunctions at the protein levels frequently correlate most accurately with diseased conditions, and most of the existing antimalarial drugs target proteins and metabolic pathways. Consequently, there is a need for proteinlevel rhythmicity analysis of the malaria-related host factors that we have shown to be involved in multiple interconnected host response pathways and immunological networks. Besides, we have investigated the rhythmicity of malaria-related host factors in circadian transcriptome datasets under healthy conditions. It is essential to study whether those host factors follow a similar cyclic pattern under Plasmodium infections. Our present in silico study provides an inclusive understanding of 24-hour oscillations in host proteins and immune factors involved in malaria pathogenesis and offers valuable directions for defining experimental strategies for further investigating circadian disruptions in malaria patient cohorts.

Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.