Plasma proteomics reveals markers of metabolic stress in HIV infected children with severe acute malnutrition

Gerard Bryan Gonzales, James M Njunge, Bonface M Gichuki, Bijun Wen, Isabel Potani, Wieger Voskuijl, Robert H J Bandsma, James A Berkley
doi: https://doi.org/10.1101/19006726

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 HIV status, although mechanisms of nutritional recovery in HIV and/or SAM are not well understood. 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 (n = 113), HIV-infected children (n = 54) had evidence (false discovery rate (FDR) 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-glycoprotein, butyrylcholinesterase, and increased levels of complement C2 resembling findings in metabolic syndrome, diabetes and other non-communicable diseases. HIV was also associated (FDR corrected p<0.05) with higher plasma levels of inflammatory chemokines. Considering evidence of biomarkers of metabolic stress, 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 suggest 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, whilst 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 clinical data and biological samples from a randomised clinical trial in Kenya and Malawi 13.

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 bias due to frequent death before testing could be undertaken or refusal of testing when a child was more severely ill.

View this table:
Table 1.

Descriptive characteristics of the study participants

Among HIV(+), 33% were already receiving 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).

HIV is associated with increased inflammatory and immune response, dysregulated lipid metabolism, and increased proteolysis in children with SAM

Among the children included in the proteomics study, 54 were HIV (+) and 113 were HIV(-) (Table 2). In this sub-population, age, sex and the presence of oedema were not significantly associated with HIV. HIV(+) children also had significantly lower MUAC and higher mortality than HIV(-) children.

View this table:
Table 2.

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 had 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).

Figure 1.
Figure 1. Univariate analysis of plasma proteome and individual plasma cytokines associated with HIV.

(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.

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), recognising our study was not powered 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, although power to detect was limited.

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.902], indicating a robust model. Pathway enrichment analysis highlighted that apart from increased immune response, HIV(+) children with SAM had 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).

Figure 2.
Figure 2. Multivariate analysis of plasma proteome associated with HIV.

(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, optimism-adjusted bootstrap validation of the generated EN model, validated AUROC = 0.90 [95% CI: 0.90 – 0.90] using 2000 iterations; (C) 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. (D) 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; (E) 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 2D), 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 2E).

Discussion

In this study, we report plasma proteomic differences associated with HIV status, suggesting 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. Circulating 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 14. 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 15–17.

On the other hand, BChE is a protein synthesized in the liver and abundant in plasma, which hydrolyses acetylcholine. Although very similar to its sister protein, acetylcholinesterase, biological functions of BChE appear to be more varied but less understood 18. 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 19. BChE is also reported to be reduced in SAM, stress and inflammation 20. In animal studies, BChE deficiency was found to strongly affect fat metabolism and promotes hepatic lipid accumulation 21. 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 22.

ZAG is a newly characterized adipokine that is involved in lipolysis, body weight regulation and may also be involved in the development of insulin resistance 23. Reduction in plasma levels of ZAG was previously reported to be implicated in dyslipidaemia in HIV(+) adults under ART treatment 23. 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 24. Serum ZAG levels have been reported lower among adults with impaired glucose tolerance and type 2 diabetes mellitus 25. Taken together, our results therefore suggest that children with both HIV and SAM manifest hallmarks of metabolic stress similar to those occurring in metabolic syndrome and other non-communicable diseases (NCD).

This study is the first proteomics investigation on the interaction between HIV and SAM. In summary, our results, which together with the previously published metabolomics study 14, strengthens 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 26–32.

