Skip to main content

Long-term survival outcomes of HIV infected children receiving antiretroviral therapy: an observational study from Zambia (2003–2015)



In 2017, 64% of children living with HIV in Zambia accessed Antiretroviral Therapy (ART). Despite expanded ART coverage, there is paucity of information on effectiveness of pediatric ART in reducing mortality. The aim of this research is to describe treatment outcomes, measure mortality rates and assess predictors of mortality among children receiving ART.


Using a retrospective cohort study design, we abstracted routinely collected clinical data from medical records of children from birth to 15 years old, who had received ART for at least 6 months at Livingstone Central Hospital in Southern Province Zambia, between January 2003 and June 2015. The primary outcome was death. Cause of death was ascertained from medical records and death certificates. Distribution of survival times according to baseline covariates were estimated using Kaplan Meier and Cox Proportional Hazards methods.


Overall, 1039 children were commenced on ART during the study period. The median age at treatment initiation was 3.6 years (IQR: 1.3–8.6) and 520 (50%) children were female. Of these, 71 (7%) died, 164 (16%) were lost to follow-up, 210 (20%) transferred and 594 (56%) were actively on treatment. After 4450 person years, mortality rate was 1.6/100 (95% CI: 1.4–1.8). Mortality was highest during the first 3 months of treatment (11.7/100 (95% CI: 7.6–16.3). In multivariable proportional hazards regression, the adjusted hazards of death were highest among children aged < 1 year (aHR = 3.1 (95% CI: 1.3–6.4), compared to those aged 6–15 years, WHO stage 4 (aHR =4.8 (95% CI: 2.3–10), compared to WHO stage 1 and 2. In the sensitivity analysis to address bias due to loss to follow-up, mortality increased 5 times when we assumed that all the children who were lost to follow up died within 90 days of their last visit.


We observed low attrition due to mortality among children on ART. Loss to follow-up was high (16%). Mortality was highest during the first 3 months of treatment. Children aged less than one year and those with advanced WHO disease stage had higher mortality. We recommend effective interventions to improve retention in care and early diagnosis of HIV in children.

Peer Review reports


The availability of Antiretroviral Therapy (ART) for children living with HIV and implementation of universal treatment of all pregnant and breastfeeding women living with HIV (Option B+) is a game changer in the global fight against HIV [1]. In 2017, about 1.8 million children were living with HIV globally and 90% of these children lived in sub-Saharan Africa [2, 3]. Although progress has been reported in the scale-up of access to treatment for children, only 52% of children living with HIV received lifesaving ART and only 51% of HIV-exposed infants were tested for HIV by the age of 2 months as recommended by the World Health Organization (WHO) guidelines [2]. These estimates fall short of the UNAIDS 90–90-90 treatment targets, a strategy to end the global HIV epidemic. This strategy aims to achieve the following by 2020: 1) 90% of people living with HIV will know their status, 2) 90% of all diagnosed people will be on ART, 3) 90% of people on ART will be virally suppressed [4]. Although some milestones have been achieved in the provision of ART, access to early HIV diagnosis and ART among infants and children remains a challenge in high HIV burden settings [2, 5].

Similarly, the pediatric HIV program in Zambia has made tremendous progress with over 64% of children living with HIV accessing ART by the end of 2017 [6]. With a population of 17 million people and an estimated HIV prevalence of 12.9%, 72,000 were children living with HIV and 8900 were newly infected in 2017 [6, 7]. This scenario creates a large pool of children living with HIV in need of treatment. An early study done in routine care settings in Zambia demonstrated that children were diagnosed at older ages with advanced WHO stage 3 or 4 disease. The same study reported that 57% of deaths occurred within the first 90 days of treatment initiation and loss to follow-up was high [8]. Early mortality was consistently associated with lower CD4 count, younger age, low weight for height and anemia at ART initiation [8,9,10]. Recent studies suggest that routine care settings in resource limited settings are still confronted with critical gaps in early infant diagnosis and pediatric HIV treatment, resulting in children initiating treatment with advanced disease and at high risk of death during the first 3–6 months of treatment [11, 12]. In addition, there are concerns around interpretation of treatment outcomes in routine care settings because of the high loss to follow-up and high numbers of transfers with no mortality ascertainment [11, 13].

