Abstract
Background:
Although there is a presence of governmental and non-governmental organizations running to provide quality HIV care services to reduce HIV-related mortality, there is rapid disease progression and death among children in developing countries including Ethiopia. Thus, this study was aimed to assess the mortality predictors of children living with HIV at Bahir Dar town public health facilities.
Method:
A facility-based retrospective follow-up study was conducted among 588 children who were enrolled in the HIV care clinic from 1 September 2010 to 30 August 2019. Data were entered into the Epi-Data entry 3.1 and then exported to STATA version 14 for analysis. Multiple imputation models were employed to handle missing data using the multivariate imputation Chained Equations technique. The Kaplan–Meier survival curve and log-rank test were used to estimate and compare the survival time of categorical variables.
Result:
About 27 (4.6%) (95% confidence interval: 2.9–6.5) deaths were observed from the 30,062.3 person-months follow-up period, and the overall incidence density rate of 0.9 per 1000 child-months (95% confidence interval: 0.6–1.3). Advanced WHO clinical stage (adjusted hazard ratio = 3.18; 95% confidence interval: 1.07–9.43), hemoglobin level less than 8 g/dL (adjusted hazard ratio = 3.54; 95% confidence interval: 1.27–8.85), children having a weight for age of <−2z (adjusted hazard ratio = 2.81; 95% confidence interval: 1.19–6.6), children with poor adherence (adjusted hazard ratio = 3.91; 95% confidence interval: 1.41–10.8), and starting the treatment beyond 1 week of being eligible (adjusted hazard ratio = 3.22; 95% confidence interval: 1.21–8.53) were predictors of HIV-related mortality among children initiated antiretroviral therapy.
Conclusion:
The hazard of mortality was higher among HIV-infected children in the early period of initiation. Enhancing antiretroviral therapy drug adherence, monitoring Hgb level, and timely initiation of antiretroviral therapy reduce HIV-related mortality.
Introduction
According to the 2018 UNAIDS report, globally there was a total of 1.8 million children living with HIV with an estimated 110,000 AIDS-related deaths. In Ethiopia, the epidemic estimate of children living with HIV and HIV-related death among children was 62,000 and 3600. 1 There were evidence from different countries regarding the survival status of children living with HIV. From the tertiary hospital of Benin, 76.4% of survival probability after 60 months of antiretroviral therapy (ART) initiation was evidenced from a 5-year survival trend and outcome among children living with HIV. 2 The median survival time for children living with HIV in Swaziland was 78 months. 3 The mortality rate from the observational study in Zambia was 1.6/100 person-years of observation. 4 In Ethiopia, mortality rates of children living with HIV after initiating ART ranged from 4.8% to 22.9%.5,6
The pre-treatment mortality rate was high among children eligible for ART.7,8 Prompt initiation of ART within the beginning week of eligibility was associated with better treatment outcome in terms of virologic outcome, 9 but its effect in mortality aspect is not addressed; this research undergoes ascertaining its effect on mortality.
The Ethiopia government has adopted UNAIDS treatment 90-90-90 targets which aimed at ending the HIV epidemics by 2030 with specific targets, 90% of individuals living with HIV know their HIV status, 90% of those eligible are initiated ART, and 90% of individuals on ART will have viral load suppression by 2020. 10
Guidelines that Ethiopia have been implementing since 2014 consider that all children living with HIV should start ART regardless of CD4 count and World Health Organization (WHO) clinical stage. 11 But the guideline that was implemented before had different eligibility criteria like WHO clinical stage 3 and 4 regardless of CD4, WHO clinical stage 1 and 2 with CD4 below the threshold of %CD4+ ⩽25, or CD4 absolute count of ⩽750 cells/mm3 all children with the age of 2–5 years. 12
While researches are conducted to look at survival status in the different areas of Ethiopia, most of them considered the previous guideline. In resource-limited settings like Ethiopia, there is a pressing need for research to identify factors and further refine HIV treatment strategies among children to have better survival after initiation of ART. The survival outcomes might be affected by different potential factors; however, those factors are poorly understood in the study area, especially after the implementation of the new guideline that was launched in 2014. The finding of this study will help the researchers to uncover the different areas of factors by going further based on the factors and outcomes identified in this study and to achieve a good treatment outcome by understanding the individual variability in the treatment outcome and acting on the modifiable factors. This study assessed the survival status and predictors of mortality among children living with HIV and who initiated ART in Bahir Dar town public health facilities.
Methods and materials
Study design and setting
A retrospective cohort study was conducted to assess survival status and predictors among children living with HIV and initiated ART in Bahirdar town public health facility. Follow-up time for each patient was calculated from the date of ART initiation to the date of file closure which was due to death, transfer out, and loss to follow-up and the last clinic visit before 30 August 2019.
The study was conducted in Bahirdar city from 1 September 2010 to 30 August 2019. Bahirdar is the administrative capital of Amhara regional state and located 565 km northeast of Addis Ababa, Ethiopia’s capital.
Study population
All medical records of children living with HIV aged less than 15 years old at initiation and who were enrolled from 1 September 2010 to 30 August 2019 at Bahir Dar city public health facilities were study participants in this follow-up study.
Inclusion criteria
All children living with HIV enrolled on ART from 1 September 2010 to 30 August 2019, and if ART was initiated in the selected health facilities.
Exclusion criteria
Patients with incomplete baseline information (at least the child’s socio-demographics, initiation date, and exit date from the cohort).
Sample size and sampling procedure
The sample size was determined using double population proportion from factors in different studies using Epi info version 7 assuming: 95% confidence interval, power 80%.
The sample size for this study was 588, CD4 percentage less than 10% as a main explanatory variable for HIV-related death gave a higher sample size from the information available in the previous literature. 5
Sampling procedure
From all public health facilities providing ART services in Bahir Dar city, a random sampling method was employed. The sampling frame (808) was constructed at selected health facilities. Then from the sampling frame constructed, participants who fulfilled the inclusion criteria were selected randomly using excel (computer-generated random sample).
The allocation of the study subjects to each health facility was computed based on the proportion of the number of medical records of HIV-positive children.
Data collection tools and procedures
In the checklist, socio-demographic characteristics, baseline clinical, laboratory, and ART information were included. Data were extracted through chart review using a structured checklist, adapted from ART intake, and follow-up form. Before the actual data collection, data collectors and the supervisor were oriented about the objectives of the study and contents of the tools. Then, all medical records of children living with HIV that fulfilled the inclusion criteria were reviewed by four trained nurses as data collectors and 1 BSc nurse as a supervisor.
Data quality assurance
Both supervisors and data collectors were closely following the data collection process. Consistency cross-checking was conducted in Bahir Dar health center for 5% of the study population 2 weeks before the actual data collection to evaluate the validity of the checklist. Every day all checklists were reviewed and checked during the data collection and any errors were corrected accordingly.
Statistical analysis
All filled checklists were entered into EpiData entry 3.1 and checked for completeness and consistency and exported to STATA version 14. The data were processed by STATA version 14 to estimate the survival time of children living with HIV.
Cox proportional hazard model assumption was checked graphically using log-log and Kaplan–Meier versus predicted survival plot and by running a global test using the Schoenfeld residual test. Bivariable Cox proportional hazard model was employed. Variables found to have a p-value less than 0.2 in the bivariate analysis were a candidate to the multi-variable Cox proportional hazards regression model. Variables having a p-value ⩽0.05 were considered statistically significant.
Missing values were identified for some of the variables. The percentage of missing values across all variables varied between 0% and 6.1%. The missing data were expressed in terms of proportion, mechanisms under which the missing data occurs, and missing data patterns, to apply appropriate procedures to handle the existing missing data. The missing mechanism of the data was found to be missing at random by using little’s test and the need to do imputation was ascertained. Multiple imputations were used to create and analyze multiply imputed data sets. Complete variables were imputed under fully conditional specification, by using multivariate imputation Chained Equations (MICE) technique. For comparison, a complete case analysis was also performed.
Operational definition
Results
Baseline socio-demographic characteristics
Among 588 patients, 308 (52.4%) of study participants were females and the majority (79.4%, n = 467) of the participants were urban residents. The mean age at the commencement of ART was 7.2 (± 3.9) years with extremes of 0.5 and 14 years (Table 1).
Baseline socio-demographic characteristics of children on ART at Bahirdar city public health facility, Amhara regional state, Northwest Ethiopia (2010–2019).
ART: antiretroviral therapy.
Baseline clinical, laboratory, and ART information
Regarding opportunistic infection, 188 (32%) of the children experience an opportunistic infection. 26.29% had recurrent upper respiratory tract infection (URTI) followed by bacterial pneumonia (19.5%) and herpes zoster (15.5%). About 60.4% of children start the treatment within 1 week of eligibility. The majority of the children have had a regimen change and treatment failure and drug side effects account for 14.4%, and 10% of the reason for regimen change (Table 2).
Baseline clinical, laboratory, and ART information of children on ART at Bahir Dar city public health facility, Amhara regional state, Northwest Ethiopia (2010–2019).
ART: antiretroviral therapy; WHO: World Health Organization; CD4: cluster of differentiation 4; n: number; INH: isonicotinic acid hydrazide/isoniazide; W/A: weight for age; W/H/L: weight for height/length; H/A: height for age.
Among 588 patients, 362 (61.6%) of study participants have initiated their treatment after 2014. From the whole cohort, more than half (60%) of the children living with HIV are live on ART (Figure 1).