Metabolic abnormalities have previously been reported to be attributed HAART use among HIV(+) patients 33. In a recent systematic review, use of 2 classes of HAART, protease inhibitors and nonnucleoside reverse transcriptase inhibitors, has been found to be associated with abnormalities in plasma lipid profiles 34. However, dysregulation in lipid metabolism has also been reported in HAART-naïve patients, which indicates that HIV infection alone cause lipid metabolism perturbations. An earlier longitudinal study of 50 men in the USA reported notable declines in serum total cholesterol after HIV infection compared to results of blood analysis from last seronegative visit. Large increases in total cholesterol and low-density lipoproteins (LDL) were detected after HAART initiation 35. However, many other studies reported increases in total cholesterol among HIV-infected patients naïve to HAART. For instance, in a study of ART-naïve HIV-infected adults in Ethiopia, malnutrition and lipid abnormalities (specifically total cholesterol) were associated with CD4+ T cell counts 36. In in vitro studies, transfection of a T-cell (RH9) with HIV led to the enhanced production of free fatty acids and LDL 37. Furthermore, monocytes isolated from HIV-infected patients both taking HAART and HAART-naïve, were found to have altered expression patters of receptors linked with lipid metabolism (i.e. FXR, PXR, PPARα, GR, RARα and RXR) compared to monocytes of HIV-uninfected controls 38. For our study however, we are unable to ascertain whether the lipid metabolism dysregulation we observed is due primarily on the viral load itself or the use of HAART due to lack of power for this sub-analysis. Majority of the participants subjected to proteomics analysis were HAART-naïve (50%), where 26% were on HAART, 6% were on Nevirapine alone and we had no data on treatment of 18% of the patients (Table 2). In all these studies cited, authors argue to need for monitoring of lipid profiles in HIV-infected populations. Hence, lipid monitoring may also inform nutritional and clinical recovery of children with SAM and HIV and could be implemented to improve clinical care for these children.

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 39. Nutritional management for in-patient children with SAM involves provision of a low-protein, low-fat milk-based food, F75, every three 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 milk/solute load, children are transitioned to F100, a higher-calorie, high-fat milk intended to boost weight gain or to Ready-to-Use Therapeutic Food (RUTF), a peanut-based calorie-dense diet. Upon discharge from in-patient care, children are referred to community based nutritional therapeutic centres where they are provided with RUTF on a 2 weekly basis.

View this table:
Table 3.

Nutritional management protocol for children with severe acute malnutrition 39

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 of F100 or RUTF modified to meet the expected 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) 40. 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 41, which has been previously reported to be associated with lipodystrophy among HIV patients 42. 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.

The European Society for Parenteral and Enteral Nutrition (ESPEN) have given 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 43. 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 44. 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.

Limitations of this study include absence of data on viral load and CD4+ counts of the patients, which could provide a deeper understanding of the results. Furthermore, in this study, we did not find association between oedematous malnutrition and HIV status, although several studies have a found higher HIV prevalence among non-oedematous children with SAM 45–47. In our study however, we found high in-patient mortality rate (16/46, 34%) among children with unknown HIV status, where 39/46 (85%) had non-oedematous SAM. Considering the high rate of mortality, these children may have been HIV(+). This highlights the need for earlier HIV screening among 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.

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 48. 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 a child’s HIV positive status was documented, HIV status was assessed by offering an antibody test at admission plus appropriate counseling. 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 49. 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 (MIP1a & 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 50. 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 using logistic regression analysis with inverse probability weighting (IPW) to correct for selection bias 51–54. Weights (w) were calculated as suggested by Samuelsen 52 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: Embedded Image 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: Embedded Image

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) 55.

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 56. 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 56,57. 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 56,57.

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 58 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. 59. 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. To test how well these proteins can discriminate HIV status, they 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 60. The y-axis represents the fold enrichment which indicates the magnitude of the enrichment, as calculated in DAVID. Fold enrichment is defined as: Embedded Image 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 60. However, the proponents of this metric warn 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 61, which is calculated as follows: Embedded Image 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 61.

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.

Funding Sources

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.

Declaration of Interests

The authors declare that we do not have any conflicts of interest.

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.

Tables

View this table:
Supplementary Table 1.

Prior use of anti-retroviral treatment among all HIV (+) children included in the clinical trial

View this table:
Supplementary table 2.

Association between individual proteins to HIV status