Although, pediatric HIV treatment programs in resource limited settings like Zambia have been treating large numbers of HIV-exposed and infected children for over 15 years, long-term survival and associated factors including barriers to long-term adherence and retention in care are not adequately elucidated and therefore the wider impact of these interventions is not clearly understood [14,15,16]. The aim of this research is to describe treatment outcomes, measure mortality rates and assess predictors of mortality among children receiving ART in a routine care setting over a 12-year period (2003–2015) at Livingstone Central Hospital in Southern Province, Zambia.


Using a retrospective cohort study design, we evaluated treatment outcomes among children living with HIV who received ART between January 2003 and June 2015.

Study setting and population

The data were collected from the Pediatric Center of Excellence clinic (PCOE) at Livingstone Central Hospital (LCH) in Southern Province, Zambia. Livingstone Central Hospital offers preventive and treatment services to nearly 1.2 million people in southern and western parts of Zambia. Pediatric HIV treatment at LCH was started in 2003 in line with national policy and in 2006 the PCOE clinic was established through a collaborative agreement between the Ministry of Health in Zambia and the Centers for Diseases Control and prevention country office (CDC) [15].

The Ministry of Health in Zambia began implementation of the WHO guidelines of treating all people who test positive for HIV in 2016. The pediatric HIV testing and treatment guidelines recommended collection of dried blood spots (DBS) cards for DNA-PCR from all HIV-exposed infants at birth, 6 weeks and 6 months of age followed by serologic testing at 9 months, 12 months, 18 months and 3 months after cessation of breastfeeding [17]. Dried blood spots cards were collected from all children who were still breastfeeding and had HIV positive serologic tests, for confirmatory testing with HIV DNA-PCR. Routine provider initiated counseling and testing was offered to guardians of all hospitalized children. Prior to 2016, pediatric HIV treatment guidelines recommended ART based on clinical and immunological criteria. First line ART regimens comprised of three antiretroviral medicines including two Nucleoside reverse transcriptase inhibitors (NRTI’s): Azido thymidine, Stavudine, Abacavir, (Tenofovir for children above 10 years old) and Lamivudine, plus one Non-Nucleoside Reverse Transcriptase Inhibitor (NNRTI), either Nevirapine or Efavirenz for ART naïve infants and children. Infants who took ART prophylaxis after delivery were commenced on Protease inhibitors (PI) based regimens (Lopinavir boosted with Ritonavir) in addition to two NRTIs. Triple NRTI based regimens were recommended for children below 3 months of age who were co-infected with Tuberculosis at baseline.

Patient data were abstracted from the medical records of children living with HIV who had received ART. The inclusion criteria were: all children aged from birth to 15 years at the time of ART initiation, children who had taken ART for 6 months or more, children who had received at least three ARVs. Children who received one antiretroviral drug for PMTCT were excluded from the study because we considered ART to be a combination of at least three antiretroviral drugs. We excluded children who started ART on or after December 31st 2014 because data collection was closed on June 30th 2015 and we included only children who took ART for at least 6 months.

Data collection

We abstracted routinely collected data from patient medical records into a Microsoft access database. We abstracted the data starting with initial patient enrollment forms which included demographics such as date of birth, gender, date of first HIV test, mode of delivery and date when they started to take ART. We assigned unique study identification numbers to each patient record. Data abstraction was done independently by two data entry clerks and their entries were verified at random intervals by a supervisor.

Study outcomes

The main outcome was death. We ascertained death from hospital medical records. We actively searched for death certificates and other documentation on the details of death in medical records. In some cases, details of death were verified by verbal reports from the child’s adult family members or care givers.

The secondary outcomes were: 1) Lost to follow-up, and 2) Transferred to another facility. Lost to follow-up was defined as no clinical contact for more than 3 months after they missed their last scheduled appointment. Patients who missed appointments for more than 90 days were visited by a clinic nurse and their status established. All the children who were lost to follow-up or transferred were censored on the day of their last visit to the clinic. Children who survived to the end of the study were censored on June 30th, 2015. Adolescents were censored on the date they transitioned into the adolescent clinic and were no longer seen in the pediatric clinic at 15 years of age.


The main exposure was initiation of ART. Since this was a retrospective cohort study, the decision to initiate ART was made by the treating physician and the child’s caregivers and not influenced by the study investigators. Covariates that were abstracted included: ART start date, the initial drug regimen and baseline laboratory measurements which included hemoglobin, CD4+ count or CD4+ percentage. Additionally, WHO clinical stage and anthropometric measurements were abstracted.

Statistical analysis

