Abstract
HIV infection affects up to 30% of children presenting with severe acute malnutrition (SAM) in Africa and is associated with increased mortality. Children with SAM are treated similarly regardless of the presence of HIV although mechanisms of nutritional recovery in HIV and/or SAM are not well understood. We aimed to test the hypothesis that metabolic and inflammatory perturbations among children with SAM differed between those with HIV (n = 54) and without HIV (n = 113) to better understand their specific nutritional and therapeutic requirements. We performed a secondary analysis of a clinical trial and plasma proteomics data among children with complicated SAM in Kenya and Malawi. Compared to children with SAM without HIV, HIV-infected children had evidence (false discovery rate corrected p<0.05) of metabolic stress, including enriched pathways related to inflammation and lipid metabolism. Moreover, we observed reduced plasma levels of zinc-α-2-glycoprotrein, butyrylcholinesterase, and increased levels of complement C2 resembling findings in metabolic syndrome, diabetes and other non-communicable diseases. HIV was also associated (false discovery rate corrected p<0.05) with higher levels of plasma chemokines: MCP1, CCL4, GCSF, IL1b, TNFa, IL-2, 5, 7, 8 and 15, IL12p40, IP-10, and IL-1RA. Considering evidence of biomarkers of metabolic stress in children with HIV and SAM, it is of potential concern that our current treatment strategy for SAM regardless of HIV status involves a high-fat therapeutic diet. The results of this study indicate a need for clinical trials of therapeutic foods that meet the specific metabolic needs of children with HIV and SAM.
Introduction
Malnutrition, specifically undernutrition in all its forms, remains a global public health burden that accounts for 45% of all death among children under 5 years old [1]. Despite careful monitoring and adherence to guidelines set by the World Health Organization, although in general, uncomplicated SAM cases treated in the community do well, up to 25% of children with complicated severe acute malnutrition (SAM) treated in a hospital environment do not survive [2-5]. Furthermore, about one in five children treated for complicated SAM and discharged alive, die in the first year after discharge in low-resource settings [6-8]. However, our understanding of the pathophysiology underlying the poor prognosis for these children is surprisingly limited.
Infection with the human immunodeficiency virus (HIV) is a common co-morbidity of SAM in sub-Saharan Africa affecting up to 30% of admissions among SAM cases [9]. HIV-infected or exposed children are significantly more likely to be stunted, wasted, and underweight [10]. They also more often present with other clinical complications and greater susceptibility to infections, thus further complicating their clinical management, which may include providing more aggressive antimicrobial therapy and higher caloric nutritional intervention [11]. Moreover, response to clinical management is also less predictable and less well-understood in HIV-infected children compared to their uninfected counterparts [12]. Although acute opportunistic infections play a key role in the outcome of these children, intestinal pathology including inflammation and malabsorption, and metabolic perturbations may also be present. However, mechanisms driving poor nutritional recovery of children with HIV even when detected co-morbidities are treated remain poorly understood [12].
We hypothesised that inflammatory, metabolic and other pathways which are likely to be involved in the response to infection, survival and nutritional recovery differ between children with SAM with and without HIV. We conducted a secondary analysis of data and biological samples from a randomised clinical trial in Kenya and Malawi [13].
Materials and Methods
Patient recruitment and study design
This is a secondary analysis of a nested case control study from a randomised controlled trial (NCT02246296), which tested the effect of a lactose-free, low-carbohydrate F75 milk to limit carbohydrate malabsorption, diarrhoea and refeeding syndrome among children hospitalized for complicated SAM at Queen Elizabeth Central Hospital in Blantyre, Malawi, Kilifi County Hospital and Coast General Hospital, Mombasa, Kenya [13]. Children aged 6 months to 13 years were eligible for enrolment into the trial at admission to hospital if they had SAM, defined as: mid-upper arm circumference (MUAC) <11.5cm or weight-for-height Z score <−3 if younger than 5 years of age, BMI Z score <−3 if older than 5 years, or oedematous malnutrition at any age and had medical complications or failing an appetite test, as defined by WHO guidelines [14]. Children were excluded if they had a known allergy to milk products and did not provide consent. Biological samples were obtained before the children received the randomised treatment irrespective of HIV status. Unless an HIV positive status was documented, HIV status was assessed by offering an antibody test at admission. For this analysis, patients that tested positive on an HIV antibody test were considered HIV(+) and children with missing or declined HIV test were excluded.
The nested case-control study was designed to investigate inpatient mortality for which proteomic, cytokine, and chemokine data was generated using plasma samples collected at admission during enrolment to the trial. To compare the proteomic profiles between HIV affected and non-affected children with SAM, we used data from a nested case-control study to investigate inpatient mortality during the trial including 79 deaths and 88 children that were clinically stabilized and discharged from the hospital matched by site of recruitment. For this analysis, all observations from the proteomics dataset were included and the analysis was designed to help overcome selection bias, as given in detail below.
Proteomics, cytokine and chemokine analysis
Untargeted proteomics and targeted cytokines and chemokines analysis of plasma samples were performed following methods described previously [15]. The targeted protein panel included: epidermal growth factor (EGF); eotaxin; granulocyte-colony stimulating factor (GCSF); granulocyte-macrophage colony-stimulating factor (GMCSF); interferon alpha-2 (IFNa2); interferon gamma (IFNg); interleukins 10, 12p40, 12p70, 13, 15, 17A, 1a, 1b, 1RA, 2 to 8; interferon gamma-induced protein 10 (IP10); monocyte chemoattractant protein 1 (MCP1), macrophage inflammatory protein 1 alpha and beta (MIP1as & b); tumour necrosis factor alpha (TNFa) and beta (TNFb); and vascular endothelial growth factor (VEGF).
Data analysis
Data analyses were performed using R v3.5 [16]. Analysis of the prevalence of HIV(+), nutritional status and their associations with inpatient mortality utilised the entire trial dataset (N=843). Analysis of categorical data was performed using Fisher’s test and generalised linear models for continuous outcomes. Logistic regression was used to analyse binary outcomes adjusting for age, sex, presence of oedema, and site of recruitment. These associations were also adjusted for MUAC. As a sensitivity analysis to address the possibility of confounding due to HIV maternal antibodies in younger children, a test of interaction between age above or below 18 months and individual proteins towards HIV status was performed.
The proteomics, cytokines and chemokines analyses were secondary analyses of data collected from a nested case-control study with inpatient mortality as its primary outcome, hence with strong selection bias. The analysis for the association between HIV status and individual proteins was therefore performed usng logistic regression analysis with inverse probability weighting (IPW) to correct for selection bias [17-20]. Weights (w) were calculated as suggested by Samuelsen [18] wherein the weight for each observation selected into the nested case-control study was computed as the inverse of the probability of being selected for the nested study from the main clinical trial. The probability of inclusion was therefore calculated as:
where p(i) is the probability of inclusion in the nested case-control study and x1, x2 … xn are HIV status, sex, age, presence of oedema, mid-upper arm circumference, and site of recruitment of the i-th observation (child) based on the entire trial population. Inverse probability weight is therefore:
Differences in individual proteins abundances were considered statistically significant when p<0.05 after adjustment for multiple comparisons using Benjamini-Hochberg false discovery rate (FDR) [21].
Multivariate analysis was undertaken in order to determine several proteins that are collectively associated with HIV status, some of which may not be significantly associated to HIV independently. This was performed using a weighted elastic net (EN) model implemented using the “glmnet” package in R [22]. EN is a penalized regression approach that was developed to help overcome problems caused by high dimensional data. It is an integration of two regularized approaches, ridge regression and least absolute shrinkage and selection operator (LASSO), wherein the contribution of each of these models to the final EN model is controlled by the α parameter [22, 23]. The strong penalization imposed by LASSO draws coefficients to zero thereby eliminating non-predictive proteins features, whereas ridge regression addresses potential multi-collinearity problems in high-dimensional data [22, 23].
Weighted EN model generation was performed with HIV status as outcome, protein profile as predictors, and w as observation weights. The penalization parameter lambda, which influences the shrinkage of variable coefficients to zero thus eliminating some non-contributing variables, was determined by estimating the area under the receiver operator curve (ROC) of the population using ten-fold cross validation. Several alpha parameter values were assessed and a final value of 0.85 was taken to achieve a compromise between predictive ability and fewer number of features extracted. The final lambda parameter was based on the value which gave the highest area under the ROC (AUROC) value.
Proteins with significant association with HIV status after correction for false discovery and those extracted by the EN model were then uploaded to The Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 Bioinformatics Resource [24] to assess the Gene ontology (GO) enriched pathways of the differentially expressed proteins.
EN model validity was judged based on the AUROC and misclassification error rate. The fitted EN model performance measured as optimism-corrected AUC was validated using bootstrap, following the procedure of Smith et al. [25]. Bootstrapping was performed on 2000 iterations using the “BootValidation” package in R. Protein features extracted at least 80% of all iterations by the bootstrap EN model were then considered to be the most relevant protein biomarkers. These proteins were then fitted on a weighted logistic regression with HIV as outcome.
Visualisation of significantly enriched GO terms
Bubble plots were used to visualise the significantly enriched pathways (p<0.05 after adjustment for FDR) obtained from DAVID. The p-values in DAVID were obtained using a modified Fisher’s exact test [26]. The y-axis represents the fold enrichment which indicates the magnitude of the enrichment, as calculated in DAVID. Fold enrichment is defined as:
where m is the number of proteins significantly associated with HIV status or proteins extracted by the EN model that belong to a particular pathway, while M is the total number of proteins belonging to the same pathway. Variable n is the number of all proteins significantly associated with HIV status or extracted by the EN model and N is the total number of all proteins in the human background. Therefore, a fold enrichment of 10 indicates that 10% of the proteins significantly associated with HIV status belong to a particular pathway, and 1% of all annotated proteins in the human background belongs to the same pathway [26]. However, the proponents of this metric warns that big fold enrichments could be obtained from a small number of proteins, which could be due to small n or pathways with fewer members.
The x-axis on the hand represents the enrichment z-score for a particular pathway [27], which is calculated as follows:
where up is the total number of proteins upregulated, down is the total number of proteins downregulated, and count is the total number of proteins in the input which belongs to a particular pathway. Variables up and down were based on the weighted logistic regression for each individual protein. Hence, if 5 proteins belonging to pathway x were upregulated and 2 were downregulated, the z-score for pathway x would be: (5-2)/√7 = 1.13. A positive z-score indicates that the particular pathway is overall upregulated in HIV(+), whereas a negative z-score indicates an overall downregulation [27].
Results
Patient characteristics
Table 1 presents the baseline characteristics of the children in the randomised trial. A total of 843 complicated SAM children were recruited for the randomised trial, of which 179 (22%) patients were HIV(+). Age was higher and MUAC was lower in HIV(+) children than HIV(-) counterparts. Most HIV cases were found in Malawi. Sex and the presence of oedema were not associated with HIV status. Mortality was more than 2 times higher among in HIV(+) compared to HIV(-) (p < 0.001). Children whose HIV status were unknown had the highest mortality of 34%, which indicates the need for more active HIV screening among children with SAM.
Descriptive characteristics of the study participants
Among HIV(+), 33% were already under an anti-retroviral treatment (ART) regime: 53/179 (30%) on highly active antiretroviral therapy (HAART), and 7/179 (4%) on Nevirapine only. About half of the children (90/179) were naïve for ART whereas HIV treatment status was unknown for 16% (29/179). Mortality was not significantly different among children on HAART, ART naïve and children with unknown HIV treatment status (Supplementary Table 1).
Prior use of anti-retroviral treatment among all HIV (+) children included in the clinical trial
HIV is associated with increased inflammatory and immune response, dysregulated lipid metabolism, and increased proteolysis in children with SAM
Table 2 shows the characteristics for children included in the proteomics study. Fifty four HIV (+) children were compared to 113 HIV(-) children with SAM. In this sub-population, age was not significantly associated with HIV, and as in the full dataset. HIV(+) children also had significantly lower MUAC and higher mortality than HIV(-) children. There were no significant associations between sex and presence of oedema with HIV status.
Patient characteristics of those subjected to proteomics analysis
A total of 204 circulating proteins were annotated and compared between children with and without HIV infection. Of these, levels of 42 proteins were found to be significantly associated with HIV status in the initial univariate analysis (Figure 1A) (Supplementary table 2). Specifically, HIV(+) was associated with higher circulating levels of immunoglobulins, inflammatory proteins such as calprotectin (S100 calcium binding protein A8 and S100 calcium binding protein A9), complement proteins, and proteins related to host response to infection (i.e. lipopolysaccharide binding protein, galectin 3 binding protein and CD5 molecule-like protein). Enrichment analysis suggested that HIV(+) children have increased complement activation and immune response, and inflammatory responses than HIV(-) children. Neutrophil aggregation and chemokine production appeared to be the pathways most highly enriched in HIV(+) compared to HIV(-) SAM children. To substantiate these results, we quantified chemokine and cytokine levels in plasma. As shown, most chemokines had the tendency to be associated with HIV infection, where elevated plasma concentration of 12 were significantly associated with HIV status in SAM children (Figure 1B), namely: monocyte chemoattractant protein 1 (MCP1), macrophage inflammatory protein 1 beta (MIP1b, CCL4), granulocyte colony-stimulating factor (GCSF), interleukin 1 beta (IL1b), tumour necrosis factor alpha (TNFa), interleukins 2,5,7, 8 and 15 (IL2, 5, 7, 8, 15), interleukin 12 subunit beta (IL12p40), interferon gamma-induced protein 10 (IP-10), and interleukin-1 receptor antagonist (IL-1RA).
Association between individual proteins to HIV status
(A) Volcano plot showing several significantly different (FDR adjusted p value < 0.05) proteins and their log2 HIV(+) versus HIV(-) fold change. Red points represent those significantly higher in plasma of HIV(-), blue points significantly enriched in plasma of HIV(+) and orange points significantly higher than 1.5 folds in HIV(+) compared to HIV(-) SAM children. Vertical lines indicate significance level at p = 0.05 and 0.01; horizontal lines indicate more than 1.5 folds enrichment. (B) Log odds plots showing association of chemokine markers analysed using Luminex platform and HIV status. Points indicate log odds ratio for every log increase in plasma protein concentration; bars indicate 95% confidence interval. (C) Gene entology (GO-terms) enrichment analysis of proteins significantly higher in plasma showing biologically relevant processes heightened among HIV(+) SAM children. X-axis represents z-scores; y-axis, fold enrichment, and size of the spheres represent the number of proteins involved in the particular pathway. Only significantly enriched pathways (p < 0.05 after FDR adjustment) are plotted. See main text for explanation of the plots. Pathways enriched are identified in the table
Out of the 43 differentially expressed proteins, three proteins were found to be negatively associated with HIV status on initial univariate analysis, namely: adiponectin, kininogen-1 and peptidase inhibitor 16. Among HIV(+) children, there were no statistically significant associations with receiving HAART (n = 53) compared to ART naïve (n = 90) children (data not shown), although our study may not be powered to observe significant difference for this comparison. Furthermore, sensitivity analysis to address the possibility of HIV maternal antibodies in younger children, showed no significant interaction of age above or below 18 months and individual proteins plasma levels towards HIV status.
The weighted EN model extracted 73 circulating proteins (Figure 2A) that are associated with HIV status with AUROC = 0.97 [95% CI: 0.95 – 0.99] (Figure 2B) and misclassification error rate of 2.4%. Optimism-adjusted validated AUROC after bootstrapping was 0.90 [95% CI: 0.90 – 0.90], indicating a robust model. Pathway enrichment analysis highlighted that apart from increased immune response, HIV(+) children with SAM have increased proteolysis and lipid mobilisation, specifically increased very low-density lipoprotein assembly, indicating metabolic dysregulation related to cholesterol and triglyceride metabolism among HIV(+) patients (Figure 2D).
(A) Elastic net (EN) regularized regression lambda parameter optimization curve, optimal lambda parameter was chosen based on the highest area under the receiver operator curve (AUROC); (B) AUROC (0.97 [95% CI: 0.95 – 0.99]) of the EN model generated using the lambda parameter, alpha parameter was set to 0.75, final model extracted 34 protein features; (C) optimism-adjusted bootstrap validation of the generated EN model, validated AUROC = 0.90 [95% CI: 0.90 – 0.90] using 2000 iterations; (D) Gene entology (GO-terms) enrichment analysis of proteins extracted by the EN model. X-axis represents z-scores; y-axis, fold enrichment, and size of the spheres represent the number of proteins involved in the particular pathway. Gold circles represent pathways enriched in HIV(+) whereas blue circles are pathways more associated with HIV(-). The grey circle indicate that there are as much proteins in this pathway that are significantly upregulated and downregulated in HIV. Only significantly enriched pathways (p < 0.05 after FDR adjustment) are plotted. See main text for explanation of the plots. Pathways enriched are identified in the table. (E) Log odds ratio plot of the 3 proteins extracted after bootstrap validation with log odds on the x-axis and bars indicating 95% confidence interval obtained using weighted logistic regression with HIV as outcome variable and the 3 proteins as covariates. Weights used were obtained by inverse probability of treatment weights; (F) Predictive ability of the weighted logistic regression model using the 3 bootstrap validated proteins with HIV as outcome variable, AUROC = 0.80 [95% CI: 0.73 – 0.87]
After 2000 bootstrap iterations during bootstrap validation, 3 proteins were consistently extracted by the EN model >80% of the time (Figure 2E), namely: butyrylcholinesterase (BChE), complement C2 and zinc-α-2-glycoprotein (ZAG), indicating that these 3 proteins are likely to be the most important features associated with HIV in children with complicated SAM. Weighted logistic regression model of these 3 proteins showed good discrimination of HIV status (AUROC = 0.80 [95% CI: 0.74 – 0.87]) (Figure 2F).
Discussion
In this study, we report plasma proteomic differences associated with HIV status, which indicates that HIV imposes additional metabolic and inflammatory insults among HIV(+) children with SAM. Our results show that pathways involved in inflammatory response, complement cascade activation and lipid metabolism dysregulation are associated with HIV status. Levels of several plasma chemokines were also found to be higher in HIV(+) among children with SAM. Greater inflammatory responses in these children could be related to the higher inpatient mortality of HIV(+) compared to HIV(-) children with SAM.
An earlier metabolomics study in Uganda reported reduced serum levels of adiponectin and leptin, whereas serum triglycerides, ketones and even-chain acylcarnitines were higher in HIV(+) children with SAM indicating perturbed lipid metabolism [28]. Our current study therefore concurs with this finding, as we also found reduced plasma levels of adiponectin in HIV(+) SAM children compared to HIV(-) SAM children, along with upregulation of pathways involved in lipid transport and metabolism, specifically very low-density lipoprotein assembly.
Using optimism-adjusted bootstrap validation of the EN model, we found 3 proteins: complement c2, BChE and ZAG robustly distinguished HIV(+) from HIV(-) in children with SAM, demonstrating the ability of multivariate analysis techniques, such as EN, to uncover underlying relationships between protein markers which would be difficult to identify when analysed individually. The activation of the complement system during HIV infection has been previously discussed at length, which is associated with the increased cellular invasion of HIV in cells [29-31]. On the other hand, in a recent study in China, low circulating BChE was found to be highly associated with HIV severity, was predictive of mortality in adults, and was proposed as a plausible strategy for severity classification among adults with HIV [32]. BChE is also reported to be reduced in protein-energy malnutrition, stress and inflammation [33]. In animal studies, BChE deficiency was found to strongly affect fat metabolism and promotes hepatic lipid accumulation [34]. Serum BChE levels have been found to have a significant negative correlation with serum total cholesterol and serum low-density-lipoprotein cholesterol among people with diabetes mellitus [35]. Furthermore, reduced plasma levels of ZAG, an adipokine with protective effect against insulin resistance, was previously reported to be implicated in dyslipidaemia in HIV(+) adults under ART treatment [36]. Reduced circulating levels of ZAG has also been found among adults with clinically diagnosed metabolic syndrome, based on guidelines of the United States National Cholesterol Education Program (NCEP) Expert Panel Adult Treatment Panel (ATP) III criteria [37]. Serum ZAG levels have been reported lower among adults with impaired glucose tolerance and type 2 diabetes mellitus [38]. Taken together, our results therefore suggest that children with both HIV and SAM children manifest hallmarks of metabolic stress similar to those occurring in metabolic syndrome and other non-communicable diseases (NCD).
The results of this study is the first proteomics investigation on the interaction between HIV and SAM, which together with the previously published metabolomics study [28], provides more evidence on the increased metabolic stress and altered metabolic response among children living with both HIV and SAM. Our results also concur with previous studies that reported elevated metabolic stress among non-malnourished adults living with HIV leading to increased prevalence or risk for metabolic syndrome, cardiovascular diseases, diabetes and other non-communicable diseases [39-45].
However, despite our knowledge that HIV-infected populations have altered metabolic requirements compared to HIV-uninfected counterparts, WHO guidelines for the nutritional management for SAM are globally the same regardless of HIV status, which is summarized in Table 3 [46]. Nutritional management for in-patient children with SAM involves provision of a low-protein, low-fat milk-based food, F75, every few hours. F75 is used during clinical stabilization occurring during the first few days after admission and is not intended for weight gain. Once the children are clinically stabilized and are able to tolerate the solute load, children are transitioned to F100, a higher-calorie, high-fat milk intended to boost weight gain. Upon discharge from in-patient care, children are referred to community based nutritional therapeutic centres where they are provided with ready-to-use therapeutic food (RUTF), a peanut-based calorie-dense diet.
Nutritional management protocol for children with severe acute malnutrition [46]
Considering evidence of biomarkers of metabolic syndrome and NCD in HIV(+) children with SAM, it is of potential concern that our current treatment strategy involves a high-fat therapeutic diet. About 50% of much needed calories during the growth catch-up phase are supplied as lipids, which HIV(+) children may not be able to efficiently assimilate. Alterations in lipid metabolism in HIV(+) children with SAM may also mean that the high amounts of dietary lipids could be deposited as ectopic fat in the liver and muscle, predisposing to insulin resistance, diabetes, cardiovascular problems and other NCDs later in life. Although long-term metabolic follow-up studies could be done for HIV(+) children previously treated for either complicated and uncomplicated SAM, significant barriers are the high mortality rate in earlier studies of HIV(+) children with SAM, cost and difficulty tracing them years later. The results of this study indicate a need for clinical trials modifying the composition of F100 or RUTF to meet the specific metabolic needs of HIV(+) children with SAM. This could initially be done in relatively small groups with outcomes that include measuring metabolic stress.
Several studies on nutritional intervention strategies among HIV-infected adults have been reported. For instance, a study in the USA showed that dietary fat intake, specifically saturated fats, was significantly associated with hypertriglyceridemia among HIV-infected adults (18 – 60 years) [47]. Moreover, in a preclinical model, high saturated fat consumption was found to accelerate immunodeficiency virus disease progression in macaques, specifically increased mortality hazard and circulating levels of pro-inflammatory cytokines, especially IL8 [48], which has been previously reported to be associated with lipodystrophy among HIV patients [49]. In our study, we also found a significant association between high plasma IL8 concentration and HIV in SAM children. Hence, modifying the saturated fat composition of the milk-based F75 and F100 could potentially lower metabolic stress.
Furthermore, the European Society for Parenteral and Enteral Nutrition (ESPEN) gave a grade A recommendation for the use of medium-chain triglyceride (MCT)-based diet on HIV(+) patients with diarrhoea and severe undernutrition in its 2006 ESPEN Guidelines on Enteral Nutrition [50]. Grade A recommendations are given to strategies based on meta-analysis or at least one randomised control trial. In this case, the recommendation was based on a prospective, randomized double-blind comparative trial on 24 adult patients with HIV and diarrhoea of more than 4-wk duration, fat malabsorption, and loss of 10–20% of ideal body weight [51]. In this study, the authors found improved outcomes from diarrhoea and fat malabsorption from MCT than long-chain triglyceride-based diet among HIV(+) adults.
Lastly, the long-term metabolic effect of nutritional intervention strategies for SAM still remains unresolved. Most specifically, the potential metabolic stress associated with the rapid weight gain during the nutritional rehabilitation phase after SAM and its implications on nutritional outcomes during adulthood demands urgent research attention, especially for HIV(+) children with SAM.
Conclusion
Plasma proteomics reveals that HIV(+) children with SAM manifest hallmarks of metabolic stress similar to those observed in non-communicable diseases. This could be related to the poor nutritional recovery and high mortality of HIV(+) children with SAM despite clinical and nutritional intervention. The results of this study indicate a need for clinical trials modifying the composition of F100 or RUTF to meet the specific metabolic needs of HIV(+) children with SAM during rehabilitation phase. This could initially be done in relatively small groups with outcomes that include measuring metabolic stress.
Data Availability
Request for data access can be sent to KEMRI-Wellcome Trust via Prof James Berkley: JBerkley@kemri-wellcome.org. Data is currently not available in public repositories but the authors are currently in the process of obtaining permission to be able to make it publicly available.
Study approvals
The secondary analyses of the trial were approved by the Kenyan national ethics committee, KEMRI-SERU (KEMRI/RES/7/3/1). The trial was registered at clinicaltrials.gov (NCT02246296).
Author contributions
GBG, JMN and JAB designed the study and data analysis. JMN and BG performed the proteomics analysis. GBG performed data analysis. BW, RB, WV, IP and JAB were involved in the clinical aspects of the study. GBG wrote the initial drafts of the manuscript and all authors contributed to editing and improving the manuscript.
Acknowledgements
The original clinical trial was supported by the Thrasher Fund (grant number: 9403), awarded to JAB. GBG is a postdoctoral fellow of the Research Foundation – Flanders (FWO). The Ghent University –VLIR-UOS Global Minds Fund supported the travel of GBG to Kenya. JAB and JMN are currently supported by the Bill & Melinda Gates Foundation within the Childhood Acute Illness and Nutrition (CHAIN) Network (grant OPP1131320). JAB is currently supported by the MRC/DfID/Wellcome Trust Joint Global Health Trials scheme (grant MR/M007367/1). The funders had no role in the design, conduct, analysis, or writing of this manuscript.
References