References

  1. ↵
    Black, R. E. et al. Maternal and child undernutrition and overweight in low-income and middle-income countries. The Lancet 382, 427–451, doi:10.1016/S0140-6736(13)60937-X (2013).
    OpenUrlCrossRefPubMed
  2. ↵
    Heikens, G. T. et al. Case management of HIV-infected severely malnourished children: challenges in the area of highest prevalence. Lancet 371, 1305–1307, doi:10.1016/s0140-6736(08)60565-6 (2008).
    OpenUrlCrossRefPubMedWeb of Science
  3. Irena, A. H., Mwambazi, M. & Mulenga, V. Diarrhea is a Major killer of Children with Severe Acute Malnutrition Admitted to Inpatient Set-up in Lusaka, Zambia. Nutr. J. 10, 110–110, doi:10.1186/1475-2891-10-110 (2011).
    OpenUrlCrossRefPubMed
  4. Maitland, K. et al. Children with Severe Malnutrition: Can Those at Highest Risk of Death Be Identified with the WHO Protocol? PLoS Med. 3, e500, doi:10.1371/journal.pmed.0030500 (2006).
    OpenUrlCrossRefPubMed
  5. ↵
    Hossain, M. et al. Efficacy of World Health Organization guideline in facility-based reduction of mortality in severely malnourished children from low and middle income countries: A systematic review and meta-analysis. J. Paediatr. Child Health 53, 474–479, doi:10.1111/jpc.13443 (2017).
    OpenUrlCrossRef
  6. ↵
    Kerac, M. et al. Follow-Up of Post-Discharge Growth and Mortality after Treatment for Severe Acute Malnutrition (FuSAM Study): A Prospective Cohort Study. PLOS ONE 9, e96030, doi:10.1371/journal.pone.0096030 (2014).
    OpenUrlCrossRefPubMed
  7. Berkley, J. A. et al. Daily co-trimoxazole prophylaxis to prevent mortality in children with complicated severe acute malnutrition: a multicentre, double-blind, randomised placebo-controlled trial. The Lancet. Global health 4, e464–473, doi:10.1016/s2214-109x(16)30096-1 (2016).
    OpenUrlCrossRef
  8. ↵
    Nemetchek, B. et al. Paediatric postdischarge mortality in developing countries: a systematic review. BMJ Open 8, e023445, doi:10.1136/bmjopen-2018-023445 (2018).
    OpenUrlAbstract/FREE Full Text
  9. ↵
    Fergusson, P. & Tomkins, A. HIV prevalence and mortality among children undergoing treatment for severe acute malnutrition in sub-Saharan Africa: a systematic review and meta-analysis. Trans. R. Soc. Trop. Med. Hyg. 103, 541–548, doi:10.1016/j.trstmh.2008.10.029 (2009).
    OpenUrlCrossRefPubMed
  10. ↵
    Ndirangu, M., Wariero, J. O., Sachs, S. E., Masibo, P. & Deckelbaum, R. J. Nutritional status of under-five children in HIV-affected households in western Kenya. Food Nutr. Bull. 32, 159–167, doi:10.1177/156482651103200208 (2011).
    OpenUrlCrossRefPubMed
  11. ↵
    Rose, A. M., Hall, C. S. & Martinez-Alier, N. Aetiology and management of malnutrition in HIV-positive children. Arch. Dis. Child. 99, 546–551, doi:10.1136/archdischild-2012-303348 (2014).
    OpenUrlAbstract/FREE Full Text
  12. ↵
    Trehan, I., O’Hare, B. A., Phiri, A. & Heikens, G. T. Challenges in the Management of HIV-Infected Malnourished Children in Sub-Saharan Africa. AIDS Res. Treat. 2012, 8, doi:10.1155/2012/790786 (2012).
    OpenUrlCrossRef
  13. ↵
    Bandsma, R. H. J. et al. A reduced-carbohydrate and lactose-free formulation for stabilization among hospitalized children with severe acute malnutrition: A double-blind, randomized controlled trial. PLoS Med. 16, e1002747, doi:10.1371/journal.pmed.1002747 (2019).
    OpenUrlCrossRef
  14. ↵
    Mody, A. et al. Effects of HIV infection on the metabolic and hormonal status of children with severe acute malnutrition. PLoS One 9, e102233, doi:10.1371/journal.pone.0102233 (2014).
    OpenUrlCrossRefPubMed
  15. ↵
    Tacnet-Delorme, P. et al. In Vitro Analysis of Complement-Dependent HIV-1 Cell Infection Using a Model System. The Journal of Immunology 162, 4088 (1999).
    OpenUrl
  16. Stoiber, H., Kacani, L., Speth, C., Würzner, R. & Dierich, M. P. The supportive role of complement in HIV pathogenesis. Immunol. Rev. 180, 168–176, doi:10.1034/j.1600-065X.2001.1800115.x (2001).
    OpenUrlCrossRefPubMedWeb of Science
  17. ↵
    Stoiber, H., Speth, C. & Dierich, M. P. Role of complement in the control of HIV dynamics and pathogenesis. Vaccine 21, S77–S82, doi: https://doi.org/10.1016/S0264-410X(03)00203-2 (2003).
    OpenUrlCrossRefPubMed
  18. ↵
    Oliveira, J. d. et al. Effect of BCHE single nucleotide polymorphisms on lipid metabolism markers in women. Genet. Mol. Biol. 40, 408–414, doi:10.1590/1678-4685-GMB-2016-0123 (2017).
    OpenUrlCrossRef
  19. ↵
    Xu, L. et al. Butyrylcholinesterase Levels on Admission Predict Severity and 12-Month Mortality in Hospitalized AIDS Patients. Mediators Inflamm. 2018, 5201652, doi:10.1155/2018/5201652 (2018).
    OpenUrlCrossRef
  20. ↵
    Santarpia, L., Grandone, I., Contaldo, F. & Pasanisi, F. Butyrylcholinesterase as a prognostic marker: a review of the literature. Journal of cachexia, sarcopenia and muscle 4, 31–39, doi:10.1007/s13539-012-0083-5 (2013).
    OpenUrlCrossRef
  21. ↵
    Chen, V. P. et al. Butyrylcholinesterase Deficiency Promotes Adipose Tissue Growth and Hepatic Lipid Accumulation in Male Mice on High-Fat Diet. Endocrinology 157, 3086–3095, doi:10.1210/en.2016-1166 (2016).
    OpenUrlCrossRef
  22. ↵
    Sridhar, G. et al. Serum butyrylcholinesterase in type 2 diabetes mellitus: a biochemical and bioinformatics approach. Lipids Health Dis. 4, 18, doi:10.1186/1476-511x-4-18 (2005).
    OpenUrlCrossRefPubMed
  23. ↵
    Ceperuelo-Mallafre, V. et al. Zinc alpha-2 glycoprotein is implicated in dyslipidaemia in HIV-1-infected patients treated with antiretroviral drugs. HIV Med. 13, 297–303, doi:10.1111/j.1468-1293.2011.00976.x (2012).
    OpenUrlCrossRefPubMedWeb of Science
  24. ↵
    Lei, L. et al. Circulating zinc-α2-glycoprotein levels are low in newly diagnosed patients with metabolic syndrome and correlate with adiponectin. Nutr. Metab. (Lond.) 14, 53, doi:10.1186/s12986-017-0210-6 (2017).
    OpenUrlCrossRef
  25. ↵
    Yang, M. et al. Zinc-α2-glycoprotein is associated with insulin resistance in humans and is regulated by hyperglycemia, hyperinsulinemia, or liraglutide administration: cross-sectional and interventional studies in normal subjects, insulin-resistant subjects, and subjects with newly diagnosed diabetes. Diabetes care 36, 1074–1082, doi:10.2337/dc12-0940 (2013).
    OpenUrlAbstract/FREE Full Text
  26. ↵
    Bastard, J. P. et al. Diabetes and dyslipidaemia are associated with oxidative stress independently of inflammation in long-term antiretroviral-treated HIV-infected patients. Diabetes Metab., doi:10.1016/j.diabet.2019.02.008 (2019).
    OpenUrlCrossRef
  27. Nix, L. M. & Tien, P. C. Metabolic syndrome, diabetes, and cardiovascular risk in HIV. Curr. HIV/AIDS Rep. 11, 271–278, doi:10.1007/s11904-014-0219-7 (2014).
    OpenUrlCrossRef
  28. Sobieszczyk, M. E. et al. Metabolic Syndrome After HIV Acquisition in South African Women. J. Acquir. Immune Defic. Syndr. 73, 438–445, doi:10.1097/qai.0000000000001123 (2016).
    OpenUrlCrossRef
  29. Policarpo, S., Rodrigues, T., Moreira, A. C. & Valadas, E. Cardiovascular risk in HIV-infected individuals: A comparison of three risk prediction algorithms. Rev. Port. Cardiol. 38, 463–470, doi:10.1016/j.repc.2019.08.002 (2019).
    OpenUrlCrossRef
  30. Lagathu, C. et al. Metabolic complications affecting adipose tissue, lipid and glucose metabolism associated with HIV antiretroviral treatment. Expert Opin. Drug Saf. 18, 829–840, doi:10.1080/14740338.2019.1644317 (2019).
    OpenUrlCrossRef
  31. Sears, S. et al. Metabolic Syndrome Among People Living with HIV Receiving Medical Care in Southern United States: Prevalence and Risk Factors. AIDS Behav., doi:10.1007/s10461-019-02487-8 (2019).
    OpenUrlCrossRef
  32. ↵
    Seth, A. & Sherman, K. E. Fatty liver disease in persons with HIV infection. Top. Antivir. Med. 27, 75–82 (2019).
    OpenUrlPubMed
  33. ↵
    Palios, J., Kadoglou, N. P. E. & Lampropoulos, S. The pathophysiology of HIV-/HAART-related metabolic syndrome leading to cardiovascular disorders: the emerging role of adipokines. Exp. Diabetes Res. 2012, 103063–103063, doi:10.1155/2012/103063 (2012).
    OpenUrlCrossRefPubMed
  34. ↵
    Silva, B., Peixoto, G., da Luz, S., de Moraes, S. & Peres, S. Adverse effects of chronic treatment with the Main subclasses of highly active antiretroviral therapy: a systematic review. HIV Med. 20, 429–438, doi:10.1111/hiv.12733 (2019).
    OpenUrlCrossRef
  35. ↵
    Riddler, S. A. et al. Impact of HIV Infection and HAART on Serum Lipids in Men. JAMA 289, 2978–2982, doi:10.1001/jama.289.22.2978 (2003).
    OpenUrlCrossRefPubMedWeb of Science
  36. ↵
    Adal, M., Howe, R., Kassa, D., Aseffa, A. & Petros, B. Malnutrition and lipid abnormalities in antiretroviral naïve HIV-infected adults in Addis Ababa: A cross-sectional study. PLOS ONE 13, e0195942, doi:10.1371/journal.pone.0195942 (2018).
    OpenUrlCrossRef
  37. ↵
    Rasheed, S., Yan, J. S., Lau, A. & Chan, A. S. HIV Replication Enhances Production of Free Fatty Acids, Low Density Lipoproteins and Many Key Proteins Involved in Lipid Metabolism: A Proteomics Study. PLOS ONE 3, e3003, doi:10.1371/journal.pone.0003003 (2008).
    OpenUrlCrossRefPubMed
  38. ↵
    Renga, B. et al. HIV-1 infection is associated with changes in nuclear receptor transcriptome, pro-inflammatory and lipid profile of monocytes. BMC Infect. Dis. 12, 274, doi:10.1186/1471-2334-12-274 (2012).
    OpenUrlCrossRefPubMed
  39. ↵
    Trehan, I. & Manary, M. J. Management of severe acute malnutrition in low-income and middle-income countries. Arch. Dis. Child. 100, 283–287, doi:10.1136/archdischild-2014-306026 (2015).
    OpenUrlAbstract/FREE Full Text
  40. ↵
    Joy, T. et al. Dietary fat intake and relationship to serum lipid levels in HIV-infected patients with metabolic abnormalities in the HAART era. AIDS (London, England) 21, 1591–1600, doi:10.1097/QAD.0b013e32823644ff (2007).
    OpenUrlCrossRefPubMedWeb of Science
  41. ↵
    Mansfield, K. G. et al. A diet high in saturated fat and cholesterol accelerates simian immunodeficiency virus disease progression. J. Infect. Dis. 196, 1202–1210, doi:10.1086/521680 (2007).
    OpenUrlCrossRefPubMedWeb of Science
  42. ↵
    Lindegaard, B. et al. Adipose tissue expression of IL-18 and HIV-associated lipodystrophy. AIDS 18, 1956–1958, doi:10.1097/00002030-200409240-00013 (2004).
    OpenUrlCrossRefPubMedWeb of Science
  43. ↵
    Ockenga, J. et al. ESPEN Guidelines on Enteral Nutrition: Wasting in HIV and other chronic infectious diseases. Clin. Nutr. 25, 319–329, doi:10.1016/j.clnu.2006.01.016 (2006).
    OpenUrlCrossRefPubMed
  44. ↵
    Wanke, C. A. et al. A medium chain triglyceride-based diet in patients with HIV and chronic diarrhea reduces diarrhea and malabsorption: A prospective, controlled trial. Nutrition 12, 766–771, doi: https://doi.org/10.1016/S0899-9007(96)00221-3 (1996).
    OpenUrlCrossRefPubMed
  45. ↵
    Amadi, B. et al. Intestinal and systemic infection, HIV, and mortality in Zambian children with persistent diarrhea and malnutrition. J. Pediatr. Gastroenterol. Nutr. 32, 550–554, doi:10.1097/00005176-200105000-00011 (2001).
    OpenUrlCrossRefPubMedWeb of Science
  46. Kessler, L., Daley, H., Malenga, G. & Graham, S. The impact of HIV infection on the clinical presentation of severe malnutrition in children at QECH. Malawi medical journal : the journal of Medical Association of Malawi 13, 30–33 (2001).
    OpenUrl
  47. ↵
    Musoke, P. M. & Fergusson, P. Severe malnutrition and metabolic complications of HIV-infected children in the antiretroviral era: clinical care and management in resource-limited settings. The American journal of clinical nutrition 94, 1716S–1720S, doi:10.3945/ajcn.111.018374 (2011).
    OpenUrlAbstract/FREE Full Text
  48. ↵
    WHO. Management of severe malnutrition: a manual for physicians and other senior health workers, <http://apps.who.int/iris/bitstream/10665/41999/1/a57361.pdf> (1999).
  49. ↵
    Njunge, J. M. et al. Biomarkers of post-discharge mortality among children with complicated severe acute malnutrition. Sci. Rep. 9, 5981, doi:10.1038/s41598-019-42436-y (2019).
    OpenUrlCrossRef
  50. ↵
    R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, <https://www.R-project.org/> (2017).
  51. ↵
    Kim, R. S. & Kaplan, R. C. Analysis of secondary outcomes in nested case-control study designs. Stat. Med. 33, 4215–4226, doi:10.1002/sim.6231 (2014).
    OpenUrlCrossRefPubMed
  52. ↵
    Samuelsen, S. O. A Pseudolikelihood Approach to Analysis of Nested Case-Control Studies. Biometrika 84, 379–394 (1997).
    OpenUrlCrossRefWeb of Science
  53. Pan, Y., Cai, J., Longnecker, M. P. & Zhou, H. Secondary outcome analysis for data from an outcome-dependent sampling design. Stat. Med. 37, 2321–2337, doi:10.1002/sim.7672 (2018).
    OpenUrlCrossRef
  54. ↵
    Jiang, Y., Scott, A. J. & Wild, C. J. Secondary analysis of case-control data. Stat. Med. 25, 1323–1339, doi:10.1002/sim.2283 (2006).
    OpenUrlCrossRefPubMedWeb of Science
  55. ↵
    Benjamini, Y. & Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society. Series B (Methodological) 57, 289–300 (1995).
    OpenUrlCrossRefPubMedWeb of Science
  56. ↵
    Friedman, J. H., Hastie, T. & Tibshirani, R. Regularization Paths for Generalized Linear Models via Coordinate Descent. 2010 33, 22, doi:10.18637/jss.v033.i01 (2010).
  57. ↵
    Zou, H. & Hastie, T. Regularization and variable selection via the elastic net. J. Roy. Stat. Soc. Ser. B. (Stat. Method.) 67, 301–320, doi:10.1111/j.1467-9868.2005.00503.x (2005).
    OpenUrlCrossRefWeb of Science
  58. ↵
    Huang, D. W., Sherman, B. T. & Lempicki, R. A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protocols 4, 44–57, doi: http://www.nature.com/nprot/journal/v4/n1/suppinfo/nprot.2008.211_S1.html (2008).
    OpenUrl
  59. ↵
    Smith, G. C. S., Seaman, S. R., Wood, A. M., Royston, P. & White, I. R. Correcting for Optimistic Prediction in Small Data Sets. Am. J. Epidemiol. 180, 318–324, doi:10.1093/aje/kwu140 (2014).
    OpenUrlCrossRefPubMedWeb of Science
  60. ↵
    Huang da, W., Sherman, B. T. & Lempicki, R. A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57, doi:10.1038/nprot.2008.211 (2009).
    OpenUrlCrossRefPubMedWeb of Science
  61. ↵
    Walter, W., Sanchez-Cabo, F. & Ricote, M. GOplot: an R package for visually combining expression data with functional analysis. Bioinformatics 31, 2912–2914, doi:10.1093/bioinformatics/btv300 (2015).
    OpenUrlCrossRefPubMed
Back to top
Posted January 27, 2020.
Download PDF
Data/Code
Email
Plasma proteomics reveals markers of metabolic stress in HIV infected children with severe acute malnutrition
Gerard Bryan Gonzales, James M Njunge, Bonface M Gichuki, Bijun Wen, Isabel Potani, Wieger Voskuijl, Robert H J Bandsma, James A Berkley
medRxiv 19006726; doi: https://doi.org/10.1101/19006726
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area