Baseline demographics and clinical features were described by estimating medians and interquartile ranges for continuous variables and frequencies and proportions for categorical variables. Observed survival was defined as the duration of time in years from the date of ART initiation to the date of death or censoring. We estimated the distribution of survival times for the baseline covariates affecting survival using the Kaplan-Meier (KM) method. The covariates were selected apriori from literature supported by a directed acyclic model (Additional file 1). Log-rank tests were used to compare the survival curves. Cox Proportional Hazards models were used to estimate the hazard ratios for death and their 95% confidence intervals to quantify the associations. We evaluated the proportional hazards assumption using log-log plots and plots of Schoenfeld’s residuals and no violations of the assumption were found. To calculate the mortality rates, we estimated the actual person time that everyone contributed to the study during the first 6 months, 1 year, 2 years, 5 years and 10 years and constructed confidence intervals.

Multivariable Cox models were used to adjust for possible confounding. All p-values were two tailed. In the analysis, age was a categorical variable with 3 levels. The age categorization was based on findings from previous studies which demonstrated that infants aged less than 12 months had a different survival experience based on disease progression [5, 18]. WHO stage was a categorical variable with 4 levels in the initial analysis. In the Cox Proportional Hazards models, we combined WHO stage 1 and 2, therefore WHO stage was a categorical variable with 3 levels.

The status of children who were lost to follow-up in this study was unknown. A meta-analysis of individual patient data from HIV treatment programs in Sub-Sahara Africa found that 20% of adults lost to follow-up were actually deceased [19]. To evaluate the “worst case scenario” effect of bias due to loss to follow-up on mortality, we performed sensitivity analysis by assuming that all the children who were lost to follow-up survived for 90 days after their last visit and then died. We estimated the mortality rates for this scenario at 3 months, 6 months, 1, 5 and 10 years and compared with the main study results.

Data analysis was done using SAS version 9.4 (SAS Institute Inc., Cary, NC) and R statistical software version 3.4 [20]. We used the Survival package to plot the KM curves and Survminer package to visualize the KM curves [21, 22].


One thousand and thirty-nine (1039) children aged less than 15 years commenced ART at LCH between January 2003 and June 2015. Overall, 71 (7%) children were confirmed to have died after commencing treatment and 594 (56%) were alive and active in care at the time of this study. A total of 164 (16%) were lost to follow-up and 210 (20%) transferred to other health care facilities (Fig. 1).

Fig. 1
figure 1

Treatment outcomes and retention in care among children on ART at LCH 2003–2015

Baseline characteristics of children started on ART

At baseline, 520 (50%) children were female and 721 (69%) were cared for by their biological mothers. A total of 179 (17%) children started treatment during their first year of life. At least 304 (29%) children were diagnosed during hospital admission and 30 (3%) were diagnosed after delivery from the delivery wards (Table 1).

Table 1 Baseline Characteristics of infants and Children who Received ART by Age Group - Livingstone Central Hospital, Zambia (2003–2015)

The median age at baseline was 3.6 years (IQR: 1.3–8.6). The immunological criteria that was used to initiate treatment for children older than 5 years was CD4 count and the baseline median CD4 count was 505 (IQR: 243–948) while CD4 percentage was used for children younger than 5 years. Children aged less than 1 year had median CD4 percentage of 19.4 (IQR: 12.4–25.6) and children aged between 1 and 5 years had median CD4 percentage of 16.7 (IQR: 11.3–21.5). At least 1002 (97%) took Cotrimoxazole at baseline and 301 (30%) had a diagnosis of clinical Tuberculosis (TB) at ART initiation. Overall, 472 (46%) were WHO stage 3 and 177 (17%) were advanced WHO stage 4 (Table 2).

Table 2 Baseline laboratory and clinical features of infants and children who received ART at Livingstone Central Hospital, Zambia (2003–2015)

The first line regimen for 907 (87%) children was two NRTIs and one NNRTI, and 54 (5%) were commenced on two NRTIs and a Protease Inhibitor (Lopinavir boosted with Ritonavir) while 78 (8%) children took triple NRTIs at baseline, (Azido thymidine or Stavudine, Lamivudine and Abacavir) (Table 2).

Mortality trends

A total of 71 deaths were documented during the study period. The highest proportion of deaths occurred between 2007 and 2011 (Additional file 2). The largest number of deaths were among young children aged less than 12 months (n = 29; 41%) (Fig. 2). Of the 71 deceased children, 44% (n = 31) were WHO stage 4, 41% (n = 29) were WHO stage 3, 8% (n = 6) were WHO stage 2 and 7% (n = 5) were WHO stage 1. We verified causes of death from 27 (38%) death certificates. Of the verified deaths 14% (n = 10) were due to diarrhea, 11% (n = 8) severe pneumonia and protein energy malnutrition respectively and 7% (n = 5) were due to TB (Table 3).

