Abstract
Introduction:
Malaria pathogenesis involves both parasite burden and host inflammatory and oxidative stress responses that contribute to haematological alterations. Full blood count (FBC)-derived indices have been explored as accessible surrogate biomarkers, reflecting systemic inflammatory changes; however, they provide limited mechanistic insight into underlying haem-driven oxidative stress pathways. Haem oxygenase-1 (HO-1) reflects cellular responses to free haem, while interleukin-6 (IL-6) mediates systemic inflammation. However, their combined role in relation to parasitaemia and haematological alterations in endemic populations remains insufficiently defined. This study evaluated plasma HO-1 and IL-6 as complementary host-response biomarkers in a hospital-based malaria-endemic population.
Methods:
This cross-sectional study was conducted at Covenant University Medical Centre (CUMC), Canaanland, Ota, Ogun State, Nigeria. A total of 647 individuals were screened, and 304 participants with complete microscopy and haematological data were retained in the analytical cohort (152 microscopy-positive and 152 microscopy-negative). An integrated biomarker subset of 40 participants with matched plasma samples and complete laboratory data was used for ELISA, correlation, and receiver operating characteristic (ROC) analyses. Group comparisons used Mann-Whitney U and Kruskal-Wallis tests, categorical variables were compared using chi-square or Fisher exact tests, and associations were assessed using Spearman correlation.
Results:
In the analytical cohort, microscopy-positive participants had lower haemoglobin, haematocrit, red-cell count, mean cell indices, and platelet counts, together with higher temperature and white-cell counts than microscopy-negative participants (all p< 0.001 for the major contrasts). Fever and anaemia were also substantially more frequent in microscopy-positive participants. In the biomarker subset, IL-6 was significantly higher in malaria-positive participants than malaria-negative participants (median 141.28 vs 106.58 pg/mL; p = 0.027), whereas HO-1 showed a non-significant upward shift (230.98 vs 193.73 ng/mL; p = 0.131). IL-6 correlated positively with HO-1 (rho = 0.702, p< 0.001), parasitaemia (rho = 0.387, p = 0.015), and temperature (rho = 0.432, p = 0.006), and negatively with haemoglobin (rho = -0.417, p = 0.008). ROC analysis showed moderate discriminatory performance for IL-6 and lower performance for HO-1. IL-6 yielded an AUC of 0.709 (95% CI 0.531–0.874), with sensitivity of 72.7% and specificity of 76.5% at the Youden-optimal threshold. HO-1 yielded an AUC of 0.641 (95% CI 0.440–0.826), while the combined IL-6 + HO-1 model yielded an AUC of 0.687 (95% CI 0.461–0.886).
Conclusions:
IL-6 showed better discriminatory performance than HO-1 in this cohort, while HO-1 demonstrated only modest and non-significant differences between groups. These findings suggest that IL-6 may be a more informative adjunct inflammatory marker in malaria, whereas the role of HO-1 requires further validation in larger studies.
1 Introduction
Malaria remains one of the most consequential infectious diseases worldwide, with the highest burden concentrated in sub-Saharan Africa. The 2025 World Malaria Report states that malaria continues to impose major morbidity and mortality, with Nigeria remaining among the countries contributing the largest share of global disease burden (, ). Yet parasite detection alone does not fully explain why some infected individuals experience greater physiological derangement than others (–).
Contemporary understanding of malaria pathogenesis has emphasised that severe or complicated disease arises from an interplay among parasite biomass, haemolysis, endothelial activation, coagulation disturbances, oxidative stress, and dysregulated inflammatory signalling (–). During blood-stage infection, parasitised erythrocytes rupture, releasing cell-free haemoglobin and free haem. Free haem is a potent pro-oxidant that can damage endothelium, amplify inflammatory cascades, and worsen tissue injury (, ). In this context, host-response biomarkers provide a mechanistic window into disease biology that is not captured by microscopy or rapid diagnostic testing alone (, ).
Haem oxygenase-1 (HO-1), encoded by HMOX1 is a stress-responsive enzyme that catalyses haem degradation to biliverdin, carbon monoxide, and free iron. In malaria, HO-1 has been interpreted as a double-edged pathway: experimental and translational studies are consistent with protective effects through detoxification of free haem and mitigation of oxidative injury, while human studies have also linked heightened HO-1 expression with severe disease states in specific clinical settings (, ). This duality makes HO-1 biologically useful as both a mechanistic readout and a candidate adjunct biomarker (, ).
IL-6 is a multifunctional inflammatory cytokine involved in fever generation, hepatic acute-phase signalling, lymphocyte activation, endothelial perturbation, and the orchestration of innate and adaptive responses (, ). Earlier clinical studies demonstrated that circulating IL-6 is elevated in complicated malaria (, ), while newer systematic reviews and meta-analyses support an overall association between higher IL-6 concentrations and malaria-associated physiological disturbance (, ).
Despite the expanding literature on host biomarkers in malaria, several important gaps remain. First, many studies have investigated either inflammatory mediators or oxidative-stress pathways independently rather than examining both biological axes within an integrated host-response framework. Second, a large proportion of biomarker studies have been conducted in narrowly defined populations such as paediatric severe malaria cohorts () or pregnancy-associated malaria (, ), where host responses may differ substantially from those observed in general endemic populations.
Other investigations have focused on specific clinical contexts, including studies evaluating inflammatory cytokines as indicators of severe malaria in hospitalised paediatric populations. For example, Obeng-Aboagye et al. examined inflammatory cytokines as potential biomarkers for early diagnosis of severe malaria in children (), while studies by Oluboyo et al. and Varo et al. investigated endothelial activation markers such as angiopoietins in relation to severe malaria pathophysiology and clinical outcomes (, ). In addition, Kortz et al. compared host biomarker profiles in malaria and non-malarial sepsis, highlighting differences in inflammatory and endothelial responses across febrile illnesses ().
Generally, only a limited number of studies have evaluated the diagnostic or discriminatory performance of host biomarkers. For instance, Foko et al. identified CRP, Ang-2, Ang-2/1 ratio, PfHRP2, and platelet count as high-value biomarkers for predicting severe malaria and cerebral malaria (). Systematic and meta-analytic assessments of inflammatory cytokines propose them as potential biomarkers for malaria disease (, ). This limitation is particularly relevant in endemic outpatient and mixed-population settings, where host biomarker signals may differ from those observed in highly selected severe-malaria cohorts. To the best of our knowledge, this study represents one of the earliest investigations to simultaneously evaluate plasma haem oxygenase-1 (HO-1) and interleukin-6 (IL-6) and assess their joint diagnostic performance in a general malaria-exposed population. Herein, we investigate plasma HO-1 and IL-6 levels in a hospital-based malaria-endemic population in southwestern Nigeria. Malaria-associated physiological disturbance was assessed using parasite density, haemoglobin concentration, platelet count, and fever. The study aimed to determine (i) whether plasma HO-1 and IL-6 are elevated in malaria infection, (ii) whether these biomarkers are associated with parasite-density and haematological indicators of malaria-associated physiological disturbance, and (iii) whether their combined diagnostic performance supports their potential utility as adjunct host biomarkers in endemic clinical settings.
2 Materials and methods
2.1 Study design and setting
This study was conducted among residents of Canaanland and its adjoining communities in Ota, Ogun State, Nigeria (Figure 1). Canaanland is a prominent religious and educational enclave located in Ado-Odo/Ota Local Government Area, approximately 40 km from Lagos, Nigeria’s commercial capital. The area hosts the international headquarters of Living Faith Church Worldwide (Winners’ Chapel), Covenant University, Faith Academy Secondary School, and Kingdom Heritage Nursery and Primary School, forming a large residential and academic community with an estimated population of approximately 11, 000 residents, including church officials, academic and support staff, and students.
Figure 1
Geographically, Canaanland lies between latitudes 6°40′–6°42′N and longitudes 3°08′–3°10′E. The region experiences a tropical wet-and-dry climate, with a prolonged rainy season from April to October and a shorter dry season from November to March. These climatic conditions favour the breeding of Anopheles mosquito vectors, contributing to stable malaria transmission in the region (
The study was conducted between March and December 2025, spanning the onset and progression of the rainy season in southern Nigeria, a period often associated with increased malaria transmission due to enhanced mosquito breeding conditions (
2.2 Study population
Participant recruitment was conducted among individuals presenting for malaria screening or routine clinical evaluation at the Medical Centre and during community health outreaches within Canaanland.
Participants were recruited from two sources: (i) individuals presenting for clinical evaluation at the health facility, and (ii) participants enrolled through community outreach screening. Individuals presenting at the health facility with fever or symptoms suggestive of malaria, or with an axillary temperature ≥ 37.5 °C, were considered clinically eligible (
A total of 647 individuals meeting these initial eligibility criteria were screened for malaria. Participants who were currently receiving antimalarial therapy or had taken antimalarial medication within the preceding three weeks were excluded to minimise potential effects on parasite detection and biomarker measurements. In addition, individuals with microscopy-confirmed mixed Plasmodium infections were excluded to ensure consistency in the interpretation of malaria-associated host responses.
An integrated biomarker subset of 40 participants was defined using a complete-case approach, including only individuals with fully matched plasma IL-6 and HO-1 measurements, parasitaemia data, and haematological profiles. Demographic variables, including age and sex, were retained (
Figure 2

Study flow diagram.
2.3 Sample collection and processing
Approximately 5 mL of venous blood was collected aseptically from each participant by venipuncture into ethylenediaminetetraacetic acid (EDTA)-anticoagulated tubes, following standard phlebotomy procedures. EDTA-anticoagulated blood was used for malaria microscopy and full blood count analysis. For plasma preparation, blood was collected in lithium-heparin tubes and was centrifuged at 3000 rpm for 10 minutes, separated and then aliquoted into sterile cryovials before being stored at −20 °C until biomarker analysis.
2.4 Malaria diagnosis and parasite quantification
Malaria screening was initially performed using the CareStart Malaria HRP2/pLDH (Pf/PAN) Combo rapid diagnostic test (Access Bio, Somerset, NJ, USA). Microscopic examination, considered the gold standard for malaria diagnosis, was used to define malaria status for the analytical comparisons in this manuscript. Thick and thin blood films were prepared on clean microscope slides using standard laboratory techniques (
2.5 Haematological measurements
Full blood counts were performed using the Sysmex XS-1000i automated haematology analyser (Sysmex Corporation, Kobe, Japan). Parameters considered in the present work included haemoglobin concentration, red blood cell count, white blood cell count, and platelet count. Fever was recorded clinically and considered alongside parasite density, haemoglobin, and platelet count as a malaria-associated physiological disturbance-related variable (
2.6 Definition of parasite burden and clinical indicators
Malaria burden in this study was assessed using parasite density together with selected haematological and clinical indicators, including haemoglobin concentration, platelet count, and fever. These parameters are widely recognised as important clinical correlations of malaria-associated haemolysis, thrombocytopenia, and disease progression in endemic settings (
2.6.1 Parasite density estimation
Parasite density was determined from Giemsa-stained thick blood smears by counting the number of malaria parasites relative to 200 white blood cells (WBCs) under oil immersion microscopy. An assumed average leukocyte count of 8000 WBCs/µL was used for the calculation. Parasite density per microlitre of blood was estimated using the formula:
For analytical purposes, parasite density was categorized as low (<1, 000 parasites/µL), moderate (1, 000–9, 999 parasites/µL), high (10, 000–99, 999 parasites/µL), and hyperparasitaemia (≥100, 000 parasites/µL) (
2.6.2 Haemoglobin measurement and anaemia classification
Haemoglobin (Hb) concentration was measured, and each participant’s Hb value was recorded in grams per decilitre (g/dL). Anaemia was classified using haemoglobin thresholds as follows: mild anaemia, 10.0–10.9 g/dL; moderate anaemia, 7.0–9.9 g/dL; and severe anaemia, <7.0 g/dL (
2.7 Measurement of plasma IL-6 and HO-1
Plasma concentrations of haem oxygenase-1 (HO-1) were measured using a Human Haem Oxygenase-1 ELISA kit (Boster Biological Technology, USA, Cat. No. EKC33968) according to the manufacturer’s instructions. The assay was performed using a 96-well microplate format designed for the quantitative detection of HO-1 in human serum and plasma samples. The kit has a detection range of 15.6–1000 ng/mL with an analytical sensitivity of 3.9 ng/mL.
Plasma interleukin-6 (IL-6) concentrations were measured using the Human IL-6 DuoSet ELISA Kit (R&D Systems, Cat. No. DY206, Minneapolis, MN, USA). Samples and standards were assayed in duplicate according to the manufacturer’s instructions.
Optical densities for both analytes were recorded at 450 nm and 630 nm, with 630 nm serving as the reference wavelength for background correction, using a BioTek ELx800 microplate reader (BioTek Instruments, USA). Concentrations were interpolated from four-parameter logistic standard curves (
2.8 Evaluation of diagnostic performance
Receiver operating characteristic (ROC) curve analysis was performed to evaluate the diagnostic performance of IL-6, HO-1, and their combined model. The area under the curve (AUC) was calculated for each biomarker. Optimal cut-off values were determined using Youden’s index, and corresponding sensitivity and specificity were reported. Diagnostic performance was assessed based on the ability of each biomarker to discriminate between malaria-positive and malaria-negative individuals, consistent with established approaches (
Table 1
| AUC range | Interpretation |
|---|---|
| 0.90–1.00 | Excellent |
| 0.80–0.89 | Very good |
| 0.70–0.79 | Good |
| 0.60–0.69 | Satisfactory |
| 0.50–0.59 | Poor/unsatisfactory |
Interpretation of area under the curve (AUC) values for diagnostic performance .
2.9 Statistical analysis
All statistical analyses were performed independently of laboratory procedures using Python version 3.14.2 with NumPy 2.4.1, Pandas 2.3.3, SciPy 1.17.0, Matplotlib 3.10.8, and Scikit-learn 1.8.0.
Continuous variables were non-normally distributed and are therefore summarized as medians and interquartile ranges (IQR), while categorical variables are presented as frequencies and percentages. Comparisons between microscopy-positive and microscopy-negative participants were performed using the Mann-Whitney U test for continuous variables and chi-square or Fisher exact tests for categorical variables, as appropriate. Comparisons across parasite-density categories were conducted using the Kruskal-Wallis test.
Associations between biomarkers and clinical variables were assessed using Spearman’s rank correlation coefficient (rho). Diagnostic performance was evaluated using receiver operating characteristic (ROC) analysis, with area under the curve (AUC), 95% confidence intervals, sensitivity, specificity, and optimal cutoffs defined by the Youden index (
3 Results
3.1 Cohort characteristics and clinical profile
The analytical cohort comprised 304 participants: 152 microscopy-positive and 152 microscopy-negative. Age and sex distributions were similar between groups, but microscopy-positive participants had substantially greater physiological disturbance. Fever was present in 68.4% of microscopy-positive participants compared with 16.4% of microscopy-negative participants, and any anaemia was present in 78.9% versus 12.5%, respectively. Summary statistics of haematological indices, fever, and anaemia status across malaria-positive and malaria-negative participants are presented in Tables 2A, B.
Table 2A
| Variable | Malaria-positive (n=152) | Malaria-negative (n=152) | P |
|---|---|---|---|
| Age (years) | 28.00 (14.00-40.00) | 27.00 (16.75-38.00) | 0.918 |
| Temperature (°C) | 38.59 (36.87-39.13) | 36.74 (36.51-37.12) | <0.001** |
| Red blood cells (×1012/L) | 4.07 (3.74-4.42) | 5.04 (4.64-5.43) | <0.001** |
| Haemoglobin (g/dL) | 10.20 (8.30-11.71) | 13.79 (12.55-14.88) | <0.001** |
| Haematocrit (%) | 33.12 (28.55-36.37) | 42.05 (39.31-44.81) | <0.001** |
| MCV (fL) | 81.75 (77.11-89.35) | 88.88 (84.79-94.06) | <0.001** |
| MCH (pg) | 27.93 (24.73-29.90) | 30.08 (28.33-31.15) | <0.001** |
| MCHC (g/dL) | 32.18 (30.98-33.59) | 33.99 (33.10-35.11) | <0.001** |
| White blood cells (×109/L) | 7.85 (6.16-9.16) | 6.44 (5.33-7.62) | <0.001** |
| Neutrophils (%) | 57.49 (52.24-64.01) | 53.36 (47.63-58.68) | <0.001** |
| Lymphocytes (%) | 32.16 (26.65-38.70) | 35.62 (31.66-39.92) | <0.001** |
| Monocytes (%) | 4.80 (3.44-6.52) | 5.22 (3.33-6.52) | 0.663 |
| Eosinophils (%) | 2.33 (1.44-3.14) | 2.94 (2.10-3.99) | <0.001** |
| Platelets (×109/L) | 142.94 (103.65-173.38) | 260.48 (222.82-308.96) | <0.001** |
Continuous clinical and haematological variables by microscopy-defined malaria status.
Values are median (IQR). P values are from Mann-Whitney U tests. All significant comparisons in this table remained significant after false-discovery-rate correction [55]. Strong significant differences are marked with **.
Table 2B
| Variable | Malaria-positive (n=152) | Malaria-negative (n=152) | P |
|---|---|---|---|
| Male sex, n (%) | 68 (44.7%) | 77 (50.7%) | 0.358 |
| Fever present, n (%) | 104 (68.4%) | 25 (16.4%) | <0.001 |
| Anaemia, n (%) | 120 (78.9%) | 19 (12.5%) | <0.001** |
| Anaemia severity: None, n (%) | 32 (21.1%) | 133 (87.5%) | <0.001** |
| Mild anaemia, n (%) | 23 (15.1%) | 18 (11.8%) | |
| Moderate anaemia, n (%) | 68 (44.7%) | 1 (0.7%) | |
| Severe anaemia, n (%) | 29 (19.1%) | 0 (0.0%) |
Categorical variables by microscopy-defined malaria status.
P values are from Fisher exact tests except anaemia-malaria-associated physiological disturbance distribution, which uses a chi-square test across categories [56]. ** p < 0.001.
Microscopy-positive participants demonstrated significantly lower red blood cell count, haematocrit, mean corpuscular volume (MCV), mean corpuscular haemoglobin (MCH), and mean corpuscular haemoglobin concentration (MCHC), alongside significantly higher body temperature and white blood cell count (all p< 0.001). Marked reductions in haemoglobin and platelet counts among microscopy-positive participants are illustrated in Figure 3.
Figure 3

Distribution of haemoglobin (A) and platelet count (B) by microscopy-confirmed malaria status.
3.2 Plasma IL-6 and HO-1 levels by malaria status
The integrated biomarker subset included 40 participants (22 malaria-positive and 18 malaria-negative). Plasma IL-6 and HO-1 concentrations stratified by malaria status are presented in Table 3, which also includes the demographic characteristics of the biomarker subset. Plasma IL-6 levels were significantly higher in malaria-positive participants compared to malaria-negative individuals (median: 141.28 vs 106.58 pg/mL; p = 0.027). In contrast, HO-1 levels showed a higher median in malaria-positive participants, but this difference was not statistically significant (median: 230.98 vs 193.73 ng/mL; p = 0.131).
Table 3
| Variable | Malaria-positive (n = 22) | Malaria-negative (n = 18) | p-value |
|---|---|---|---|
| Age (years) | 29 (21–38) | 22 (12–38) | 0.18 |
| Sex (Male/Female) | 12/10 | 7/11 | 0.34 |
| IL-6 (pg/mL) | 141.28 (113.60–173.44) | 106.58 (10.00–118.48) | 0.027* |
| HO-1 (ng/mL) | 230.98 (205.33–295.73) | 193.73 (168.55–278.77) | 0.131 |
Plasma biomarker concentrations and demographic characteristics by malaria status in the integrated subset.
*Statistically significant at p< 0.05. Values are presented as median (IQR). p-values for continuous variables were obtained using the Mann–Whitney U test, and categorical variables using the Chi-square test. IL-6 values were available for 39 participants because one malaria-negative sample had missing IL-6 data.
3.3 Variation of IL-6 and HO-1 across parasite-density categories
Among the 22 microscopy-positive participants included in the biomarker subset, parasite density was categorized as low (n = 7), moderate (n = 4), high (n = 7), or hyperparasitaemia (n = 4) (
Figure 4

Plasma IL-6 (A) and HO-1 (B) concentrations across parasite-density categories among microscopy-positive participants.
3.4 Correlation analysis of biomarkers and clinical parameters
Table 4 presents Spearman correlation coefficients between biomarkers and key clinical variables in the integrated biomarker subset. IL-6 was positively correlated with HO-1, parasitaemia, and temperature, and negatively correlated with haemoglobin. HO-1 was positively correlated with temperature and negatively correlated with haemoglobin, whereas its correlation with parasitaemia did not reach statistical significance. The strongest association was between IL-6 and HO-1 (rho = 0.702, p< 0.001). IL-6 correlated with parasitaemia (rho = 0.387, p = 0.015), temperature (rho = 0.432, p = 0.006), and lower haemoglobin (rho = -0.417, p = 0.008). HO-1 correlated with temperature (rho = 0.412, p = 0.008) and lower haemoglobin (rho = -0.367, p = 0.020).
Table 4
| Biomarker | Variable | N | Spearman rho | p | FDR q |
|---|---|---|---|---|---|
| IL-6 (pg/mL) | HO-1 (ng/mL) | 39 | 0.702 | <0.001** | <0.001 |
| Parasitaemia (/µL) | 39 | 0.387 | 0.015* | 0.025 | |
| Temperature (°C) | 39 | 0.432 | 0.006* | 0.017 | |
| Haemoglobin (g/dL) | 39 | -0.417 | 0.008* | 0.017 | |
| Platelets (×109/L) | 39 | -0.261 | 0.109 | 0.121 | |
| HO-1 (ng/mL) | Parasitaemia (/µL) | 40 | 0.267 | 0.096 | 0.120 |
| Temperature (°C) | 40 | 0.412 | 0.008* | 0.017 | |
| Haemoglobin (g/dL) | 40 | -0.367 | 0.020* | 0.028 | |
| Platelets (×109/L) | 40 | -0.216 | 0.181 | 0.181 |
Key Spearman correlations in the integrated biomarker subset.
Values shown are Spearman’s rank correlation coefficients (rho). FDR q-values represent p-values adjusted for multiple comparisons using the false discovery rate approach. Associations with q< 0.05 were considered robust after correction for multiple testing. *p< 0.05; **p< 0.001.
3.5 Diagnostic performance of biomarkers (ROC analysis)
Receiver operating characteristic (ROC) curve analysis was performed to evaluate the diagnostic performance of IL-6, HO-1, and their combined model in distinguishing malaria-positive from malaria-negative participants. As shown in Figure 5, IL-6 demonstrated the highest discriminative ability, followed by the combined IL-6 + HO-1 model, while HO-1 alone showed lower performance.
Figure 5

Receiver operating characteristic (ROC) curves for IL-6, HO-1, and their combined model.
ROC curves illustrating the diagnostic performance of IL-6 (AUC = 0.709), HO-1 (AUC = 0.641), and the combined IL-6 + HO-1 model (AUC = 0.687) in discriminating malaria-positive from malaria-negative participants. The dashed red line represents the line of no discrimination (AUC = 0.5).
IL-6 demonstrated acceptable discriminatory performance, with an area under the curve (AUC) of 0.709, sensitivity of 72.7%, and specificity of 76.5%. In contrast, HO-1 showed more modest discrimination, with an AUC of 0.641, sensitivity of 77.3%, and specificity of 66.7%. The combined IL-6 + HO-1 logistic regression model yielded an AUC of 0.687 (95% CI 0.456–0.882), with improved sensitivity (86.4%) but reduced specificity (58.8%).
4 Discussion
This study evaluated plasma biomarkers haem oxygenase-1 (HO-1) and interleukin-6 (IL-6) as potential adjunct indicators of malaria disease burden in a hospital-based, endemic population in southwestern Nigeria.
The findings show that plasma IL-6 was significantly elevated in malaria-positive participants, whereas HO-1 showed a higher median concentration that did not reach statistical significance. This suggests that IL-6 more clearly distinguished malaria status in this cohort, while HO-1 may reflect a more variable host-response pattern. Our IL-6 findings are consistent with previous studies showing increased circulating IL-6 in malaria infection (
Taken together, these findings support IL-6 as a more responsive marker of malaria-associated inflammation in this dataset (
More recently, Wilairatana et al. synthesised available studies in a systematic review and meta-analysis and concluded that IL-6 is significantly increased in severe compared with non-severe malaria (
At the same time, our study contributes nuance. That is, current findings extend previous observations by providing additional evidence on the relationship between these biomarkers and malaria-associated physiological changes. While IL-6 in our integrated subset showed acceptable discrimination for malaria positivity, its AUC did not approach the level required for standalone clinical replacement of parasite-based tests. This is consistent with the broader severe-malaria biomarker literature, which increasingly recognizes that host biomarkers are complementary rather than substitutive (
The HO-1 findings are particularly informative when interpreted in the context of haemolysis-driven oxidative stress during malaria infection. Experimental and translational work has established a strong mechanistic basis for HO-1 involvement in malaria pathogenesis. Ferreira et al. proposed that free haem is a central driver of severe malaria pathophysiology and demonstrated that induction of HO-1 and carbon monoxide generation can be protective in experimental models (
However, human data also suggest that the role of HO-1 in malaria is complex. Walther et al. reported that higher HO-1 levels and promoter variants associated with greater inducibility were linked to severe malaria in Gambian children (
In comparison with more recent biomarker studies, the current findings situate the novelty of this work, which brings additional perspective by examining both inflammatory (IL-6) and oxidative stress (HO-1) pathways within the same host-response framework. Nahm (62) emphasised that ROC analysis is useful for comparing diagnostic performance and selecting cut-off values, but that AUC should be interpreted carefully within the clinical context rather than overstated as proof of strong diagnostic utility. Çorbacıoğlu and Aksel (63) further noted that modest AUC values indicate limited discrimination, supporting cautious interpretation when biomarker performance does not clearly separate groups. Consistent with this, IL-6 showed better discriminatory performance than HO-1, whereas HO-1 demonstrated only limited discrimination and no significant difference between groups. The combined model did not meaningfully improve performance beyond IL-6 alone (Figure 5).
Contemporary studies have explored endothelial or multi-protein signatures in severe paediatric malaria or malaria in pregnancy. Kanoi et al. identified HO-1 among differentially expressed plasma proteins in pregnant women with falciparum malaria using high-throughput proteomics (
An important methodological consideration in this study relates to how malaria-associated physiological disturbance was conceptualised. Because the cohort was not assembled as a classical severe-malaria case series, formal classification using WHO-defined severe malaria criteria was not feasible. Instead, malaria-associated physiological disturbance was inferred from parasite density, interpreted according to WHO-recommended thresholds (
5 Study limitations
This study has several limitations. First, the cross-sectional design does not permit temporal inference regarding biomarker dynamics during malaria infection. Second, IL-6 is a non-specific inflammatory marker, and alternative febrile comorbidities were not systematically excluded (67). Third, the integrated biomarker subset was relatively small and should be regarded as exploratory. Accordingly, these findings should be interpreted cautiously and validated in larger prospective cohorts. Future studies should include larger prospective cohorts, explicit WHO severe-malaria phenotyping, and comparison with other host-response candidates such as angiopoietins, C-reactive protein (CRP), and endothelial activation markers (
6 Conclusion
Plasma IL-6 was significantly associated with malaria positivity, whereas HO-1 showed a non-significant upward trend and limited discriminatory ability. These findings support IL-6 as a clearer inflammatory marker in this cohort, while HO-1 may reflect oxidative or haemolytic host stress but should not be interpreted as a strong standalone discriminator based on the present data. Larger studies are needed to clarify the clinical relevance of HO-1 and the potential added value of combined host biomarkers. These findings support IL-6 as a clearer inflammatory marker in this cohort and may provide additional biological insight beyond parasite detection and conventional haematological indices in understanding malaria-associated physiological disturbance in endemic settings. Larger prospective studies are needed to validate these observations and clarify their clinical relevance.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by Covenant University Health Research Ethics Committee (CHREC/432/2024). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.
Author contributions
TW: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Writing – original draft, Writing – review & editing. IK: Data curation, Formal analysis, Methodology, Visualization, Writing – original draft, Writing – review & editing. TO: Data curation, Formal analysis, Software, Visualization, Writing – original draft, Writing – review & editing. SC: Conceptualization, Supervision, Validation, Writing – original draft, Writing – review & editing. IA: Conceptualization, Project administration, Supervision, Validation, Writing – original draft, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This work was supported by Covenant Applied Informatics and Communication–Africa Centre of Excellence (CApIC-ACE), Covenant University, Nigeria, through the World Bank IMPACT Project scholarship awarded to TW, which supported the laboratory materials and analyses. In addition, the Covenant University Centre for Research, Innovation and Discovery (CUCRID) has offered up to USD 1, 000 toward the article processing charge (APC).
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Abbreviations
AUC, area under the curve; CUMC, Covenant University Medical Centre; ELISA, enzyme-linked immunosorbent assay; FDR, false-discovery rate; Hb, haemoglobin; HO-1, haem oxygenase-1; IL-6, interleukin-6; IQR, interquartile range; ROC, receiver operating characteristic; RDT, rapid diagnostic test.
References
1
BamikoleTH. Confronting malaria – addressing a critical health crisis among vulnerable groups in Nigeria. In: Neglected tropical diseases - unsolved debts for the one health approach. Geneva: IntechOpen (2024). doi: 10.5772/intechopen.1002999
2
World Health Organization. World malaria report 2025. Geneva: World Health Organization (2025). Available online at: https://www.who.int/publications/i/item/world-malaria-report-2025 (Accessed February 11, 2026).
3
CunningtonAJ. The importance of pathogen load. PloS Pathog. (2015) 11:e1004563. doi: 10.1371/journal.ppat.1004563. PMID:
4
SeneISCostaDLZacariasDAdos SantosJCFerreiraGRAndradeDRet al. Immune–pathological correlates of disease severity in new-world kala-azar: the role of parasite load and cytokine profiles. Pathogens. (2025) 14:615. doi: 10.3390/pathogens14070615. PMID:
5
DokunmuTMOlasehindeGIOladejoDOAdjekukorCUAkinbohunAEOnileereOAet al. Evaluation of Plasmodium falciparum K13 gene polymorphism and susceptibility to dihydroartemisinin in an endemic area. Biomed Res Ther. (2018) 5:2651–7. doi: 10.15419/bmrat.v5i9.474
6
VaroRSitoeAMadridLAidePCuambaICossaA. Host biomarkers and parasite biomass are associated with severe malaria in Mozambican children: a case-control study. Sci Rep. (2025) 15:14262. doi: 10.1038/s41598-025-98154-1. PMID:
7
Garcia-SantosDChiesJAB. HO-1 polymorphism as a genetic determinant behind the malaria resistance afforded by haemolytic disorders. Med Hypotheses. (2010) 74:807–13. doi: 10.1016/j.mehy.2009.12.010. PMID:
8
BoachieJAhiableDAjabuinLAAmissahRAsmah-BrownAApalebilahSet al. Diagnostic value of full blood count derived systemic inflammatory biomarkers in malaria infection. Pract Lab Med. (2025) 46:e00494. doi: 10.1016/j.plabm.2025.e00494. PMID:
9
FerreiraABallaJJeneyVBallaGSoaresMP. A central role for free heme in the pathogenesis of severe malaria: the missing link? J Mol Med (Berl). (2008) 86:1097–111. doi: 10.1007/s00109-008-0368-5. PMID:
10
VasquezMZunigaMRodriguezA. Oxidative stress and pathogenesis in malaria. Front Cell Infect Microbiol. (2021) 11:768182. doi: 10.3389/fcimb.2021.768182. PMID:
11
OjeniyiFDAyoolaAOIbitoyeOOpaleyeOOOloweOAEhigieLOet al. Performance and challenges of malaria rapid diagnostic tests in endemic regions of Africa. (2025) 768182. doi: 10.1101/2025.07.30.25331638
12
ParijaSCPoddarA. Parasitic diagnosis: a journey from basic microscopy to cutting-edge technology. Trop Parasitol. (2025) 15:71–80. doi: 10.4103/tp.tp_61_24. PMID:
13
PáduaTASouzaMC. Heme on pulmonary malaria: friend or foe? Front Immunol. (2020) 11:1835. doi: 10.3389/fimmu.2020.01835. PMID:
14
TagamiYHaraYMurohashiKNagasawaRFujiiHIzawaAet al. Serum heme oxygenase-1 as a prognostic biomarker in patients with acute exacerbation of interstitial lung disease. Sci Rep. (2023) 13:22639. doi: 10.1038/s41598-023-49342-4. PMID:
15
WakaiTChineduSAfolabiI. Heme oxygenase-1 and malaria pathogenesis. Front Immunol. (2026) 17:1777520. doi: 10.3389/fimmu.2026.1777520. PMID:
16
WilairatanaP. Increased interleukin-6 levels associated with malaria infection and disease severity: a systematic review and meta-analysis. Sci Rep. (2022) 12:5983. doi: 10.1038/s41598-022-09848-9. PMID:
17
TanakaTNarazakiMKishimotoT. IL-6 in inflammation, immunity, and disease. Cold Spring Harbor Perspect Biol. (2014) 6:a016295. doi: 10.1101/cshperspect.a016295. PMID:
18
MshanaRNBoulandiJMshanaNMMayomboJMendomeG. Cytokines in the pathogenesis of malaria: levels of IL-1 beta, IL-4, IL-6, TNF-alpha and IFN-gamma in plasma of healthy individuals and malaria patients in a holoendemic area. J Clin Lab Immunol. (1991) 34:131–9.
19
WenischCLinnauKFLooareesuwanSRumpoldH. Plasma levels of the interleukin-6 cytokine family in persons with severe Plasmodium falciparum malaria. J Infect Dis. (1999) 179:747–50. doi: 10.1086/314630. PMID:
20
HashmiFAqeelSZuberiUFKhanW. A systematic review and meta-analysis of inflammatory biomarkers associated with malaria infection and disease severity. Cytokine. (2023) 169:156305. doi: 10.1016/j.cyto.2023.156305. PMID:
21
FokoLPKNarangGTamangSHawadakJJakhanJSharmaAet al. The spectrum of clinical biomarkers in severe malaria and new avenues for exploration. Virulence. (2022) 13:634–53. doi: 10.1080/21505594.2022.2056966. PMID:
22
WaltherMCaulAAkaPNjieMAmambua-NgwaAWaltherB. HMOX1 gene promoter alleles and high HO-1 levels are associated with severe malaria in Gambian children. PloS One. (2012) 8(3):e1002579. doi: 10.1371/journal.ppat.1002579. PMID:
23
AubouyAOlagnierDBertinGEzinmegnonSMajorelCMimarS. Nrf2-driven CD36 and HO-1 gene expression in circulating monocytes correlates with favourable clinical outcome in pregnancy-associated malaria. Malar J. (2015) 14:358. doi: 10.1186/s12936-015-0888-8. PMID:
24
KanoiBNWaweruHKobiaFMMukalaJKiriraPMogereDet al. Plasma protein biomarkers of Plasmodium falciparum infection in pregnant women: a high-throughput proteomics study. Front Cell Infect Microbiol. (2025) 15:1594088. doi: 10.3389/fcimb.2025.1594088. PMID:
25
Obeng-AboagyeEFrimpongAAmponsahJADansoSEOwusuEDAOforiMF. Inflammatory cytokines as potential biomarkers for early diagnosis of severe malaria in children in Ghana. Malar J. (2023) 22(1):220. doi: 10.1186/s12936-023-04652-w. PMID:
26
OluboyoAOChukwuSIOluboyoBOOdewusiOO. Evaluation of angiopoietins 1 and 2 in malaria-infested children. J Environ Public Health. (2020) 2020:2169763. doi: 10.1155/2020/2169763. PMID:
27
KortzTBNyirendaJTemboDElfvingKBaltzellKBandaweG. Distinct biomarker profiles distinguish Malawian children with malarial and non-malarial sepsis. Am J Trop Med Hyg. (2019) 101:1424–33. doi: 10.4269/ajtmh.18-0635. PMID:
28
AkinbileCOOlatunjiASOjoAO. Effects of industrialization on groundwater quality in Shagamu and Ota industrial areas of Ogun state, Nigeria. Heliyon. (2020) 6(7):e04353. doi: 10.1016/j.heliyon.2020.e04353
29
BelloABHassanAAOladejoSOAwoniranDRPopoolaKOK. Influence of meteorological variables on the abundance and distribution patterns of Anopheles mosquitoes in Epe and Orimedu, Lagos, Nigeria. Discover Anim. (2025) 2:28. doi: 10.1007/s44338-025-00076-9. PMID:
30
OforkaCLOmotayoAIAdelekeMA. Seasonal diversity in mosquito larval ecology and its public health implications in urban slums of Lagos, Nigeria. Am J Trop Med Hyg. (2024) 110:448. doi: 10.4269/ajtmh.23-0192. PMID:
31
OranusiSUMamehEOOyegbadeSABalogunDOAririguzohV-GAtokoloA. Prevalence of malaria infection among symptomatic patients of selected healthcare centers in Ota, Ogun State Nigeria: proactive prevention strategies. In: IsiborPOAkindutiPA, editors.Harnessing biotechnology tools for product development: selected papers from the international biotechnology conference exhibition and workshop (IBCEW 2024). Springer Nature Switzerland, Cham (2025). p. 333–44. doi: 10.1007/978-3-031-86002-7_19
32
WakaiTFiamitiaCKintungIJohngweMChineduSAfolabiI. Knowledge, practices, and perceptions towards malaria prevention and control among residents of Canaanland and surrounding areas in Ota, Ogun State, Nigeria: a cross-sectional study. FrontTrop Dis. (2025) 6:1686197. doi: 10.3389/fitd.2025.1686197
33
OlasehindeGI. Studies on incidence of malaria and comparative efficacy of diagnostic test methods for plasmodium falciparum and p. vivax. Int J Infect Dis. (2016) 45:154. doi: 10.1016/j.ijid.2016.02.371. PMID:
34
KiemdeFTahitaMCLompoPRouambaTSomeAMTintoHet al. Treatable causes of fever among children under five years in a seasonal malaria transmission area in Burkina Faso. Infect Dis Poverty. (2018) 7:60. doi: 10.1186/s40249-018-0442-3. PMID:
35
OkoroCIIhenetuFCDungaKEAchigbuKObasiCCOdinakaKKet al. Fever cases associated with Plasmodium falciparum malaria infection among children attending a tertiary health facility in Imo State, Nigeria. Int J Trop Dis Health. (2020) 6:7–14. doi: 10.9734/ijtdh/2020/v41i530274. PMID:
36
RanganathanPHunsbergerS. Handling missing data in research. Perspect Clin Res. (2024) 15:99–101. doi: 10.4103/picr.picr_38_24. PMID:
37
KhogaliMZachariahRKeiluhuAVan den BrandeKTayler-SmithKAyadaLet al. Detection of malaria in relation to fever and grade of malnutrition among malnourished children in Ethiopia. Public Health Action. (2011) 1:16–8. doi: 10.5588/pha.11.0008. PMID:
38
SinghRK. Malaria microscopy: a challenge in malaria elimination. Trans R Soc Trop Med Hygiene. (2025) 119:e8–9. doi: 10.1093/trstmh/traf014. PMID:
39
CheesbroughM. District laboratory practice in tropical countries, part 2. Cambridge: Cambridge University Press (2006). doi: 10.1017/CBO9780511543470
40
OyiboWLathamVOladipoONtadomGUhomoibhiPOgbulaforNet al. Malaria parasite density and detailed qualitative microscopy enhances large-scale profiling of infection endemicity in Nigeria. Sci Rep. (2023) 13:1599. doi: 10.1038/s41598-023-27535-1. PMID:
41
SumbeleIUNTehRNNkeudemGASandieSMMoyehMNSheyRAet al. Asymptomatic and sub-microscopic Plasmodium falciparum infection in children in the Mount Cameroon area: a cross-sectional study on altitudinal influence, haematological parameters and risk factors. Malar J. (2021) 20:382. doi: 10.1186/s12936-021-03916-7. PMID:
42
AchameMSGedefieADebashHTesfayeATirunehKTKassawAB. Prevalence of thrombocytopenia among patients with malaria in Ethiopia: a systematic review and metanalysis. Malar J. (2025) 24:61. doi: 10.1186/s12936-025-05296-8. PMID:
43
EyongJEEEbaiCBAwahPKadaAS. Prevalence, awareness, treatment-seeking behaviours and its implications in the control of malaria in Dumbu community, Donga-Mantung Division, North West Region, Cameroon. Malar J. (2025) 24:161. doi: 10.1186/s12936-025-05249-1. PMID:
44
WhiteNJ. Severe malaria. Malar J. (2022) 21:284. doi: 10.1186/s12936-022-04301-8. PMID:
45
KonaduEAObeseAODanso-CoffieCKAddisonTKyei-BaafourEObiriDet al. Reevaluating malaria severity metrics: insights from PfHRP2-derived biomass estimates in Ghanaian children. Malar J. (2026) 25:119. doi: 10.1186/s12936-026-05804-4. PMID:
46
HermanRASchererPNShanG. Evaluation of logistic and polynomial models for fitting sandwich-ELISA calibration curves. J Immunol Methods. (2008) 339:245–58. doi: 10.1016/j.jim.2008.09.001. PMID:
47
YangJLiuAPerkinsNChenZ. Youden index estimation based on group-tested data. Stat Methods Med Res. (2025) 34:45–54. doi: 10.1177/09622802241295319. PMID:
48
RuoppMDPerkinsNJWhitcombBWSchistermanEF. Youden index and optimal cut-point estimated from observations affected by a lower limit of detection. Biometrical J Biometrische Z. (2008) 50:419–30. doi: 10.1002/bimj.200710415. PMID:
49
LaiY. A statistical method for the conservative adjustment of false discovery rate (q-value). BMC Bioinf. (2017) 18:69. doi: 10.1186/s12859-017-1474-6. PMID:
50
MortazaviSELugaajjuAIvarssonA-CKarlsson SöbirkSNorrgrenHPerssonKEM. High rate of false positive malaria rapid diagnostic tests in a district hospital in Uganda. Front Malar. (2025) 3:1545825. doi: 10.3389/fmala.2025.1545825
51
WroczyńskaANahorskiWBakowskaAPietkiewiczH. Cytokines and clinical manifestations of malaria in adults with severe and uncomplicated disease. Int Marit Health. (2005) 56(1–4):103–14.
52
VallelianFBuehlerPWSchaerDJ. Hemolysis, free hemoglobin toxicity, and scavenger protein therapeutics. Blood. (2022) 140:1837–44. doi: 10.1182/blood.2022015596. PMID:
53
KernPHemmerCJDammeJGrussHJDietrichM. Elevated tumor necrosis factor alpha and interleukin-6 serum levels as markers for complicated Plasmodium falciparum malaria. Am J Med. (1989) 87:139–43. doi: 10.1016/S0002-9343(89)80688-6. PMID:
54
DayNPJHienTTSchollaardtTLocPPChuongLVChauTTH. The prognostic and pathophysiologic role of pro- and antiinflammatory cytokines in severe malaria. J Infect Dis. (1999) 180:1288–97. doi: 10.1086/315016. PMID:
55
AlansariANZaazoueeMSElshanbaryAAManiSMessaoudM. Diagnostic accuracy of interleukin-6 (IL-6) as a significant biomarker in late-onset neonatal sepsis: an updated systematic review and meta-analysis. Eur J Pediatr. (2025) 184:587. doi: 10.1007/s00431-025-06409-w. PMID:
56
DeySMazumderSSiddiquiAAIqbalMSBanerjeeCSarkarSet al. Association of heme oxygenase 1 with the restoration of liver function after damage in murine malaria by Plasmodium yoelii. Infect Immun. (2014) 82(8):3113–26. doi: 10.1128/IAI.01598-14
57
PereiraMLM. Could heme oxygenase-1 be a new target for therapeutic intervention in malaria-associated acute lung injury/acute respiratory distress syndrome? Front Cell Infect Microbiol. (2018) 8:161. doi: 10.3389/fcimb.2018.00161. PMID:
58
SinghNAhmadZBaidNKumarA. Host heme oxygenase-1: friend or foe in tackling pathogens? IUBMB Life. (2018) 70:869–80. doi: 10.1002/iub.1868. PMID:
59
Ong’echaJMDavenportGCVululeJMHittnerJBPerkinsDJ. Identification of inflammatory biomarkers for pediatric malarial anemia severity using novel statistical methods. Infect Immun. (2011) 79(11):4674–80. doi: 10.1128/IAI.05161-11
60
AwosoluOBYahayaZSFarah HaziqahMT. Prevalence, parasite density and determinants of falciparum malaria among febrile children in some peri-urban communities in southwestern Nigeria: a cross-sectional study. Infection Drug Resistance. (2021) 14:3219–32. doi: 10.2147/IDR.S312519. PMID:
61
PunnathKDayanandKKChandrashekarVNAchurRNKakkilayaSBGhoshSK. Association between inflammatory cytokine levels and thrombocytopenia during Plasmodium falciparum and Plasmodium vivax infections in South-Western Coastal Region of India. Malar Res Treat. (2019) 2019:4296523. doi: 10.4103/tp.tp_66_18. PMID:
62
NahmFS. Receiver operating characteristic curve: overview and practical use for clinicians. Korean J Anesthesiol. (2022) 75:25–36. doi: 10.4097/kja.21209. PMID:
63
ÇorbacıoğluŞKAkselG. Receiver operating characteristic curve analysis in diagnostic accuracy studies: a guide to interpreting the area under the curve value. Turk J Emerg Med. (2023) 24:1–9. doi: 10.4103/tjem.tjem_182_23. PMID:
64
SturgeonCHillRHortinGLThompsonD. Taking a new biomarker into routine use – a perspective from the routine clinical biochemistry laboratory. Proteomics Clin Appl. (2010) 4:892–903. doi: 10.1002/prca.201000073. PMID:
65
OjeniyiFDAyoolaAOIbitoyeOOpaleyeOOOloweOAEhigieLOet al. Performance and challenges of malaria rapid diagnostic tests in endemic regions of Africa. Sci Rep. (2025) 15:43663. doi: 10.1038/s41598-025-97259-x. PMID:
66
BohliJKAnsahPBKwamenaTEAlloteyEDuneehRVYeboahEBet al. Variations in haematological and inflammatory biomarkers and their association with Plasmodium falciparum malaria: a cross-sectional comparative study at a clinic in Ghana. Malar J. (2025) 24:220. doi: 10.1186/s12936-025-05474-8. PMID:
67
HarmonisJAKusumaSAFRukayadiYHasanahAN. Exploring biomarkers for malaria: advances in early detection and asymptomatic diagnosis. Biosensors. (2025) 15:106. doi: 10.3390/bios15020106. PMID:
Summary
Keywords
biomarkers, haem oxygenase-1, inflammation, interleukin-6, malaria, oxidative stress, parasite density, thrombocytopenia
Citation
Wakai T, Kintung I, Ogundimu T, Chinedu S and Afolabi I (2026) Plasma haem oxygenase-1 and interleukin-6 as adjunct host biomarkers associated with malaria. Front. Immunol. 17:1835939. doi: 10.3389/fimmu.2026.1835939
Received
21 March 2026
Revised
22 April 2026
Accepted
27 April 2026
Published
15 May 2026
Volume
17 - 2026
Edited by
Patricia Maria Lourenço Dutra, Universidade do Estado do Rio de Janeiro, Brazil
Reviewed by
Joel Karikari Nyarkoh, University of Cape Coast, Ghana
Raquel Matavele Chissumba, Centro de Investigacao e Desenvolvimento em Etnobotanica, Mozambique
Updates

Check for updates
Copyright
© 2026 Wakai, Kintung, Ogundimu, Chinedu and Afolabi.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Theophilus Wakai, theophilus.wakaipgs@stu.cu.edu.ng
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.