Outcomes of HIV-positive children initiated ART at Bahir Dar city public health facility, Amhara regional state, Northwest Ethiopia (2010–2019).
Survival status after initiation of ART
A total of 588 children living with HIV enrolled in care were followed and given 30,062.3 person-months observation. When the survival was stratified by the turning point for guideline change happened in 2014 (before and after 2014), there was no significant difference, but the proportion of deaths before and after 2014 was 5.7% and 3.9%, respectively. At the end of the follow-up, more than half (60%) of the patients were retained in care in those public health facilities, and they initiate their treatment. The median follow-up time was 51 months. In this study, 4.6% (95% confidence interval, CI: 2.9–6.5) of the study participant died during the follow-up period. The death incidence rate was calculated using person-months of follow-up. The overall incidence density rate (IDR) of the cohort was 0.9 per 1000 child-month observations (95% CI: 0.6–1.3). About half (48.2%) of the deaths occurred within the first 6 months of ART initiation which gives a mortality incidence of 3.9 per 1000 child months. The probability of death was estimated to be 2.3%, 2.7%, and 3.1% at 6, 12, and 24 months after initiation, respectively. The cumulative probabilities of survival at 6, 12, and 24 months of ART initiation were found to be 97.7%, 97.3%, and 97.1%, respectively (Figure 2).

Cumulative Kaplan–Meier survival curve with 95% confidence intervals of children on ART at Bahir Dar city public health facility (2010–2019).
Predictors of mortality after ART initiation
In the multivariable Cox proportional hazard adjusted model, five factors were identified (Table 3). The hazard of child mortality in children with advanced WHO clinical stage was 3.2 times higher when compared with those who were on WHO clinical stage 1 and 2 at initiation (adjusted hazard ratio (AHR): 3.2; CI: 1.1–9.4) (Figure 3).
Bivariable and multivariable Cox regression among HIV-positive children on ART at Bahir Dar public health facility Amhara regional state, North West Ethiopia (N = 588) (2010–2019).
ART: antiretroviral therapy; CHR: crude hazard ratio; AHR: adjusted hazard ratio; CI: confidence interval; Others in the initial regimen: (ABC-3TC-LPV, TDF-3TC-ATV, and ABC-3TC-LPV); W/A: weight for age; n: number; WHO: World Health Organization; NVP: Nevirapine; EFV: Efavirenz.
Significant (p-value < 0.05).

Kaplan–Meier survival estimate of mortality-free survival proportion based on baseline WHO clinical stage at enrollment at Bahir Dar public health facility (2010–2019).
The hazard of child mortality in children with poor adherence was 3.9 times higher among children when compared with those with good adherence ((AHR: 3.9; CI: 1.4–10.8). The hazard of child mortality was 2.8 times higher among children who had a low W/A (<−2z) when compared with those who had ⩾−2z (AHR: 2.8; 95% CI: 1.2–6.6). The hazard of mortality in children who had a hemoglobin level less than 8 g/dL was 3.5 times higher when compared with those who had hemoglobin levels greater than 8 g/dL (AHR: 3.5; 95% CI: 1.3–8.9). Furthermore, the hazard of mortality in children who started the treatment after 1 week of being eligible was 3.2 times higher when compared with those who started within 1 week (AHR: 3.2; 95% CI: 1.2–8.5) (Figure 4).

Kaplan–Meier survival estimate of mortality-free survival proportion based on the length of time from eligibility to initiation at Bahirdar public health facility (2010–2019).
Discussion
The findings of this study indicated that the overall incidence of death was 0.9 per 1000 child-months observation and cumulative incidence of 4.6% (95% CI: 2.9–6.5). The cumulative incidence finding was consistent with the comparative analysis of outcomes of ART in children in East Africa (4.3%), Asia (5.4%), and South Africa (5.7%). 14 But this finding is lower than study findings from West Africa (7.4%), 14 Dessie referral hospital, and Debre tabor general hospital (22.9%), 6 and studies conducted in Addis Ababa (10.4%). 15 This difference may be due to the study period variability since there was a change in treatment modality and care of children on ART starting from eligibility criteria to medication changes.
The overall incidence rate of death was found to be 0.9 per 1000 child-months. This finding is lower than studies conducted in China, 16 Asia, 17 South Africa, 18 Zambia, 19 Wolaita zone health facilities, Ethiopia. 20 The possible reason may be the variety in study populations; some studies consider only younger age groups, which tend to have a high incidence rate of mortality.
About half of the deaths occurred in the first 6 months of ART. This is consistent with surveillance data for children less than 5 years of age receiving ART in 48 HIV/AIDS treatment programs in Africa and Asia. 21 It is also similar to the finding from Addis Ababa, in which 70% of deaths occurred in the first 6 months of ART. 15
Related to the probability of death, the probability of death was estimated to be 2.3%, 2.7%, and 3.1% at 6, 12, and 24 months after initiation, respectively. This finding was consistent with a systematic review report from 29 articles in African countries (5%) at 6 months and(6%) at 12 months of initiation 22 and from 70 health facilities, Ethiopia, which was 3.4% at 6 months, 4.1% at 12 months, and 4.8% at 24 months. 23
In this study, the overall cumulative survival probability was found to be 92.5% at 138 months after ART initiation. The survival probability was also found to be 97.7% and 97.3%, at the first, 6th, and 12th months of initiation, this finding was higher than studies done in Adama referral hospital (93.9%) at 6 months and (92.8%) at 12 months 24 and study conducted in Debre Tabor referral hospital and Dessie referral hospital (85%) at 12 months. 6 This difference might be due to the characteristics of the population included in our follow-up. In this study, most of the study participants initiated ART regardless of any criteria; they started the treatment as soon as their status is known. It is well known that those children who have initiated ART at an early phase had good clinical and immunological improvement and the risk of severe morbidity is lower in those children.25,26
Baseline advanced WHO clinical stages at the initiation of ART was found to be a significant predictor of mortality among HIV-positive children on ART. The hazard of child mortality was higher among children with advanced WHO clinical stage compared with WHO clinical stages 1 and 2 at the initiation. This finding is supported with other studies conducted in China (HR = 2.4), 16 sub-Saharan Africa, 22 Adama, 24 Mekelle, 5 and Zambia (WHO stage 4 (AHR = 4.8), compared to WHO stage 1 and 2). 4 The mortality is high when the disease is advanced more because the viral load is high and it destroyed body defense mechanisms and exposes to several opportunistic infections. It finally results in death.
Mortality was closely linked to the children’s level of adherence to ART. The hazard of mortality among HIV-positive children who had end-line poor treatment adherence was higher than those children who had good treatment adherence.
This finding is supported by several studies. A study in Benin revealed an increased hazard of death among patients with poor adherence to therapy. 2 Also, a study in Addis Ababa, Ethiopia, provides evidence that suboptimal ART adherence increases the risk of mortality. 15 It is due to the development of drug resistance and an increase in the risk of virologic failure and subsequent disease progression.27,28
In this study, the hazard of mortality was higher among HIV-positive children who had a baseline (hemoglobin level of less than 8 g/dL) than those children who had greater than 8 g/dL during initiation of ART. This finding is supported by other studies that were conducted in Mekelle, Ethiopia, 5 and Addis Ababa, Ethiopia. 29 This could be due to the impact of anemia in HIV-infected patients on their physical functioning and quality of life. 30
The hazard of mortality in children who did not start the treatment within 1 week of being eligible was higher when compared with those who started within 1 week. This finding is supported by different studies done in resource-limited settings as mortality was associated with delays between clinic entry and ART initiation there was evidence of an increase in mortality from 11.0% to 14.7% with a 10-week delay in ART for patients entering care. 31 It could be because of the benefits of immediate initiation in terms of good treatment response as there is evidence of low risk of treatment failure in infants. 32 Prompt initiation might help in decreasing the burden of anemia and its subsequent complications. 33 Only the discussed factors were identified in our study, but there was evidence that children with young age (age < 1.5 years), 34 age group of <1 year, 3 age less than 18 months and children with chronic diarrhea 5 had a high chance of mortality and children who took isoniazide prophylaxis and Cotrimoxazole prophylactic therapy 35 and participants who had psychosocial support during follow-up 6 had higher survival probability. The survival status of HIV-positive children needs to be studied with a prospective study design to include different areas of factors, including the reasons for delay after being eligible for ART.
However, the follow-up period was long, which increases the patient’s period of observation and enabled us to know the long-term impact of ART on mortality and appropriate measures taken to handle the missing value; the present study does have some limitations. As the data were collected from medical records, socioeconomic factors were not possible to collect for the whole cohort, as they might be associated with mortality; therefore, this variable was not included in analyses.
Conclusion and recommendation
The survival rate of children living with HIV was higher as compared to other previous studies. The hazard of mortality was relatively higher among HIV-positive children in the early period of initiation. The risk of mortality is increased if the child was at a young age at the initiation of ART, lower baseline CD4 count, advanced WHO Clinical stage, lower hemoglobin values, and who have had poor adherence. However, prompt initiation of ART within 1 week of being eligible may minimize mortality. There was a considerable delay that can affect the long-term survival outcome, so health facilities that have been giving ART service better initiate ART as soon as the case is confirmed followed by strengthening careful and regular monitoring.
Footnotes
Acknowledgements
Our gratitude goes to supervisors, the data collectors, and the staff at the ART clinic of Felege hiwot comprehensive and specialized hospital, Bahir Dar, Abay, and Han health centers. We would also like to thank Mr Getnet Dessie (MSc, Adult health nursing) who supported us a lot in model development and the technical issues related to the software.
Author contributions
B.C. prepared the proposal, provided training for data collectors, supervised data collection, analyzed and interpreted the data, and wrote the main manuscript. A.B. supervised and evaluated the overall research process and reviewed the whole manuscript. E.A. prepared the proposal and reviewed the manuscript. A.G. reviewed the manuscript. A.T. prepared graphs and tables.
Availability of data materials
The datasets analyzed during the current study are available from the corresponding author on reasonable request.
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Ethical approval
Ethical approval for this study was obtained from Research, Ethical Clearance Committee of Bahirdar University with a protocol number of 0046/2020. Following the approval, the official letter was obtained from the College of Medicine and Health Sciences of Bahir Dar University and given to the hospital managers, medical directors, and ART focal persons. Permission was obtained from the hospital’s administration and ART focal persons before data collection.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research has been financially supported by Wolkite University for data collection.
Informed consent
Written informed consent was not sought from the legally authorized representatives of the patients directly since we used a medical record.