Fig. 2
figure 2

Age at time of Death Among 71 Deceased Children on ART: Livingstone Central Hospital, Zambia (2003–2015)

Table 3 Cause of death as recorded on death certificate

For an overall person-year of observation (PYO) of 4450 PY the mortality rate was 1.6/100 PYO (95% CI:1.4–1.8). Mortality was highest within the first 3 months of treatment initiation with estimated mortality rate of 11.9/100 PYO (95% CI:7.6–16.3), which accounted for 41% (n = 29) of the deaths (Table 4). The 3-month survival probability was 0.97 (Fig. 3). After 10 years of treatment, the observed mortality rate was 1.6/100 PYO (95% CI: (1.4–1.8).

Table 4 Mortality among children receiving ART at Livingstone Central Hospital, (2003–2015)
Fig. 3
figure 3

KM Survival plot of time to death after ART initiation: Livingstone Central Hospital, Zambia (2003–2015)

In a sensitivity analyses that assumed the “worst case scenario” that all children who were lost to follow-up died 90 days after their last clinic visit, the mortality rate after 3 months of treatment was 14.5/100 person years (95% CI: 9.7–19.2). After 10 years of treatment, the mortality rate was 5.4/100 person years (95% CI: 4.7–6.1) which is 5 times higher than actual observed mortality rate (Table 5).

Table 5 Sensitivity analysis - effect of bias due to loss to follow-up – worst case scenario

Children who took triple NRTI regimens had the shortest survival time (86% at 1 year) compared to children who took 2NRTI+ 1NNRTI and protease inhibitor based regimens (log rank test, p < 0.0001, Fig. 4 d). Other baseline factors associated with shorter survival were: children who had anemia at ART initiation (Hemoglobin < 8 g/dl) (p = 0.032), children of mothers who did not take antiretroviral drugs during pregnancy (p = 0.012), severe wasting as estimated by WAZ score < -3SD (p = 0.051) and infants who were given ARVs for prophylaxis (p = 0.056) (Additional file 3). Tuberculosis at baseline, calendar year of ART initiation and whether the infant took antiretroviral drugs for prophylaxis at enrollment were not statistically significant in the univariable Cox proportional hazards model and we did not include them in the multivariable model (Table 6).

Fig. 4
figure 4

a-4d: 4A-Age at enrollment, 4b, WHO stage at enrollment, 4c, HIV disclosure at enrollment, 4d: ART regimen at Baseline

Table 6 Factors Associated with Mortality Among Children on ART at Livingstone Central Hospital (2003–2015)

Infants who started ART within the first year of life experienced shorter survival time compared to those who started treatment between 1 to 5 years and between 6 to 15 years of age (Fig. 4 a, log-rank test, p-value < 0.0001). The survival probability for infants aged less than 1 year was 0.93 after 3 months of treatment compared to 0.96 for those aged 1 to 5 years and 0.99 for those aged 6–15 years. In the multivariable Proportional Hazards regression, the adjusted hazards of death among children aged less than 1 year were 3.1 (95% CI: 1.3–6.4), and children aged 1–5 years aHR = 1.1 (95% CI: 0.6–2.3) when compared to the referent group 6–15 years (Table 6).

Children with severe immune-suppression (WHO stage 4) at baseline had lower survival probability of 93, 97% for the WHO stage 3 and 99% for WHO stage 1 and 2, (Fig. 4 b, log-rank test, p-value < 0.0001). In the multivariable Proportional Hazards regression, children with advanced WHO stage 4 had the highest hazards of mortality (aHR = 4.8 (95% CI: 2.3–10) (p < .0001), WHO stage 3 (aHR = 1.8 (0.9–3.6) p = 0.1188), compared to the referent group WHO stage 1 and 2, (Table 6).


In this study carried out in a routine pediatric HIV treatment setting, we found that after 12 years of treatment (2003–2015) attrition due to death was low (7%; n = 71). Death was highest within the first 3 months of starting ART (mortality rate: 11.7/100 PYO; 95% CI:7.6–16.3). The hazard of death was highest among children aged less than 12 months (aHR: 3.1: 95% CI: 1.3–6.4), compared to the referent group 6–15 years. Loss to follow-up was high (16%; n = 164).

The mortality rate observed in this study was much lower than that observed in other parts of Zambia [8, 23]. In a study conducted in an urban area in Zambia, the mortality rate was 6.6/100 PYO over 3018 PYO [8] and another study in a rural area in Zambia found a high mortality of 14.4% after 6 months of treatment [23]. These 2 studies were characterized by high mortality within the first 3 months (17.4/100 PYO) and the associated risk factors were similar to what we observed [8]. Mortality rates in our study population were also similar to a study from a rural setting in Malawi where 12 months mortality rate was 6.6/100 PYO [95% CI: 5.5–7.9] and overall mortality was 3.4/100 PYO [95% CI:2.9–4] after 5 years of follow-up [24]. Our observed mortality rate is however higher than findings from developed countries where 10-year mortality was as low as 0.3/100 PYO in the Dutch cohort [25].

The observation that hazards of death were highest within the first few months after ART initiation and then declines after the first 6 months of treatment is consistent with studies done in other parts of Zambia and in similar settings [8, 26,27,28]. Other studies have attributed the high early mortality to late presentation to care [29]. Studies done in routine clinic settings in resource limited settings found that the cumulative incidence of mortality during the first year of treatment among older children between 5 and 10 years of age was less than 2% to more than 45% among infants aged less than 12 months with severe disease [10]. Once children get past the initial 6 months of ART, their risk of mortality declines to very low rates especially when adherence to treatment and follow-up is optimized [11, 27, 30].

Our finding that the hazards of mortality were highest among infants and children aged less than 12 months of age is consistent with predictive models that have been done in both developed and developing countries and in fact motivated the universal ART policy by WHO to test and treat all HIV-exposed infants and children [31]. Results of a modelling study done in six countries in Sub-Sahara Africa suggested that mortality is higher among perinatally HIV-infected children than those infected through breastfeeding [32]. Perinatally HIV-infected infants have a higher risk of mortality and disease progression [5]. Early diagnosis of HIV infection and early initiation of ART improves the outcomes of infants and children [33]. Although early infant diagnosis of HIV is critical, there are still challenges in resource limited settings to diagnose children. Children who are missed by the PMTCT programs are at highest risk of mortality because they usually present to the hospital after an illness when they are already immunocompromised with high risk of mortality [34].

We assessed the causes of death in our study but only found death certificates for 27 children and no postmortem or thorough investigations had been done to ascertain cause of death. However, causes of death recorded on the few death certificates were similar to the causes of childhood deaths in Zambia, which are diarrhea, severe pneumonia and protein energy malnutrition [35]. This setting has a very high infant mortality rate (65/1000 in 2015) [36]. In the pre-ART era, AIDS related deaths led to an increase in infant mortality in high HIV burden countries. Antiretroviral therapy improved child survival and evidence is now showing that early initiation of ART leads to better clinical outcomes [33]. In a study that pooled results of clinical trial data from Zimbabwe and Uganda, the mortality risk was attributed to pre-ART risks that persist until the antiretroviral drugs reach their maximal effectiveness [27]. These early deaths have been attributed to suboptimal management of malnutrition, Tuberculosis and other related medical conditions during the early treatment stage including Immune Reconstitution Syndrome (IRIS) [11]. However, the role of IRIS in mortality during the early months of treating HIV seropositive children with ART is poorly understood [37]. Our study had insufficient evidence to support the role of IRIS as a risk factor of early mortality.

The major strength of our study is that it was done in a routine clinic setting among children commencing ART in a high HIV burden setting. This provides a real-world effectiveness of pediatric HIV treatment outcomes and provides information on important predictors of mortality. Our sample size was large and we had enough outcomes to make reasonable conclusions.

A major limitation is that there was a high proportion of loss to follow-up (16%) and many transfers (20%). This made ascertainment of mortality very challenging because some of the children who were lost to follow-up may be deceased and therefore misclassified and the mortality rates that we observed may be an underestimate. The estimated mortality rates under the “worst case scenario” assumption that all children who were lost to follow-up had died was 7 times higher after 1 year of follow-up and 5 times higher after 10 years of follow-up. A study from Malawi found that only 11% of the children who were lost to follow-up had died [13]. This study was done in an urban setting unlike our study setting which is both urban and rural. Our sensitivity analysis approach allowed us to estimate maximum mortality rates which was more informative in our setting.

Transfers are high in this clinic, which is motivated by Ministry of Health policy that encourages people to seek care at the health facility nearest to them. Studies have showed that this is an effective strategy as patients can additionally benefit from support groups and other programs that may be available at their nearest health facility [38, 39]. However, this makes it difficult to study long-term treatment outcomes because patient tracking is impossible once they transfer to another facility.


We observed low attrition due to mortality among children on ART in Zambia. Loss to follow-up was high (16%) and could underestimate the mortality. Mortality was highest during the first 3 months of treatment. Children aged less than one year and children with advanced WHO disease stage had higher mortality. We recommend to program managers and clinicians to develop effective interventions to improve retention in care and strengthen early infant diagnosis of HIV. Pretreatment screening and treatment of opportunistic infections among children commencing ART needs to be strengthened to reduce early mortality.



Acquired Immunodeficiency Syndrome


Antiretroviral Therapy


Human Immunodeficiency Virus


Livingstone Central Hospital


Non-Nucleoside Reverse Transcriptase Inhibitors


Nucleoside Reverse Transcriptase Inhibitors


Pediatric Center of Excellence


Prevention of Mother to Child Transmission of HIV


World Health Organization


  1. Joint United Nations Programme on HIV/AIDS U. On the fast track to an AIDS free generation; the incredible journey of the global plan towards the elimination of new HIV infections by 2015 among children and keeping their mothers alive. Geneva Switzerland: UNAIDS; 2016.

  2. UNICEF. Paediatric care and treatment New York 2018 [Available from: Accessed 17 Nov 2018.

  3. Joint United Nations Programme on HIV/AIDS U. Fact Sheet: Global HIV & AIDS statictics - 2018. 2018 [Available from: Accessed 17 Nov 2018.

  4. Joint United Nations Programme on HIV/AIDS U. 90–90-90 an ambitious target to help end the AIDS epidemic. Geneva Switzerland: UNAIDS; 2014.

  5. Violari A, Cotton M, Gibb D, Babiker A, Steyn J, Madhi S, et al. Early antiretroviral therapy and mortality among HIV-infected infants. N Engl J Med. 2008;359(21):2233–44.

    Article  CAS  Google Scholar 

  6. MOH/CSO/PEPFAR/ICAP/CDC/TDRC/UTH/UNZA. Zambia Population-Based HIV Impact Assessment ZAMPHIA 2015-2016. 2016.

    Google Scholar 

  7. Joint United Nations Programme on HIV/AIDS U. Fact Sheets: Zambia 2017. Geneva, Switzerland2017 [Available from: Accessed 18 Nov 2018.

  8. Bolton-Moore C, Mubiana-Mbewe M, Cantrell RA, et al. Clinical outcomes and cd4 cell response in children receiving antiretroviral therapy at primary health care facilities in Zambia. JAMA. 2007;298(16):1888–99.

    Article  CAS  Google Scholar 

  9. Dijk JH, Sutcliffe CG, Munsanje B, Sinywimaanzi P, Hamangaba F, Thuma PE, et al. HIV-infected children in rural Zambia achieve good immunologic and virologic outcomes two years after initiating antiretroviral therapy. PLoS One. 2011;6.

  10. Davies M-A, May M, Bolton-Moore C, Chimbetete C, Eley B, Garone D, et al. Prognosis of children with HIV-1 infection starting antiretroviral therapy in southern Africa: a collaborative analysis of treatment programs. Pediatr Infect Dis J. 2014;33(6):608–16.

    Article  Google Scholar 

  11. Davies M-A, Gibb D, Turkova A. Survival of HIV-1 vertically infected children. Curr Opin HIV AIDS. 2016;11(5):455–64 Web.

    Article  CAS  Google Scholar 

  12. Koller M, Patel K, Chi BH, Wools-Kaloustian K, Dicko F, Chokephaibulkit K, et al. Immunodeficiency in Children Starting Antiretroviral Therapy in Low-, Middle-, and High-income countries. J Acquir Immune Defic Syndr (1999). 2015;68(1):62–72.

  13. Ardura-Garcia C, Feldacker C, Tweya H, Chaweza T, Kalulu M, Phiri S, et al. Implementation and Operational Research: Early Tracing of Children Lost to Follow-Up From Antiretroviral Treatment: True Outcomes and Future Risks. J Acquir Immune Defic Syndr (1999). 2015;70(5):e160-e1e7.

  14. Mutanga J, Raymond J, Towle M, Mutembo S, Fubisha R, Lule F, et al. Institutionalizing provider-initiated HIV testing and counselling for children: an observational case study from Zambia. PLoS One. 2012;7(4):e29656.

    Article  CAS  Google Scholar 

  15. Kankasa C, Carter RJ, Briggs N, Bulterys M, Chama E, Cooper ER, et al. Routine Offering of HIV Testing to Hospitalized Pediatric Patients at University Teaching Hospital, Lusaka, Zambia: Acceptability and Feasibility. J Acquir Immune Defic Syndr (1999). 2009;51(2):202–208.

  16. Stringer JA, Zulu I, Levy J, et al. Rapid scale-up of antiretroviral therapy at primary care sites in Zambia: feasibility and early outcomes. JAMA. 2006;296(7):782–93.

    Article  CAS  Google Scholar 

  17. Ministry of Health Z. Zambia Consolidated Guidelines for Prevention and Treatment of HIV Infection. 2018.

    Google Scholar 

  18. Newell M-L, Coovadia H, Cortina Borja M, Rollins N, Gaillard P, Dabis F. Mortality of infected and uninfected infants born to HIV-infected mothers in Africa: a pooled analysis. Lancet. 2004;364(9441):1236–43.

    Article  Google Scholar 

  19. Chammartin F, Zürcher K, Keiser O, Weigel R, Chu K, Kiragga AN, et al. Outcomes of patients lost to follow-up in African antiretroviral therapy programs: individual patient data meta-analysis. Clin Infect Dis. 2018;67(11):1643–52.

    Article  Google Scholar 

  20. Team RC. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2017.

  21. Alboukadel Kassambara MK. survminer: Drawing Survival Curves using 'ggplot2'. R package version 040999. 2017; Accessed 2 May 2017.

  22. Therneau T. _A Package for Survival Analysis in S_. version 2.38. 2015. Accessed 27 Aug 2017.

  23. van Dijk JH, Sutcliffe CG, Munsanje B, Sinywimaanzi P, Hamangaba F, Thuma PE, et al. HIV-infected children in rural Zambia achieve good immunologic and Virologic outcomes two years after initiating antiretroviral therapy. PLoS One. 2011;6(4):e19006.

    Article  Google Scholar 

  24. Brophy JC, Hawkes MT, Mwinjiwa E, Mateyu G, Sodhi SK, Chan AK. Survival outcomes in a pediatric antiretroviral treatment cohort in southern Malawi. PLoS One. 2016;11(11):e0165772.

    Article  Google Scholar 

  25. Cohen S, Smit C, van Rossum AMC, Fraaij PLA, Wolfs TFW, Geelen SPM, et al. Long-term response to combination antiretroviral therapy in HIV-infected children in the Netherlands registered from 1996 to 2012. AIDS. 2013;27(16):2567–75.

    Article  CAS  Google Scholar 

  26. Wamalwa DC, Obimbo EM, Farquhar C, Richardson BA, Mbori-Ngacha DA, Inwani I, et al. Predictors of mortality in HIV-1 infected children on antiretroviral therapy in Kenya: a prospective cohort. BMC Pediatr. 2010;10:33.

    Article  Google Scholar 

  27. Walker AS, Prendergast AJ, Mugyenyi P, Munderi P, Hakim J, Kekitiinwa A, et al. Mortality in the year following antiretroviral therapy initiation in HIV-infected adults and children in Uganda and Zimbabwe. Clin Infect Dis. 2012;55(12):1707–18.

    Article  CAS  Google Scholar 

  28. Njom Nlend AE, Loussikila AB. Predictors of mortality among HIV-infected children receiving highly active antiretroviral therapy. Med Mal Infect. 2017;47(1):32–7.

    Article  CAS  Google Scholar 

  29. Braitstein PBM, Dabis F. Mortality of HIV-1-infected patients in the first year of antiretroviral therapy: comparison between low-income and high-income countries. Lancet. 2006;367(9513):817–24.

    Article  Google Scholar 

  30. Cotton M, Violari A, Otwombe K, Panchia R, Dobbels E, Rabie H, et al. Early time-limited antiretroviral therapy versus deferred therapy in south African infants infected with HIV: results from the children with HIV early antiretroviral (CHER) randomised trial. Lancet. 2013;382(9904):1555–63.

    Article  Google Scholar 

  31. Ministry of Health Z, National AIDS Council, Zambia. Zambia Country Report: Monitoring the declaration of commitment on HIV/AIDS and the Universal Access. Lusaka, Zambia; 2015 5/2/2017.

  32. Marston M, Becquet R, Zaba B, Moulton L, Gray G, Coovadia H, et al. Net survival of perinatally and postnatally HIV-infected children: a pooled analysis of individual data from sub-Saharan Africa. Int J Epidemiol. 2011;40(2):385–96.

    Article  Google Scholar 

  33. Luzuriaga K. Early combination antiretroviral therapy limits HIV-1 persistence in children. Annu Rev Med. 2016;67(1):201–13.

    Article  CAS  Google Scholar 

  34. Wagner A, Slyker J, Langat A, Inwani I, Adhiambo J, Benki-Nugent S, et al. High mortality in HIV-infected children diagnosed in hospital underscores need for faster diagnostic turnaround time in prevention of mother-to-child transmission of HIV (PMTCT) programs. BMC Pediatr. 2015;15:10.

    Article  Google Scholar 

  35. Central Statistics Office MoH, Tropical Disease Research Centre University of Zambia and Macro International Inc. Zambia Demographic and Health Survey Maryland, USA: ICF International Rockville; 2014. p. 13–4.; 2013-2014.

  36. Mundi I. Zambia-Mortality rate 2015 [Available from: Accessed 21 Aug 2017.

  37. Boulware DR, Callens S, Pahwa S. Pediatric HIV immune reconstitution inflammatory syndrome. Curr Opin HIV AIDS. 2008;3(4):461–7.

    Article  Google Scholar 

  38. Escamilla V, Chibwesha CJ, Gartland M, Chintu N, Mubiana-Mbewe M, Musokotwane K, et al. Distance from household to clinic and its association with the uptake of prevention of mother-to-child HIV transmission regimens in rural Zambia. J Acquir Immune Defic Syndr (1999). 2015;70(3):e94-e101.

  39. Sam-Agudu NA, Ramadhani HO, Isah C, Erekaha S, Fan-Osuala C, Anaba U, et al. The Impact of Structured Mentor Mother Programs on Presentation for Early Infant Diagnosis Testing in Rural North-Central Nigeria: A Prospective Paired Cohort Study. J Acquir Immune Defic Syndr (1999). 2017;75:S182-S1S9.

Download references


We sincerely thank the management and staff of Livingstone Central Hospital. Special thank you to the wonderful clinical team and caregivers at Livingstone Pediatric Center of Excellence. We acknowledge the contribution of Dr. Nathan Hansen of the University of Georgia, for his valuable contribution to this work. We wish to thank our colleagues: Dr. Juliet N. Sekandi, Allan K. Nkwata, Maria Eugenia Castellanos, Dr. Leonardo Martinez and Dr. Robert Kakaire for their contributions during the Epidemiology in Action research group meetings at the Department of Epidemiology, University of Georgia.


This work was supported in part by a fellowship award to J.N.M from the Schlumberger Foundation Faculty for the Future Fellowship.

Availability of data and materials

The data that support the findings of this study are available from Livingstone Central Hospital but restrictions apply to the availability of these data, which were used under license for the current study and are currently not publicly available. Data are however available from the authors upon reasonable request and with the permission of the Zambia National Research Authority.

Author information

Authors and Affiliations



JNM (study concept and data analysis). JNM, S.M, CW, PET, AEE, XS (study design, development of data collection tools, data interpretation, manuscript writing and manuscript review). JNM, SM, DS, BS, KMK, RCF and JC (study implementation, data collection, data cleaning, manuscript writing and critical revision the manuscript). AEE, XS, SM, CW (critical revision of data analysis and interpretation). All the authors read and approved the manuscript for publication.

Corresponding author

Correspondence to Jane N. Mutanga.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Zambia National Health Research Authority and the institutional review boards at Macha Research Trust and the University of Georgia. We conducted analysis of anonymized routinely collected program data, hence all the institutional review boards waived informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. Directed Acyclic Graph (DAG) showing adjusted covariates of ART Initiation among Children on ART. (DOCX 108 kb)

Additional file 2:

KM plots for 1) Figure S4E: Anemia, 2) Figure S4F: Mom took ART, 3) Figure S4G: Wasting (low weight at baseline, 4) Fig. 4h: Baby took NVP for Prophylaxis after delivery. (DOCX 66 kb)

Additional file 3:

Trends in treatment outcomes among 1039 children on ART at Livingstone Central Hospital, Zambia: 2003–2015. (DOCX 45 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mutanga, J.N., Mutembo, S., Ezeamama, A.E. et al. Long-term survival outcomes of HIV infected children receiving antiretroviral therapy: an observational study from Zambia (2003–2015). BMC Public Health 19, 115 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: